# Electromigration

Electromigration is the transport of material caused by the gradual movement of the ions in a conductor due to the momentum transfer between conducting electrons and diffusing metal atoms. The effect is important in applications where high direct current densities are used, such as in microelectronics and related structures. As the structure size in electronics such as integrated circuits (ICs) decreases, the practical significance of this effect increases.

## History

The phenomenon of electromigration has been known for over 100 years, having been discovered by the French scientist Gerardin. The topic first became of practical interest during the late 1960s when packaged ICs first appeared. The earliest commercially available ICs failed in a mere three weeks of use from runaway electromigration, which led to a major industry effort to correct this problem. The first observation of electromigration in thin films was made by I. Blech. Research in this field was pioneered by a number of investigators throughout the fledgling semiconductor industry. One of the most important engineering studies was performed by Jim Black of Motorola, after whom Black's equation is named. At the time, the metal interconnects in ICs were still about 10 micrometres wide. Currently interconnects are only hundreds to tens of nanometers in width, making research in electromigration increasingly important.

## Practical implications of electromigration

Electromigration decreases the reliability of integrated circuits (ICs). It can cause the eventual loss of connections or failure of a circuit. Since reliability is critically important for space travel, military purposes, anti-lock braking systems, medical equipment like Automated External Defibrillators and is even important for personal computers or home entertainment systems, the reliability of chips (ICs) is a major focus of research efforts.
Due to the difficulty of testing under real-world conditions, Black's equation is used to predict the life span of integrated circuits. To use Black's equation, the component is put through high temperature operating life (HTOL) testing. The component's expected life span under real conditions is extrapolated from data gathered during this testing.
Although damage from electromigration ultimately results in the failure of the affected IC, the first symptoms are intermittent glitches, which are quite challenging to diagnose. As some interconnects fail before others, the circuit exhibits seemingly random errors, which may be indistinguishable from other failure mechanisms (such as electrostatic discharge damage). In a laboratory setting, electromigration failure is readily imaged with an electron microscope, as interconnect erosion leaves telltale visual markers on the metal layers of the IC.
With increasing miniaturization, the probability of failure due to electromigration increases in VLSI and ULSI circuits because both the power density and the current density increase. Specifically, line widths will continue to decrease over time, as will wire cross-sectional areas. Currents are also reduced due to lower supply voltages and shrinking gate capacitances. However, as current reduction is constrained by increasing frequencies, the more marked decrease in cross-sectional areas (compared to current reduction) will give rise to increased current densities in ICs going forward.
In advanced semiconductor manufacturing processes, copper has replaced aluminium as the interconnect material of choice. Despite its greater fragility in the fabrication process, copper is preferred for its superior conductivity. It is also intrinsically less susceptible to electromigration. However, electromigration (EM) continues to be an ever-present challenge to device fabrication, and therefore the EM research for copper interconnects is ongoing (though a relatively new field).
In modern consumer electronic devices, ICs rarely fail due to electromigration effects. This is because proper semiconductor design practices incorporate the effects of electromigration into the IC's layout. Nearly all IC design houses use automated EDA tools to check and correct electromigration problems at the transistor layout-level. When operated within the manufacturer's specified temperature and voltage range, a properly designed IC device is more likely to fail from other (environmental) causes, such as cumulative damage from gamma-ray bombardment.
Nevertheless, there have been documented cases of product failures due to electromigration. In the late 1980s, one line of Western Digital's desktop drives suffered widespread, predictable failure after 12–18 months of field usage. Using forensic analysis of the returned bad units, engineers identified improper design-rules in a third-party supplier's IC controller. By replacing the bad component with that of a different supplier, WD was able to correct the flaw, but not before significant damage was done to the company's reputation.
Electromigration can be a cause of degradation in some power semiconductor devices such as low voltage power MOSFETs, in which the lateral current through the source contact metallisation (often aluminium) can reach the critical current densities during overload conditions. The degradation of the aluminium layer causes an increase in on-state resistance, and can eventually lead to complete failure.

## Fundamentals

The material properties of the metal interconnects have a strong influence on their life span. The characteristics are predominantly the composition of the metal alloy and the dimensions of the conductor. The shape of the conductor, the crystallographic orientation of the grains in the metal, procedures for the layer deposition, heat treatment or annealing, characteristics of the passivation, and the interface to other materials also affect the durability of the interconnects. There are also important differences with time dependent current: direct current or different alternating current waveforms cause different effects.

### Forces on ions in an electrical field

