Contact Properties of a Microelectromechanical (MEMS) switch with gold-on-gold contacts. Sumit Majumder ECE, Northeastern University

Similar documents
Chapter 3 Contact Resistance Model with Adhesion between Contact

Chapter 2 A Simple, Clean-Metal Contact Resistance Model

Chapter 1 Introduction

CONTACT MODELING OF AN RF MEMS SWITCH

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT

Unloading of an elastic plastic loaded spherical contact

ADHESION OF AN AXISYMMETRIC ELASTIC BODY: RANGES OF VALIDITY OF MONOMIAL APPROXIMATIONS AND A TRANSITION MODEL

A Micro Model for Elasto-Plastic Adhesive-Contact in Micro-Switches: Application to cyclic loading

A General Equation for Fitting Contact Area and Friction vs Load Measurements

Contact Modeling of Rough Surfaces. Robert L. Jackson Mechanical Engineering Department Auburn University

ADHESIVE ROUGH SURFACE CONTACT. by Clint Adam Morrow B.S., University of Kentucky, 1996 M.S., University of Kentucky, 1999

! Importance of Particle Adhesion! History of Particle Adhesion! Method of measurement of Adhesion! Adhesion Induced Deformation

A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius

A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS

Module-4. Mechanical Properties of Metals

Contact, Adhesion and Rupture of Elastic Solids

Supporting Information. Interfacial Shear Strength of Multilayer Graphene Oxide Films

Nano-Scale Effect in Adhesive Friction of Sliding Rough Surfaces

EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY

Deterministic repeated contact of rough surfaces

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study

Review of Thermal Joint Resistance Models for Non-Conforming Rough Surfaces in a Vacuum

Figure 43. Some common mechanical systems involving contact.

Contents. Preface XI Symbols and Abbreviations XIII. 1 Introduction 1

Roughness picture of friction in dry nanoscale contacts

A Short Course on Nanoscale Mechanical Characterization: The Theory and Practice of Contact Probe Techniques

Numerical modeling of sliding contact

CONTACT MODEL FOR A ROUGH SURFACE

Contact Mechanics Modeling of Homogeneous and Layered Elastic-Plastic Media: Surface Roughness and Adhesion Effects

Normal contact and friction of rubber with model randomly rough surfaces

The Frictional Regime

An analysis of elasto-plastic sliding spherical asperity interaction

Elastic-plastic Contact of a Deformable Sphere Against a Rigid Flat for Varying Material Properties Under Full Stick Contact Condition

Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

A Finite Element Study of the Residual Stress and Deformation in Hemispherical Contacts

Abstract. 1 Introduction

Supplementary Figures

A Finite Element Study of an Elastic- Plastic Axisymmetric Sinusoidal Surface Asperity in Contact Against a Rigid Flat with Strain Hardening

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

Effect of Tabor parameter on hysteresis losses during adhesive contact

FCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

Chapter 4 Normal Contact with Adhesion

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

Experimental Investigation of Fully Plastic Contact of a Sphere Against a Hard Flat

Chapter 7. Highlights:

Introduction to Engineering Materials ENGR2000. Dr. Coates

An Accurate Model for Pull-in Voltage of Circular Diaphragm Capacitive Micromachined Ultrasonic Transducers (CMUT)

Mechanical Properties of Materials

4.MECHANICAL PROPERTIES OF MATERIALS

A van der Waals Force-Based Adhesion Model for Micromanipulation

Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.

Experimental research on electrical resistance of microcontacts

MATERIALS. Why do things break? Why are some materials stronger than others? Why is steel tough? Why is glass brittle?

Module17: Intermolecular Force between Surfaces and Particles. Lecture 23: Intermolecular Force between Surfaces and Particles

INTERACTIONS BETWEEN FINE PARTICLES

For an imposed stress history consisting of a rapidly applied step-function jump in

Johns Hopkins University What is Engineering? M. Karweit MATERIALS

