Next Article in Journal
A Low-Noise Interface ASIC for MEMS Disk Resonator Gyroscope
Previous Article in Journal
A Review of Fingerprint Sensors: Mechanism, Characteristics, and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Integrated Fan-Out Wafer-Level Package Micro-Bump Electromigration Study

1
Hangzhou Institute of Technology, Xidian University, Hangzhou 311231, China
2
School of Electro-Mechnical Engineering, Xidian University, Xi’an 710071, China
3
Zhongkexin Integrated Circuit Co., Wuxi 214035, China
4
Yangzhou Yangjie Electronic Technology Co., Ltd., Yangzhou 225008, China
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(6), 1255; https://doi.org/10.3390/mi14061255
Submission received: 19 April 2023 / Revised: 30 May 2023 / Accepted: 12 June 2023 / Published: 15 June 2023
(This article belongs to the Section A:Physics)

Abstract

:
To meet the demands for miniaturization and multi-functional and high-performance electronics applications, the semiconductor industry has shifted its packaging approach to multi-chip vertical stacking. Among the advanced packaging technologies for high-density interconnects, the most persistent factor affecting their reliability is the electromigration (EM) problem on the micro-bump. The operating temperature and the operating current density are the main factors affecting the EM phenomenon. Therefore, when a micro-bump structure is in the electrothermal environment, the EM failure mechanism of the high-density integrated packaging structure must be studied. To investigate the relationship between loading conditions and EM failure time in micro-bump structures, this study established an equivalent model of the vertical stacking structure of fan-out wafer-level packages. Then, the electrothermal interaction theory was used to carry out numerical simulations in an electrothermal environment. Finally, the MTTF equation was invoked, with Sn63Pb37 as the bump material, and the relationship between the operating environment and EM lifetime was investigated. The results showed that the current aggregation was the location where the bump structure was most susceptible to EM failure. The accelerating effect of the temperature on the EM failure time was more obvious at a current density of 3.5 A/cm2, which was 27.51% shorter than 4.5 A/cm2 at the same temperature difference. When the current density exceeded 4.5 A/cm2, the change in the failure time was not obvious, and the maximum critical value of the micro-bump failure was 4 A/cm2~4.5 A/cm2.

1. Introduction

In recent years, in order to meet the needs of miniaturization and multi-functional and high-performance applications for electronic products, the requirements for microelectronic packaging have increased. When improving the performance of microelectronic systems, it is necessary to reduce the size and cost, which brings great challenges to traditional wafer-level packaging. As the traditional fan-in package is no longer able to complete the multilayer rewiring and bump array arrangement in its chip area, fan-out wafer-level packaging (FOWLP) has emerged. By increasing the area of a single package through wafer reconfiguration, the fan-out package breaks the limitation of the number of I/O terminals. Then, the advanced manufacturing process of wafer-level packaging is applied to complete the multilayer rewiring and bump preparation, in addition to cutting and separating, to obtain a package that can interconnect with external electrical properties [1,2,3].
Fan-out wafer-level packaging (FOWLP), an advanced packaging technology developed in recent years, eliminates the lap wire and flip chip bumps in traditional packaging structures and replaces the traditional PCB substrate with the redistribution layer (RDL). It has the advantages of high integration, high bandwidth, and short interconnection and is gradually developing into the mainstream technology of packaging [4,5,6,7]. On the other hand, FOWLP faces many critical issues to be solved in terms of reliability, such as multilayer-interconnect migration, large residual stresses, bump voids, bonding failures, and other failure behaviors. In advanced packaging technologies for high-density interconnects, EM problems on micro-bumps are the most persistent cause of reliability. Micro-bumps function as electrical and mechanical connections between FOWLP components, and their integrity and reliability are critical factors for the stable operation of electrical components [8]. While achieving high-density interconnections, current density aggregation and increased Joule heating have led to current congestion and rising temperature problems on micro-bumps [9,10,11,12]. Therefore, the operating temperature and operating current density are the main factors affecting the EM phenomenon on micro-bumps. Electromigration is a phenomenon in which atoms migrate rapidly with electrons under the action of high-density currents, creating voids or bumps and leading to component failure within the micro-bumps [13,14]. With the development and application of three-dimensional integration technology, the electromigration phenomenon in micro-bumps has attracted much attention [15,16]. Compared to other packaging technologies, FOWLP has a smaller package size and, therefore, is in a denser arrangement of micro-bumps. In the high-density interconnections of FOWLP, the electromigration failure phenomenon is highlighted due to the interaction of heat and electricity [17,18,19,20,21]. At present, the analysis of EM failure mechanisms for micro-bump structures is mostly based on experimental phenomena. In the literature [22], an interconnection model was developed for GaInSn/Cu materials. By applying different current densities to the model, the EM phenomenon under the current flux was obtained, which verified the migration of atoms from the cathode to the anode through the GaInSn path. In addition, the room temperature EM critical current density of Cu/GaInSn/Cu was two orders of magnitude lower compared to other EM systems. In multi-physical field conditions, the literature [23] studied the EM behavior caused by current crowding in Cu/Sn3.0Ag0.5Cu/Cu ball grid array (BGA) solder joints. By loading the current density and temperature conditions, it predicted the cavity formation location and critical value. The literature combines experimental phenomena to analyze the critical value of cavity formation from the perspective of atomic concentration. It does not analyze the relationship between loading conditions and lifetime in depth.
In this study, in order to obtain the relationship between the loading conditions and lifetime, the failure mechanism of the micro-bump EM phenomenon was investigated in a multi-physical field environment. Firstly, developing an equivalent model for the micro-bump structure of a high-density integrated package and its failure behavior was numerically simulated in an electrical-thermal environment. Then, the data were substituted into the MTTF equation, with Sn63Pb37 as the bump materials to predict the law, which is about the EM life of electronic products in relation to the loading conditions. Due to the adoption of the equivalent model of the advanced package form and using a high-density micro-bump interconnection structure as the research object, the research content will provide an important reference for extending the lifetime and improving the reliability of advanced packaging components.

