Passive Intermodulation at Contacts of Rough Conductors: Comparison
Please note this is a comparison between Version 2 by Lindsay Dong and Version 9 by Lindsay Dong.

Passive intermodulation (PIM) is a nonlinear phenomenon causing frequency mixing and harmonic generation. PIM products, generated by weak nonlinearities of passive devices, interfere with information signals and degrade the performance of communications, navigation and astronomical systems. The passive nonlinearities of different types: electrical, thermal and mechanical are coupled intrinsically despite their different time scales. Roughness of the contact surfaces notably enhances the PIM effect at conductor joints. 

  • passive intermodulation (PIM)
  • nonlinearity
  • multiphysics effects
  • electrical contacts
  • surface roughness

1. Introduction

The continuously growing volume and speed of data transmission pose major challenges to the performance of existing and future wireless and satellite communications and navigation systems [1][2][3][4][5][6][7]. The stringent requirements of the integrity of information signals push the limits of radio frequency (RF) hardware. Weak nonlinearities of passive devices such as antennas, filters, couplers and multiplexers at the RF front-end of the smart multi-radio base stations and user terminals produce spurious emission, corrupt information signals and debilitate the system performance [8][9][10]. Therefore, efficient mitigation of nonlinear signal distortions is a major requirement to reconfigurable antennas, dynamically adjustable RF front-ends and their individual passive components [11][12][13][14][15][16].

Constituent materials and contacts in RF devices proved to be the main sources of passive nonlinearities [17][18][19]. The state-of-the-art RF materials have low loss, high thermal conductivity and good mechanical properties. But when exposed to the high power of RF signals, they exhibit weakly nonlinear behaviour and generate frequency harmonics and passive intermodulation (PIM) products. Figure 1 illustrates that the 3rd order PIM products, PIM3L­ and PIM3U, are close to the two carrier frequencies f1 and f2, and PIM3L­ falls in the receive band Rx. In “linear” passive devices, PIM is usually the result of mixing high-power electromagnetic signals by weak nonlinearities of good conductors, their contacts and surface finish. The PIM effects manifest themselves in spurious emission, extra noise and nonlinear distortion of the information signals. PIM products are particularly harmful to radars, wireless and space communications systems and radio astronomical instruments [5][6][7][20][21].

Figure 1. The two high power carriers with frequencies f1 and f2 in the transmit band Tx generate PIM products PIM3L and PIM3U with frequencies 2f1-f and 2f2-f1, respectively. The frequency of PIM3L falls in the receive band Rx and causes nonlinear distortions of the received weak information signals.

PIM phenomenology has been studied for more than 50 years but still remains a nagging problem due to its complex multiphysics nature. The basic mechanisms of nonlinearities and PIM generation were explored in metal contacts [22][23][24][25][26][27][28][29], printed RF transmission lines [30][31][32][33][34][35][36], cable assemblies [17][18][31][37][38][39][40], coaxial connectors [41][42][43]  and antennas [21][44][45][46][47][48]. The main sources of passive nonlinearities in conductor joints and their interlinks are summarised in Figure 2. They include Metal–Insulator–Metal (MIM) junctions [22][23][24][26][27][28][29][22,23,24,26,27,28,29], electro-thermal effects  [37][40][37,40,[49][50], surface roughness  [32][32][51][52] and mechanical deformations [52][52][53][54][55][56]. The electrical contacts with rough surfaces and their RF performance have been recently studied in micro-electro-mechanical systems (MEMS). Whilst the mechanical properties of MEMS have been explored in great detail (see [53][54][55][56] [53,54,55,56] and references therein) the existing models remain limited to the linear devices exposed to weak RF signals. The effects of high RF power on the contacts and junctions of conductors with rough surfaces were studied only in coaxial connectors [38][39][40] [38,39,40] and waveguide flanges [57]. But the practical means of mitigating nonlinear distortions and PIM in passive RF circuits remain mostly semi-empirical [10][10][58].


Figure 2. Main types of contact nonlinearities of rough conductor surfaces. Links between different mechanisms of nonlinearities are colour coded in relation to the primary source.

2. Electrical Nonlinearities of Conductor Joints

PIM at contacts of conductors with rough surfaces is a nonlinear multiphysics process with several distinct time scales. The fastest nonlinearities in contacts of good conductors are associated with the electrical effects of charge tunnelling and current constriction at the MIM junctions [17]. The charges are funnelled through the contact spots of rough surfaces as shown in Figure 3(a). The current magnitude depends on the size of a contact area and thickness of an insulating layer that are determined by the applied pressure, local temperature and deformations of the surface asperities. The contact nonlinearities mix and modulate the high-power carriers due to variations of the contact size, temperature and resistance [27][38][39][40][27,38,39,40]. The thermal effects and mechanical deformations develop much slower than the tunnelling current and the oscillations of the RF carriers. However, both fast and slow nonlinearities remain intrinsically coupled to each other despite their different time scales.

Current flow through the contact asperities of rough surfaces is described by the equivalent circuit shown in Figure 3(b), where Rv1 and Rv2 are the resistances of conductors outside of the asperities. The resistance and capacitance of a pair of compressed or non-touching asperities are represented by Rcm and Ccm. The linear resistance Rcm of the thin oxide film is very high even in contacts of good conductors. For example, Al2O3 film of thickness δs = 1 nm, area 10 μm2 and the macroscopic resistivity ρAl2O3 = 1 × 1014 Ω·cm has its resistance Rcm = 1014 Ω. Such a high value of Rcm is typical for good insulators and suggests that the charge tunnelling must be responsible for the measured resistance of conductor junctions with thin insulating films. It is also necessary to note that the capacitive reactance 1/jωCcm of the conductor joint can noticeably affect the contact impedance. Indeed, the reactance of this junction is 1/jωCcm = −j89.9 Ω alters the impedance of the contact joint.

Figure 3. (a) Sketch of contacts between conductors with rough surfaces. Tips of the collided asperities are deformed and the constriction current is funnelled through the contact areas. Only the displacement current (broken lines) is shown at the non-contact asperities. (b) Equivalent circuit of an asperity pair: Rv1 and Rv2 are the resistances of bulk conductors outside asperities; Rcm and Ccm are the resistance and capacitance of either the collided asperities or a pair of isolated asperities.

2.1. Contact and Constriction Resistances

Contacts of conductors with rough surfaces contain the MIM and metal-to-metal junctions. Their resistance depends on the number of touching asperities, size of each contact spot and thickness of an insulating layer in MIM junctions [59]. The mechanical and thermal deformations of individual asperities determine an overall size of the contact area and its contact resistance. The resistance Rc of a contact spot with an equivalent radius a is usually approximated as [60].

R c = ρ c 2 a [ f ( λ a ) + 8 λ 3 π a ]



is an average electrical resistivity of a pair of contact asperities,


is an electron-free path and






) is an interpolation function describing a contribution of Maxwell resistance