Princeton University. Adhsion and Interfacial Failure in Drug Eluting Stents

A detailed 3D finite element analysis of the peeling behavior of a gecko spatula

Characterization of MEMS Devices

Lecture #2: Split Hopkinson Bar Systems

Extended hard-sphere model and collisions of cohesive particles

Modelling particle-particle and particle-fluid interactions in DEM

Basic Laboratory. Materials Science and Engineering. Atomic Force Microscopy (AFM)

Forces Acting on Particle

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

Outline. 4 Mechanical Sensors Introduction General Mechanical properties Piezoresistivity Piezoresistive Sensors Capacitive sensors Applications

Measurement of the Deformation and Adhesion of Rough Solids in Contact

Chapter 6: Mechanical Properties of Metals. Dr. Feras Fraige

How materials work. Compression Tension Bending Torsion

ME 2570 MECHANICS OF MATERIALS

Background and Overarching Goals of the PRISM Center Jayathi Murthy Purdue University. Annual Review October 25 and 26, 2010

Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment. P. Krulevitch G. C. Johnson

SUPPLEMENTARY INFORMATION

Nanomechanics Measurements and Standards at NIST

Mechanics of Earthquakes and Faulting

A COMPACT MODEL FOR SPHERICAL ROUGH CONTACTS

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture,

Strain Gages. Approximate Elastic Constants (from University Physics, Sears Zemansky, and Young, Reading, MA, Shear Modulus, (S) N/m 2

Supplementary Figures

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

ScienceDirect. A Numerical Model for Wet-Adhesive Line Contact

MECHANICS OF MATERIALS

The Mechanics of CMP and Post-CMP Cleaning

Sinan Müftü Associate Professor Department of Mechanical Engineering Northeastern University, 334 SN Boston, MA

arxiv: v1 [cond-mat.soft] 10 Oct 2018

Adhesive elastic contacts JKR and more

A statistical model of elasto-plastic asperity contact between rough surfaces

Mechanics of Earthquakes and Faulting

A finite element study of the deformations, forces, stress formations, and energy losses in sliding cylindrical contacts

PARTICLE ADHESION AND REMOVAL IN POST-CMP APPLICATIONS

ME 383S Bryant February 17, 2006 CONTACT. Mechanical interaction of bodies via surfaces

Study of electrostatic potential induced by friction

ELASTICITY (MDM 10203)

Finite-Element Analysis of Stress Concentration in ASTM D 638 Tension Specimens

Thermal Contact Resistance of Nonconforming Rough Surfaces, Part 1: Contact Mechanics Model

MAS.836 PROBLEM SET THREE

SCALING OF THE ADHESION BETWEEN PARTICLES AND SURFACES FROM MICRON-SCALE TO THE NANOMETER SCALE FOR PHOTOMASK CLEANING APPLICATIONS

Transcription:

Contact Properties of a Microelectromechanical (MEMS) switch with gold-on-gold contacts Sumit Majumder ECE, Northeastern University

MEMS switches - background Advantages size, flexible configuration, integration capability zero steady-state power consumption low leakage, high isolation compared to FET switches Applications Pin diode replacement in RF applications (no steady-state power, easy scalability to multi-pole, multi-throw configurations) Reed relay replacement in various applications (ultrasound, ATE)

MEMS switches - actuation mechanisms Electrostatic Vs electromagnetic Electrostatic switches are easier to build, smaller and faster, do not draw steady-state power Electromagnetic switches have higher actuation force (more reliable contact), easier to build 4-terminal relay

Why study contact adhesion? F = εav x( V ) x(v) V A = 0.005-0.05 mm, x(v) = 0.5-1 µm, V less than 100 V Total actuator force typically ~ 1 mn Contact force 100 µn - 1 mn, restoring force 100 µn - 1 mn

Why study contact adhesion? To obtain low and stable contact resistance, require a minimum clean, metallic contact area => minimum contact force minimum restoring force to overcome contact adhesion Actuator force 1 order higher in MEMS electromagnetic relays, orders higher in reed relays Much more difficult to meet contact force and restoring force constraints Objective of this thesis: increase our understanding of adhesion in electrostatically actuated microswitch contacts

Evolution of measured characteristics 1.4 1.4 Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. Initial Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. 10 cycles 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) Contact Resistance (Ohms) 1.4 1. 1.0 0.8 0.6 0.4 0. Contact Resistance (Ohms) 1.4 1000 cycles 1. 10000 cycles 1.0 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn)

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Northeastern University microswitch Predecessor: Zavracky & Morrison, Foxboro (1984). Devices used in this work fabricated in Microfabrication Laboratory (MFL) at NU in 1999. Since 1999, further development and commercialization at Analog Devices and Radant MEMS. Fabricated using all-metal surface micromachining process developed at NU.