2. Theoretical Analysis and Model Design

2.1. Electrothermal Interaction Theory

Under single environmental stresses, such as temperature cycling, high temperatures, EM, and mechanical stresses, the study of failure behavior and the development of life models of micro-bumps are gradually being advanced and improved [24]. However, under complex practical working conditions, it is difficult to singularize the effects of various environmental stresses by primary and secondary effects. From the perspective of close to the actual working conditions, it is necessary to explore the reliability of the micro-bump structure under multi-physical fields. When loading electronic products with both current and temperature conditions, the electromagnetic heat generated by the current will increase the temperature. Because of this, electrical–thermal conditions can have an impact on the heat transfer environment. In studying the relationship between loading conditions and lifetime, a theoretical model of the electrical–thermal interaction is required.
When an operating current is applied, the current density j follows the current conservation equation:
j = ( σ + ε 0 ε r t ) E + j e ,
E = V ,
σ = σ 0 1 + T T 0 ,
where σ is the material resistivity, j e is the external current density, ε 0 is the vacuum dielectric constant, ε r is the relative dielectric constant, V is the potential gradient, σ 0 is the resistivity at the temperature T 0 , and is the material temperature coefficient.
During the numerical simulation, the temperature field follows the heat transfer equation:
ρ C p u T + q = Q + Q t e d ,
q = k T ,
where ρ is the material density, Cp is the constant pressure heat capacity, q is the heat flux density, Q is the heat, Qted is the thermo-elastic damping, and k is the thermal conductivity.
In the electrothermal environment, current flow generates Joule heat, which interacts with the temperature field. The current and temperature field produce electromagnetic–thermal effects, following the electromagnetic heat equation for a multi-physics field:
ρ C p u T = ( k T ) + Q e ,
Q e = J E ,
where ρ is the material density, Cp is the constant-pressure heat capacity, k is the thermal conductivity, and Qe is the heat generated by the electric field.

2.2. Numerical Simulation Model Construction

FOWLP maximizes the interconnect density and enables the seamless connection of semiconductor devices through two-to-five-micron spaced lines. At the same time, achieving high-bandwidth data transmission and removing the substrate results in significant cost savings. In addition, in the field of advanced packaging, high-density integrated packages with built-in vertical stacks are more prone to micro-bump EM failure problems. For this reason, a typical top-and-bottom vertically stacked micro-assembly structure was chosen as the experimental model. The actual model comes from DDR3 micro-assemblies produced by the 58th Research Institute of China Electronics Technology Group Corporation. DDR3 has a higher external data transfer rate and a more advanced topology architecture. On the other hand, it further reduces power consumption while ensuring performance. It used two 2 Gb DDR3 bare cores inside the overall micro-assembly model. The package structure was integrated with TSV technology on both sides, with vertical stacking of the top and bottom of the chip. Figure 1 shows the physical diagram of the micro-assembly.
The two DDR3s died in the micro-assembly, were pin-processed using fan-out technology and wired based on the redistribution layer (RDL) process. The model had a theoretical total capacity of 4 Gb and a bit width of 32 bits. Without affecting the simulation results, simplifying the model structure can reduce the workload and improve the calculation accuracy. Since the main object of this numerical simulation was the micro-bump structure, the TSV structure on both sides was simplified. The specific model was constructed as shown in Figure 2.
For the overall model, the micro-bump was the focus and difficulty, which needed to be refined. The size of the micro-bump changed in a thermo-electric coupled environment. Therefore, adjustments were made to the modeled size of the bump, as shown in Figure 2. In order to obtain the current density distribution that meets the engineering reality, a daisy chain structure was added to the micro-bump part. The daisy chain structure is made of Cu material. For forming a pathway between the micro-bumps, a wrap-around type complete connection was constructed between the entire micro-bump structure. Between the micro-bumps, the conductive properties of the material and the daisy chain structure design can realize the effective current flow. The simplified micro-assembly structure consisted of two chips (Si) and a fan-out package structure (polyimide, PI) stacked vertically, two RDL layers equivalent to a thin plate structure [25,26], and two 16 × 22 micro-bump structures (Sn63Pb37) on the top and bottom. The specific dimensions are shown in Table 1.
According to the dimensional parameters and the refined structural design of the key parts, the equivalent numerical simulation model was constructed. The specific model is shown in Figure 3.

3. Electrothermal Environmental Loading Conditions

3.1. Material Parameters

