Open access

Using Numerical Simulations to Study and Design Semiconductors Devices in Micro and Nanoelectronics

Written By

Salvador Pinillos Gimenez and Marcello Bellodi

Submitted: 22 November 2011 Published: 19 September 2012

DOI: 10.5772/48334

From the Edited Volume

Numerical Simulation - From Theory to Industry

Edited by Mykhaylo Andriychuk

Chapter metrics overview

2,393 Chapter Downloads

View Full Metrics

1. Introduction

The survival of the integrated circuits(ICs) industries is directly linked to its ability to innovate in terms of the developing new semiconductor devices (planar and three-dimensional) and new manufacturing process technologies (Silicon-On-Insulator, etc), in theircapability to reduce the dimensions of semiconductor devices focusing on the increasing of integration capability of the devices into a single integrated circuit (IC) and in their efficiency to develop new digital and analog ICs applications to attend the market in different forms of consumption. In this context, the numerical simulators represent a very important role for the ICs industries and very high investments have been done in this area of the process and device simulation tools (Technology Computer-Aided Design, TCAD) in order to use the computer simulations to develop and optimize semiconductor processing technologies and devices [1-2].

Significant costs reductions can be achieved through the use of numerical simulators for developing new semiconductor manufacturing technologies and devices that are associated with ICs industries. Thinking on developing new devices, the numerical simulations are able to supply a lot of information about the devices in terms of theirs electrical characteristics curves (drain current as a function of gate voltage for the MOSFETs, etc) and mainly in terms of the internal electrical and physical behaviors of the semiconductor devices (current density, electric field density, potential distribution along of the device structure, mobile carriers density along of the channel width and length of a MOSFET, etc) previously to the physical implementation, which is impossible to be obtained, visualized and understood with only the electrical characterization of the physical devices [1-2]. The involved costs with numerical simulations are associated mainly with human resources and infra-structure (hardware and software) while that, if new devices developing were associated only with the manufacturing process use, would involve very high production costs, very high lead time and many uncertainties about the efficiency of the proposed devices [1-2].

