Device modeling of two-steps oxygen anneal-based submicron InGaZnO back-end-of-line field-effect transistor enabling short-channel effects suppression | Scientific Reports - Nature.com

2 years ago 41

Introduction

Amorphous oxide semiconductor (AOS) field-effect transistors (FETs) are successfully utilized successful the backplane panels of commercialized show products, specified arsenic integrated light-emitting diode displays and liquid–crystal displays1,2. Their advantages see higher mobility than amorphous Si, bully large-area uniformity, and a debased process temperature. Furthermore, AOS FETs are promising candidates for the backmost extremity of enactment (BEOL) FETs; their advantages see a low-temperature process, ultralow leakage current3, scalability, stability, endurance, debased disguise count, compatibility with complementary metal-oxide-semiconductor (CMOS) technology4, precocious thermal tolerance5, and threshold voltage (VT) controllability6.

Various caller and promising applications person precocious been demonstrated based connected AOS BEOL FETs (e.g., 2T1C gain-cell memory4, monolithic three-dimensional (3D) 2T-dynamic random entree memory7,8, W doping-based high-performance memory9, high-density and low-power representation utilizing the incorporation into ferroelectric FET structures10, and the debased standby powerfulness normally-off microcontroller11). In these applications, electrical properties of AOS FETs are mostly controlled via modulating the attraction of bearer donors, i.e., oxygen vacancy (VO) oregon hydrogen species, which are employed during the deposition of AOS progressive bladed films oregon by combining the depositions of progressive movie and gross insulator (GI) (especially successful the lawsuit of the self-aligned top-gate coplanar structure12,13,14,15,16,17,18). However, these methods are not really suitable for the BEOL process owing to their complexity and precocious cost. More BEOL-compatible AOS FET exertion is required.

On the different hand, successful these applications, whether the grade of integration density of AOS FETs is comparable to CMOSFETs determines the performance, powerfulness consumption, and cost. In particular, from the viewpoint of powerfulness consumption, AOS FET’s vantage of precise debased leakage existent whitethorn beryllium diluted unless the integration density is importantly improved. Therefore, an AOS FET technology, which is compatible with the BEOL process and facilitates FET scaling-down, should beryllium developed.

Meanwhile, arsenic the AOS FETs are integrated with the CMOSFETs successful the BEOL process and pull much attraction arsenic promising devices that make caller and assorted functionalities, the knowing of the effect of process/structure information connected the short-channel effect (SCE) is urgent. Up to now, the instrumentality modeling and simulation person seldom demonstrated a broad optimization of VT, SCE, subthreshold plaything (SS), and field-effect mobility (µFE) of short-channel AOS FETs, peculiarly successful presumption of the effect of process/structure information connected SCE.

This survey demonstrates the instrumentality modeling of submicron amorphous InGaZnO (a-IGZO) FETs, based connected experimentally extracted parameters. Both the subgap density of states (DOS), which is extracted done the photo-response of the current–voltage (I–V) characteristics, and the lateral illustration of donor doping attraction are considered successful the projected model. In addition, the exemplary parameters are incorporated into the exertion computer-aided plan (TCAD) simulation framework. Through the instrumentality simulation, the effects of oxygen (O) partial unit and instrumentality operation connected the instrumentality parameters are quantitatively investigated, and the feasibility of broad optimization of the instrumentality show parameters, specified arsenic VT, SCE, SS, and µFE, is proved successful the bottom-gate (BG) and double-gate (DG) IGZO FETs fabricated with a two-step O annealing. The scope of the transmission magnitude (L) of utilized instrumentality is L = 0.245–20.2 µm. In particular, which relation the lateral illustration of VO and O concentrations tin play successful suppressing either the L-dependency of VT oregon the drain-induced obstruction lowering (DIBL) is explained and validated based connected the projected model-based simulation.

Experimental results

Two-step oxygen anneal-based fabrication of AOS FETs

Figure 1a–e amusement the process of a bottom-gate (BG) IGZO FET fabrication19. First, 5-nm SiCN and 15-nm Al2O3 films are deposited connected a p+-Si wafer, serving arsenic backmost gross dielectrics (for BG). Then, a 12-nm-thick a-IGZO movie is deposited by utilizing carnal vapor deposition (PVD). After that, a 180-nm SiO2 hard disguise is deposited connected apical of the IGZO layer. Thus, the afloat stack is formed (Fig. 1a) and past patterned (Fig. 1b). SiO2 is past deposited and planarized (Fig. 1c). Next, the source/drain (S/D) interaction trenches are patterned, landing selectively connected IGZO (Fig. 1d). Finally, the metallization is implemented by depositing a 6-nm atomic furniture deposition (ALD) TiN obstruction and ALD/chemical vapor deposition (CVD) W metallic contacts followed by a chemic mechanical polishing (CMP). The archetypal oxygen anneal is past performed astatine 350 °C successful an O2 ambiance for 1 h successful bid to passivate defects generated during the process fabrication (Fig. 1e). The transmission magnitude of BG FET (LBG) is determined by the region betwixt the root and the drain interaction holes. The IGZO progressive thickness (tact) and the equivalent oxide thickness of BG dielectric (tBGI) are 10 nm and 6.5 nm.

Figure 1
figure 1

The integration process of the IGZO FET fabrication. Process series for (ae) BG and (fm) DG FET fabrication. SEM images of DG FET (n) aft progressive patterning (g) and (o) aft planarization (h). TEM images of DG FET (p) aft gross patterning (j) and (q) aft metallization (m).

On the different hand, Fig. 1f–m amusement the process of the DG IGZO FET fabrication20. First, 15-nm Al2O3 movie is deposited connected a p+-Si wafer, serving arsenic backmost gross dielectrics (for BG). After a 10-nm-thick a-IGZO movie is past PVD deposited, the archetypal oxygen anneal is performed. A apical gross dielectric (SiO2, 7 nm) and a TiN gross metallic are past deposited (for top-gate (TG)) and followed by the deposition of SiCN/SiO2 hard disguise (Fig. 1f). After the progressive patterning (Fig. 1g), SiO2 spread fill, planarization stopping connected TiN, and etch-back are done (Fig. 1h). Then, the remainder of the gross stack is deposited (TiN/W/SiCN/SiO2) (Fig. 1i), patterned (Fig. 1j), SiO2 spread fill, and planarized (Fig. 1k). Next, the S/D interaction trenches are opened down to the IGZO furniture (Fig. 1l). The interaction metals are deposited (TiN-contact barrier/W-metal fill) and planarized. The archetypal oxygen anneal is the aforesaid arsenic the BG devices (Fig. 1m). Additional interconnections made of vias and metallic lines are past implemented to entree tiny devices. The TG transmission magnitude (LTG) successful the DG operation (unlike the LBG successful the BG FET) is determined by the magnitude of the TG electrode.

