A Two Dimensional Tunneling Resistance Transmission Line Model for Nanoscale Parallel Electrical Contacts

Contact resistance and current crowding are important to nanoscale electrical contacts. In this paper, we present a self-consistent model to characterize partially overlapped parallel contacts with varying specific contact resistivity along the contact length. For parallel tunneling contacts formed between contacting members separated by a thin insulating gap, we examine the local voltage-dependent variation of potential barrier height and tunneling current along the contact length, by solving the lumped circuit transmission line model (TLM) equations coupled with the tunneling current self consistently. The current and voltage distribution along the parallel tunneling contacts and their overall contact resistance are analyzed in detail, for various input voltage, electrical contact dimension, and material properties (i.e. work function, sheet resistance of the contact members, and permittivity of the insulating layer). It is found the existing one-dimensional (1D) tunneling junction models become less reliable when the tunneling layer thickness becomes smaller or the applied voltage becomes larger. In these regimes, the proposed self-consistent model may provide a more accurate evaluation of the parallel tunneling contacts. For the special case of constant ohmic specific contact resistivity along the contact length, our theory has been spot-checked with finite element method (FEM) based numerical simulations. This work provides insights on the design, and potential engineering, of nanoscale electrical contacts with controlled current distribution and contact resistance via engineered spatially varying contact layer properties and geometry.

