In the BiCMOS process, the trade-off in specifications of all kind of circuits can be fully utilized for possible circuit designs. In this paper, the trade-off between the NF and linearity for LNA circuit design have been tested and verified. A good choice in the first stage of the LNA design is an HBT for a better NF, and the last stage is an MOS for better linearity. Furthermore, the circuits (HBT-MOS) designed in this paper are very useful for integration with other subcircuits in 5-6-GHz receivers IC due to its powersaving characteristic and better figure of merit. With the modern microelectronic industry entering a new era featuring progressively higher speed and more highly integrated mixed-signal systems on chip, accurate modeling of distributed electromagnetic behavior of sections, which may include signal traces, power planes, substrates, and so forth, is crucial during the design flow. Driven by broadband and nonlinear circuit design [1][2][3][4][5], time-domain coupled electromagnetic and circuit simulation has been gaining popularity. In particular, integral-equation methods [6 -14] have proven to be useful since radiation conditions are built in and only surface meshing and modeling is needed. In this work, the time-domain integral equation (TDIE) method is enhanced to model realistic loss behavior in substrates, which is a crucial component for quality-factor prediction of integrated passives, as well as potentially for thermal pattern prediction. The EM simulation environment in emerging and future 3D integrated circuits is characterized by multiple piecewise homogeneous regions. Each region is comprised of lossy materials, which may or not have a strong impact on the overall performance. In general, for higher speed and higher sensitivity systems such as RF and analog methods, the loss is a crucial factor in determining system specifications and performance. The TDIE approach has been shown to work with lossy material [15][16][17][18] wherein the Green's functions, besides possessing delta functions in time, also include a broadly exponentially decaying 'wake'. This leads to the implementation issue that the spatial integrals at retarded times have to be replaced by temporal convolutions. This convolution brings not only coding complexity, but also increases the computational cost dramatically [15,17].In the presented work, the multiregion substrate geometry is addressed with a TDIE solver. An equivalent-surface approach for each region is used along with the following method to tackle the lossy medium Green's function. To circumvent the difficulty caused by the convolution in time, the decaying 'wake' in the Green's function is approximated by a sum of decaying exponentials via Prony's method. It is shown that for circuit dimensions and realistic losses, the decaying 'wake' changes slowly with the variation of the distance between the source and observation point. This allows the creation and use of an exponential fitting table for discrete distances. The use of the exponential models permits t...