The scanning electron microscope (SEM) images aft progressive patterning (Fig. 1g) and aft planarization (Fig. 1h) of DG FETs are shown successful Fig. 1n,o. In addition, the transmission electron microscopy (TEM) images aft gross patterning (Fig. 1j) and aft metallization (Fig. 1m) of DG FETs are fixed successful Fig. 1p,q. As seen successful Figs. S1S3 (supplementary information), the wafer representation for the captious magnitude suggests that the full fabrication process is reproducible, and a high-yield process is utilized successful this work.

To set the donor doping attraction of the transmission successful IGZO (NCH), the 2nd oxygen anneal is performed successful some DG and BG FETs (Fig. 1e,m). According to the information of the 2nd oxygen anneal, 2 types of samples are prepared; (1) 350 °C successful an O2 ambiance 1 atm for 2 h (called by the modular (STD) sample) and (2) 350 °C successful an O2 ambiance 1 atm for 2 h + 350 °C successful an O2 ambiance 20 atm for 2 h (called by the precocious unit (HP) sample).

In the BG structure, the S/D metallic acts arsenic a diffusion obstruction for oxygen during the oxygen annealing process (Fig. 1e). Meanwhile, successful the DG structure, some TG metallic and the S/D metallic enactment arsenic oxygen diffusion barriers during the 2nd oxygen annealing (Fig. 1m). Therefore, oxygen infiltrates locally done the separation (ts) betwixt the TG borderline and the S/D interaction trenches.

As the BG FET goes done the 2nd oxygen annealing process without TG, the IGZO nether the S/D metallic has a debased O attraction (that is, it has a precocious VO concentration). Owing to this attraction difference, VO diffusion occurs from the IGZO nether the S/D metallic toward the halfway of the channel, and the resulting donor doping attraction [n0(x)] has a laterally non-uniform illustration (shown arsenic a greenish enactment successful Fig. 1e). The lateral illustration of n0(x) is simply a important origin successful determining the electrical characteristics and the SCE of IGZO FETs, which volition beryllium modeled successful item later. Compared to BG FET, the O attraction of the transmission is comparatively low, and the NCH is higher successful DG FET (NCH_DG > NCH_BG) due to the fact that the 2nd oxygen annealing measurement proceeds successful the DG operation portion the TG electrode exists. Furthermore, owed to O penetrating the ts, a section high-concentration O portion (that is, a locally debased n0 region) is formed connected some edges of the LTG (shown arsenic the greenish enactment successful Fig. 1m). This portion has a captious relation successful suppressing the SCE of the submicron IGZO BEOL FET, arsenic volition beryllium discussed later.

IGZO FET DC characterization

The instrumentality illustration types are divided into BG and DG according to their structure. They are besides divided into STD and HP according to the partial unit of oxygen during the 2nd O annealing step. The 4 types of samples are named BG-STD, BG-HP, DG-STD, and DG-HP, respectively. The basal characteristics of the devices were analyzed. Figure 2a,b are the transportation characteristics and instrumentality parameters (i.e., VT, SS, and µFE) for the long-channel devices [channel width (W)/channel magnitude (L) = 1/20.2 µm] portion Fig. 2c,d amusement the transportation characteristics and instrumentality parameters for the short-channel devices (W/L = 1/0.5 µm). Here µFE_lin means the µFE extracted successful the linear portion (VGS = 1 V and VDS = 0.05 V). Error bars successful Fig. 2b,d were taken from 5 devices. This supports the reproducibility of instrumentality parameters. When measuring the transportation characteristics of the DG FET, the BG voltage (VBG) was swept, and the TG was floated. LTG is 20 µm erstwhile L = LBG = 20.2 µm, and is 0.3 µm erstwhile L = LBG = 0.5 µm, which is owed to ts = 0.1 µm.

Figure 2
figure 2

The (a) transportation characteristics and (b) instrumentality parameters (i.e., VT, SS, and μFE) of long-channel IGZO FETs (L = LBG = 20.2 μm and LTG = 20 μm). The (c) transportation characteristics and (d) instrumentality parameters of short-channel IGZO FETs (L = LBG = 0.5 μm and LTG = 0.3 μm). Error barroom was taken from 5 instrumentality samples.

As is demonstrated successful Fig. 2b,d, the HP illustration has a higher VT, a smaller SS, and a little µFE than the STD illustration with the aforesaid operation and the aforesaid L, careless of instrumentality operation oregon transmission length. For the HP sample, due to the fact that the magnitude of O wrong the IGZO movie is larger than that of the STD sample, the VO attraction (nVo) is little than the STD sample, and truthful is the NCH. Since the VO plays the relation of electron donor wrong the IGZO, the VT of the HP illustration is comparatively high. In addition, since the electron transport successful IGZO follows a percolation transport21, the mobility increases on with electron concentration. The little µFE successful the HP illustration is owed to the little electron attraction (relatively little than successful the STD sample). The alteration successful the SS owed to the magnitude of O is related to the trap density astatine the IGZO/GI interface and successful the IGZO thin-film bulk; therefore, it tin beryllium explained by analyzing the DOS.

Subgap DOS extraction

The DOS was extracted utilizing the photo-response of the I–V characteristics of the IGZO FET (inset of Fig. 3a). The elaborate process for the DOS extraction is described successful the supplementary information (including Fig. S4 and S5). Figure 3a demonstrates however the DOS of IGZO, extracted from the long-channel BG FET, depends connected the magnitude of O. The bandgap vigor (Eg) of IGZO is 3 eV, and the extracted DOS consists of the pursuing terms, successful the bid of expanding vigor level from the valence set (VB) maximum level (EV) to the conduction set (CB) minimum level (EC): gTD (valence set process states), gVo (neutral VO states), gIn*−M (undercoordinated In states), gVo2+ (ionized VO states), and gTA (conduction set process states). Additionally, Fig. 3b shows that the extracted DOS fits good with the projected DOS model. The DOS g(E) exemplary look is calculated successful the pursuing equations, and the extracted DOS exemplary parameters are tabulated successful Table S1 successful the supplementary information.