Two forces affect ionized atoms in a conductor: 1) The direct electrostatic force F$_e$, as a result of the electric field $E$, which has the same direction as the electric field, and 2) The force from the exchange of momentum with other charge carriers F$_p$, toward the flow of charge carriers, is in the opposite direction of the electric field. In metallic conductors F$_p$ is caused by a so-called "electron wind" or "ion wind".
The resulting force F$_{res}$ on an activated ion in the electrical field can be written as
    $F_{res}=F_{e}-F_{p}=q\cdot (Z_{e}-Z_{p})\cdot E=q\cdot Z^{*}\cdot E=q\cdot Z^{*}\cdot j\cdot \rho$
where $q$ is the electric charge of the ions,$Z_{e}$ and $Z_{p}$ the valences corresponding to the electrostatic and wind force respectively, $Z^{*}$ the so-called effective valence of the material, $j$ the current density, and $\rho$ the resistivity of the material . Electromigration occurs when some of the momentum of a moving electron is transferred to a nearby activated ion. This causes the ion to move from its original position. Over time this force knocks a significant number of atoms far from their original positions. A break or gap can develop in the conducting material, preventing the flow of electricity. In narrow interconnect conductors, such as those linking transistors and other components in integrated circuits, this is known as a void or internal failure (open circuit). Electromigration can also cause the atoms of a conductor to pile up and drift toward other nearby conductors, creating an unintended electrical connection known as a hillock failure or whisker failure (short circuit). Both of these situations can lead to a malfunction of the circuit.

## Step Bunching due to Electromigration

Step bunching on DC-heated sublimating vicinal crystal surfaces of Si(111) was observed by A. Latyshev et al. in 1989. Soon after, Stoyan Stoyanov advanced a model in which as the reason for step bunching is identified the biased diffusion of the adatoms. In 1998, Stoyanov and Tonchev extended Stoyanov’s model by incorporating step-step repulsions and derived a scaling relation for the minimal step-step distance in a bunch under diffusion-limited sublimation, non-transparent steps, and step-down current conditions:
$l_{\min }\sim N^{-2/3}$
where $N$ is the number of steps in the bunch, and the proportionality coefficient has the dimension of length. This scaling law has been confirmed by numerous experimental studies. In 2018, Toktarbaiuly et al. reported electromigration-induced step bunching on vicinal W(110) surfaces. Their study revealed that step bunching occurred for both step-up and step-down current directions at the same temperature, T = 1500°C, with distinct size-scaling exponents depending on the current direction.
More recently, Usov et al. (2020) demonstrated that electromigration-induced step bunching is not limited to silicon surfaces but can also occur on dielectric surfaces, such as sapphire (Al₂O₃(0001)).
This study suggests that the fundamental mechanism of step bunching on W(110), Al₂O₃(0001), and Si(111) follows similar principles. Moreover, annealing W(110) offcut in the [001] direction with an up-step current produced a morphology where the bunch edges formed zigzag segments meeting at right angles.

## Failure mechanisms
### Diffusion mechanisms

In a homogeneous crystalline structure, because of the uniform lattice structure of the metal ions, there is hardly any momentum transfer between the conduction electrons and the metal ions. However, this symmetry does not exist at the grain boundaries and material interfaces, and so here momentum is transferred much more vigorously. Since the metal ions in these regions are bonded more weakly than in a regular crystal lattice, once the electron wind has reached a certain strength, atoms become separated from the grain boundaries and are transported in the direction of the current. This direction is also influenced by the grain boundary itself, because atoms tend to move along grain boundaries.
Diffusion processes caused by electromigration can be divided into grain boundary diffusion, bulk diffusion and surface diffusion. In general, grain boundary diffusion is the major electromigration process in aluminum wires, whereas surface diffusion is dominant in copper interconnects.

### Thermal effects

In an ideal conductor, where atoms are arranged in a perfect lattice structure, the electrons moving through it would experience no collisions and electromigration would not occur. In real conductors, defects in the lattice structure and the random thermal vibration of the atoms about their positions causes electrons to collide with the atoms and scatter, which is the source of electrical resistance (at least in metals; see electrical conduction). Normally, the amount of momentum imparted by the relatively low-mass electrons is not enough to permanently displace the atoms. However, in high-power situations (such as with the increasing current draw and decreasing wire sizes in modern VLSI microprocessors), if many electrons bombard the atoms with enough force to become significant, this will accelerate the process of electromigration by causing the atoms of the conductor to vibrate further from their ideal lattice positions, increasing the amount of electron scattering. High current density increases the number of electrons scattering against the atoms of the conductor, and hence the rate at which those atoms are displaced.
In integrated circuits, electromigration does not occur in semiconductors directly, but in the metal interconnects deposited onto them (see semiconductor device fabrication).
Electromigration is exacerbated by high current densities and the Joule heating of the conductor (see electrical resistance), and can lead to eventual failure of electrical components. Localized increase of current density is known as current crowding.