The chip layer, fan-out layer, and micro-bump were composed of a single material, and the material parameters used are shown in Table 2. The RDL structure consists of multilayer wiring, whose components mainly include PI and Cu. When assigning materials, the error of using a single material property is large. In order to ensure the accuracy of the electrical and thermal conductivity properties of the RDL layer, the material properties were obtained using the volume percentage method [27,28]. The volume percentage method is the most common method to determine the properties of equivalent materials, which can be obtained with isotropic and homogeneous properties. Based on the composition of the DDR micro-assembly rewiring layer, the PI was calculated to account for 57% and Cu for 43%. According to this ratio, the RDL layer material properties were obtained. The material parameters used in the model are shown in Table 2.
Since Sn63Pb37 is the key material to be studied in the electrothermal environment, its thermal conductivity properties are particularly important. The thermal conductivity of Sn63Pb37 as a bump material at different temperatures is labeled in Table 3.

3.2. Mesh Division and Boundary Conditions

Referring to the experimental case of DDR micro-assembly, the overall temperature constraint (temperature = 398.15 K (corresponding to 125 °C)) and thermal convection constraint (natural convection heat transfer coefficient of 20 W/(m2·°C)) were imposed. A current density of 4 A/cm2 was applied at the chip end, which flowed into the micro-bump end through the RDL layer. In addition to this, a ground constraint was applied below the micro-bump. The overall meshing is shown in Figure 4, with 950,778 mesh vertices and 4,649,295 cells. The average quality of the mesh is around 0.7. The closer the mesh quality is to 1, the more accurate the calculation result.
After completing the meshing, the simulation calculation was carried out on the model. The overall current density distribution and temperature distribution are shown in Figure 5 and Figure 6. From the result cloud diagram, it can be seen that the local current dense point location appeared in the middle area of the chip at the RDL layer and the bump connection, which is consistent with the actual current flow direction. The current density was 593.60 A/cm2. The highest local temperature point location appeared at the edge of the chip, and the temperature at the hot spot was 461.36 K.
The overall DDR model simulation was roughly calculated for the current density and temperature distribution. Since the heat generation at the chip was not negligible, the location distribution of current gathering points and hot spots was in accordance with the electric–thermal interaction principle. In order to improve the mesh quality, more accurate calculation results of key bumps can be obtained and the workload streamlined, while the 1/4 micro-component model can be used for subsequent work. The daisy chain structure was reset to the wraparound mutual circulation form. Symmetry constraints were added to the sub-model in order to make the results fit the DDR micro-component operating context. After fine mesh division, 325,534 mesh vertices, 1,629,738 cells, and an average mesh quality of 0.8 were available. The 1/4 micro-component model construction and the mesh division are shown in Figure 7.
This study loaded the model with an electrical-thermal environment using a numerical simulation. The simulation software used COMSOL, which could work on the principle of electrical–thermal interaction and imposed boundary conditions in conjunction with the actual operating conditions. The operating temperature of the micro-assembly was −40 °C~85 °C, and the operating voltage was in the range of −0.4 V~1.975 V. Therefore, 55 °C (328.15 K) and 125 °C (398.15 K) were selected as the normal operating temperature and accelerated experimental temperature, respectively. In the selection of current density, there were differences between experimental and numerical simulations. Based on the micro-assembly actual measurement dataand combined with the normal operating voltage range, 3.5 A/cm2~5 A/cm2 interval 0.5 was selected as the current density [22,29]. In addition to this, when the current critical value is reached, the change in MTTF will be small as the current density changes. The predicted normal operating current maximum was within this range. To simulate the natural heat dissipation environment, thermal convection constraints (natural convection heat transfer coefficient of 20 W/(m2·°C)) were still added. The specific loading conditions are shown in Table 4.
At the selected electrothermal environment, the micro-bump simulation results at both the operating and accelerating temperatures can be observed. Moreover, these experimental phenomena varied greatly with the current density at the same interval.

3.3. Theoretical Model for Electromigration Lifetime Prediction