Currently there is two-pronged research in the world focusing on developing new semiconductor devices [3]. The first one is based on designing of the new tridimensional (3D) devices, such as the multi-gates MOSFETs (FinFETs, Three-gates, Four-gates, Surrounding-gate) that involve high cost and investments related to the new techniques of ICs manufacturing process andnew materials. The other tries to explore the many real improvement opportunities of the current planar ICs manufacture technology. One of them, it is to use the “drain/source and channel region interface engineering” (D/SCRIE) approach [4-6], where it can implement new and innovative layout styles with non-conventional geometries, focusing on enhancing the device performance and consequently the performance of analog and digital ICs applications [4-7]. An example is very promising of this layout approach is the “Diamond MOSFET” that it was carefully and specially designed to improve the longitudinal (parallel) electric field (ε//) in order to increase the mobile carriers drift velocity in the channel region (ν//) and consequently the drain current (IDS), transconductance (gm) and to reduce the on-state drain/source series resistance (RDS_on) by a simple changing of the gate geometry from the conventional (rectangular) to hexagonal, as illustrated inFigure 1, showing an example of a SOI nMOSFET3D structure [4-7].

Figure 1.

Example of a SOI nMOSFET 3D structure.

InFigure 1, the shorter and longer channel dimensions are represented by b and B, respectively, W is the channel width and Leffis the effective channel length, α angle is the angle of the drain/channel and the channel/source interfaces due to hexagonal gate geometry, tox, tSIand tBOXare gate oxide, silicon film and buried oxide thickness, respectively.

Figure 2presents the Diamond SOI nMOSFET layout style, where the longitudinal corner effect (LCE) acting in the point P in the channel is shown and the resultant longitudinal electric field (ε//) is given by the vector sum of the longitudinal electric field components (ε//1andε//2). Last is higher than the one observed in the conventional SOI nMOSFET counterpart. These longitudinal electric field components (ε//1andε//2) are generated by the drain bias (VDS) and are perpendiculars to the metallurgical junctions of the drain/silicon film (channel region) [4-7].

The effective channel length (Leff) of the Diamond SOI nMOSFET is given by theequation (1).

Leff=B+b2E1

Figure 2.

Example of a SOI nMOSFET 3D layout showing the LCE acting in the point P.

In this structure, the resultant longitudinal electric field (ε//) in the center of the gate region is higher than the one observed in the edges of the gate region, due to the higher interaction between the longitudinal electric field components (ε//1andε//2) in the center of the device (LCE) [4-7] than the one found in the edges of the gate region, and therefore the drain current tends to flow more in the center of the gate region than in the edges of the gate region. Besides that, the Diamond structure can be electrically represented by the parallel association of the N transistors, where the channel width of each transistor of this parallel association is given by W/N and the channel length of these transistors varies from b to B. As a consequence of this last feature, the drain current tends to flow more by the edges than by the center of the gate region. This effect is defined here by “Variable Channel Lengths of Transistors Associated in Parallel Effect” (VACLETAPE). So, in the Diamond structure occurs two effects at the same time, i.e., the LCE and the VACLETAPE generating a Diamond drain current higher than the one observed in the conventional (rectangular) SOI nMOSFET counterpart, regarding the same gate area (AG) and bias conditions.

Advertisement

2. Devices used to perform the study

So, in order to qualify and quantify the benefits of Diamond layout style to improve the transistor performance, it were implemented three pairs of SOI nMOSFETs, being three Diamond transistors with different α angles and other three conventional SOI nMOSFETs counterparts, regarding the same gate area (AG). Table I presents the devices dimensions used to implement these devices in the DevEdidt3D from (TCAD/Silvaco Inc) [8]. A similar TCAD to implement devices is the Sentaurus Structure Editor of the Synopsys, Inc [9].

Conventional SOI nMOSFETDiamond SOI nMOSFET
#W (µm)L (µm)W/L#b (µm)B (µm)(o)
CSM16.07.00.86DSM11.013.053.1
CSM26.04.01.50DSM21.07.090.0
CSM36.02.52.40DSM31.04.0126.9

Table 1.

Devices Dimensions used to the study.

The constructive characteristics of these devices are: gate-oxide (tox), silicon-film (tSI) and buried-oxide (tBOX) thickness are 2 nm, 100 nm and 400 nm, respectively, and drain/source and channel doping concentrations are 5.5x1017cm-3and 1x1020cm-3, respectively.

Advertisement

3. Three-dimensional simulations results

The ATLAS (Device Simulation Software) [9] was used to perform the 3D numerical simulations of the Diamond and Conventional SOI nMOSFETs. The curves of IDSversus gate voltage (VGS) and IDSversus drain voltage (VDS) were obtained. Based on these curves, the main parameters of these devices were extracted and compared in order to verify the advantages and disadvantages of the Diamond layout style in relation to the conventional one counterpart.

Basically, the ATLAS device simulator uses two files inputs: the first one is the text file that contains commands for the ATLAS. Usually, these commands are to include a file that define the devices structure (usually built by using DevEdit3D [8]), to define the dimensions, the terminals and the meshof the structure of the device. The second one is responsible to define the physic phenomenon models (mobility, band width, impact ionizations, etc), the solution methods and the bias commands of the terminals that must be taking in account in the simulation of the device [1].

Besides that, ATLAS generates three types of output files. The first one is the “run-time output file” that provides all the error and warning messages. The second one is the “log files”, which store allterminal voltages and currents of the device [direct current (DC), small signal alternate current (AC) and transient] generated the command SOLVE. This last file can be seen with the TonyPlot (visualization tool) to visualize the electrical behavior of the device. The last one is the “solution files”, and it is responsible to stores the two-dimensional and three-dimensionaldata of the solution variables within the device for a single bias point (electric field, potential, etc.). This last file can be visualized by the use of TonyPlot too [1].

The following models are used to perform the 3D numerical simulations of the devices: I. Lombardi’s vertical and horizontal electric-field-dependent mobility model, which includes effects of low (surface acoustic phonon scattering, optical inter-valley scattering and surface roughness) and high (Thonber’s drift velocity model to calculate the velocity saturation) electric field; II. The Aurora model takes into account mobility degradation due to the lattice temperature; III. The Fowler-Nordheim Tunneling (electron and holes) model considersthe tunneling of electrons from the semiconductor Fermi level into the insulator conduction band, when high electric field across the gate is applied; IV. The Lucky-Electron Hot Carrier Injection model, considers the electrons emissions into the oxide by first gaining enough energy from the electric field in the channel to surmount the insulator/semiconductor barrier [1].

Figure 3presents the three-dimensional numerical (3D) simulations results of the DSMs and CSMs IDSas a function of VGS,for VDSequal to 10 mV (Triode Region).

Figure 3.

Curves of three-dimensional numerical simulations of the IDSas a function of VGSfor VDS=10mV.

InFigure 3, IDS_=127o, IDS_=90oand IDS_=53oare respectively the differences between the IDSof the DSMs with different angles (127o, 90oand 53o) and the respective CSMs counterparts.

By analyzing theFigure 3, considering the same AGand bias conditions in the Triode Region, it can be observe that DSMs IDSwith angles equals to 127o, 90oand 53oare respectively 20.3 %, 37.2 % and 126.5 % (2.26 times) higher than those observed in the CSMs counterparts and these differences decreaseswhile the DSM angle increases, due the two effects existents in the Diamond SOI nMOSFET (LCE and VACLETAPE). So, this innovative layout style can be used to increase the fun-out and fun-in capability of the logic gates.

Figure 4presents the three-dimensional numerical (3D) simulations results of the DSMs and CSMs gmas a function of gate voltage (VGS),for drain voltage (VDS) equal to 10 mV.

Figure 4.

Curves of three-dimensional numerical simulations of the gmas a function of VGSfor VDS=10 mV (Triode Region).

InFigure 4, gm_=127o, gm_=90oand gm_=53oare respectively the differences between the gmof DSMs with different angles (127o, 90oand 53o) and the respective CSMs counterparts.

Basing onFigure 4, we can see that, DSMs gmwith angles equals to 127o,90oand 53oare respectively 20.5 %, 47.6 % and 122.9 % (around 2.23 times) higher than those verified in the CSM counterparts and these differences also decreaseswhile DSM angle increases, due the LCE and VACLETAPE,considering the same AGand bias conditions.So, this Diamond structures can be used to improve significantly the voltage gain (AV) and unit voltage gain frequency (fT) of the amplifiers used in analog ICs.

Figure 5presents the three-dimensional numerical (3D) simulations results of DSM (=53o) and CSM counterpart Log(IDS) as a function of gate voltage (VGS),for drain voltage (VDS) equal to 10 mV.

Figure 5.

Curves of three-dimensional numerical simulations of DSM (α=53o) and CSM Log(IDS) as a function of VGSfor VDS=10 mV.

AnalyzingFigure 5, we can see that the subthreshold slopes (S) from of the both devices are practically similar around 60 mV/decade. Similar results are found to the others devices studied.

Figure 6presents IDSas a function of VDSwith VGSequal to 0.4 V in the saturation region to DSM with angles equals to 127o(Figure 6.a), 90o(Figure 6.b) and 53o(Figure 6.c) respectively.

AnalyzingFigure 6, it is possible toobserve that always the IDSin saturation region (IDSsat) of all DSMs with different angles are higher than those found in the CSMs counterpart, considering the same AG, geometric factor and bias conditions. Besides that, IDSof DSMs with angles equals to 127o, 90oand 53ois higher 27%, 52% and 284% (1.84 times) than those observed in CSMs counterparts. This can be justified because of LCE and VACLETAPE effects present in the Diamond structure. Additionally, itcan see that DSMs Early voltages (VEA) are strongly degraded in relation to CSMs counterparts, due to the high resultant longitudinal electric field in the drain region of the Diamond structures, as for example, DSM (=53o) VEAis approximately equal to 5V while CSM counterpart VEAis approximately equal to 36V (around seven times smaller) and therefore can limit the use of this device to the output stages of analog ICs.

Figure 6.

Curves of three-dimensional numerical simulations of IDSas a function of VDSfor VGS=0.4 V (saturation region) for DSMs with different α angles [127o(Figure 6.a), 90o(Figure 6.b) and 53o(Figure 6.c)and CSMs counterparts.

The on-state series resistance (RDS_on) can be extracted from theFigure 6in the linear region of the triode region. So, notice that, DSMs with different angles RDS_onare always higher than those observed in CSMs counterparts, as for example, DSM (=53o) RDS_onis around 82% (6 K) smaller than the one found in CSM counterpart (110 K) and therefore, by using DSM counterpart instead CSM, we can improve significantly the processing clock velocity of digital ICs.

Advertisement

4. Electrical and physical analysis

In order to understand the electrical and physical behavior of the Diamond structure in relation the to conventional one counterpart, the longitudinal electrical field vectors picture was extracted of DSM ( =90o) and CSM counterpart structures, regarding VGSand VDSequals to 0.4 V and 1.2 V, respectively, in saturation region, as illustrated inFigure 7.

Figure 7.

Vectors of the resultant longitudinal electric field of the DSM (=90o) (Figure 7.a) and CSM counterpart (Figure 7.b) channel regions for VGSand VDSequal to 0.4 V and 1.2 V, respectively (saturation region).

With the use of TonyPlot3D (Silvaco Inc.) [1] and analyzingFigure 7, it is possible to see how is the behavior of DSM ( =90o) resultant longitudinal electric field (ε//) along of the channel length. Notice that, DSM ( =90o) resultant longitudinal electrical field (ε//) in the channel region edges (1.5x106V/cm) is smaller than the one found in the center of the channel from of the middle of the channel length (2.0x106V/cm), due to the smaller interaction between the longitudinal electric field components (ε//1andε//2) next to DSM edges regions. Besides that, it can beverify that DSM ( =90o) average resultant longitudinal electric field (1.9x106V/cm) is higher than the one observed in CSM counterpart (1.2.104V/cm), due to the presence of LCE and VACLETAPE in the Diamond structure.

Additionally, it is plotted a picture of the DSM ( =53o) total drain current density, considering VGSand VDSequals to 0.4 V and 1.2 V, respectively, in saturation region, as illustrated inFigure 8.

Figure 8.

Curves of the total drain current density of DSM (=53o) (Figure 8.a) and CSM counterpart (Figure 8.b) channel regions for VGSand VDSequal to 0.4 V and 1.2 V, respectively (saturation region).

Figure 8shows that DSM ( =90o) total drain current density is higher in the center of the channel region than in the vertices (source/drain and channel regions interfaces) of the device, indicating that this layout style can be evolved to the octagonal gate geometric [10].

Advertisement

5. SOI MOSFEts investigationat high temperatures

Today, there are much electronic equipments that operate under high-temperature environment (from room temperature up to 300oC) such as avionics, automotive, aircrafts, spacecrafts, ships, train, among other [11-13], where conventional (bulk) CMOS integrated circuits can operate satisfactorily at moderate temperatures (up to 150oC), but when the temperature increases beyond this, the devices present failures arising the threshold voltage reduction and excessive drain leakage current increases [14-18].

Fortunately, thanks to the advantages of SOI CMOS technology, especially at high temperatures, their applications at high temperatures can be extended up to 300oC, where its electrical performance is less impaired than the one found in the conventional MOSFETs [18-20].

In this context, the drain leakage current (IDLeak) of the Diamond SOI nMOSFET (DSM) is analyzed from room temperature up to 300oC by the three-dimensional numerical simulator ATLAS [1]. The DSM IDLeakresults are then compared tothe one found in the SOI MOSFETs (CSM) counterpart, taking into account the same die area, geometric factor and bias and temperature conditions.

Advertisement

6. Three-dimensional numerical simulations results at high temperatures

To analyze DSM operating at high temperatures, from room temperature up to 300oC, CSM and DSM with different angles are implemented by using DevEdit3D [8].

In order to investigate IDLeakbehavior, it is necessary to extract the IDSas a function of VGScurves, at high temperatures. Once IDLeakis obtained in the subthreshold region. In this study, IDLeakis extracted considering VGSequal to -0.5V, as it can be seen inFigure 9, for different temperatures (27oC, 100oC, 200oC and 300oC) and for all devices under evaluation in this work.

Figure 9.

DSM (=90o) Log(IDS) as a function of VGSwith VDSequal to 100 mV, considering different temperatures.

Figure 10shows some results concerning IDLeakas a function of the temperature, considering DSM with different angles and the CSMs counterparts, operating at same bias and temperature conditions.

Figure 10.

DSM and CSM counterpart IDLeakas a function of the temperature.

FromFigure 10, it is possible to observe that, as the temperature increases, IDLeakalso increases, as expected [11]. Besides that, IDLeakis higher in DSM when compared to CSM devices, independently of the angle and the temperature. It happens since the longitudinal electric field is higher in DSM when compared to CSMs and consequently, its drain leakage current (IDLeak) is also larger. On the other hand, it worthwhile to describe that as reduces, IDLeakdifference observed between both devices becomes higher. As for example, when is equal to 127o, DSM IDLeakis 20% higher than the one observed in the CSM counterpart (L= 2.5µm) for VGS=-0.5 V, but providing an IDSalso 20% more than rectangular gate transistor, for VGS=1.2V. These increasing proportions between IDLeakand drain currents are maintained for all equivalent devices, i.e., for = 53oand 127o.

Figure 11shows DSM IDLeakbehavior as a function of the angle α. It is important to comment that this comparison is possible once IDLeakis normalized with the geometric factor (W/L), thus providing a more realistic comparison of IDLeakas α changes. Then, analyzing the results it is possible to observe that, for the same temperature, as α increase, IDLeakbecomes smaller. Otherwise, for the same α, IDLeakalso increases according to the temperature rises.

Figure 11.

Diamond SOI nMOSFET IDLeakbehavior as a function of α, at high temperatures.

Anyway, for both CSM and DSM, IDLeakincreases when the temperature increases in the same proportion,mainly due to the thermal generation process [15], but in DSM IDLeakis higher because the longitudinal electric field along the DSM channel length is higher than the one present in its CSM counterpart.

Advertisement

7. Conclusions

Comparative studies between Diamond SOI MOSFETs and Conventional ones counterparts by numerical simulations demonstrate that Diamond layout style can improve current drive, transconductance, on-state series resistance, frequency response (voltage gain and unit voltage gain frequency) when it is compared to the Conventional counterpart, considering the same gate area, geometry factor (aspect ratio) and bias conditions.

By using Diamond layout style, instead rectangular gate geometry, the die area of analog integrated circuits can be significantly reduced since the drain current of the Diamond structure is higher than the one observed in the Conventional counterpart, considering the same gate area, geometry factor (aspect ratio) and bias conditions.

Therefore, the Diamond layout style becomes an important device alternative to improve the performance of the transistors and consequently the performance of the analog and digital ICs.

This temperature studyin the Diamond layout structure shows the drain leakage current behavior in the Diamond SOI MOSFETs, in comparison to the conventional one counterpart, regarding these devices operating since room temperature (27oC) up to 300oC. From this investigation it is observed that IDLeakin Diamond SOI MOSFETs is higher than the one observed in Conventional devices operating at same bias and temperatures conditions, because the longitudinal electric field along the Diamond channel length is higher than the one found in the conventional counterpart. On the other hand, the DSM drain current is increased in the same proportion of IDLeak, regarding to the same die area, geometric factor and bias and temperature conditions. So, these results show that also it is possible to analyze and investigate the behavior of devices at high temperatures, through three-dimensional numerical simulations, in order to predict its performance in harsh environments before submit them in such real operation conditions, and thus, even get additional information to improve the design/performance of the ICs in harsh temperature conditions. Therefore, Diamond layout style IDLeak(35 pA and 228 pA for equals to 127oand 53o, respectively) presents practically the same magnitude order than the one observed in the CSM counterpart (67 pA and 28 pA for equals to 127oand 53o, respectively) at 300oC and therefore when it is used DSM instead of the conventional one, Diamond structure does not degrade the performance the battery time life of the portable electronics equipment operating in harsh (high temperature) environments.

Advertisement

Acknowledgement

The author would really acknowledge CNPq and FAPESP for the financial support.

References

  1. 1. Atlas User’s Manual2011Device Simulation Software, Silvaco Inc.
  2. 2. Sentaurus Device User Guide (2009) Synopsys Inc., v. C-2009.06.
  3. 3. Colinge J-P2008FinFETs and Other MOSFET Multi-Gate Transistors, Springer.
  4. 4. Gimenez S P, Bellodi M (2010) Diamond SOI MOSFET. PI0802745-5 A2, Publication Date: 03/23/2010, (RPI 2046).
  5. 5. GimenezS. P.BellodiM.(200DiamondS. O. I. M. O. S. F. E. T. A.NewDrain.SourceChannel.InterfaceLayout.toImprove.DrainCurrent.Int. EUROSOI2009, Goteborg, Sweden.18788
  6. 6. Gimenez S P2010Diamond MOSFET: An innovative layout to improve performance of ICs. Int. solid-state electron. j.5416901699
  7. 7. BellodiM.GimenezS. P.2009Drain Leakage Current Evaluation in the Diamond SOI nMOSFET at High Temperatures. Int. ecstransactions.25243253
  8. 8. DeviEdit User’s Manual2011Silvaco Inc.
  9. 9. Sentaurus Structure Editor User Guide2008Synopsys, Inc.
  10. 10. GimenezS. P.(200layoutO. C. T. O.styleP.I0903005A.PublicationDate.10/05/20. R. P.RPI: 2105).
  11. 11. McDonnald J A1996Int. hte rep., III-Vs Review.95759
  12. 12. GrabinskiW.StrickerA.FichtnerW.(199VersatileA.Setupfor.SemiconductorTesting.upto.550ºC”.Int. HITEC 1998. 1: 143.
  13. 13. OsmanM. A.(199Investigationof.hightemperature.effectson. M. O. S. F. E. T.transconductance(gm.Int. HITEC 1998. 1: 301.
  14. 14. EmamM.TinocoJ. C.Vanhoennacker-JanvierD.RaskinJ. P.(200High-temperatureD. C.behaviorsR. F.ofpartially-depleted. S. O. I. M. O. S. F. E. T.transistorsInt. solid-state Electron. j.5219241932
  15. 15. Colinge J-P2004Silicon-on-Insulator Technology: Materials to VLSI. 3rded., Kluwer Acad. Pub., Massachusetts.
  16. 16. FlandreD.FrancisP.GentinneB.ColingeJ. P.1993Demonstration of the potential of accumulation-mode MOS transistors on SOI substrates for high-temperature operation (150-300C). In. ieee trans. electon. dev.141014
  17. 17. FrancisP.TeraoA.GentinneB.FlandreD.Colinge-PJ.1992SO1 Tecnology for high-temperature applications. Int. iedm techn. Dig.1353356
  18. 18. Jeon-SuDeok.BurkD. E.(199temperature-dependentA.modelS. O. I. M. O. S. F. E. T.forhigh-temperature.Application(27o.C-300oC.Int. ieee transact. on elect. dev., j.3821012111
  19. 19. KaamouchE.DelatteP.MoussaM.Si-PRaskin. J.Vanhoenacker-JanvierD.(200TemperatureBehavior.ofSpiral.Inductorson.highresistive.substratein. S. O. I. C. M. O. S.technologyInt. solid-state Electron. J.5219151923
  20. 20. Eggermont-PJ.CeusterD.De FlandreD.GentinneB.JespersG. A.Colinge-PJ.1996Design of SOI CMOS Operational Amplifiers for Applications up to 300oC. Int. ieee journal of solid-state circ. j.,31179186

Written By

Salvador Pinillos Gimenez and Marcello Bellodi

Submitted: 22 November 2011 Published: 19 September 2012