. A good approximation of






) proposed in [60] has a maximum error less than 1% at any




f ( λ / a ) = 1 + 0.83 · λ / a 1 + 1.33 · λ / a

The values of f(λ/a) vary in a relatively narrow range between f(λ/a) ≈ 0.624 at λ >> a and f(λ/a) ≈ 1 at λ << a. The second term in (1) represents Sharvin resistance [61]: R S = 4 ρ c λ 3 π a 2 that is associated with the collision-free motion of charges. It plays a notable role when the size a of a contact spot is smaller than λ. Then Sharvin resistance RS can exceed Maxwell resistance RM and its contribution to the constriction current becomes significant.

Current constriction by asperities is an inherent feature of the contacts of conductors with rough surfaces [53,57][53][57][62]. The constriction resistance is determined by the number of the compressed asperities, sizes of their contact spots and thicknesses of the insulating layer. The constriction current decreases when the size of the contact spot is larger than the mean free path of electrons, a > λ, as evident from (2). Then the charge transport is predominantly diffusive and is dominated by Maxwell resistance RM, defined in (1). In high-quality conductor contacts, e.g., in MEMS, the constriction current is much smaller than the conduction current and is usually combined with the conduction or tunnelling currents [63].

The effective resistance Ru of the whole junction surface containing M asperities can be averaged and approximated as suggested in [53].[53]

R u = ( m = 1 M R c m 1 ) 1 ρ a v 2 a e f f [ f ( λ a e f f ) + 8 λ 3 π a e f f ]

where M is the number of compressed asperity pairs, and Rcm is defined by (1) and describes the contact resistance of mth pair of asperities. Thus, Ru is evaluated with (1) where a and ρ are replaced by an effective radius aeff and average resistivity ρav [53]. Such an approximation proved to be fairly accurate when the contact area is large and aeff > λ. Then the contribution of Sharvin resistance is small and deviations of resistivities ρcm from their average value ρav become negligible.

The contact resistances Rcm of individual contact pairs and the whole ensemble, Ru, are the main electrical parameters used for characterising the performance of MEMS switches. They include the effects of both the electrical properties of the contact surfaces and the mechanical deformations. The electro–mechanical models of the contacts in the MEMS switches have been discussed in detail in [53][54][55][63][53,54,55,63].

2.2. Nonlinearity of MIM Contacts

The tunnelling current in MIM contacts subjected to RF power is the main source of their fast nonlinearity [22][23][24][22,23,24]. When a thin insulating film of a few nanometre thick separates asperities of rough conductors, the tunnelling current flows through. An efficient interaction of the high-power RF carriers with free charges is the result of the small timescale of the tunnelling effects in the conductor junctions. However, when the thickness of the insulating layer is commensurate with the free path of electrons, charges cannot penetrate the potential barrier. The tunnelling current rapidly decays and becomes practically negligible as the insulator between asperities is thicker than 5-10 nm.