The EM failure mechanism of micro-bump structures is the mass migration of atoms under the action of current stress. In order to obtain more accurate EM failure times, a thermo-electric reliability model needs to be constructed for Sn63Pb37. The conventional MTTF theory was proposed by Black in 1969, which analyzes the EM lifetime prediction [30].
M T T F = A 1 ( j ) n exp E a k T ,
where MTTF is related to the current density j and temperature T through three parameters, Ea is the activation energy of atomic diffusion in EM, A is the coefficient value, and n is the current density power factor. In general, both A and n are empirical valuations under practical applications.
The current flow will have Joule heating, and the entropy generated by this process is waste heat. The accumulation of entropy is carried out through the EM of atoms. Therefore, to revisit Black’s equation, the equation is considered in terms of entropy production in irreversible processes. Consider the link between entropy production and microstructural changes under mechanical damage. The theory of Onsager’s entropy yield in irreversible processes is [31]:
T d S V d t = J X e ,
where T is the temperature, d S d t is the entropy yield, V is the volume of the specimen, and J and Xe are the electron wind flux and driving force, respectively.
In EM, we consider the MTTF as the accumulation time of the threshold entropy (Sthreshold). According to Equation (9), the total entropy production up to the failure is derived as follows:
t f a i l u r e = M T T F = T S t h r e s h o l d V J e X e ,
where Je is the electron wind flux, Xe is the driving force, and D is the atomic diffusivity. The MTTF is derived as follows:
J e = c D k T Z e   ρ   j ,
X e = Z e   ρ   j ,
D = D 0 exp ( E a k T ) ,
M T T F t f a i l u r e = T S V J e X e = A 1 j 2 1 D = A j 2 exp ( E a k T ) .
From the derivation process above, the empirical value of Black’s equation can be defined as n = 2. One of the most widely used bump materials is Sn63Pb37. Finding the general A value of this material can quickly and accurately predict the EM lifetime of electronic products. The EM parameters of Sn63Pb37 are shown in Table 5.
The value of coefficient A of the MTTF equation has a significant impact on the electrical migration failure time. In the electric-thermal environment, concluding the first failure time of electronic products under different working conditions and the A value are obtained by substituting n = 2. According to the literature [34], the Weibull distribution of the MTTF data was obtained. By reading the article’s logic, it was found that there were data that clearly did not match the engineering reality. Following the article [35], the data were optimized within a reasonable range, and the experimental data of Sn63Pb37 were obtained. Next, the MTTF data from the Weibull distribution were quoted, and the A-value was calculated according to Equation (14).
It was calculated that the five sets of A-value data for this material in the literature [20] were calculated within the same order of magnitude. When the data error is not significant, the average value can be used for theoretical guidance. The average value is of the theoretical guidance. As a consequence, the A-value in the MTTF equation for Sn63Pb37, as a bump material, is 3.64 × 10−3. Substituting the A-value into Equation (14) yields Black’s equation for Sn63Pb37 as:
M T T F = 3.64 × 10 3 j 2 exp ( E a k T ) .
Equation (15) is the established reliability model for life calculation, which provides the theoretical basis for finding the average failure time of the micro-bump structure. Finally, the data are used to analyze the current density and temperature distribution in relation to the lifetime of the micro-bump.

4. Discussion

By observing the simulation results, the dangerous bumps that are most prone to EM failure were identified. The current density distribution and temperature distribution of current aggregation points and hot spots were counted, and the MTTF was calculated using Equation (15). After the mesh quality was improved, the simulation data of the overall model and the sub-model under the same working conditions were compared. The ambient temperature was 398.15 K, the natural convection heat transfer coefficient was 20 W/(m2·°C), and the current density applied at the chip end was 4 A/cm2. Figure 8 and Figure 9 show the resultant clouds at the current aggregation points and hot spots.
By observing the cloud diagram, it can be directly seen that in the simulation calculation of the 1/4 model, the maximum current density was 717.03 A/cm2, and the hot spot temperature was 483.15 K compared to 593.60 A/cm2 and 461.36 K of the overall model in the same environment, respectively. The calculation results were more accurate after the grid quality was improved. In order to specifically describe the location of the hazard bumps, the coordinates of the bumps were set for the sub-model as in Figure 10. The highest local temperature points could be obtained at points H2, H3, and the local current gathering point at point G11 after the electrical–thermal simulation.
In the sub-model, the location of the current gathering point remained the same, however, the location of the highest local temperature point was shifted. According to the existing research, the failure phenomenon first occurs at the edges and corners when micro-convex structures are affected by the temperature [36,37]. The location of the highest local temperature point changed from the inside edge of the chip in the overall model to the outside edge. However, the maximum point of the temperature was always at the edge of the chip with a small offset. This indicates that when the simulation reached the steady state with increased computational accuracy, the effect of the temperature was more pronounced. The hot spot location did not exactly follow the current gathering point, and the effect of the temperature on EM failure came to the fore. The location of the current density aggregation point appeared to be consistent with the existing experimental phenomena [22].
Under different operating conditions, numerical simulations were carried out to obtain the current density distribution and temperature distribution. From the available studies in the literature, it can be concluded that dangerous solder joints exist at the maximum current and at the hot spots [38]. The locations of the dangerous bumps in the simulation results were found, and the values were counted. Then, the MTTF results were calculated using Equation (15). The simulation results and calculation results are shown in Table 6. For the same operating conditions, the first row shows the results at the current density, and the second row shows the results at the hot spots.
The data in the table above illustrate that under the same operating conditions, the calculated MTTF values at the current aggregation point in a steady state were smaller than those at the hot spot. Therefore, the influence of current density on the EM lifetime always dominates. The hot spot coordinates were shifted when the current densities were both 3.5 A/cm2 and the temperatures were 328.15 K and 398.15 K. The hot spot coordinates are H2 and H3, respectively. The results are shown in Figure 11 and Figure 12. This indicates that the hot spot locations were not fixed but moved within a certain range. When current density was applied to the chip to generate heat and add ambient temperature, the hot spot coordinates appeared at the location of the chip edge in the sub-model. When the current density dominates, the effect on the EM is greater as the temperature rises. Therefore, as the temperature changes, the hot spot location will move to a small area at the edge of the chip.
The normal operating temperatures of DDR components were −40 °C–85 °C (233.15 K–358.15 K). Therefore, the temperature for accelerated experimental conditions was 398.15 K. The MTTF was most obvious when the current density was 3.5 A/cm2, and the MTTF was shortened rapidly with the increase in temperature. Figure 13 shows the comparison of MTTF with the change in current under different temperatures.
By observing the trend in the comparison curves, it is found that the accelerating effect of the temperature on EM failure becomes less and less as the circuit density increases. When the temperature increased from 328.15 K to 398.15 K, and the current density was 3.5 A/cm2, the MTTF at the current aggregation changed from 7239.30 h to 1553.34 h. Comparing the MTTF at the two temperatures shows that the MTTF at 398.15 K was shortened by 21.46% at 328.15 K for a current density of 3.5 A/cm2. At the current density of 4.5 A/cm2, the MTTF at the current aggregation changed from 64.51 h to 31.59 h, and the MTTF at 398.15 K is shortened to 48.97% of that at 328.15 K. Therefore, at the same temperature difference, the MTTF decreased at a current density of 3.5 A/cm2 at a faster rate than at 4.5 A/cm2, with a 27.51% increase in the percentage decrease. Compared to the normal BGA package, the temperature variation and current density variation were more pronounced under the fan-out package [39]. From the comparison curves, it can be seen that the decreasing trend in MTTF changed the same when the temperature was different. The decreasing rate was the fastest at 3.5 A/cm2 to 4 A/cm2 and then gradually slowed down. It can be seen that the fundamental reason for the change in MTTF at different temperatures is that the increase in temperature intensifies the atomic migration rate, leading to faster EM failure.
When the temperature was 328.15 K, and the current density changed from 3.5 A/cm2 to 4 A/cm2, the MTTF at the current aggregation was changed from 7239.30 h to 643.84 h, which was shortened to 8.89% of the original. The data shows that the current density was the main influencing factor for the occurrence of EM failure. The experimental results are logically identical to those of the literature [34]. From the trend of Figure 13, it can be seen that as the current density increases to a certain value, its MTTF does not change significantly. At this time, the maximum threshold of the working current density has been reached. The maximum working current density of DDR micro-assembly is between 4 A/cm2–4.5 A/cm2. The MTTF decreased to about 97% of the original when the current density changed from 3.5 A/cm2 to 5 A/cm2, confirming that 5 A/cm2 has become an accelerated experimental condition. In addition, existing experimental studies have verified the correctness of the threshold current density [40].