Northeastern University microswitch Contact Detail Gate G Drain Beam Source S DRAIN GATE SOURCE D 100 µm DRAIN GATE SOURCE Gold cantilever beam, sputtered gold contacts.

Northeastern University microswitch with gold-gold contacts Φ= 0 75 15 35 Threshold voltage 65-70 V. Restoring force 70 µn per contact. Actuation voltage 90-10 V => contact force 70-180 µn per contact. Operated in dry nitrogen. 0.3 0.6 DRAIN GATE 6 SOURCE Contact resistance initially 0.5-1 Ω, after 10 3-10 4 cycles 0.1-0. Ω. Lifetime ~10 6 cycles with 10 ma per contact, limited by permanent stuck-closed failures.

MEMS switch contacts: Related work Majumder et al. (1997,1998,001) Hyman et al. (1998) thermal model Yan et al.(001) thermal model Kruglick (1999), Lafontan (001) simple contact models from device perspective AFM-based work in MEMS switch force range Beale and Pease (199) Schimkat (1998), Bromley et al. (001) MEMS switches First switch Petersen (1979). Many switches since 1995 (majority electrostatic, some electromagnetic and thermal) Large body of work on contact mechanics and adhesion (cited while developing models)

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Model of contact surface d z 4 z z 3 z 1 a 1 a Contact bump surface much rougher than drain surface - drain assumed to be a smooth surface. Greenwood and Williamson (1966) - surface represented by asperities of known shape and varying heights n asperities of radius R, heights z 1,z,... Deformation of asperities under load => contact spots are formed. Asperities assumed to deform independently, total contact force F n = i= 1 F i

GW (asperity-based) model versus fractal model Real surfaces have roughness at different length scales, GW model represents roughness at a single length scale (represented by asperity radius) Fractal-based contact models capture roughness at different length-scales (Bhushan and Majumdar (1990)), but: more complicated than GW model requires more accurate surface profile than GW model published models do not include adhesion

Selecting parameter values for GW model d z 4 z z 3 z 1 a 1 a Surface described by Radius of curvature (R) Number of asperities (n) distribution of heights - Gaussian distribution => standard deviation σ SEM of drain shows 10-50 distinct spots, ~10 nm size. R=0.01 µm too small for this spot size. Assume R=0.1 µm, later vary R between 0.1-1 µm. Distribution of heights: assume 100 asperities, σ=0.1 µm (few contact spots) and σ=0.01 µm (10 s of contact spots). Later vary σ between 0.003 µm and 0.1 µm.

Mechanical properties Nanoindentation measurements at Hysitron, MN Elastic modulus E=91 GPa Hardness H=. GPa Poisson s ratio ν=0.5 (assumed; reported as 0.44 for bulk gold)

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Single asperity without adhesion - elastic/perfectly plastic R F=0 F=F 1 Effective elastic modulus Elastic deformation (Hertz): contact radius a and vertical deformation α Plastic deformation: a and α α modified to obtain geometric continuity with elastic model (Chang, Etsion and Bogy (1988)) Elastic to perfectly plastic transition 1 K = 3 1 ν1 1 ν ( + ) 4 E E 1 FR 1/3 a = ( ) K F = πa a = H RY a p = 11 E α = α = α p Rα( ) α a R a R