In the widely accepted Simmons model, the current density J(Vg) in MIM junction is represented as

J ( V g ) = J 0 [ φ ¯ e - A φ ¯ - ( φ ¯ + e V g ) e - A φ ¯ + e V g ]



J 0 = e 2 π h ( β Δ s ) 2

and is

A = 4 π β Δ s h 2 m


φ ¯

an average height of potential barrier inside an insulator layer,


is the voltage between the contact conductors,


is the electron charge,


is the Planck constant,


is the electron mass,


is the correction factor used in the approximation of Δ


– an effective thickness of insulator film which is usually smaller than the actual thickness. The value of the parameter


is typically close to 1, and at the insulator thicknesses Δ


~4–5 nm, the error of the approximation

β = 1 is less than a few percent [23].

= 1 is less than a few percent [23].

2.3. The Effect of Surface Roughness

Surface roughness significantly affects the performance of electrical contacts in high frequency applications [52][53][54][55][57][52,53,54,55,57]. The statistical models have been used for simulating the contacts of rough surfaces. Gaussian distribution of contact asperity heights with the standard deviation up to 20% from an average value was examined in [64], taking into account the effect of an insulator thickness on the contact resistance of rough surfaces. Using the Brinkman–Dynes–Rowell model [65], it has been shown that a single thickness model reasonably estimates the surface conductance when the standard deviation of asperity heights remains within 5% of the average value.

2.4. Thermionic Emission

Thermionic emission is the process of electron discharge from the free metal surface. In good conductors, it occurs at very high temperatures, above 1000°C [66], when hot electrons can gain enough energy to break their bonding and escape from the conductor surface.

At operational temperatures of electronic devices, the thermionic emission current in joints of good conductors is practically negligible. It also decreases exponentially with the thickness of the insulating layer and becomes undetectable as the insulator layer is thicker than 5 nm. Thus, the thermionic emission current can be neglected in contacts of good conductors.

3. Thermal Nonlinearities at Contacts of Good Conductors

3.1. Self-Heating Effect and Nonlinearity of Contact Resistivity

Heat generation is an inherent feature of electromagnetic (EM) wave interactions with conductors and their contact joints. This multiphysics process couples the electric and thermal domains as illustrated in Figure 4. Namely, the dissipative losses of high-power RF signals generate heat which, in turn, alters the resistance of conductors and causes their thermal expansion and mechanical deformations [67]. These nonlinear effects are interlinked despite the time scales of the thermal and mechanical processes notably differ.

Figure 4. Compact model of a pair of the compressed asperities where the contact resistance Rcm is coupled to a thermal node with the thermal resistance Rth and thermal capacitance Cth.

3.2. Electro-Thermal Effect in Contact Joints and PIM Generation

Electrothermal PIM (ET–PIM) is caused by self-heating due to conductor and dielectric losses [37]. ET–PIM has a distinct signature of the nonlinear coupling of electrical and thermal domains as illustrated in Figure 4. The basic mechanism of this nonlinearity is the thermal modulation of resistivity. In essence, the heat due to RF losses alters the resistance of a conductor, which is heated by the high-power carriers, and this results in the generation of PIM products. The effects of ET–PIM have been studied in the termination resistors [37][50][37,50], printed transmission lines [36][3667][68] and thin-film coplanar waveguides with spatially inhomogeneous current distributions [35]. The developed theory of the ET–PIM and supporting experiments [36][37][68][36,37,68] have revealed that the baseband resistivity of conductors is modulated by the heat oscillations. The effect of resistivity variation on the skin depth due to modulation of the RF carriers was analysed in [45][46] [45,46]. These studies have enabled the development of the qualitative analytical model, which sheds light on the principal mechanisms of ET–PIM generation and provides an accurate assessment of signal distortion by ET–PIM in the full-wave EM simulations [33][34][35][36][68][33,34,35,36,68].

Evaluations of the dissipative losses and the rate of self-heating are the critical steps in the ET–PIM analysis. While the skin effect is routinely accounted for in the EM simulations when calculating the losses of imperfect conductors, surface roughness is normally ignored despite its proven major impact on the performance of printed circuits [51][52][51,52] and MEMS switches [53][54][55][53,54,55].

3.3. Thermal Expansion of the Contact Area

The heat, generated by the high RF power, is dissipated in conductors and their joints and causes the thermal expansion and deformation of contact surfaces. As the result, the contact areas and the conductor resistivity ρc(x) increase with temperature. The contact resistance Rcm of asperities also grows with temperature, whilst the pace of its growth depends on the relation between the temperature coefficient of resistivity α, the rate of the thermal expansion of the contact areas and asperity deformations. These thermal effects cause nonlinearity of the contact resistance that inflict [69]nonlinear signal distortions.

4. Mechanical Nonlinearities at Conductor Contacts

4.1. Effect of Asperity Deformations on Conductor Resistivity