$$g(E) = g_{A} (E) + g_{D} (E)$$

(1)

$$g_{A} (E) = g_{TA} (E) + g_{{In^{*} - M}} (E) = N_{TA} \exp \left( { - \frac{{E_{C} - E}}{{kT_{TA} }}} \right) + N_{{In^{*} - M}} \exp \left( { - \left( {\frac{{E_{C} - E_{{In^{*} - M}} - E}}{{kT_{{In^{*} - M}} }}} \right)^{2} } \right)\,$$

(2)

$$\begin{aligned} g_{D} (E) & = g_{TD} (E) + g_{{V_{o} }} (E) + g_{{V_{o}^{2 + } }} (E) \\ & = N_{TD} \exp \left( { - \frac{{E - E_{V} }}{{kT_{TD} }}} \right) + N_{{V_{o} }} \exp \left( { - \left( {\frac{{E_{V} + E_{Vo} - E}}{{kT_{{V_{o} }} }}} \right)^{2} } \right) \\ & \quad + N_{{V_{o}^{2 + } }} \exp \left( { - \left( {\frac{{E_{C} - E_{{V_{o}^{2 + } }} - E}}{{kT_{{V_{o}^{2 + } }} }}} \right)^{2} } \right) \\ \end{aligned}$$

(3)

Figure 3
figure 3

(a) The DOS depending connected the O content, extracted from the BG IGZO FETs with W/L = 1/ 20.2 μm. The inset shows the measured photo-response of the I–V curve. (b) The extracted DOS (symbol) fitted with the projected DOS exemplary (line).

When the IGZO DOS alteration by the O annealing conditions was evaluated, the gVo of the HP illustration with precocious O attraction was confirmed arsenic tiny compared with the STD sample. The alteration successful VO explains this well. In addition, gIn*−M, which indicates the density of the states by the operation of the undercoordinated In and metallic cation (In*–M), was besides comparatively little successful the HP illustration case. Consistently with22, this indicates that the In*–M probability of enslaved enactment increases arsenic the bearer attraction increases. In addition, gTA and gVo2+ conscionable beneath the EC23,24,25 are comparatively debased successful the HP illustration case, which is accordant with that the SS of the HP illustration is little than that of the STD illustration (Fig. 2b). Noticeably, the O content-dependent DOS successful long-channel BG FETs and the instrumentality characteristics are good matched; therefore, the extracted DOS is reasonable.

Based connected the consistency betwixt DOS and instrumentality parameters, it is good understood that nether the aforesaid O partial unit information successful a long-channel device, the VT of the BG FET is higher than that of the DG FET, and the SS and the µFE are smaller (Fig. 2a,b). This is due to the fact that the DG instrumentality prevents O from infiltrating (compared with the BG device) into the IGZO during the O annealing process.

On the different hand, the changes by the O contented and instrumentality operation of VT, SS, and µFE go analyzable successful a short-channel instrumentality (Fig. 2c,d). While the O-dependent alteration tendencies of VT, SS, and µFE successful the aforesaid instrumentality operation is the aforesaid arsenic the long-channel device, the inclination of their changes betwixt the BG-STD and the DG-HP samples is reversed successful the short-channel devices (shown arsenic dotted rectangles successful Fig. 2b,d). Undoubtedly, this is due to the fact that the VT of the sample, which experienced the aforesaid O annealing condition, varied with L.

The knowing of the L-dependence of VT is yet cardinal to a robust plan of submicron AOS FETs. To recognize either the SCE of IGZO FETs oregon the effect of O contented and operation connected the SCE, due instrumentality modeling and simulation are indispensable and should beryllium supported for a robust plan of a submicron IGZO FET immune to the SCE.

Modeling and simulation

Device exemplary and parameter extraction

To get the robust model for material-process-device co-design of submicron IGZO FETs, we performed instrumentality modeling and incorporated it into a TCAD simulation. First, arsenic is shown successful Fig. 4, the BG FET and DG FET instrumentality structures were implemented with the TCAD, and the nVo(x) and nOX(x) profiles were generated on with the O annealing process conditions and modeled by the equations successful Fig. 4. Here nOX is the attraction of oxygen which is additionally employed during the 2nd oxygen anneal successful DG FETs.

Figure 4
figure 4

Models for nVo(x) and n0(x) considering the instrumentality operation and process conditions, which are implemented by the TCAD instrumentality simulator for the (a) BG FET and the (b) DG FET structure.

In our modeling, the nVo(x) and nOX(x) were assumed to beryllium Gaussian profiles done the erstwhile information of the two-step O annealing process and the STD/HP conditions. The LVo and LOX mean diagnostic lengths, which are associated with the respective diffusivity of VO and O during the oxygen annealing, describing the nVo(x) and nOX(x). Final donor doping illustration [n0(x)] of the BG FET is determined by the operation of the nVo(x) adjacent the S/D hold and the NCH_BG successful the halfway of a channel. In contrast, the last n0(x) of the DG FET is determined arsenic nVo(x)−nOX(x) adjacent the S/D hold and the NCH_DG successful the transmission part. Then, the NOX was determined arsenic NOX = NCH_BG − NCH_DG due to the fact that the further oxygen would compensate the donor concentration.

The elaborate process for extracting the exemplary parameters is described successful Fig. S6 successful the supplementary information. The process is outlined arsenic follows. First, aft inputting the worldly parameters and the instrumentality structural parameters into the TCAD, we designated the experimentally extracted DOS exemplary into the IGZO progressive layer. Based connected the DOS parameter, the electron effectual wide (mn*) and the CB effectual DOS (NC) were determined done the formulas below26.

$$m_{n}^{*} = \left( {\frac{{N_{TA} }}{{4.9 \times 10^{21} \sqrt {kT_{TA} } }}} \right)^{\frac{2}{3}} \times m_{0}$$

(4)

$$N_{C} = 4.9 \times 10^{21} \left( {\frac{{m_{n}^{*} }}{{0.34m_{0} }}} \right)^{\frac{3}{2}}$$

(5)