Single asperity without adhesion - elasto-plastic model Physically realistic transition from elastic to perfectly plastic Elastic (Hertz) FR 1/3 a = ( ) K α = a R Fully plastic F = πa H Elastic to elastoplastic - from von Mises yield criterion Elastoplastic to fully plastic Semi-empirical formula in elastoplastic (Studman (1976), modified after Maugis and Pollock (1984)) F 1.1Y. = πa a = a p 60 F πa RY E = Y (1.1 + 0.7 ln Ea ) 3.9YR

Single asperity without adhesion - contact force-contact radius 10-7 Contact radius (m) 10-8 elastic-fully plastic model, elastic elastic-fully plastic model, plastic modified Studman model, elastic modified Studman model, elasto-plastic modified Studman model, plastic R=0.1 µm 10-9 10-8 10-7 10-6 10-5 10-4 Force (N) Elastic to perfectly plastic: F=0.55 µn, a=8.9 nm Elastic to elasto-plastic: F=0.05 µn, a=3.1 nm, elastoplastic to perfectly plastic F=16 µn, a=49 nm

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Single asperity model with adhesion Physical origins of contact adhesion Elastic contact models with adhesion JKR DMT Elasto-plastic and perfectly plastic models with adhesion Unloading and separation of contact surfaces; adherence force

Physical origins of adhesion Various surface forces act between two surfaces in contact or proximity - van der Waals, ionic, covalent, metallic. Interaction characterized by a potential U(z). When surfaces are in contact, U(z) reaches a (negative) minimum. Defined as adhesion energy w (energy required to separate unit areas of two surfaces in contact). A surface is characterized by surface energy ν (energy required to create unit area of the surface). w=ν for two identical surfaces

Van der Waals vs metallic bonds van der Waals forces long-range (several nm) w~50 mj/m A U ( z) = 1πZ 0 Z ( z Lennard-Jones 0 ) 1 4 Z ( z 0 ) 8 Metallic bonds short range (characteristic length ~ 0.5 A 0 ) w=-.5 J/m for gold U ( z ) = U (1 + βz ) e * * βz 0 Ferrante et al. (198) * Metallic bonds can be screened by adsorbed molecules, films.

Adhesion and elastic stresses are part of the same continuum 1.5 1.5 1.0 1.0 U/U 0 0.5 σ/σ th 0.5 0.0 0.0-0.5-4 - 0 4 6 8 10 1 14 z* -0.5-0 4 6 8 10 1 z* Equilibrium at z*=0. Stress-strain is linear only when z*~0. Compressive stress when z*<0 Tensile stress/adhesion when z*>0

Adhesive elastic contact and fracture mechanics Elastic loading and unloading of sphereflat contact is analogous to retreat and advance of a crack ( external circular crack ) Crack advances by creating new surface area - requires energy wda. Crack advances if -du E >wda. In equilibrium Contact area under load, and contact adherence force can be calculated by invoking this equilibrium ( U / A) w G E = = δ crack

Single asperity with adhesion - JKR model σ σ Hertz model Compressive Force F 1 > F a r Flat punch model Force F -F 1 Boussinesq (1885) Sneddon (1946) a r σ JKR model Force F Johnson, Kendall and Roberts (1971) a r Physically unrealistic solution at contact boundary, due to flat-punch model assumption of abrupt contact boundary

Single asperity with adhesion - JKR model Contact radius a: 3 a K = R ( 3π ) F + 3π wr ± 6πwRF + wr Vertical deformation α: α = a R 8πaw 3K Contact adherence force is the tensile force at which the equilibrium G=w becomes unstable 3 F adh = πwr 3 πwr amin = K 1/.3