For decades, the basic models of tunneling current between electrodes separated by thin insulating films have been those of Simmons [36][37][38] in 1960s. Simmon's formula have since been used for evaluating tunneling current in tunneling junctions 24,29,39 . Though there have been attempts to extend Simmons' models to the field emission and space-charge-limited regimes 11,40,41 , it is always assumed that the tunneling junctions are one-dimensional (1D), i.e. there is no variations on the voltages drops along the length of the tunneling junction and the insulating film thickness is uniform. Thus, these existing models of tunneling junctions give no hint on the variation of tunneling current along the contact length and the importance of current crowding near the contact area when the two contacting members are partially overlapping (cf. Fig. 1). On the other hand, the widely used transmission line models (TLM) for electrical contacts typically assume the contact resistivity of the interface layers are constant 10, [42][43][44] . It is questionable to apply these models to study the tunneling contacts, as the tunneling resistance depends on the junction voltage that varies spatially along the contact length.
In this paper, we propose a two-dimensional (2D) transmission line model for partially overlapped parallel contacts with spatially varying specific contact resistivity. Spatial dependence of specific contact resistivity of the contact interface may be introduced by many factors, such as nonuniform distribution of the resistive contaminants, oxides, or foreign objects at the contact interface, formation of contact interfaces with spatially varying thickness, or the presence of tunneling contacts between contact members. In the latter case, because of the nonlinear current-voltage characteristics of the tunneling junctions 11,36 , the specific resistivity along the contact length will become spatially dependent, even for a tunneling layer with uniform thickness (Fig. 1). For the tunneling-type contacts, the model considers the variation of potential barrier height and tunneling current along the contact length, by solving the TLM equations coupled with the tunneling current self consistently. We provide comprehensive analysis of the effects of contact geometry (i.e. dimension of the contact, and distance between the contact electrodes), and material properties (i.e. work function, sheet resistance of the contact members, and permittivity of the insulating layer) on the spatial distributions of currents and voltages across these contacts, and the overall contact resistance of parallel contacts.
The formulation of our 2D contact resistance TLM model is given in the next section. We would like to point out that, albeit an application of the standard transmission line theory based on the Kirchhoff 's laws, the TLM has been used extensively with great success to characterize mesoscale and nanoscale electrical contacts 10, 30,42,43 .
Here we further extend the TLM model with the effects of spatially dependent contact resistivity. In the results and Discussion section, we consider three cases of parallel contacts: (1) constant specific contact resistivity, (2) linearly varying specific contact resistivity, and (3) tunneling contact resistivity depending on local junction voltages along the contact length. The first case of uniform specific contact resistivity along the contact length has been verified with COMSOL 45 2D simulations. For the last case, for simplicity, we use the Simmons' model [36][37][38] to determine the local current-voltage characteristics across the tunneling junction. Though full scale quantum mechanical calculations may have to be used to accurately evaluate the nanoscale circuits, our model based on Simmons formula reveals the fundamental scalings and parametric dependence of current and voltage profiles, as well as electric contact resistance of tunneling contacts. Summary and suggestions for future research are given in the last section.
Note that, although this paper is focused on the normal Schrödinger tunneling type electrical contacts, the proposed TLM with spatially varying contact resistivity can be used for many other types of electric contacts, such as nanoscale Schottky contacts based on 2D materials heterostructure 46,47 , and Klein tunneling junctions 48 . formulation Consider a parallel contact formed between two nanowires or nanotubes or between two conducting thin films or layers, as shown in Fig. 1(a,b), respectively. The distance between the two contact members is D, and the contact length is L. A thin resistive interface layer is sandwiched between them. Both contacts in Fig. 1(a,b) can be described by a two-dimensional (2D) model, as shown in Fig. 1(c). Note that the proposed formulation is generally applicable to parallel Cartesian nanojunctions with different shape of the electrodes, for example, electrical contact between a nanowire and a thin film. In the 2D model, the effects of the transverse dimension (perpendicular to the paper) can be included in the effective sheet resistances R sh1 and R sh2 for conductor 1 and 2, respectively, such that there is no variation along the width w in the transverse dimension. The spatial dependent specific interfacial resistivity (also termed specific contact resistivity) is ρ c (x), which is either predefined, or calculated from the local tunneling current in case of insulating tunneling layer [36][37][38] . We use the DC equivalent lump circuit transmission line model (TLM) 10,42-44 , as shown in Fig. 1(d), to model the 2D parallel contact in Fig. 1(c).
In the contact region PQNM in Fig. 1(c,d), using Kirchoff 's laws for current and voltage, we get the following equations, where I 1 (x) and I 2 (x) represent the current flowing at x through the lower contact member, MN and the upper contact member, PQ respectively, and V 1 (x) and V 2 (x) the local voltage at x along MN and PQ, respectively, and w is the effective transverse dimension of the contacts. When Δx → 0, Equation (1) becomes, where J c (x) = V g (x)/ρ c (x) and V g (x) = V 1 (x) − V 2 (x) are the local current density and the local voltage drop across the contact interface at x, respectively. Note that, from Eqs (2a) and (2c), I 1 (x) + I 2 (x) = I tot = constant, where I tot is the total current in the circuit, to be determined from the boundary conditions. The boundary conditions for Eq. (2) are, 2 where, without loss of generality, we assume the voltage of the upper contact member at x = L is 0, and the externally applied voltage at x = 0 of the lower contact member is V 0 . Note that I 1 (x = 0) = I tot . From Eqs (2) and (3), it is easy to show , where a prime denotes a derivative with respect to x. For the contact model in Fig. 1(d), the contact resistance is defined as, It is convenient to introduce non-dimensional quantities (2) can be recast into the following second order differential equations, www.nature.com/scientificreports www.nature.com/scientificreports/ . The corresponding boundary conditions to Eq. 5(a-c) are respectively, and (  1) ( 1), (6a) 1 1 where the unknown constant α = I tot /I o is the normalized total current in the circuit, and prime denotes a derivative with respect to x . Note that integrating Eq. (5b) subject to Eq. (6b) gives ∫ , which means that the total current is conserved across the contact interface.
Equations (5) and (6) are solved to give the voltage and current distribution along and across the contact interface as well as the total contact resistance, for a given electrical contact ( Fig. 1) with spatially dependent interface specific contact resistivity ρ x ( ) c . An example of the procedure to solve Eqs (5) and (6) numerically is as follows. For an initially guess on α, Eq. (5b) is solved using the shooting method, subject to Eq. (6b). Next, Eq. (5a) is solved with the initial values of V(0) 1 and is satisfied. Finally, Eq. (5c) is solved to get I 1 (and I 2 ).
In principle, Eqs (5) and (6) can be solved numerically for arbitrary spatial dependence of specific contact Here, we focus on a few special cases of practical importance. We first consider the case of constant ρ c , where analytical solutions can be obtained, which also serve to validate our numerical approach. We then consider the effects of spatially dependent ρ x ( ) c on the parallel electrical contacts. We focus on two situations: linearly varying specific contact resistivity along x, and thin tunneling junction with uniform thickness, where analytical solutions to the TLM current and voltage equations are no longer available, and Eqs (5) and (6) are solved numerically.