5. Conclusions

In this study, the equivalence model was based on the DDR micro-assembly technology manual. Using the theory of electrothermal environment interaction, numerical simulations were performed in the electrothermal environment. According to the resultant plots, the current density distribution and temperature distribution of the hazardous bumps were obtained. The data were substituted into the MTTF equation for Sn63Pb37 to calculate the EM lifetime of the micro-bump structure. Finally, the relationship between the working environment and EM failure time was explored using statistical data.
The results showed that the most likely location for the EM failure of the bump structures was at the current aggregation. In the electrothermal environment, the current dominated the EM failure. The accelerating effect of the temperature on EM failure time was more obvious when the current density was 3.5 A/cm2, which was 27.51% shorter than 4.5 A/cm2 at the same temperature difference. In addition, when the current density exceeded 4.5 A/cm2, the change in failure time was not obvious, and the maximum critical value of the working current density of micro-bump structures was at 4 A/cm2–4.5 A/cm2.
Therefore, the study of the EM failure time of electronic products under different working conditions is extremely useful for exploring the relationship between the working environment and EM life. The conclusion helps to extend the EM life of the product represented by DDR3.

Author Contributions

W.T. designed and guided the complete related simulation and experiment. R.G. wrote this paper. L.G. offered simulation method guidance. H.J. provided physical model information. L.Z. offered the help to modify the grammar. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by China Postdoctoral Science Foundation (No. 2022M712975).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tian, W.; Li, B.; Li, Z.; Cui, H.; Shi, J.; Wang, Y.; Zhao, J. Using Chiplet Encapsulation Technology to Achieve Processing-in-Memory Functions. Micromachines 2022, 13, 1790. [Google Scholar] [CrossRef]
  2. Chang, K.T.; Huang, C.; Kuo, H.; Jhong, M.; Hsieh, T.; Hung, M.; Wang, C. Ultra High Density IO Fan-Out Design Optimizationwith Signal Integrity and Power Integrity. In Proceedings of the 2019 IEEE 69th Electronic Components and Technology Conference (ECTC), Las Vegas, NV, USA, 28–31 May 2019; pp. 41–46. [Google Scholar]
  3. Chong, C.T.; Guan, L.T.; Ho, D.; Choong, C.S.; Faxing, C.; Bhattacharya, S. High Density Interconnection for Heterogeneous Integration on FOWLP Platform. In Proceedings of the 2018 IEEE 20th Electronics Packaging Technology Conference (EPTC), Singapore, 4–7 December 2018; pp. 186–190. [Google Scholar]
  4. Zhu, J.; Ming, X.; Yao, X. Research on Key Process Technology of RDL-first Fan-out Wafer Level Packaging. In Proceedings of the 2018 19th International Conference on Electronic Packaging Technology (ICEPT), Shanghai, China, 8–11 August 2018. [Google Scholar]
  5. Beth, K.; Kroehnert, S. (Eds.) Advances in Embedded and Fan-Out Wafer Level Packaging Technologies; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  6. Ghaffarian, R. Microelectronics packaging technology roadmaps, assembly reliability, and prognostics. Facta Univ. Ser. Electron. Energ. 2016, 29, 543–611. [Google Scholar] [CrossRef] [Green Version]
  7. Yoon, S.W. Challenges of advanced wafer level packaging technology: Cost-effectiveness, integration and scalability. In Proceedings of the Electronics Packaging Technology Conference (EPTC), Singapore, 5–7 December 2012. [Google Scholar]
  8. Tian, R.; Hang, C.; Tian, Y.; Xu, J. Brittle fracture of Sn-37Pb solder joints in-duced by enhanced intermetallic compound growth under extreme temperature changes. J. Mater. Process. Technol. 2019, 268, 1–9. [Google Scholar] [CrossRef]
  9. Shen, Y.-A.; Zhou, S.; Li, J.; Tu, K.; Nishikawa, H. Thermomigration induced microstructure and property changes in Sn-58Bi solders. Mater. Des. 2019, 166, 107619. [Google Scholar] [CrossRef]
  10. Tu, K.N.; Hsiao, H.Y.; Chen, C. Transition from flip chip solder joint to 3D IC microbump: Its effect on microstructure anisotropy. Microelectron. Reliab. 2013, 53, 2–6. [Google Scholar] [CrossRef]
  11. Baek, S.-M.; Park, Y.; Oh, C.; Chun, E.-J.; Kang, N. Modeling and Experimental Verification of Intermetallic Compounds Grown by Electromigration and Thermomigration for Sn-0.7Cu Solders. J. Electron. Mater. 2018, 48, 142–151. [Google Scholar] [CrossRef]
  12. Feng, J.; Hang, C.; Tian, Y.; Wang, C.; Liu, B. Effect of electric current on grain orientation and mechanical properties of Cu-Sn intermetallic compounds joints. J. Alloys Compd. 2018, 753, 203–211. [Google Scholar] [CrossRef]
  13. Chen, W.-J.; Lee, Y.-L.; Wu, T.-Y.; Chen, T.-C.; Hsu, C.-H.; Lin, M.-T. Effects of Electrical Current and External Stress on the Electromigration of Intermetallic Compounds Between the Flip-Chip Solder and Copper Substrate. J. Electron. Mater. 2017, 47, 35–48. [Google Scholar] [CrossRef]
  14. Chen, C.; Tong, H.; Tu, K. Electromigration and Thermomigration in Pb-Free Flip-Chip Solder Joints. Annu. Rev. Mater. Res. 2010, 40, 531–555. [Google Scholar] [CrossRef]
  15. Chen, C.; Hsiao, H.-Y.; Chang, Y.-W.; Ouyang, F.; Tu, K. Thermomigration in solder joints. Mater. Sci. Eng. R Rep. 2012, 73, 85–100. [Google Scholar] [CrossRef]
  16. Zhang, Z.H.; Cao, H.J.; Chen, H.T. Formation mechanism of a cathodic serrated inter-face and voids under high current density. Mater. Lett. 2018, 211, 191–194. [Google Scholar] [CrossRef]
  17. Cui, Z.; Fan, X.; Zhang, Y.; Vollebregt, S.; Fan, J.; Zhang, G. Coupling model of electromigration and experimental verification-Part I: Effect of atomic concentration gradient. J. Mech. Phys. Solids 2023, 174, 105257. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Zhang, J.; Wang, Y.; Fang, Y. Effect of Grain Structure and Ni/Au-UBM Layer on Electromigration-Induced Failure Mechanism in Sn-3.0Ag-0.5Cu Solder Joints. Micromachines 2022, 13, 953. [Google Scholar] [CrossRef] [PubMed]
  19. Kim, J.M.; Jeong, M.H.; Yoo, S.; Park, Y.B. Effects of surface finishes and current stressing on interfacial reaction characteristics of Sn-3.0Ag-0.5Cu solder bumps. J. Electron. Mater. 2012, 41, 791–799. [Google Scholar] [CrossRef]
  20. Zhao, Y.H.; Zhang, Y.X.; Xv, Y.J.; Yang, L.H. Simulation of electromigration failure and optimization analysis of Sn0.7Cu lead-free solder joints. Electron. Compon. Mater. 2021, 35, 2133–2138. [Google Scholar]
  21. Hsu, H.C.; Ju, S.W.; Lu, J.R.; Chang, H.S.; Wu, H.H. Electromigration analysis and electro-thermo-mechanical design for semiconductor package. In Proceedings of the International Conference on Electronic Packaging Technology & High Density Packaging, Beijing, China, 10–13 August 2009. [Google Scholar]
  22. Gao, Z.; Dong, C.; Shang, S.; Huang, M.; Ma, H.; Wang, Y. Pronounced electromigration of GaInSn/Cu interconnects under super low critical current density. Mater. Lett. 2021, 300, 130137. [Google Scholar] [CrossRef]
  23. Jiang, Y.; Li, H.; Chen, G.; Mei, Y.; Wang, M. Electromigration behavior of Cu/Sn3.0Ag0.5Cu/Cu ball grid array solder joints. J. Mater. Sci. Mater. Electron. 2019, 30, 6224–6233. [Google Scholar] [CrossRef]
  24. Passage, J.M.; Azhari, N.; Lloyd, J.R. Stress Migration Followed by Electromigration Reliability Testing. In Proceedings of the 2019 IEEE International Reliability Physics Symposium (IRPS), Monterey, CA, USA, 31 March 2019–4 April 2019. [Google Scholar]
  25. Su, M.; Gao, R.; Li, J.; Cao, L. Simulation analysis of system-level fan-out package based thermal management. Power Electron. Technol. 2018, 52, 4. [Google Scholar]
  26. Wang, J.; Yuan, Y.; Zhu, Y.; Zhang, Z. Study on the fatigue life of solder balls in fan-out BGA packages. Electron. Prod. Reliab. Environ. Test. 2022, 040. [Google Scholar]
  27. Cheng, T.; Lv, G.; Li, Y.; Yun, H.; Zhang, L.; Deng, Y.; Nan, J. Low Dielectric Polyimide/Fluorinated Ethylene Propylene (PI/FEP) Nanocomposite Film for High-Frequency Flexible Circuit Board Application. Macromol. Mater. Eng. 2021, 306, 2100086. [Google Scholar] [CrossRef]
  28. Lee, C.C.; Wang, C.W.; Chen, C.Y. Comparison of Mechanical Modeling to Warpage Estimation of RDL-First Fan-Out Panel-Level Packaging. IEEE Trans. Compon. Packag. Manuf. 2022, 12, 1100–1108. [Google Scholar] [CrossRef]
  29. Chen, F.; Chen, S.; Fu, Z.; Huang, Y.; Qin, F.; An, T. Prediction of electromigration lifetime of copper pillar bumps in ceramic packaging device. In Proceedings of the 2020 21st International Conference on Electronic Packaging Technology (ICEPT), Guangzhou, China, 12–15 August 2020. [Google Scholar]
  30. Black, J.R. Electromigration—A brief survey and some recent results. IEEE Trans. Electron Devices 1969, 16, 338–347. [Google Scholar] [CrossRef] [Green Version]
  31. Prigogine, I. Introduction to Thermodynamics of Irreversible Processes; Interscience: New York, NY, USA, 1967. [Google Scholar]
  32. Yeo, A.; Lee, C.; Pang, J.H. Flip chip solder joint reliability analysis using viscoplastic and elastic-plastic-creep constitutive models. IEEE Trans. Compon. Packag. Technol. 2006, 29, 355–363. [Google Scholar] [CrossRef]
  33. Basaran, C.; Lin, M. Electromigration induced strain field simulations for nanoelectronics lead-free solder joints. Int. J. Solids Struct. 2007, 44, 4909–4924. [Google Scholar] [CrossRef] [Green Version]
  34. Su, P.; Li, L.; Lai, Y.S.; Chiu, Y.T.; Kao, C.L. A comparison study of electromigration performance of Pb-free flip chip solder bumps. In Proceedings of the Electronic Components & Technology Conference, San Diego, CA, USA, 26–29 May 2009. [Google Scholar]
  35. Liu, Y.X.; Gusak, A. Fast prediction of electromigration lifetime with modified mean-time-to-failure equation. Mater. Lett. 2022, 325, 132880. [Google Scholar] [CrossRef]
  36. Rahim, M.K.; Zhou, T.; Fan, X.; Rupp, G. Board Level Temperature Cycling Study of Large Array Wafer Level Package. In Proceedings of the Electronic Components & Technology Conference, San Diego, CA, USA, 26–29 May 2009. [Google Scholar]
  37. Dudek, R.; Rzepka, S.; Dobritz, S.; Doring, R.; Keybig, K.; Wiese, S.; Michel, B. Fatigue Life Prediction and Analysis of Wafer Level Packages with SnAgCu Solder Balls. In Proceedings of the Electronics System-integration Technology Conference, Dresden, Germany, 5–7 September 2006. [Google Scholar]
  38. Cui, H.; Tian, W.; Zhang, Y.; Chen, Z. The Study of the Reliability of Complex Components during the Electromigration Process. Micromachines 2023, 14, 499. [Google Scholar] [CrossRef]
  39. Liu, B.; Tian, Y.; Qin, J.; An, R.; Zhang, R.; Wang, C. Degradation behaviors of micro ball grid array (μY), solder joints under the coupled effects of electromigration and thermal stress. J. Mater. Sci. Mater. Electron. 2016, 27, 11583–11592. [Google Scholar] [CrossRef]
  40. Yeh, Y.T.; Chou, C.K.; Hsu, Y.C.; Chen, C.; Tu, K.N. Threshold current density of electromigration in eutectic SnPb solder. Appl. Phys. Lett. 2005, 86, 5451. [Google Scholar] [CrossRef] [Green Version]