Single asperity with adhesion - DMT model Derjaguin, Muller and Toporov (1975). Adhesion forces assumed to exist in a ring surrounding the contact spot, but does not change Hertz solution for deformed surface profiles, stress within contact spot. Solution is not self-consistent. Contact adherence force: = πwr F adh

Validity of JKR and DMT models JKR DMT Height of neck in JKR model at small loads: Characteristic length for adhesive forces λ~0.5 A 0 R=0.1-1 µm, w=.5 J/m, K=81 GPa => h>> λ Therefore forces outside contact spot can be neglected => JKR is valid. h 1/3 Rw K

Single asperity with adhesion - elasto-plastic / plastic models Roy Chowdhury and Pollock (1981), Maugis (1984) Adhesion forces assumed to exist in a ring surrounding the contact spot, but does not change adhesionless (elasto-plastic/perfectly plastic) solution for deformed surface profiles, stress within contact spot. Solution is not self-consistent. Contact radius in elasto-plastic regime: F = πa pm πrw Ea p m = Y (1.1 + 0.7 ln ) 3.9YR Contact radius in plastic regime: F = πa H πrw Vertical deformation modified to obtain geometric continuity with JKR model a α = R αe + 1 8πa ew 3K

Single asperity with adhesion - contact force-contact radius Contact radius (m) 10-7 No adhesion, elastic No adhesion, elasto-plastic No adhesion, fully plastic w=.5 J/m, elasto-plastic w=.5 J/m, fully plastic 10-8 10-9 10-9 10-8 10-7 10-6 10-5 10-4 Force (N)

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Single asperity with adhesion - recovered asperity radius after elasto-plastic / plastic deformation Elasto-plastic/plastic deformation after loading to F f (no adhesion): flattening of asperity (radius R to R 1 ), concavity in drain radius=-r. Equivalent to contact between flat and sphere of radius R eff, 1 R eff 1 = R 1 1 R Re-loading the contact to F f forms contact spot of same radius a f as given by earlier plastic loading: R eff = Ka F f 3 f Elasto-plastic/plastic deformation with adhesion: R eff, is assumed to be the same as in the absence of adhesion (Johnson (1977)). R eff a 3 f = πa f K p m = a f πp K m R eff a 3 f = πa f K H = a f K πh

Single asperity with adhesion - recovered asperity radius after elasto-plastic / plastic deformation Asperity radius (micron) 1.0 0.8 0.6 0.4 0. No adhesion, elastic No adhesion, elasto-plastic No adhesion, fully plastic w=.5 J/m, elasto-plastic w=.5 J/m, fully plastic 0.0 10-8 10-7 10-6 10-5 10-4 Contact Force (N)

Unloading and separation of single asperity Elastic loading: unloadinaccording to JKR model. Equilibrium becomes unstable at 3 F adh = πwr a min Elasto-plastic loading: 3 πwr = K Initial unloading follows elastic unloading of a flat punch with constant contact radius α α f 3( F Ff ) = Ka f 1/.3 Followed by one of three possible modes.

Unloading and separation of single asperity Ductile (F d ) mode: axial stress becomes tensile and equal to H - (yield point) plastic yield followed by cohesive separation a f a f 6Kw < π ( H + p 3Kw < πh m ) Brittle (F m ) mode: contact force F is such that JKR equilibrium point is reached, and a<a min - unstable equilibrium - adhesive separation. Brittle (F b ) mode: contact force F is such that JKR equilibrium point is reached, and a>a min - elastic (JKR) unloading until F=F adh and a=a min - adhesive separation. 6Kw π ( H + p a f a f m ) 3Kw > πp m 3Kw > πh < a f < 3Kw πp m

Unloading and separation of single asperity If contact interface is metallic it becomes indistinguishable from the bulk. Even if the contact separation is predicted to be brittle, the contacts may not separate along the interface (material transfer).