### Balance of atom concentration

A governing equation which describes the atom concentration evolution throughout some interconnect segment, is the conventional mass balance (continuity) equation
    ${\frac {\partial N}{\partial t}}+\nabla \cdot {\vec {J}}=0$
where $N({\vec {x}},t)$ is the atom concentration at the point with a coordinates ${\vec {x}}=(x,y,z)$ at the moment of time $t$, and $J$ is the total atomic flux at this location. The total atomic flux $J$ is a combination of the fluxes caused by the different atom migration forces. The major forces are induced by the electric current, and by the gradients of temperature, mechanical stress and concentration. ${\vec {J}}={\vec {J}}_{c}+{\vec {J}}_{T}+{\vec {J}}_{\sigma }+{\vec {J}}_{N}$.
To define the fluxes mentioned above:
    ${\vec {J}}_{c}={\frac {NeZD\rho }{kT}}\,{\vec {j}}$.
Here $e$ is the electron charge, $eZ$ is the effective charge of the migrating atom, $\rho$ the resistivity of the conductor where atom migration takes place, ${\vec {j}}$ is the local current density, $k$ is the Boltzmann constant, $T$ is the absolute temperature. $D({\vec {x}},t)$ is the time and position dependent atom diffusivity.
    ${\vec {J}}_{T}=-{\frac {NDQ}{kT^{2}}}\nabla T$. We use $Q$ the heat of thermal diffusion.
    ${\vec {J}}_{\sigma }={\frac {ND\Omega }{kT}}\nabla \!H$,
here $\Omega =1/N_{0}$ is the atomic volume and $N_{0}$ is initial atomic concentration, $H=(\sigma _{11}+\sigma _{22}+\sigma _{33})/3$ is the hydrostatic stress and $\sigma _{11},\sigma _{22},\sigma _{33}$ are the components of principal stress.
    ${\vec {J}}_{N}=-D\,\nabla \!N$.
Assuming a vacancy mechanism for atom diffusion we can express $D$ as a function of the hydrostatic stress $D=D_{0}\exp \left({\tfrac {\Omega H-E_{A}}{kT}}\right)$ where $E_{A}$ is the effective activation energy of the thermal diffusion of metal atoms. The vacancy concentration represents availability of empty lattice sites, which might be occupied by a migrating atom.

## Electromigration-aware design
### Electromigration reliability of a wire (Black's equation)

At the end of the 1960s J. R. Black developed an empirical model to estimate the MTTF (mean time to failure) of a wire, taking electromigration into consideration. Since then, the formula has gained popularity in the semiconductor industry:
    ${\text{MTTF}}={\frac {A}{J^{n}}}\exp {\left({\frac {E_{\text{a}}}{kT}}\right)}$.
Here $A$ is a constant based on the cross-sectional area of the interconnect, $J$ is the current density, $E_{\text{a}}$ is the activation energy (e.g. 0.7 eV for grain boundary diffusion in aluminum), $k$ is the Boltzmann constant, $T$ is the temperature in kelvins, and $n$ a scaling factor (usually set to 2 according to Black). The temperature of the conductor appears in the exponent, i.e. it strongly affects the MTTF of the interconnect. For an interconnect of a given construction to remain reliable as the temperature rises, the current density within the conductor must be reduced. However, as interconnect technology advances at the nanometer scale, the validity of Black's equation becomes increasingly questionable.

### Wire material

Historically, aluminium has been used as conductor in integrated circuits, due to its good adherence to substrate, good conductivity, and ability to form ohmic contacts with silicon. However, pure aluminium is susceptible to electromigration. Research shows that adding 2-4% of copper to aluminium increases resistance to electromigration about 50 times. The effect is attributed to the grain boundary segregation of copper, which greatly inhibits the diffusion of aluminium atoms across grain boundaries.
Pure copper wires can withstand approximately five times more current density than aluminum wires while maintaining similar reliability requirements. This is mainly due to the higher electromigration activation energy levels of copper, caused by its superior electrical and thermal conductivity as well as its higher melting point. Further improvements can be achieved by alloying copper with about 1% palladium which inhibits diffusion of copper atoms along grain boundaries in the same way as the addition of copper to aluminium interconnect.

### Bamboo structure and metal slotting