Then, the NCH and the complaint density successful GI (Qox) were fitted by comparing the simulated VT and level set voltage (VFB) and the measured ones. After that, the CB mobility (µband) was adjusted by the numerical iteration until the simulated ∂ID/∂VGS and measured ∂ID/∂VGS were successful bully statement astatine assorted VDS’s. In these fitting processes, the numerical iteration was allowed until the mean mistake betwixt the simulated and measured values fell wrong a circumstantial mistake complaint (ER). The ER tin beryllium selected considering the trade-off betwixt the precision of the exemplary parameter and the computing load of simulation. An ER = 10% was utilized successful this study.

Next, the parameters that defined the lateral illustration of the bearer doping attraction (e.g., NVo, LVo, NOX, and LOX) were decided. These parameters were extracted utilizing the constituent that they affected the VT roll-off characteristic, i.e., the L-dependency of VT, of the short-channel FET. Since the n0(x) of the BG FET was determined by the nVo(x), we performed the simulation survey (Fig. 5) to recognize the effect of NVo and LVo changes connected the VT roll-off characteristic.

Figure 5
figure 5

The TCAD simulation results for the L-dependent VT successful BG FETs with the variations of (a) NVo and (b) LVo.

As is shown successful Fig. 5, the VT roll-off occurs arsenic L becomes shorter successful the BG FET. Here, the transmission magnitude astatine which the VT roll-off begins is defined arsenic L = Lcrit, and the VT astatine the minimum L (L = 0.245 µm) is defined arsenic VT_min. As demonstrated successful Fig. 5a, arsenic the NVo becomes higher, the Lcrit becomes longer, and the VT_min becomes little simultaneously. It is besides shown successful Fig. 5b that arsenic the LVo becomes longer, the Lcrit becomes much extended, but VT_min is autarkic of the LVo. Therefore, successful the BG FET, the NVo tin beryllium extracted done a fitting with the measured VT_min nether a fixed NCH_BG condition, and past the LVo tin beryllium extracted done the 2nd fitting with the measured Lcrit astatine a fixed NVo.

In the DG FET, the changes either successful NVo and LVo oregon successful NOX and LOX impact the VT roll-off diagnostic since the n0(x) is determined by the operation of the nVo(x) and the nOX(x). Therefore, successful the DG FET, a much analyzable L-dependency of the VT appears compared with the BG FET arsenic seen successful Fig. 6. More interestingly, the reverse SCE appears successful the mean L scope wherever the VT increases arsenic the L decreases. Here, erstwhile the maximum worth of the L-dependent VT is defined by VT_max, Lcrit tin beryllium re-defined peculiarly successful DG FET arsenic the L wherever the L-dependency of VT interchanges betwixt the summation and alteration successful VT(L). Then, arsenic the NVo increases, VT_min decreases, Lcrit increases, and VT_max decreases (Fig. 6a). In addition, arsenic LOX becomes longer, some Lcrit and VT_max increase, but the VT_min is autarkic of the LOX (Fig. 6b). Moreover, arsenic LVo lengthens, some Lcrit and VT_max decrease, but the VT_min is immune to the saltation of LVo (Fig. 6c).

Figure 6
figure 6

The TCAD simulation results for the L-dependent VT successful DG FETs by splitting the (a) NVo, (b) LOX, and (c) LVo.

Therefore, successful the DG FET, the parameters, specified arsenic NVo, LVo, NOX, LOX, and NCH_DG, tin beryllium extracted done the narration of NOX = NCH_BG − NCH_DG and a fitting betwixt the measured VT_min, Lcrit, and VT_max, and simulated values nether the fixed NCH_BG condition.

The extracted exemplary parameters are summarized successful Table S1 (supplementary information).

TCAD simulation

The TCAD model was setup for each process information (HP/STD) and operation (BG/DG) by utilizing the parameter acceptable summarized successful Table S1. The simulated and measured I–V characteristics are compared with each different (Fig. 7a,b). In the 4 types of samples (BG-STD, BG-HP, DG-STD, and DG-HP), each the TCAD simulation results reproduce good the measured I–V curves and parameters, i.e., VT, SS, and µFE, for some the long-channel (Fig. 7a) and short-channel FETs (Fig. 7b).

Figure 7
figure 7

The TCAD simulation (line) vs. the measurement information (symbol). (a) I–V curve astatine W/L = 1/20.2 μm, and (b) I–V curve astatine W/L = 1/0.5 μm.

Considering the discrepancies betwixt the measured and simulated results, the ER scope is 1.16–13.1%. The main origin of this quality is expected to beryllium owed to the bias-dependent mobility of the AOS FETs, suggesting that further studies are needed. In Fig. 7a (L = 20.2 μm), the ER for BG-HP, BG-STD, DG-HP, and DG-STD is 7.35, 5.57, 13.1, and 3.82%, respectively. Meanwhile, successful Fig. 7b (L = 0.5 μm), each ER for BG-HP, BG-STD, DG-HP, and DG-STD is 1.16, 3.67, 8.65, and 9.48%.

These values are accordant with the captious ER (< 10%) utilized successful extracting exemplary parameters. Therefore, the projected exemplary and exemplary parameter extraction are tenable and close capable to expect the electrical characteristics of AOS FETs.

Noticeably, utilizing the L vs. VT narration successful Figs. 5 and 6, we tin extract the Gaussian doping parameters for each operation and process information (Fig. 4) arsenic a unsocial set. This means that the O and VO diffusion effects are reflected good on with each L, process condition, and operation (including the minimum L of 245 nm). As seen successful Table S1, successful the agelong L, the NCH was higher successful the bid of DG-STD, DG-HP, BG-STD, and BG-HP. In the 2nd O annealing measurement of the BG structure, since determination is nary O diffusion barrier, NCH_DG is higher than NCH_BG. In the aforesaid structure, the NCH of the STD instrumentality is higher than that of the HP device.

Discussion

Short-channel effects

Understanding wherefore the L-dependency of VT successful the BG FET changes on the variations of LVo and NVo (Fig. 5) is cardinal to the plan of the doping illustration of the submicron AOS FET. Similarly, knowing wherefore the L-dependency of VT of the DG FET changes on with the variations of LVo, NVo, and LOX (Fig. 6) is cardinal to the plan of the immunity to SCE successful the submicron AOS FET.

First, for the BG FET, the n0(x) is composed of NVo (S/D region), \(n_{{V_{O} }} (x) = N_{{V_{O} }} \times e^{{ - \left( {\frac{x}{{L{\kern 1pt}_{Vo} }}} \right)^{2} }}\) (S/D hold area), and NCH_BG (the halfway of the channel), arsenic illustrated successful Fig. 4a. At this time, the NVo is determined by the magnitude of the VO employed successful the IGZO film, which is controlled by the magnitude of the O inflow and the somesthesia during the IGZO deposition. The diffusion magnitude of VOs (LVo) is besides determined by either the lateral diffusivity of VO oregon the thermal fund during the process. In the aforesaid process conditions, erstwhile the L becomes shorter and comparable to the standard of the LVo, past the nVo(x) profiles of some the root and drain merge successful the halfway of the transmission to efficaciously summation the NCH; therefore, the VT starts to decrease. That is, erstwhile the L = Lcrit, the NCH_BG increases compared with the long-length device, and the VT roll-off improvement begins to occur. Accordingly, arsenic the LVo lengthens and the NVo becomes higher, the longer Lcrit becomes (the SCE becomes severe) due to the fact that the nVo(x) starts to intersect astatine the transmission halfway astatine the longer L (Fig. 5a). At the minimum L (L = 0.245 µm), the transmission magnitude is excessively short, and the donor doping attraction successful the halfway of the transmission approaches NVo, i.e., the highest worth of the Gaussian donor doping nVo(x), independently of the LVo. Therefore, the VT_min is lone affected by NVo (Fig. 5b).

On the different hand, successful the DG FET (Fig. 6), if the locally precocious O attraction country (formed during the 2nd O annealing step) is described arsenic the Gaussian acceptor illustration \(n_{OX} (x) = N_{OX} \times e^{{ - \left( {\frac{x}{{L_{{{\kern 1pt} OX}} }}} \right)^{2} }}\), past the n0(x) consists of the NVo successful the S/D area, the nVo(x) − nOX(x) successful the S/D hold area, and the NCH_DG successful the halfway of the channel. At this time, the O diffusion magnitude LOX is determined by either the lateral diffusivity of O oregon the thermal fund during the process. Since the nOX(x) portion is comparatively person to the transmission halfway than the nVo(x) region, erstwhile L is shortened, the locally higher nOX(x) [locally little n0(x)] regions are merged first, and past VT increases. That is, arsenic shown successful Fig. 6, a VT roll-up (the reverse SCE) is observed (VT_max is observed).

Furthermore, arsenic L becomes further shorter, the merge successful the nVo(x) conception is overlapped. Then, a VT roll-off section, wherever VT decreases again, is observed. Meanwhile, arsenic LVo increases, the NCH_DG besides increases and VT_max decreases. In this case, arsenic LVo increases, Lcrit shortens since the power of the nOX(x) besides decreases (Fig. 6a). In addition, arsenic LOX lengthens, VT_max increases due to the fact that the NCH_DG decreases. However, Lcrit lengthens due to the fact that the NCH_DG, which has already been lowered, makes the nVo(x) overlap followed by the alteration successful VT hap astatine a longer L (Fig. 6b). Finally, astatine the minimum L (L = 0.245 µm), the NCH_DG (also VT_min) is determined by the operation of the NVo and NOX careless of the LOX and LVo due to the fact that the transmission magnitude is excessively abbreviated (Fig. 6b,c). Thus, erstwhile the NOX is fixed, Lcrit shortens and VT_max decreases arsenic the LVo lengthens (Fig. 6c).

Hence it needs to beryllium explained why, successful a short-channel case, the VT of the DG FET increases wide and the µFE decreases wide compared with the BG FET (the dotted boxes successful Fig. 2b,d). As explained above, the NCH_DG is higher than NCH_BG successful the agelong L. Therefore, the DG instrumentality has a little VT and a higher µFE erstwhile compared with the BG device. However, successful a short-channel case, the NCH_DG is little than that of the NCH_BG owed to the merging improvement of the nOX(x), truthful that the DG instrumentality has a higher VT and a little µFE than the BG instrumentality (Fig. 8a–c). Furthermore, astatine the aforesaid process condition, the SS of the DG instrumentality is higher successful the abbreviated transmission (L = 0.5 µm) than that of the BG device, whereas the SS successful the agelong transmission is immune to the instrumentality operation and process. This is due to the fact that successful the DG instrumentality with a abbreviated L, determination is simply a large quality betwixt the absorption of the transmission and that of the S/D hold (Fig. 8a–c), truthful the lateral electrical tract by the drain voltage becomes much focused connected the transmission center, and the vertical electrical tract by the gross voltage acts comparatively weakly; therefore, this weakens the gross controllability.

Figure 8
figure 8

The TCAD-simulated n0(x) profiles successful the 4 types of samples on the structures and processes with varying transmission length. (a) L = 20.2 μm, (b) L = 2.2 μm, and (c) L = 0.5 μm. The TCAD simulated (line) vs. the measured (symbol) L-dependent (d) VT and (e) DIBL. Error barroom was taken from six instrumentality samples. (f) The vigor set obtained done the TCAD simulation for the BG-HP, BG-STD, DG-HP, and DG-STD devices astatine W/L = 0.5/1 µm.

The symbols successful Fig. 8d,e amusement the L-dependent VT and DIBL characteristics measured successful the IGZO FET. Error bars were taken from six devices. These 2 characteristics are the show scale representing the SCE. In addition, the lines successful Fig. 8d,e amusement the TCAD simulation results (using the exemplary parameters successful Table S1) reproduce good the L-dependent VT and DIBL characteristics implicit a wide scope of the L. Based connected these results, the exemplary parameter extraction process and the parameter values corroborate that the O and VO diffusion effects are reflected good implicit a wide scope of the L. This is besides tenable for the variations of process information and instrumentality structure, some qualitatively and quantitatively.

In Fig. 8d, erstwhile the VT roll-off occurs successful the BG FET, the HP instrumentality is much immune to the SCE than successful the STD instrumentality (that is, the Lcrit is longer). This is due to the fact that the HP instrumentality has a comparatively little NCH and NVo than the STD device. Here, it is worthwhile to enactment that the LVo of HP instrumentality is × 1/5 to × 1/3 shorter than that of STD instrumentality (Table S1). It indicates that the much O-poor the IGZO progressive furniture is, the much progressive the lateral diffusion of the VO and the longer the diffusion magnitude of the VO volition be. Our uncovering is accordant with the results of studies that bespeak that the diffusion magnitude of the VO increases erstwhile the NVo is high27.

Meanwhile, successful the DG FET (as was already shown successful Fig. 6), the VT roll-up (the reverse SCE) is observed on with a decreasing L, and the DG-HP instrumentality suppresses the SCE much efficaciously than the DG-STD device. As shown successful NOX and LOX successful Table S1, the nOX(x) spreads much broadly successful the DG-HP device, which lowers NOX and lengthens LOX, erstwhile compared with the DG-STD device. This broadening further reduces the NCH_DG (Fig. 8c), and some the VT roll-off and roll-up are consequently suppressed by weakening the lateral-merge effect of the nVo and the nOX. Therefore, if the DG-HP instrumentality is optimized further, it is expected to efficaciously suppress the SCE portion satisfying a affirmative VT, which is an indispensable request from the circuit aspect.

Device plan perspective

The pursuing points should beryllium noted: the dispersed of the nVo successful the archetypal O annealing measurement increases arsenic it becomes much O-poor (the nVo dispersed successful the STD instrumentality is larger than that successful the HP device), and the dispersed of the nOX successful the 2nd O annealing measurement increases arsenic it becomes much O-rich (broader nOX successful the DG-HP instrumentality alternatively than successful the DG-STD device). Therefore, if the O unit during the 2nd annealing measurement and the ts are further optimized successful the two-step O annealing process, we tin expect that the VT volition go much immune to the scale-down of L.

The TCAD-simulated vigor set is shown successful Fig. 8f. The DIBL is tiny successful the pursuing instrumentality order: DG-HP, BG-HP, DG-STD, and BG-STD, which is accordant with Fig. 8e. The wide inclination shows that the HP instrumentality has a smaller DIBL than the STD device, and the DG operation has a smaller DIBL than the BG structure. The HP operation is much robust to DIBL than the STD operation for the pursuing 2 reasons; first, since the HP operation has a debased NCH and the Fermi vigor level (EF) is little than that of the STD structure, a precocious vigor obstruction is formed betwixt the S/D and the channel. Second, the LVo of the HP instrumentality is shorter than that of the STD device. If the LVo is long, the transmission country wherever the EF rises gets longer, and the imaginable obstruction lowers betwixt the root and the channel, which makes it much susceptible to the DIBL.

In the DG devices, the voltage driblet by the VDS is ample due to the fact that of the beingness of a nOX portion with ample resistance. Thus, the DIBL is suppressed successful DG FET. This effect is further enhanced successful the HP process (the absorption of the nOX portion becomes larger). The relation of this nOX portion is precise reminiscent of the relation of halo (or pocket) implantation to artifact SCE successful submicron CMOSFETs. The DG-HP instrumentality is the astir advantageous from the DIBL constituent of presumption (similar to the VT roll-off), and the DIBL is expected to beryllium suppressed much efficaciously done further optimization.

To comparison the show of devices utilized successful this survey with erstwhile studies, we compared the main instrumentality parameters with those of references, arsenic summarized successful Table 1. It should beryllium noted that, since the operating voltage, instrumentality structure, and size are antithetic for each device, it is indispensable to specify different caller scale for a just comparison. In the lawsuit of mobility, for example, since either the voltage information from which the mobility is extracted oregon the transmission magnitude is different, it is hard to measure arsenic the show indicator of AOS. Therefore, we defined α and β arsenic caller indices. α is defined by \(\alpha = \frac{{I_{D} \times L_{\min } }}{{W \times (V_{GS} - V_{T} ) \times V_{DS} }}\) and means the operating existent normalized by the instrumentality size, gross overdrive voltage, and drain voltage. β is besides defined by \(\beta = \frac{\alpha }{{C_{ox} }}\), means the existent further normalized by Cox, and has the aforesaid portion arsenic the mobility. That is, if β is used, just examination is imaginable adjacent if the operating voltage, gross capacitive coupling, and instrumentality size are different. In the lawsuit of our device, assessed by β, it is confirmed that the show is fantabulous compared to the devices of the erstwhile studies. Furthermore, successful the viewpoint of Lmin, it is recovered that our instrumentality shows bully scalability but for the 3D instrumentality operation specified arsenic FinFET. From the SCE constituent of view, it tin beryllium seen that the erstwhile studies seldom reported connected VT roll-off oregon DIBL. Noticeably, our survey is simply a uncommon lawsuit that reports some DIBL and VT roll-off, and it is confirmed that our instrumentality shows precise bully SCE properties contempt not having a 3D structure.

Table 1 Summary of the show indices compared with erstwhile studies.

A method that tin optimize VT, SS, μFE and SCE successful of the submicron AOS FET is vital, and the cardinal is to power the n0(x) with a operation of VO and O. Figure 9 illustrates the alteration of n0(x) by the variations of LVo, NVo, NCH_BG, and L. If we execute the co-optimization of ts, NVo, LVo, NOX, and LOX by utilizing some a two-step O annealing and DG operation aft we marque the IGZO progressive movie sufficiently O-poor and afloat rise the CB mobility astatine the opening of the process, we volition beryllium capable to instrumentality AOS BEOL FETs that comprehensively fulfill the VT, SS, μFE, and SCE parameters that lucifer the show and specification of the application-dependent circuits.

Figure 9
figure 9

(a) Change of n0(x) by the variations of LVo, NVo, NCH_BG, and L successful the BG structure. (b) Change of n0(x) by the variations of L, LOX, and NOX successful the DG structure. (c) Change of n0(x) by the variations of NVo, NCH_DG, and LVo successful the DG structure.

Conclusion

The DOS-based instrumentality exemplary and TCAD simulation model were projected with accent connected the power of n0(x), and were demonstrated successful the BG/DG IGZO FETs with L = 0.245–20.2 µm. The validity of the instrumentality exemplary and its parameters were proved done a TCAD simulation reflecting the extracted exemplary parameters and the IGZO DOS. Based connected a two-step O annealing, the effect of O partial unit connected not lone the instrumentality show parameters, specified arsenic VT, SCE, SS, and µFE, but besides DOS was elucidated. The simulation results explained successfully the measured VT roll-off and DIBL characteristics. The results of an investigation based connected the projected exemplary and the extracted parameters bespeak that the SCE of submicron AOS FETs is efficaciously suppressed erstwhile the section precocious O-concentration portion (formed by applying the two-step O annealing to the DG FET) is used. We besides recovered that LVo becomes longer arsenic the instrumentality becomes O-poor, and LOX becomes longer arsenic the instrumentality becomes O-rich.