Unloading and separation of single asperity 10-4 Adherence Force (N) 10-5 10-6 Ductile separation Brittle separation (F m mode) Brittle separation (F b mode) 10-9 10-8 10-7 10-6 10-5 10-4 Contact force (N) Contact radius (nm) 34 3 30 8 6 4 0 18 Load Unload from Ff=4 µn, brittle separation Unload from Ff=1.9 µn, brittle separation Unload from Ff=0.7 µn, ductile separation 16 14-6 -4-0 4 6 Contact Force (µn)

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Multiple asperity model N asperities, heights z 1 > z > > z N First n asperities in contact, z n > d > z n+1. For each contacting asperity i,vertical deformation: α i = z i d. d z 4 z z 3 z 1 a 1 a Using elastic (JKR), elasto-plastic or plastic model, determine F i (α i ) and a i (α i ). n Total contact force F = F i i= 1

Asperities in contact 100 80 60 40 0 0 Multiple asperity model-loading All asperities in contact Elasto-plastic asperities Plastic asperities Contact spot radius (m) 10-6 Asperity 1 Asperity Asperity 5 10-7 Asperity 10 Asperity 50 10-8 10-7 10-6 10-5 10-4 10-3 Contact Force (N) 10-9 σ=0.01 µm 10-6 10-5 10-4 10-3 Contact Force (N) Asperities in contact 0 15 10 5 0 All asperities in contact Elasto-plastic asperities Plastic asperities σ=0.1 µm Contact spot radius (m) 10-6 Asperity 1 Asperity Asperity 5 Asperity 10 10-7 10-8 10-7 10-6 10-5 10-4 10-3 Contact Force (N) 10-9 10-6 10-5 10-4 10-3 Contact Force (N)

Multiple asperity model-unloading Number of asperities 50 40 30 0 10 Total asperities Ductile separation Brittle (F m ) separation Brittle (F b ) separation Number of asperities 10 8 6 4 0-100 -50 0 50 100 150 00 Contact Force (µn) σ=0.01 µm 0-100 -50 0 50 100 150 00 Contact Force (µn) Number of asperities 5 0 15 10 5 Total asperities Ductile separation Brittle (F m ) separation Brittle (F b ) separation σ=0.1 µm Number of asperities 10 8 6 4 0-40 -0 0 0 40 60 80 100 Contact Force (µn) 0-40 -0 0 0 40 60 80 100 Contact Force (µn)

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Contact resistance - single contact spot R ρ Sharvin resistance (contact radius < electron mean free path length) 4 l = e + 3 r r π ν ρ Constriction resistance (contact radius > electron mean free path length) (Maxwell) ρ = resistivity r = radius of contact spot l e = electron mean free path length 0.69<ν<1 (interpolation factor: Wexler) Measured ρ = 6.75x10-8 Ω-m (drain) 3.85x10-8 Ω-m (beam). Used average.

Contact resistance - multiple contact spots Lower bound: contact spots are far apart, conduct independently 1/ R con, lb / Rcon, i i = 1 Upper bound: spots are merged into a single large spot a eff = a i R con, ub 4ρl = 3πa e eff + υ( l e / a eff ) ρ a eff

Contact resistance model - with and without adhesion Contact Resistance (Ω) 3.0.5.0 1.5 1.0 0.5 w=.5j/m,lower bound w=.5j/m,upper bound no adhesion, lower bound no adhesion, upper bound 0.0 10-7 10-6 10-5 10-4 10-3 Contact Force (N) N=100, R=0.1µm and σ=0.01 µm

Contact resistance model - load and unload characteristics Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 0.8 0.6 0.4 0. Contact Resistance (Ω) 1.4 Load, lower bound 1. Load, upper bound Unload, lower bound 1.0 Unload, upper bound 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-100 0 100 00 Contact Force (µn) N=100, R=0.1µm and σ=0.01 µm

Contact resistance model - effect of σ and R on number of contacting asperities 50 50 Asperities 40 30 0 R=0.1 µm R=0.3 µm R=1 µm Asperities 40 30 0 10 σ=0.003 µm σ=0.01 µm σ=0.1 µm 10 0 0 0.001 0.01 0.1 σ (µm) 0.0 0. 0.4 0.6 0.8 1.0 1. R (µm)

Contact resistance model - effect of σ and R on contact adherence force -100-10 Adherence force (µn) -80-60 -40-0 R=0.1 µm R=0.3 µm R=1 µm Adherence force (µn) -0-30 -40-50 -60-70 σ=0.003 µm σ=0.01 µm σ=0.1 µm 0 0.001 0.01 0.1 σ (µm) -80 0.0 0. 0.4 0.6 0.8 1.0 1. R (µm)

Contact resistance model - effect of σ and R on contact resistance bounds Contact resistance (Ω) 0.4 0.3 0. 0.1 R=0.1 µm, upper bound R=0.1 µm, lower bound R=0.3 µm, upper bound R=0.3 µm, lower bound R=1 µm, upper bound R=1 µm, lower bound Resistance (Ω) 0.4 0.3 0. 0.1 σ=0.003 µm, upper bound σ=0.003 µm, lower bound σ=0.01 µm, upper bound σ=0.01 µm, lower bound σ=0.1 µm, upper bound σ=0.1 µm, lower bound 0.0 0.001 0.01 0.1 σ (µm) 0.0 0.0 0. 0.4 0.6 0.8 1.0 1. R (µm)

Effect of σ on sensitivity of resistance to force while loading Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 0.8 0.6 0.4 0. Contact Resistance (Ω) 1.4 Load, lower bound 1. Load, upper bound Unload, lower bound 1.0 Unload, upper bound 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) Contact Force (µn) N=100, R=0.1µm and σ=0.1 µm N=100, R=0.1µm and σ=0.03 µm

Effect of R on sensitivity of resistance to force while loading Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 0.8 0.6 0.4 0. Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) Contact Force (µn) N=100, R=1µm and σ=0.01 µm N=100, R=0.1µm and σ=0.01 µm

Model of film-covered surface Asperities with load < F th have zero surface energy, are electrically non-conductive N=100, R=0.1µm and σ=0.01 µm Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 0.8 0.6 0.4 0. Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) Contact Force (µn) F th = 5 µn F th = 0

Outline Introduction to Northeastern University switch Related work Modeling of contact surface with asperities (Greenwood & Williamson) Deformation model for single asperity loading without adhesion loading with adhesion unloading and separation Multi-asperity contact model Contact resistance model Measurements and discussion Conclusions

Measurement of contact resistance Contact bump Source V Drain 1 Drain Source 1 1.4 Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. 0.0 0 0 40 60 80 100 Actuation Voltage (V)

Translating actuation voltage into contact force 1.4 Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. 0.0 0 0 40 60 80 100 M R 1 F E Gate Contact Force Actuation Voltage (V) Contact Resistance (Ohms) 1.4 1. 1.0 0.8 0.6 0.4 0. Contact Force (µn) 400 300 00 100 0-100 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) -00 0 0 40 60 80 100 10 140 Actuation Voltage (V)

Evolution of measured characteristics 1.4 1.4 Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. Initial Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. 10 cycles 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) Contact Resistance (Ohms) 1.4 1. 1.0 0.8 0.6 0.4 0. Contact Resistance (Ohms) 1.4 1000 cycles 1. 10000 cycles 1.0 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn)

Can measured evolution be explained by change in σ? R=0.1µm, σ=0.1 µm, meas. 10 cycles Contact Resistance (Ω) Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 Measured after 10 cycles 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 1.4 Load, lower bound 1. Load, upper bound Unload, lower bound 1.0 Unload, upper bound Measured after 10000 cycles 0.8 0.6 0.4 0. Contact Force (µn) R=0.1µm, σ=0.003 µm, meas. 10 4 cycles R=0.1µm, σ=0.01 µm, meas. 1000 cycles Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 Measured after 1000 cycles 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn)