Figure 1. DDR3 micro-assembly physical diagram: (a) top component; (b) bottom component.
Figure 1. DDR3 micro-assembly physical diagram: (a) top component; (b) bottom component.
Micromachines 14 01255 g001
Figure 2. Schematic diagram of DDR simplified model.
Figure 2. Schematic diagram of DDR simplified model.
Micromachines 14 01255 g002
Figure 3. Micro-component Equivalent Model.
Figure 3. Micro-component Equivalent Model.
Micromachines 14 01255 g003
Figure 4. Overall meshing.
Figure 4. Overall meshing.
Micromachines 14 01255 g004
Figure 5. Current density at current aggregation: 593.60 A/cm2.
Figure 5. Current density at current aggregation: 593.60 A/cm2.
Micromachines 14 01255 g005
Figure 6. Temperature at hot spot: 461.36 K.
Figure 6. Temperature at hot spot: 461.36 K.
Micromachines 14 01255 g006
Figure 7. The 1/4 micro-component model: (a) three-dimensional diagram; (b) divided grid diagram.
Figure 7. The 1/4 micro-component model: (a) three-dimensional diagram; (b) divided grid diagram.
Micromachines 14 01255 g007
Figure 8. Cloud plot of the results at the current aggregation.
Figure 8. Cloud plot of the results at the current aggregation.
Micromachines 14 01255 g008
Figure 9. Resulting cloud at hotspot.
Figure 9. Resulting cloud at hotspot.
Micromachines 14 01255 g009
Figure 10. Sub-model bump coordinates.
Figure 10. Sub-model bump coordinates.
Micromachines 14 01255 g010
Figure 11. Hot spot temperature at 328.15 K and 3.5 A/cm2 is 424.26 K.
Figure 11. Hot spot temperature at 328.15 K and 3.5 A/cm2 is 424.26 K.
Micromachines 14 01255 g011
Figure 12. Hot spot temperature at 398.15 K and 3.5 A/cm2 is 485.50 K.
Figure 12. Hot spot temperature at 398.15 K and 3.5 A/cm2 is 485.50 K.
Micromachines 14 01255 g012
Figure 13. MTTF results comparison curve.
Figure 13. MTTF results comparison curve.
Micromachines 14 01255 g013
Table 1. Dimensional parameters.
Table 1. Dimensional parameters.
ComponentsSize (mm3)
Si7.805 × 5.105 × 0.2
PI11.3 × 8.95 × 0.3
RDL11.3 × 8.95 × 0.063
Table 2. Material parameters.
Table 2. Material parameters.
MaterialsSiCuPolyimideSn63Pb37
Density (kg/m3)2330894015508425
Modulus of elasticity (Pa)130 × 109255 × 10922 × 10934 × 109
Coefficient of thermal expansion (ppm/K)3.611711.52.39
Poisson’s ratio0.280.220.280.38
Thermal conductivity (W/m·K)1194000.25In Table 3
Resistivity (Ω·m)2.52 × 10−41.72 × 10−81014.51.46 × 10−4
Specific heat capacity (J/(kg·°C))7123851150227.2
Relative dielectric constant11.713.2\
Table 3. Thermal conductivity of Sn63Pb37.
Table 3. Thermal conductivity of Sn63Pb37.
Temperature (K)293307.45321.9336.35350.8365.25379.75394.15408.65423.05
Thermal conductivity (W/mK)52.0251.8951.7651.6351.551.3751.2451.1150.9850.85
Table 4. Electrothermal loading conditions.
Table 4. Electrothermal loading conditions.
Current Density (A/cm3)3.53.5444.54.555
Temperature (K)328.15398.15328.15398.15328.15398.15328.15398.15
Table 5. Electromigration parameters of 63Sn37Pb [32,33].
Table 5. Electromigration parameters of 63Sn37Pb [32,33].
Electromigration ParametersNumerical ValueSymbolsUnit
activation energy1EaeV
Effective charge number−33Z*/
Effective self-diffusion coefficient3.14 × 10−5D0m2/s
Heat transfer0.0094Q*eV
Atomic volume 2.48 × 10−29Ωm3
Initial resistivity1.55 × 10−7R0Ω-m
Boltzmann’s constant1.38 × 10−23KBJ/K
Table 6. Experimental results and coefficient calculation for different working conditions.
Table 6. Experimental results and coefficient calculation for different working conditions.
Current Density (A/cm2)Temperature (K)Current Density Max (A/cm2)Temperature Max (K)MTTF (h)
3.5328.15627.37423.747239.30
394.31424.2617,721.10
3.5398.15627.4448.971553.34
457.38449.242877.76
4328.15716.99459.89643.84
450.64460.621565.95
4398.15717.03485.01174.23
522.72485.49320.18
4.5328.15806.62500.8864.51
506.87501.82156.43
4.5398.15806.66516.831.59
507517.9376.16
5328.15896.24546.567.54
563.30547.7818.20
5398.15896.29569.453.21
563.34570.847.73
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, W.; Gao, R.; Gu, L.; Ji, H.; Zhou, L. Three-Dimensional Integrated Fan-Out Wafer-Level Package Micro-Bump Electromigration Study. Micromachines 2023, 14, 1255. https://doi.org/10.3390/mi14061255

AMA Style

Tian W, Gao R, Gu L, Ji H, Zhou L. Three-Dimensional Integrated Fan-Out Wafer-Level Package Micro-Bump Electromigration Study. Micromachines. 2023; 14(6):1255. https://doi.org/10.3390/mi14061255

Chicago/Turabian Style

Tian, Wenchao, Ran Gao, Lin Gu, Haoyue Ji, and Liming Zhou. 2023. "Three-Dimensional Integrated Fan-Out Wafer-Level Package Micro-Bump Electromigration Study" Micromachines 14, no. 6: 1255. https://doi.org/10.3390/mi14061255

APA Style

Tian, W., Gao, R., Gu, L., Ji, H., & Zhou, L. (2023). Three-Dimensional Integrated Fan-Out Wafer-Level Package Micro-Bump Electromigration Study. Micromachines, 14(6), 1255. https://doi.org/10.3390/mi14061255

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop