FACTA UNIVERSITATIS Series: Electronics and Energetics Vol. 30, N o 2, June 2017, pp. 257 - 265 DOI: 10.2298/FUEE1702257S ANALYTICAL TEST STRUCTURE MODEL FOR DETERMINING LATERAL EFFECTS OF TRI-LAYER OHMIC CONTACT BEYOND THE CONTACT EDGE  Neelu Shrestha, Geoffrey K. Reeves, Patrick W. Leech, Yue Pan, Anthony S. Holland School of Engineering, RMIT University, Melbourne, Victoria, Australia Abstract. Contact test structures where there is more than one non-metal layer, are significantly more complex to analyse compared to when there is only one such layer like active silicon on an insulating substrate. Here, we use analytical models for complex test structures in a two contact test structure and compare the results obtained with those from Finite Element Models (FEM) of the same test structures. The analytical models are based on the transmission line model and the tri-layer transmission line model in particular, and do not include vertical voltage drops except for the interfaces. The comparison shows that analytical models for tri-layer contacts to dual active layers agree well with FEM when the Specific Contact Resistances (SCR) of the contact interfaces is a significant part of the total resistance. Overall, there is a broad range of typical dual-layer-to-TLTLM contacts where the analytical model works. The insight (and quantifying) that the analytical model gives on the effect of the presence of the contact, on the distribution of current away from the contact is shown. Key words: Ohmic contact, specific contact resistance, transfer length, Transmission Line Model 1. INTRODUCTION The study of specific contact resistivity (c, [Ω.cm 2 ]), one of the most important parameters to investigate metal semiconductor interfacial properties, has been reported using several test structures [e.g. 1-6]. The effect of this parameter on the resistance of device components is also an area of investigation and an example is the effect of interfaces in via liners [7]. The transmission line model (TLM) is the commonly used among them and was first applied for analysing/determining the resistance components of semiconductor ohmic contacts, by Shockley in 1964 [5]. For many years the TLM test structure was considered accurate for contacts such as Aluminium to Silicon. A more complex TLM model for contacts such as alloyed Ni/Ge/Au to GaAs was introduced by Received September 26, 2016; received in revised form November 16, 2016 Corresponding author: Neelu Shrestha School of Engineering, RMIT University, Melbourne, Victoria, Australia (E-mail: neelu_shrestha@hotmail.com) 258 N. SHRESTHA, G. K. REEVES, P. W. LEECH, Y. PAN, A. S. HOLLAND Reeves [6]. The example given has a metal, an alloyed and a semiconductor layer. Another example of a tri-layer contact structure is Aluminium-silicide–Silicon. The analytical models due to the dual-layer and tri-layer structures considered the contact as beginning at the leading edge of the metal. Yao Li [8] developed an analytical model to show (and quantify) how the presence of the contact affects the current (and hence voltage distribution) away from the leading edge. This was also investigated by Reeves [9]. The works in these references [6, 8, 9] were combined in the investigation for this paper and their utilisation and further demonstration of accuracy is demonstrated. The analytical model demonstrated in this paper can provide a solution to the 2L-TLM, contact structure (a metal contact to a dual semiconductor layer (see Fig. 1(a))), considered intractable in 1994 [10]. For most semiconductor devices the sheet resistance of the semiconductor active layer is sufficient to quantify its resistive effect even for typical planar contacts such as source/drain contacts. The principle of the TLM models for ohmic contacts relies on this being so. Only in cases where the specific contact resistance of the metal/semiconductor interface is extremely small, there will be significant vertical voltage drop in the semiconductor layer. TLM models (analytical type models) do not account for vertical voltage drop in the semiconductor layer but FEM models do. This paper investigates the difference. FEM is a powerful technique and is used here to demonstrate the accuracy of the analytical models which are also expanded here to demonstrate their usefulness as one approach to quantifying the total resistance encountered by current between two Tri-layer contacts to a dual-layer (parallel sheet resistances) structure as shown in Fig 1(a) and the corresponding resistor network model shown in Fig 1(b). Fig. 1 Test structure for investigating resistance effects of contacts to dual-active layers. (a) Schematic of test structure and (b) resistor network model showing the TLTLM [6] contact model connected to the dual active layer model [9]. The one resistor that is common to the resistor network of both analytical models is indicated Analytical Test Structure Model for Determining Lateral Effects of Tri-Layer Ohmic Contact... 259 2. METHODOLOGY The general tri-layer contact structure investigated in this work is illustrated diagrammatically with its three layers, namely a metal layer, an intermediary active layer A and a bottom active layer B in Fig.1. The metal layer is an ideal metal layer for TLM models where the metal is at an equipotential. The FEM model considers this by making the metal in the FEM metal have the same effect by having extremely low resistivity, so it is effectively at an equipotential. Two interfaces with associated specific contact resistances (SCR) are included in the model and these occur between the metal and layer A and layer A and B with their SCR given as ca and cu respectively. The sheet resistance of the metal at the metal-semiconductor is considered to be zero and, Rsa and Rsu are the sheet resistance of layer A and B respectively. (The subscripts -sa, -su, -ca and -cu are used to maintain uniformity in the expressions with Refs [6] and [9].) The total current given to the structure is i0. i1 and i2 represent the current in layer B and A respectively. Here, i3 is the total current exiting through the contact. The length of the contact is d, the length between two contacts is l and w is the width of the structure. The analytical expressions for the calculation of total resistance are derived combining the Reeves’ Tri-layered Transmission Line Model (TLTLM) [6] and Yao Li’s Model [8]. Because the models are based on the transmission line model (TLM), (which can at best be regarded as a 2-D model) the vertical voltage dropped is not included except for the interfaces. As defined by Berger [11], the  parameter gives the ratio of voltage drop across the contact interface to the vertical voltage drop in the semiconductor layer and thus determines whether a metal and a semiconductor contact is in 3D circumstances or not (FEM is utilised in this work to investigate 3-D effects).  = ρc / (ρb * t) (1) where, ρb and t are the resistivity and thickness of the individual semiconductor layer respectively. Consequently, the model presented in this work best suits to modelling contacts when  is typically greater than 1 because the model will be 2-D. When the two models, dual active layer [8] and TLTLM [6] are combined to describe the test structure, the voltage drop across the common resistor (see insert in Fig.1) between two models must be the same and thus the current division factor, f can be determined by solving for voltage using this boundary condition. Reeves [9] in 1997 reported a similar structure applied to the source/drain region of a MOSFET (with a short extended dual (silicide/silicon) active layer before the contact) and considered all current entering the dual layer through layer B. In this work, we focus on the more general case of a test structure in which current enters the dual layer through both the layers as shown in Fig.1. The assumption made is that the two contacts are geometrically and electrically identical. In the following equations, f1 is the current division factor at the intersection where current is leaving one TLTLM contact and entering the dual active layer. Factor f is the modified current division factor at the intersection of a TLTLM contact and the dual layer while current is entering the TLTLM contact from the dual active layer. Thus the boundary condition used is i1(x)=iof1 at x=l and i1(x)=iof at x=0. The expression of modified current division, f is given by. ( ) ( ) (2) 260 N. SHRESTHA, G. K. REEVES, P. W. LEECH, Y. PAN, A. S. HOLLAND where, Transfer Length for the dual active layer, √ (3) ( ) (4) ( ( ) ( ( ) ) ( ) ( ) ) (5) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (6) ( ( ) ( )) ( ) ( ) ( ) (7) Reeves [6] demonstrated that the current (i1(x)) through layer B and contact resistance (Rc) in the TLTLM are given by ( ) ( ( ) ) ( ( ) ( ) ( ) ( ) ) (8) { ( ) ( ) } (9) where, ( ) ( ) (10) ( ) ( ) (11) ( ) (12) The equation for a, b and c are shown in the appendix A1, A2 and A3. Using the work reported by Li et al. [8], the current (i1(x)) through the lower layer B of the dual-layer structure and total resistance, Rtot(Li) is given by ( ) ( ( ) ( ) ( ) ( ) ( ) ( ) ) (13) ( ) ( ) (14) where, (15) ( ) (16) For two contacts (in a TLM type test structure) with identical geometrical and electrical features, the total resistance, Rtot(Std) between probes connected to each contact is usually given by ( ) (17) The numerical evaluation of total resistance and current flow in all layers with different contact parameters giving various  values are shown using the table and graphs obtained from MATLAB. Also, models with similar contact parameters were simulated using Finite Element Modelling (FEM) and the results for Rtot(FEM), Rtot(Li) and Rtot(Std) are compared in this work. Analytical Test Structure Model for Determining Lateral Effects of Tri-Layer Ohmic Contact... 261 3. RESULTS AND DISCUSSION The table below gives sets of contact material parameters used for analysis of the test structure of Fig. (1). The values a and u for two layers, (=ρc/(ρb*t), where ρb and t are the semiconductor layer resistivity and thickness, differs in sets due to different contact parameters. The contact length is 10m and the length between the contacts (dual active layer) is 20m for set 1, 2 and 3. But, the length for set 4 is 5.2m. The thickness of layers A and B are considered as 0.2 and 0.6m respectively. The distances were chosen as being typical of a TLM type test structure. The value came about because of the mesh density used, which was the same for all models. Table 1 Contact Parameters for the test structure Set Rsa(Ω/sq) Rsu(Ω/sq) ρca (Ω.cm 2 ) ρcu (Ω.cm 2 ) l (m) a u Rtot (FEM) Rtot (Li) Rtot (Std) 1 40 60 8.0E-7 1.6E-6 20 50 7 576.1 575.7 587.8 2 30 40 8.0E-7 9.0E-6 20 66.67 62.50 445.9 445.4 482.9 3 40 60 8.0E-9 1.0E-8 20 0.50 0.04 487.4 507.8 507.7 4 40 60 8.0E-7 1.6E-6 5.2 50 7 220.8 220.1 232.4 The analytical and FEM results are in excellent agreement for sets when the  value is greater than 1. Moreover, the error percentage is less between Rtot(FEM) and Rtot(Li) results compared to Rtot(Std). This clearly shows that the redistribution of current near the intersection of TLTLM and the dual layer affects the total resistance between the contacts. In sets 1 and 2, Rtot(FEM) and Rtot(Li) are practically the same and the error percentage is less than 1 in both cases(considering the FEM model to give the most realistic result). However, the error is 2% between Rtot(FEM) and Rtot(Std) for set 1 and it increased to 8% for set 2. This is because the transfer length in the dual active layer is longer for set 2 than set 1. The distance before at the leading edge of the contact that is affected by the redistribution of current due to the contact is directly proportional to the transfer length of the dual active layer. Consequently, the total resistance in the dual layer of the test structure will be different to Rsh*l/w (where Rsh is the effective sheet resistance of the two parallel active layers). In set 3, both Rtot(Std) and Rtot(Li) have approximately 4% disagreement with the FEM results. In set 4, all the parameters are the same as set 1 except that the length between the contacts is 5.2m. As the dual active layer is shorter, and the transfer lengths do not change, the effect of the contacts on current distribution in the dual layer is relatively bigger. The error is ~5% between Rtot(FEM) and Rtot(Std) and less than 1% with Rtot(Li). Thus, the expressions described in this work can give an accurate insight into the electrical behaviour of the test structure, provided that >1. The graph plotted below in Fig. 2(a) is for set 1 and the point of contact for TLTLM and dual layer is considered to be point 0. The length of the TLTLM contact region is 10 µm and the dual layered region is 20µm. It shows the flow of current starting from the point at which current exits one TLTLM contact and enters the dual layer. Here, i1 (solid line) is the current flowing in layer B, i2 (dotted line) in layer A and i3 (fine dotted line) represents the total current in the TLTLM contact. The graph illustrates that i1 distributes itself and tends to remain constant at a ratio of Rsa:Rsu (i.e. i1:i2=4:6) in the dual layered structure. But, the flow is disturbed near the point when current enters the TLTLM structure from the dual layer. 262 N. SHRESTHA, G. K. REEVES, P. W. LEECH, Y. PAN, A. S. HOLLAND The FEM result confirms the same as shown in Fig. 2(b). The voltage contour is uniform at the middle of the dual layer, and the contour is disturbed near the TLTLM contact region. Fig. 2 (a) Current distribution in one TLTLM contact and dual layer contact regions for a test structure with Rsa=40 Ω/sq, Rsu= 60Ω/sq, ρca=8e-7Ω.cm 2 , ρcu= 1.6e-6Ω.cm 2 with dual active layer of 20m (b) distribution of corresponding voltage contours, determined by FEM for the test structure Fig. 3(a) shows results for set 4 when all the other parameters are same as set 1 except the length between the contacts is reduced to 5.2m. Because of this, the current flow with ratio Rsa:Rsu is over a shorter length. This leads to greater variation in total resistance. The higher error value of ~5% between Rtot(Std) and Rtot(FEM) shows that the total resistance between contacts is not given by Rsh*l/w. Fig. 3 (a) Current distribution in a TLTLM contact and dual layer contact regions for a test structure with Rsa=40 Ω/sq, Rsu= 60Ω/sq, ρca=8e-7Ω.cm 2 , ρcu= 1.6e-6Ω.cm 2 with a dual active layer of 5.2m in length (b) distribution of corresponding voltage contours, determined by FEM for the test structure Analytical Test Structure Model for Determining Lateral Effects of Tri-Layer Ohmic Contact... 263 4. CONCLUSION Two variations of Transmission Line Model networks, namely Tri-Layer TLM and a dual layer network, for modelling current in semiconductor contact regions, were combined to model a test structure with multiple layers. A comparison of the mathematical analysis and a two-dimensional finite element model of the test structure with two metal contacts to a dual- active layer, show that the combination of TLTLM and the dual-layer network expressions provides accurate analysis for these test structures. The limitations on the accuracy of expressions have been presented in terms of the  parameter. The distribution of current through the dual-layer and TLTLM contact region is discussed in detail to understand its influence on the total resistance of the test structure. This distribution is accurately represented by the combined TLTLM-dual-active layer model investigated which is an improvement on models where the current distribution and sheet resistance is considered uniform between contacts. REFERENCES [1] A. M. Collins, Y. Pan, A. S. Holland, “Using a Two-contact circular test structure to determine the specific contact resistivity of contacts to bulk semiconductors”, Facta Universitatis, Electronics and Energetics, vol. 28, no. 3, September 2015, pp. 457-464. [2] Y. Pan, A. M. Collins and A. S. Holland, "Determining Specific Contact Resistivity to Bulk Semiconductor Using a Two-Contact Circular Test Structure", In Proceedings of the IEEE International Conference on MIEL, May 2014, pp. 257-260. [3] V. Gudmundsson, P. Hellstrom, and M. Ostling, “Error propagation in contact resistivity extraction using cross-bridge Kelvin resistors,” IEEE Trans. Electron Devices, vol. 59, no. 6, pp. 1585–1591, June 2012. [4] A. S. Holland, G. K. Reeves, “New Challenges to the Modelling and Electrical Characterisation of ohmic contacts for ULSI devices”, In Proc. of the 22nd International Conference on Microelectronics (MIEL 2000), vol. 2, Niš, Serbia, pp. 461-464, 2000. [5] W. Shockley, “Research and Investigation of inverse epitaxial UHF power transistors”, Air Force Atomic Laboratory, Wright-Patterson Air Force Base, Rep. No. AL-TDR-64-207, Sept. 1964. [6] G. K. Reeves, H. B. Harrison, “An analytical model for alloyed ohmic contacts using a Tri-layer Transmission Line Model”, IEEE Trans. Electron Devices, vol. 42, no. 8, p. 1536. 1995. [7] G. K. Reeves, A. S. Holland, P.W. Leech, “Influence of Via Liner Properties on the Current Density and Resistance of Vias”, In Proc. of the 23rd International Conference on Microelectronics (MIEL 2002), vol. 2, Niš, Yugoslavia, pp. 535-538, 2002. [8] Y. Li, G. K. Reeves, H. B. Harrison, “Correcting separating errors related to contact resistance measurement”, Microelectronics Journal, vol. 29, 1996. [9] G. K. Reeves, H. B. Harrison, “Using TLM principles to determine MOSFET contact and parasitic resistance”, Solid-State Electronics, vol. 41, no.8, 1997. [10] Y. Shiraishi, N. Furuhata, A. Okhamoto, “Influence of metal/n-InAs/interlayer/n-GaAs structure on nonalloyed ohmic contact resistance”, Journal of Applied Physics, vol. 76, p. 5099, 1994. [11] H. H Berger, “Models for contacts to planar devices”, Solid State Electronics, vol. 12, 1972. APPENDIX The expressions for the TLTLM √* { √( ( )⁄ )} + (A1) √ { √( ( )⁄ )} (A2) (A3) 264 N. SHRESTHA, G. K. REEVES, P. W. LEECH, Y. PAN, A. S. HOLLAND MATLAB CODE %Calculation of f factor Rs=Rsa+Rsu; Rsh=(Rsu*Rsa)/(Rsu+Rsa); al=sqrt(Rs/pcu); c=(Rs/pcu)+(Rsa/pca); z=sqrt((c*c)-(4*Rsu*Rsa/(pcu*pca))); a=sqrt((c-z)/2); b=sqrt((c+z)/2); D=al*pcu*coth(al*l); E=al*pcu*(Rsa*cosh(al*l)+(f1*Rs-Rsa))/(Rs*sinh(al*l)); H=((b*((Rsu-(pcu*a*a))*tanh(a*d)))-(a*((Rsu-(pcu*b*b))*tanh(b*d))))/((b*b- a*a)*tanh(b*d)*tanh(a*d)); G=Rsa*(b*tanh(a*d)-a*tanh(b*d))/((b*b-a*a)*tanh(b*d)*tanh(a*d)); f2=(E+G)/(D+H+G); % Current Flow Dual Active Layer A=Rsh/Rsu; Y1=sqrt(pcu/(Rsu+Rsa)); i11=i0*(A+((f2-A)*sinh((l-x)/Y1)/sinh(l/Y1))+((f1-A)*sinh(x/Y1)/sinh(l/Y1))); i21=i0-i11; % Current Flow TLTLM Contact c=((Rsa+Rsu)/pcu)+(Rsa/pca); z= (c*c)-(4*Rsu*Rsa)/(pcu*pca); a=sqrt((c-sqrt(z))/2); b=sqrt((c+sqrt(z))/2); P= f2*(Rsu-pcu*a*a)-(1-f2)*Rsa; Q= f2*(Rsu-pcu*b*b)-(1-f2)*Rsa; i12=(i0/(pcu*(b*b-a*a)))*((P*sinh(b*(d+y))/sinh(b*d))-(Q*sinh(a*(d+y))/sinh(a*d))); i23=(i0/(Rsa*pcu*(b*b-a*a)))*((P*(Rsu-pcu*b*b)*sinh(b*(d+y))/sinh(b*d))-(Q*(Rsu- pcu*a*a)*sinh(a*(d+y))/sinh(a*d))); itot=i0-(i12+i23); plot(x,i11); hold on; plot(y,i12); hold on; plot(x,i21); hold on, plot(y,i23); hold on; plot(y,itot); % Contact Resistance c1=((Rsa+Rsu)./pcu)+(Rsa./pca); z1= (c1.*c1)-((4.*Rsu.*Rsa)./(pcu.*pca)); a1=sqrt((c1-sqrt(z1))/2); b1=sqrt((c1+sqrt(z1))/2); K1=Rsu./(pcu.*w.*(b1.*b1-a1.*a1)); Analytical Test Structure Model for Determining Lateral Effects of Tri-Layer Ohmic Contact... 265 X1=tanh(b1*d); Y1=tanh(a1*d); P1=f2.*(Rsu-(pcu.*a1.*a1)); Q1=f2.*(Rsu-(pcu.*b1.*b1)); k1=(1-f2).*Rsa; Rc=K1.*(((P1-k1)./(b1.*X1))-((Q1-k1)./(a1.*Y1))); % Total Resistance using standard Formula Rtot=2*Rc+(Rsh*l/w) % Total Resistance using Yao et al. Formula when dual layer is longer beta=2*((((f1+f2)*Rsu)/(2*Rsh))-1); bcor=(Rsh*beta)/(w*al); Rtotli=2*(Rc+(bcor/2))+(Rsh*l)/w;