Does the surface film model explain measured evolution? σ=0.003 µm, F th = 3 µn meas. initial Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 Measured - initial 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) σ=0.003 µm, F th = 0 meas. 10000 cycles Contact Resistance (Ω) 1.4 Load, lower bound 1. Load, upper bound Unload, lower bound 1.0 Unload, upper bound Measured after 10000 cycles 0.8 0.6 0.4 0. σ=0.003 µm, F th = 1 µn meas. 1000 cycles Contact Resistance (Ω) 1.4 Load, lower bound Load, upper bound 1. Unload, lower bound Unload, upper bound 1.0 Measured after 1000 cycles 0.8 0.6 0.4 0. 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn) 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn)

Evolution towards flat punch geometry? Contact resistance load-half cycle extremely forceinsensitive after a high number of cycles, similar to unload-half cycle Suggests asperities may tend towards a flat punch-like geometry (very flat tops), and deformation becomes completely elastic. Completely force-insensitive contact while loading and unloading, all asperities separate simultaneously

Sample calculation for a simple flat-punch model: assume n contacting asperities, cylinders of equal heights and radii total contact area under 70 µn ~ 3x10-14 m (assuming plastic deformation), effective spot radius a eff = 10-7 m, each spot a = aeff / Total contact adherence force = n n 1/ 4 3 6πa Kw or π H eff a eff 1.4 F adh =-86 µn n=10=> contact resistance 0. Ω - 0.06 Ω n=50=> contact resistance 0. Ω - 0.03 Ω Contact Resistance (Ohms) 1. 1.0 0.8 0.6 0.4 0. 10000 cycles 0.0-80 -60-40 -0 0 0 40 60 80 100 Contact Force (µn)

Conclusions Initial contact not completely metallic (measured contact resistance too high) Metal-to-metal contact is obtained after cycling (typically 100-1000 cycles) measured contact resistance measured contact adherence evidence of material transfer Model predicts some plastic deformation of asperities even at very small loads. With a large number of contact spots (0-50) - ductile separation of asperities likely, makes material transfer more likely. Asperities unload like flat punches - explains forceinsensitive contact resistance while unloading

Conclusions - changes with cycling With cycling, measured contact resistance decreases, contact adherence increases, load half-cycle of characteristic becomes force-insensitive Decrease in contact resistance: removal of surface film, decrease in surface roughness. Increase in contact adherence: removal of surface film, decrease in surface roughness, blunting of asperities. Force-insensitive contact resistance while loading: blunting of asperities.

Conclusions - mechanisms of contact surface evolution Surface film removed by repeated contact, possibly plastic deformation and material transfer. Increasing metallic contact area may accelerate process further. Blunting of asperities by repeated contact, removal of surface film (surface forces cause necking of asperity), ductile separation? Decrease of surface roughness: ductile separation could cause material to be transferred from asperity to flat. Other mechanisms that can cause surface changes with time/cycling: electromigration, visco-elasticity /creep.

Conclusions - device design Central goal: obtain low and stable contact resistance and a favorable balance between restoring force and contact adherence force. To obtain better device performance: change actuator design to obtain greater actuator force at a given voltage; changes choice or deposition methods of the contact materials, or treatment of the contact surface, to reduce the contact adherence force.

Conclusions - device design Actuator design: force can be increased by increasing area, reducing gap, making beam more rigid in closed position trade-off with beam-gate capacitance, switching speed, manufacturability Contact material/surface contact surface should be clean and inert => low and stable contact resistance with small contact force. (gold) minimize material transfer for stable contact resistance - hard metal (Pt group) (mate hard metal with soft metal?) surface should have low surface energy; some Pt group metals form conducting oxide, but cleanliness hard to control, form frictional polymer.

Acknowledgments Prof. McGruer & Prof. Adams Prof. Zavracky Professors, staff and students at MFL: Rick Morrison, Weilin Hu, Keith Warner, Greg Jenkins, Michael Miller, Xiaomin Yan, Prof. Hopwood. Analog Devices and Radant MEMS.