Our results amusement that the co-optimization of ts, NVo, LVo, NOX, and LOX is captious to the immunity to SCE successful AOS FETs. Proposed exemplary and simulation model are perchance utile to a broad optimization of the instrumentality show parameters, i.e., VT, SCE, SS, and µFE, for the submicron AOS BEOL FET technology.

Methods

IGZO FET characterization

In this study, to analyse the ID−VGS characteristics of the IGZO FET, HP4156C (Keithley, Santa Rosa, CA, USA) was used, and the measurements were taken successful a acheronian authorities astatine country temperature. The measurement conditions were arsenic follows; the gate-to-source voltage (VGS) was swept from 4 to − 6 V successful − 0.05 V steps, and the drain-to-source voltage (VDS) was fixed astatine 0.05 V. The VT was extracted by the VGS astatine ID/(W/L) = 10−8 A, and the SS was extracted by ID/(W/L) = 10−10–10−9 A. In addition, the µFE_lin was calculated astatine VGSVT = 3 V by utilizing the pursuing equation.

$$\mu_{FE\_lin} = \frac{{dI_{D} }}{{dV_{GS} }} \times \frac{L}{{W \times V_{DS} \times C_{ox} }}$$

The DIBL was calculated from the quality betwixt the VT astatine VDS = 0.05 V and the VT astatine VDS = 1.05 V.

Error evaluation

In the lawsuit of measured data, the mean worth and mistake barroom were calculated by utilizing the basal mean quadrate (RMS) worth arsenic follows (x = measurement value, m = mean, n = number of data):

$${\text{Error}}\;{\text{bar }} = \sqrt {\frac{{\sum\nolimits_{{}}^{n} {\left( {x - m} \right)^{2} } }}{n}} \quad {\text{and}}\quad m = \frac{{\sum\nolimits_{{}}^{n} {\left| x \right|} }}{n}.$$

These mean and mistake barroom were utilized successful Fig. 2b,d, and successful Fig. 8d,e.

In addition, for comparing the measured and simulated data, the mistake complaint (ER) was calculated arsenic follows (y = simulation data, mabs = average of implicit values of measured data).

$${\text{RMS }} = \sqrt {\frac{{\sum\nolimits_{{}}^{n} {\left( {y - x} \right)^{2} } }}{n}} ,\quad m_{{{\text{abs}}}} = \frac{{\sum\nolimits_{{}}^{n} {\left| x \right|} }}{n},\quad {\text{ and}}\quad {\text{ER }} = \frac{RMS}{{m_{abs} }} \times 100\;(\% )$$

This ER was utilized successful Fig. 7.

Data availability

The datasets utilized oregon analyzed during the existent survey are disposable from the corresponding writer connected tenable request.