A wider wire results in smaller current density and, hence, less likelihood of electromigration. Also, the metal grain size has influence; the smaller grains, the more grain boundaries and the higher likelihood of electromigration effects. However, if you reduce wire width to below the average grain size of the wire material, grain boundaries become "crosswise", more or less perpendicular to the length of the wire. The resulting structure resembles the joints in a stalk of bamboo. With such a structure, the resistance to electromigration increases, despite an increase in current density. This apparent contradiction is caused by the perpendicular position of the grain boundaries; the boundary diffusion factor is excluded, and material transport is correspondingly reduced.
However, the maximum wire width possible for a bamboo structure is usually too narrow for signal lines of large-magnitude currents in analog circuits or for power supply lines. In these circumstances, slotted wires are often used, whereby rectangular holes are carved in the wires. Here, the widths of the individual metal structures in between the slots lie within the area of a bamboo structure, while the resulting total width of all the metal structures meets power requirements.

### Blech length

There is a lower limit for the length of the interconnect that will allow higher current carrying capability. It is known as "Blech length". Any wire that has a length below this limit will have a stretched limit for electromigration. Here, a mechanical stress buildup causes an atom back flow process which reduces or even compensates the effective material flow towards the anode. The Blech length must be considered when designing test structures to evaluate electromigration. This minimum length is typically some tens of microns for chip traces, and interconnections shorter than this are sometimes referred to as 'electromigration immortal'.

### Via arrangements and corner bends

Particular attention must be paid to vias and contact holes. The current carrying capacity of a via is much less than a metallic wire of same length. Hence multiple vias are often used, whereby the geometry of the via array is very significant: multiple vias must be organized such that the resulting current is distributed as evenly as possible through all the vias.
Attention must also be paid to bends in interconnects. In particular, 90-degree corner bends must be avoided, since the current density in such bends is significantly higher than that in oblique angles (e.g., 135 degrees).

### Electromigration in solder joints

The typical current density at which electromigration occurs in Cu or Al interconnects is 10$^6$ to 10$^7$ A/cm$^2$. For solder joints (SnPb or SnAgCu lead-free) used in IC chips, however, electromigration occurs at much lower current densities, e.g. 10$^4$ A/cm$^2$. It causes a net atom transport along the direction of electron flow. The atoms accumulate at the anode, while voids are generated at the cathode and back stress is induced during electromigration. The typical failure of a solder joint due to electromigration will occur at the cathode side. Due to the current crowding effect, voids form first at the corners of the solder joint. Then the voids extend and join to cause a failure. Electromigration also influences formation of intermetallic compounds, as the migration rates are a function of atomic mass.

### Electromigration and technology computer aided design

The complete mathematical model describing electromigration consists of several partial differential equations (PDEs) which need to be solved for three-dimensional geometrical domains representing segments of an interconnect structure. Such a mathematical model forms the basis for simulation of electromigration in modern technology computer aided design (TCAD) tools. Use of TCAD tools for detailed investigations of electromigration induced interconnect degradation is gaining importance. Results of TCAD studies in combination with reliability tests lead to modification of design rules improving the interconnect resistance to electromigration.

### Electromigration due to IR drop noise of the on-chip power grid network/interconnect

The electromigration degradation of the on-chip power grid network/interconnect depends on the IR drop noise of the power grid interconnect. The electromigration-aware lifetime of the power grid interconnects as well as the chip decreases if the chip suffers from a high value of the IR drop noise.

### Machine Learning Model for Electromigration-aware MTTF Prediction

Recent work demonstrates MTTF prediction using a machine learning model. The work uses a neural network-based supervised learning approach with current density, interconnect length, interconnect temperature as input features to the model.

## Electromigrated nanogaps

Electromigrated nanogaps are gaps formed in metallic bridges formed by the process of electromigration. A nanosized contact formed by electromigration acts like a waveguide for electrons. The nanocontact essentially acts like a one-dimensional wire with a conductance of $G=2\,e^{2}\!/h$. The current in a wire is the velocity of the electrons multiplied by the charge and number per unit length, $\,I=veN/L\$ or $\ G=veN/LV$. This gives a conductance of $G=ve^{2}\!N/LE$. In nanoscale bridges the conductance falls in discrete steps of multiples of the quantum conductance $G=2\,e^{2}\!/h$.
Electromigrated nanogaps have been proposed for use as electrodes in molecular scale electronics and as quantum tunneling sensors. Researchers have used feedback controlled electromigration to investigate the magnetoresistance of a quantum spin valve.

## Reference standards

* EIA/JEDEC Standard EIA/JESD61: Isothermal Electromigration Test Procedure.
* EIA/JEDEC Standard EIA/JESD63: Standard method for calculating the electromigration model parameters for current density and temperature.
* Fundamentals of electromigration, Chapter 2