Results and Discussion
Case 1: Constant specific contact resistivity ρ c along the contact length L. For the special case of constant specific contact resistivity ρ c , the TLM equations, Eqs (5) and (6), can be solved analytically to give, . Figure 2 shows the current and voltage distributions along the contact length and across the contact interface for various specific contact resistivity ρ c , for a parallel contact formed between similar contact members, . The voltage along both contact members V 1 and V 2 decrease with x , as shown in Fig. 2(b,c), respectively. The current I 1 in contact member 1 decreases with x ( Fig. 2(e)), whereas I 2 in contact member 2 increases with x ( Fig. 2 www.nature.com/scientificreports www.nature.com/scientificreports/ length, with the minimum at the center of the contact structure x = 0.5 and the maximum at the contact edges, as shown in Fig. 2(a,d), respectively. The current crowding effects near the contact edges are well-known phenomena 42,49 , as the current density is distributed to follow the least resistive path (i.e. minimum overall resistance). It is important to note that as the specific contact resistivity ρ c decreases, the interface current density J c becomes more crowded towards the contact edges, as shown in Fig. 2(d). In other words, the less resistive the contact interface layer, the more severe of the current crowding effects, which is in agreement with previous studies using both TLM 42,44 and field theory 10, 49,50 . Figure 3 shows the current and voltage distributions along the contact length and across the contact interface for various parallel contacts formed between dissimilar materials, R sh2 = R sh2 /R sh1 , with fixed specific contact resistivity ρ c = 1. The voltage V 1,2 and the current I 1,2 along the two contact members show similar behaviors as those in Fig. 2. However, the voltage drop across the interface layer V x ( ) g and the contact current density J x ( ) C are no longer symmetric, as shown in Fig. 3(a,d), respectively. When R sh2 < 1, the maximum of V x ( ) occurs at x = 1. This current crowding effect can again be explained by the fact that current flows are self-arranged to take the least resistive path in the circuit by adjusting the current distribution according to the local resistance.
The normalized contact resistance, R c calculated from Eq. (8) is plotted in Fig. 4 for various ρ c and R sh2 . It is clear that R c increases with both ρ c and R sh2 . In general, the contact resistance R c depends more strongly on the the specific contact resistivity of the interfacial layer ρ c than on the sheet resistance ratio of the contact members R sh2 . For the special case of R sh2 = 0, Eq. (8) becomes, , which is also plotted in Fig. 4. Note that Eq. (9) is identical to the expression typically used for metal-semiconductor contact 10,42 .
To verify the results obtained from our analytical solution, we performed numerical simulations using the COMSOL multiphysics software 45 , for various combinations of R sh2 and ρ c on the geometry shown in Fig. 1(c). The finite-element-method (FEM) based COMSOL 2D simulation results are included in Fig. 4 (cross symbols), showing excellent agreement with our theory. The convergence iteration error in the simulation was less than 10 −9 for each point.
Case 2: Specific resistivity ρ c varies linearly along the contact length L. We assume the specific resistivity varies linearly along the contact length ( Fig. 1) By solving Eqs (5) and (6) numeri- , current along (e) contact member 1, I x ( ) 1 , and (f) contact member 2, I x ( ) 2 ,for different values of specific contact resistivity ρ c , for R sh2 = R sh2 /R sh1 = 1. All the quantities are in their normalized forms defined in the Formulation section. www.nature.com/scientificreports www.nature.com/scientificreports/ cally, we obtain the current and voltage distributions along the contact interface, as shown in Fig. 5. As A increases, the overall contact interface becomes more resistive, therefore, the voltage drop V x ( ) g across the interface layer increases (Fig. 5a), whereas the current density J x ( ) C across the interface layer decreases in general (Fig. 5d). The maximum V g occurs at the contact edge with the highest specific resistivity ρ c (i.e., at x = 0 when A < 0, and at x = 1 when A > 0), while the maximum interface current J C occurs at the contact edge with the lowest   (Fig. 1). R c as a function of (a) normalized specific contact resistivity, ρ c and (b) normalized sheet resistance of contacting member 2, R sh2 . Dashed lines are for Eq. (9), the limiting case of R sh2 → 0. The cross symbols are from COMSOL 45 2D simulations. In the simulation, the length and height of both upper and lower contacting members are assumed to be 20 nm and 10 nm respectively, and the thickness of the resistive interfacial layer is assumed to be 0. . ρ c (i.e., at x = 1 when A < 0, and at x = 0 when A > 0). The effects of A on the voltage V 1,2 and the current I 1,2 along the two contact members are also shown in Fig. 5(b,c,e,f), respectively. The normalized contact resistance, R c calculated from Eq. (4) for linear specific contact resistivity c is plotted in Fig. 6. As A increases, R c increases, since the contact interface becomes more resistive. As R sh2 increases, the contact resistance R c depends more strongly on the linear constant A. Case 3: Tunneling contact resistance. Here, we assume the parallel contacts are formed through a tunneling interface layer between the two contact members. For simplicity, we have made the following assumptions: www.nature.com/scientificreports www.nature.com/scientificreports/ (1) the thickness of interfacial insulating film in the contact area is uniform, and (2) the insulating film is sufficiently thin (in the nano-or subnano-meter scale) so that charge trappings are ignored 51,52 .
For dissimilar contact members, the (normalized) current density at any location along the contact from contact member 1 to contact member 2 is calculated using Simmons' formula 37 , is the local voltage drop across the contact interface at x , = .
, ϕ 1 = W 1 − χ and ϕ 2 = W 2 − χ. W 1 and W 2 are the work functions of contacting member 1 and 2 respectively, χ is the electron affinity of the insulating layer, which is 0 for vacuum.
On the other hand, the definitions of ϕ I , y 1 and y 2 for reverse bias (when higher work function contacting member is given positive bias) are, ; For the special case of the same material for contact members 1 and 2 36 Note that we use Simmon's formula, Eq. (10) here for simplicity, which is reliable only when the barrier height is relative high and the gap voltage is low in the direct tunneling regime 11 . More accurate results for the tunneling current may be calculated using quantum models by solving the coupled Schrodinger equation and Poisson equation with the inclusion of space charge and exchange-correlation effects 11,25 .
We keep the normalization consistent with our previous calculations in the Formulation section. For a given parallel tunneling contact (Fig. 1), the inputs of our model are the applied voltage V 0 , sheet resistance (R sh1 , R sh2 ) and work function (W 1 , W 2 ) of contacting members 1 and 2, permittivity (ε r ), thickness (D), and electron affinity (χ) of the interfacial insulator layer, and the contact length L. Using Eq. (10), the specific contact resistivity is , which is inserted into the TLM equations, Eqs (5) and (6), to give a self-consistent solution to the voltage and current profiles, as well as the contact resistance for the parallel tunneling contact.
We consider CNT-vacuum-CNT parallel contact as an example. Both contact members are made of the same single-walled CNTs. Using the typical value of linear resistivity of single-walled CNT ρ L = 20 kΩ/μm 53,54 , and diameter (or the width w) of 3 nm, an equivalent sheet resistance for both CNT contact members are estimated as where the unit of the sheet resistance  Ω/ means "ohm per square" 10, 42,49 . The work function of CNTs is W 1 = W 2 = 4.5 eV 55 . The interfacial layer is assumed to be vacuum (relative permittivity ε r = 1.0, and electron affinity χ = 0). The voltage drop V g (x) across and the tunneling current density J c (x) through the contact interface are shown in Fig. 7 for various contact length L, vacuum gap distance D, and applied voltage V 0 . The profiles of both V g (x) and J c (x) are symmetric about the center of the contact, as expected for similar contact members (similar to Fig. 2a,d above). As the contact length L increases, the local voltage drop V g (x) across the contact interface decreases, so does the tunneling current density J c (x), as shown in Fig. 7a,b. However, the total current in the contact structure, increases with L, since the total contact resistance of the tunneling junction decreases as the contact length increases (cf. Fig. 8a below). As shown in Fig. 7c,d, when the gap distance D increases, the voltage drop V g (x) increases, but the current density J c (x) decreases, which is because the tunneling junction becomes more resistive 11,36 . Figure 7e,f shows both voltage drop V g (x) and current density J c (x) increase when the applied voltage V o increases. More importantly, both V g (x) and J c (x) exhibit a stronger spatial dependence as V 0 increases. This strong voltage dependence of electrical properties of the tunneling junction is in sharp contrast with those of ohmic contacts, where the shape of the profiles of V g (x) and J c (x), and the total contact resistance are independent of the applied voltage, and the current density scales linearly with the voltage drop, as discussed in Cases 1 and 2 above.
Also plotted in Fig. 7 are the analytical results from Eq. (7), by assuming constant tunneling contact resistivity across the contact length L (i.e. the typically assumed one-dimensional tunneling junction 24 ), by (a), setting V g = V 0 and using Eq. (10) (dashed lines) and (b), using ohmic approximations for the tunneling junction, in the limit of V g → 0 (dotted lines) 36,37 . In the latter case, the tunneling current density is a linear function of V g .
[Å] for similar contacting members. λ , ϕ Δ , ϕ 2 , ϕ 1 , ϕ 0 have the same definition as in Eq. (10). It is found that both assumptions of constant contact resistivity are not sufficiently reliable, especially when the tunneling thickness D decreases or the applied voltage V o increases. As the tunneling junction resistance becomes nonlinear in these cases, it is necessary to use the coupled TLM equations, Eqs (5) and (6), and the localized tunneling equation, Eq. (10), to provide more accurate predictions.
The total contact resistance R c of the CNT-vacuum-CNT parallel contact is shown in Fig. 8, as functions of contact length L, vacuum gap distance D, and applied voltage V 0 . The total contact resistance R c increases very rapidly with increasing insulating layer thickness, D, and decreases with contact length, L. For the low applied voltage regime (V 0 < 0.3 V), R c is almost independent of V 0 , as shown in Fig. 8c,d. When the applied voltage V 0 > 0.3 V, R c decreases sharply with V 0 . This is because the junction is no longer ohmic and the tunneling resistivity ρ c decreases nonlinearly with the junction voltage, as a function of position along the contact length. Ohmic approximations (Eqs 8, 11) fail to give accurate results in the latter case and it is necessary to use the self-consistent numerical model. As L increases, the dependence of contact resistance on L becomes less significant. Similar profiles of contact resistance with L were observed in other experimental and theoretical works 10, 22,31,42,44 . The contact resistance lies between 5 kΩ to 10 MΩ for the cases shown in Fig. 8, which agrees with previously reported experimental and theoretical works 24,26,29 . The existing 1D models give an inaccurate estimation of the contact resistance because they do not consider the variation of tunneling current density along the contact length.
Next, we extend our calculations for contacts of CNT with different metals -calcium (Ca), aluminum (Al), copper (Cu) and gold (Au). The work functions of Ca, Al, Cu and Au are taken as 2.9, 4.08, 4.7 and 5.1 eV respectively 56 . The work functions and dimensions of the CNT are kept same as before. In addition, the dimensions of the CNT and contacting-metal-2 are assumed to be same (width of 3 nm, thickness of 3 nm) for the simplicity of calculations. The resistivity of Ca, Al, Cu and Au are known to be 3.36 × 10 −8 Ωm, 2.7 × 10 −8 Ωm, 1.68 × 10 −8 Ωm and 2.2 × 10 −8 Ωm respectively 56,57 . All the material properties are specified in the main text. Solid lines are for self-consistent numerical calculations using Eqs 5, 6 and 10, dashed and dotted lines are for analytical calculations from Eq. 7, with ρ c calculated using V g = V 0 in Eq. 10 and using ohmic approximations for the tunneling junction, Eq. 11, in the limit of V g → 0, respectively. Figure 9 shows the effects of the work function of contacting member 2 (W 2 ) and the permittivity of the thin insulating layer (ε r ), on the current and voltage characteristics in CNT-insulator-metal contacts. As the two contact members are different, the voltage drop V g (x) and the tunneling current density J c (x) are no longer symmetric along the contact length L. Figure 9(a,b) show that the voltage drop increases and the tunneling current density decreases with increasing W 2 . Figure 9(c,d) show that the voltage drop increases and the tunneling current density reduces significantly when the relative permittivity of the insulating layer increases from 1 to 3.9. Analytical solutions obtained by assuming constant tunneling resistivity along the contact length are also included, similar to the previous cases of Fig. 7. In general, for the chosen value of D = 0.5 nm, the ohmic approximations using Eq. 11 do not yield accurate results. The constant tunneling resistivity approximation using Eq. 10 by setting V g = V 0 could be a good approximation for the self-consistent TLM model (Eqs 5, 6 and 10), for tunneling layers with higher permittivity ε r . Figure 10 shows the contact resistance (in Ω) for various contact metals and tunneling films for CNT-insulator-metal contacts. Contact resistance increases with insulating layer thickness D, insulating layer permittivity ε r and work function of contacting member W 2 . It decreases with contact length L, as in the similar contacts in Fig. 8. The potential barrier in the insulating layer increases with the increase of work function of the contact metal, resulting in lower tunneling current and higher contact resistance.

Summary
In this paper, we proposed a self-consistent model to characterize partially overlapped parallel contacts. Our model considers the spatial variation of contact resistivity along the contact structure. We solved the TLM equations for three cases: (1) constant specific contact resistivity, (2) linearly varying specific contact resistivity, and (3) spatial dependent specific contact resistivity along the contact length due to current tunneling. The analytical solutions for the first case have been verified with FEM based numerical simulations. Our study provides a Solid lines are for self-consistent numerical calculations using Eqs 5, 6 and 10, dashed and dotted lines are for analytical calculations from Eq. 8, with ρ c calculated using V g = V 0 in Eq. 10 and using ohmic approximations for the tunneling junction, Eq. 11, in the limit of V g → 0, respectively.
www.nature.com/scientificreports www.nature.com/scientificreports/ thorough understanding of the contact tunneling resistance, current and voltage distributions across nano and sub-nano scale metal-insulator-metal (MIM) junctions in parallel electrical contacts. The effects of contact geometry (i.e. dimension of the contact, and distance between the contact electrodes), and material properties (i.e. work function, sheet resistance of the contact members, and permittivity of the insulating layer) on the spatial distributions of currents and voltages across these contacts, and the overall contact resistance are studied in detail. While predominately classical in nature, the inclusion of tunneling current starts to address quantum effects in these small scale objects.
It is found that in general the ohmic approximation of tunneling junctions (Eq. 11) is not reliable for predicting the contact resistance of parallel tunneling contacts. The one-dimensional (1D) tunneling junction models (Eq. 10 with constant voltage across the whole junction) are good approximations of the parallel contacts only when the thickness D or the permittivity ε r of the tunneling film is relatively large, or the applied voltage across the contact V 0 is relatively small. When the 1D models become unreliable for small D or ε r , or large V 0 , the self-consistent TLM equations coupled with the tunneling current (Eqs 5, 6 and 10) need to be used to accurately characterize the parallel tunneling contacts.
The parallel tunneling contact in this work may be considered as the basic building block to better understand the macroscopic electrical conductivity of CNT fibers, which contains a very large number of such parallel contacts between individual CNTs. Furthermore, our study elucidates key parameters for parallel electrical contacts over a wide range of spatially dependent contact resistivity, which paves the way to strategically design of contact structures with controlled current distribution profiles and contact resistance, by spatially varying the contact layer properties and geometry.
In this formulation, we have ignored the effects of space charge and exchange-correlation inside the tunneling gap 11,40 . We have also ignored possible charge trapping inside contact junctions. The model is assumed two-dimensional, where the effects of the transverse dimension are neglected. These issues will be the subjects of future studies. It is important to note that the transmission line model (TLM) is only a simplified approximation in CNT-insulator-Al contacts, for different insulating layer permittivity ε r , with fixed D = 0.5 nm, L =50 nm, V 0 = 3V. Solid lines are for self-consistent numerical calculations using Eqs 5, 6 and 10, dashed and dotted lines are for analytical calculations from Eq. 7, with ρ c calculated using V g = V 0 in Eq. 10 and using ohmic approximations for the tunneling junction, Eq. 11, in the limit of V g → 0, respectively.
www.nature.com/scientificreports www.nature.com/scientificreports/ of the 2D electrical contacts, where the current crowding and the fringing fields near the contact corners cannot be fully accounted for 10, 49 . In order to accurately evaluate these effects as well as the impact of finite thickness in the contact members and the contact junction, and the possible parallel component of current flows in the interface layer, field solution methods need to be used 10, 49,50,58 . Future studies may also include the effects of various contact geometry, insulator layer non-uniformities and AC response on the electrical properties of tunneling type contacts. An additional feature might include the role of capacitance and inductance in nano-and micro-scale structures, especially when a large contact resistance coupled with these reactive effects might introduce new time scales into time-dependent dynamic problems.  Figure 10. The total contact resistance R c of the CNT-insulator-metal parallel contact. Contact resistance is plotted as a function of (a) contact length L, (b) insulator layer thickness D and (c) insulator layer permittivity ε r , for CNT-insulator-metal contacts for different contacting metals (Ca, Al, Cu, Au). (d) Contact resistance as a function of work function of contacting member 2 (W 2 ). The material properties and dimensions for (a-c) are specified in the text (the same as in Fig. 9). For (d), the resistivity of contacting member 2 is assumed to be 2.0 × 10 −8 Ωm. The results are from the self-consistent numerical calculations using Eqs 5, 6 and 10.