References

  1. Ha, C. et al. Self-aligned coplanar operation for large-sized ultrahigh-definition OLED TV. SID Int. Symp. Dig. Tech. Paper 46, 1020–1022 (2015).

    Article  CAS  Google Scholar 

  2. Yeh, B. & Lin, C. High-performance 4K × 2K 65-in. TV with BCE-type oxide TFTs. SID Int. Symp. Dig. Tech. Paper 46, 943–945 (2015).

    Article  CAS  Google Scholar 

  3. Oota, M. et al. 3D-Stacked CAAC-In-Ga-Zn oxide FETs with gross magnitude of 72 nm. IEEE International Electron Devices Meeting (IEDM), 3.2.1–3.2.4 (2019).

  4. Ishizu, T. et al. A 140 MHz 1 Mbit 2T1C gain-cell representation with 60-nm indium-gallium-zinc oxide transistor embedded into 65-nm CMOS logic process technology. Symposium connected VLSI Circuits, C162–C163 (2017).

  5. Kunitake, H. et al. High thermal tolerance of 25-nm c-axis aligned crystalline In-Ga-Zn oxide FET. IEEE International Electron Devices Meeting (IEDM), 3.6.1–13.6.4 (2018).

  6. Kunitake, H. et al. A c-axis-aligned crystalline In-Ga-Zn oxide FET with a gross magnitude of 21 nm suitable for representation applications. IEEE J. Electron Devices Soc. 7, 495–502 (2019).

    Article  CAS  Google Scholar 

  7. Belmonte, A. et al. Capacitor-less, long-retention (>400s) DRAM compartment paving the mode towards low-power and high-density monolithic 3D DRAM. 2020 IEEE International Electron Devices Meeting (IEDM), 28.2.1–28.2.4 (2020).

  8. Ye, H. et al. Double-gate W-doped amorphous indium oxide transistors for monolithic 3D capacitorless summation compartment Edram. IEEE International Electron Devices Meeting (IEDM), 28.3.1–28.3.4 (2020).

  9. Chakraborty, W. et al. BEOL compatible dual-gate ultra thin-body W-doped indium-oxide transistor with Ion = 370 μA/μm, SS = 73 mV/dec and Ion/Ioff ratio > 4 × 109. Symp. VLSI Technol. 6, 10–11 (2020).

    Google Scholar 

  10. Mo, F. et al. Experimental objection of ferroelectric HfO2 FET with ultrathin-body IGZO for high-density and low-power representation application. 2019 Symposium connected VLSI Technology, T42–T43 (2019).

  11. Ishizu, T. et al. A 48 MHz 880-nW standby powerfulness normally-off MCU with 1 timepiece afloat backup and 4.69-μs wakeup featuring 60-nm crystalline In–Ga–Zn oxide BEOL-FETs. Symposium connected VLSI Circuits, C48 (2019).

  12. Oh, S. et al. Comparison of top-gate and bottom-gate amorphous InGaZnO thin-film transistors with the aforesaid SiO2/a-InGaZnO/SiO2 stack. IEEE Electron Device Lett. 35, 1037–1039 (2014).

    Article  ADS  CAS  Google Scholar 

  13. Oh, S., Baeck, J. H. & Bae, J. U. Effect of interfacial excess oxygen connected positive-bias somesthesia accent instability of self-aligned coplanar InGaZnO thin-film transistors. Appl. Phys. Lett. 108, 41604 (2016).

    Google Scholar 

  14. Choi, S. et al. Systematic decomposition of the affirmative bias accent instability successful self-aligned coplanar InGaZnO thin-film transistors. IEEE Electron Device Lett. 38, 580–583 (2017).

    Article  ADS  CAS  Google Scholar 

  15. Kim, D. H. et al. Experimental decomposition of the affirmative bias somesthesia stress-induced instability successful self-aligned coplanar InGaZnO thin-film transistors and its modeling based connected the aggregate stretched-exponential functions. J. Soc. Inf. Disp. 48, 298–301 (2017).

    Google Scholar 

  16. Hong, S.-Y. et al. Study connected the lateral bearer diffusion and source-drain bid absorption successful self-aligned top-gate coplanar InGaZnO thin-film transistors. Sci. Rep. 9, 6588 (2019).

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  17. Jang, J. T. et al. Cation composition-dependent instrumentality show and affirmative bias instability of self-aligned oxide semiconductor thin-film transistors: Including oxygen and hydrogen effect. ACS Appl. Mater. Interfaces 14, 1389–1396 (2022).

    Article  CAS  PubMed  Google Scholar 

  18. Choi, S. et al. Excessive oxygen peroxide model-based investigation of positive-bias-stress and negative-bias-illumination-stress instabilities successful self-aligned top-gate coplanar In–Ga–Zn–O thin-film transistors. Adv. Electron. Mater. 8, 2101062 (2022).

    Article  CAS  Google Scholar 

  19. Kljucar, L. et al. 300 mm IGZO nFETs with low-T Ru contacts for localized doping and accrued BEOL compatibility. Solid State Devices and Materials (SSDM), J-6-03 (2020).

  20. Rassoul, N. et al. IGZO front-gated TFTs for 3D DRAMs: Process and instrumentality advancement. Solid State Devices and Materials (SSDM), J-6-03 (2021).

  21. Lee, D. H., Nomura, K., Kamiya, T. & Hosono, H. Diffusion-limited a-IGZO/Pt Schottky junction fabricated astatine 200 °C connected a flexible substrate. IEEE Electron Device Lett. 32, 1695–1697 (2011).

    Article  ADS  CAS  Google Scholar 

  22. Nahm, H. & Kim, Y. Undercoordinated indium arsenic an intrinsic electron-trap halfway successful amorphous InGaZnO4. NPG Asia Mater. 6, e143 (2014).

    Article  Google Scholar 

  23. Yang, G. W. et al. Total subgap scope density of states-based investigation of the effect of oxygen travel complaint connected the bias accent instabilities successful a-IGZO TFTs. IEEE Trans. Electron Devices 69, 166–173 (2022).

    Article  ADS  CAS  Google Scholar 

  24. Yen, C. et al. Oxygen-related reliability of amorphous InGaZnO bladed movie transistors. IEEE J. Electron Devices Soc. 8, 540–544 (2020).

    CAS  Google Scholar 

  25. Kamiya, T., Nomura, K., Hirano, M. & Hosono, H. Electronic operation of oxygen deficient amorphous oxide semiconductor a-InGaZnO4-x: Optical investigation and first-principle calculations. Phys. Stat. Solidi (c) 5, 3098–3100 (2008).

    Article  CAS  Google Scholar 

  26. Wager, J. F. Real- and reciprocal-space attributes of set process states. AIP Adv. 7, 125321 (2017).

    Article  ADS  Google Scholar 

  27. Liu, L. et al. Oxygen vacancies: The root of n-type conductivity successful ZnO. Phys. Rev. B Condens. Matter 93, 235305 (2016).

    Article  ADS  Google Scholar 

Download references

Acknowledgements

This enactment was supported successful portion by Samsung Electronics Co., Ltd nether Grant IO200424-07306-01, successful portion by the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (MSIT) of Korea Government nether Grant 2016R1A5A1012966 and 2020R1A2B5B01001979, and successful portion by the Institute of Information and Communications Technology Planning and Evaluation (IITP) funded by the Korea authorities (MSIT) nether assistance 2021-0-01764. Submicron IGZO FET samples was fabricated by imec. TCAD simulation was supported by SILVACO. The EDA instrumentality was supported by the IC Design Education Center (IDEC), southbound Korea.

Author information

Author notes

  1. These authors contributed equally: Donguk Kim and Je-Hyuk Kim.

Authors and Affiliations

  1. School of Electrical Engineering, Kookmin University, Seoul, 02707, Republic of Korea

    Donguk Kim, Je-Hyuk Kim, Woo Sik Choi, Tae Jun Yang, Jun Tae Jang & Dae Hwan Kim

  2. imec, Kapeldreef 75, 3001, Leuven, Belgium

    Attilio Belmonte, Nouredine Rassoul, Subhali Subhechha, Romain Delhougne & Gouri Sankar Kar

  3. Advanced Device Research Lab, Semiconductor R&D Center, Samsung Electronics Co., Ltd., Gyeonggi-do, Hwaseong-si, 18448, Republic of Korea

    Wonsok Lee, Min Hee Cho & Daewon Ha

Contributions

D.K. and D.H.K. designed this enactment and wrote the main portion of manuscript. D.K., J.-H.K., W.S.C., T.J.Y., and J.T.J. performed the instrumentality characterization, modeling and TCAD simulation. A.B., N.R., S.S., R.D., and G.S.K. contributed to the plan and fabrication of the submicron IGZO FETs. D.K., A.B., S.S., R.D., W.L., M.H.C., D.H., and D.H.K. discussed the results. The manuscript was written done contributions of each authors. All authors person fixed support to the last mentation of the manuscript.

Corresponding author

Correspondence to Dae Hwan Kim.

Ethics declarations

Competing interests

The authors state nary competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with respect to jurisdictional claims successful published maps and organization affiliations.

Supplementary Information

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kim, D., Kim, JH., Choi, W.S. et al. Device modeling of two-steps oxygen anneal-based submicron InGaZnO back-end-of-line field-effect transistor enabling short-channel effects suppression. Sci Rep 12, 19380 (2022). https://doi.org/10.1038/s41598-022-23951-x

Download citation

  • Received: 18 August 2022

  • Accepted: 08 November 2022

  • Published: 12 November 2022

  • DOI: https://doi.org/10.1038/s41598-022-23951-x

Read Entire Article