Thermo-mechanical simulation demands a temperature distribution in the structure. It was found out that the electrical characteristics of the complete system do not considerably change with stress during standard operation. Therefore, the simulation can be separated into an electro-thermal and a thermo-mechanical part within small time periods as long as there is no void nucleation in the interconnect lines or the passivation is not broken.
In the simulation sequence as displayed in Fig. 8.1, the first part is the three-dimensional transient electro-thermal simulation of the interconnect structure in order to calculate the temperature distribution. Additionally this simulation delivers the potential and the current density. With the temperature distribution from the first part, the three-dimensional thermo-mechanical simulation can be performed subsequently. With the electro-thermal and the thermo-mechanical simulation all necessary capabilities for the rigorous simulation of electromigration are available.
Although STAP is also based on FEM, it is not appropriate to use it for thermo-mechanical simulations, because STAP is specialized and optimized for fast and accurate electro-thermal simulations. An extension of STAP to a more universal tool which can also handle mechanical problems would reduce its performance significantly. So even with the necessary data exchange the decoupled simulations with STAP and FEDOS are more efficient than a coupled simulation only performed by STAP.
The electro-thermal simulation is performed with the simulator STAP which uses also the finite element method for the calculation of the electric potential and temperature distribution. For the numerical calculation of Joule self-heating effects, caused by the electric current flow through the wire, two partial differential equations have to be solved [125,126]. The first one describes the electric subproblem
The next step is to compute the power loss density described by
The temperature dependence of the thermal and electrical conductivities is modeled with second order approximations:
High tensile stresses in the copper interconnects can cause break-up of the material and development of voids [127]. On the other side compressive stresses can induce the generation of extrusions. In case of temperature changes thermo-mechanical stress is build up because of (significant) different thermal expansion coefficients of adjacent materials.
The thermo-mechanical stress simulation is carried out with the program package FEDOS. The modeling of thermo-mechanical stress is similar to the stress calculation during thermal oxidation as described in Section 7.1. For the assumed elastic materials the stress tensor can be written in the form
The temperature change loads the mechanical problem on every finite element, because the internal force vector is
With the stress in principle also the hydrostatic pressure is given with the formula