key: cord-0432125-yle49u42 authors: Guo, Zhenhua; Ren, Xintong; Jiang, Jiahua; Liu, Haoyang; Huo, Yongjun; Zhao, Xiuchen; Tu, K. N.; Liu, Yingxia title: The design of quaternary eutectic solder by machine learning date: 2021-10-09 journal: nan DOI: nan sha: 54f96b7b248926f0881daae212fc5a2571c28faf doc_id: 432125 cord_uid: yle49u42 In this paper, we obtain a Sn-Bi-In-Pb quaternary near eutectic alloy composition from machine learning model. The eutectic points and the alloy composition were evaluated and continuously improved by experimental input. The actual composition is near the result given by machine learning. We conclude that the application of machine learning in solder design has shown the potential to overcome the challenge in searching for the next generation eutectic solders, which will have a broad impact on the industry. Due to the COVID-19 virus pandemic, the trend of distance teaching, distance medicine, home office, and on-line meeting has increased greatly the need of advanced consumer electronic products, demanding more function, smaller form factor, larger memory, faster and greater rate of data transmission, cheaper cost, and superb reliability. To sustain the trend, advanced packaging technologies with high complexity and integration, such as vertically stacked packaging-on-packaging (POP) and threedimensional integrated circuits (3D IC) have already shown their impact to our society [1] [2] [3] . In these new packaging structures, multiple reflows of solder joints will be required to achieve the vertical stacking. To avoid the remelting of the solders, there is a need for solders with hierarchical melting points. In addition, as a variety of devices are developing, there is a need for a variety of different melting point solders in corresponding packaging strategies. For example, for bio-medical devices and neuromorphic devices, we need to develop low melting point interconnect materials, 3 since these devices can only sustain a low temperature during manufacturing processes. So far, we haven't known any appropriate solders for these applications. At present, we have the high temperature solder, Sn70Au alloy, with a melting point of 280 ℃ and a reflow temperature of about 300 ℃. We also have the medium temperature solder SAC 305 alloy with a melting point of 217 ℃ and a reflow temperature of about 250 ℃. For the low temperature solder, we have SnBi eutectic solder with a melting point of 139 ℃ and a reflow temperature of about 150 ℃ [4] [5] [6] [7] . However, the brittleness of SnBi eutectic solder limits its application in most mobile electronic products. Therefore, it is crucial to develop a solder with a low melting point and appropriate mechanical properties. The major reason that we need eutectic alloys for solder is to have a single melting point. Thus, thousands of solder joints can melt and solidify at the same time. Therefore, our goal is to search for a low melting point quaternary, or even quintenary eutectic solder. Because the pool of candidates is too big and it's impossible to do experiments of each composition one by one. This is the reason we use machine learning (ML) to solve the problem. Machine learning can build models and find correlations to predict material properties, including composition, microstructure, physical and chemical properties. Recently, several applications have been successfully applied ML to find new HEA. Huang et al. developed three ML algorithms to predict HEA phase selection and found that atomic size difference and valence electron concentration were two key parameters [8] . Wen [10] . In this study, ML was applied to quantitatively determine the composition of QEA. With experimental data, we report a QEA with excellent mechanical properties. The method of finding this QEA can have a wide range of applications. II. Experimental To build up the database, we introduce the concept of eutectic degree, which means the deviation to the eutectic composition. The smaller the deviation, the higher the eutectic degree. Eutectic degree can be evaluated by primary phase mole fraction and melting point range. The primary phase molar fraction is the percentage of the solid phase solidified before the eutectic reaction. A small primary phase molar fraction means a high eutectic degree. The melting point range is the temperature interval between the beginning and the end of eutectic reaction. We consider the alloy to be eutectic when the melting point range is less than 3 ℃. Primary phase molar fraction and melting range can be obtained from phase diagram. Although the phase diagrams of multi-component alloys with four or more elements are too complicated to be obtained, we can have a "pseudo phase diagram" by the phase diagram calculation (CALPHAD) software JMAT. To apply ML in searching for the proper solder alloy, a database for the Sn-Bi-In-Pb alloy system was established by the CALPHAD calculation of JMAT. The reason that we choose Sn-Bi-In-Pb system is because these four elements are the most common low melting point elements. To ensure the accuracy of ML, a data base containing at least several hundred data sets is needed. In this study, 300 data sets were obtained using JMAT by randomly changing the content of each element in the Sn-Bi-In-Pb system. Each data set contains the elemental composition, the primary phase mole fraction, and the melting point range. The data in the database were split into 70% as training, 15% as testing and 15% as validation. Then the ML models used for predicting QEA with near-eutectic composition were trained. The learning model was constructed based on convolutional neural network (CNN). The architecture of the entire CNN model includes 1D convolution layers to extract features, dropout layers to avoid overfitting of the data, flatten layer and standard fully connected layers. Fig. 1 shows the pictorial elucidation of the CNN model. The input is the percentage of elements for each combination. The output is the melting range and melting temperature. Furthermore, to ensure the accuracy of the composition, the result computed by the CNN model was verified and improved by the experimental data. In order to verify the results of ML model, we did experiments. We first prepared the samples according to the composition given by ML results. SnBiInPb solder was prepared using high purity (>99.9%) Sn, Bi, In and Pb. The ingots melted completely at around 300 o C in a vacuum induction furnace under argon atmosphere. After melting, the samples were cooled down to room temperature in air. Pieces about 5-10 mg were cut from the bulk solder alloy. These pieces were analyzed by differential scanning calorimetry (DSC). If the result given by ML is not a perfect eutectic, we will adjust the composition until finding a perfect eutectic with the DSC endothermic and exothermic 7 peaks both having only one peak, which we define it as QEA. The solder alloy bulk was cut into tiny pieces. These pieces were placed on a hotplate and reflowed on Cu under flux at 100 o C for different length of time. After soldering, the samples were removed from the hotplate and cooled to room temperature. The samples were mounted using epoxy resin and then cross-sectioned and polished. Besides, original samples without reflow process were made as comparison. The crosssectional interfaces of the polished samples were observed using scanning electron microscope (SEM). The elemental composition of solder and intermetallic compounds (IMC) were analyzed using energy dispersive X-ray spectroscope (EDX). After we obtained the QEA from both ML output result and the verifying experiment, we carried out shear test and tensile test to investigate the mechanical properties of the alloy. For the tensile test, QEA solder bulk was processed into tensile samples with length of 10 mm and thickness of 1.5 mm. The tensile samples were tested by Instron 5966 tensile test machine with an extension rate of 1ⅹ10 -3 s -1 . For the shear test, 4.5 ± 0.5 mg of diced QEA solder, Sn58Bi and pure Sn pieces were reflowed on 1 mm diameter circular Cu substrate. Subsequently the shear tests were performed using PTR-1100 shear test machine at room temperature with an elongation rate of 0.05 mm/s. Based on the ML model, a near eutectic composition of sample 1 (Sn-35Bi-47.5In-8 4Pb) was obtained. Three samples (sample 2, 3, 4) whose composition are Sn-30Bi-52.5In-4Pb, Sn-40Bi-42.5In-4Pb and Sn-37.5Bi-45In-4Pb, respectively, were made to search for a perfect eutectic composition. Fig. 2 (a) shows the strategy in composition modification. On the basis of sample 1, the content of Bi decreased by 5% and the content of In increased by 5%, we obtained sample 2. Then the content of In decreased by 5% and the content of Bi increased by 5%, we obtained sample 3. To confirm the results given by ML, we did the DSC tests for the samples, as shown in Fig. 2(b) to (e). Fig. 2(b) shows that the ML sample has only one endothermic peak, and two exothermic peaks in the DSC curve, indicating that the alloy is almost eutectic, but there is phase transformation during the solidification. In Fig. 2(c) , the endothermic peak of sample 2 is irregular, and there are more exothermic peaks than sample 1, showing that the composition adjustment in this direction does not meet our requirements. In Fig. 2(d) , the endothermic peak of sample 3 is multiple, but the exothermic peak becomes single, just contrary to sample 1. Therefore, we believe the eutectic composition should be between the composition of sample 1 and sample 3. As a result, we obtained sample 4. As shown in Fig. 2 (e), sample 4 has only one peak for both exothermic and endothermic process, so it can be considered as an QEA alloy. The melting point of this alloy is 60.81 o C. Sample 2 is far from the desired eutectic composition, so we didn't do further investigation on that sample. Fig. 3 shows the SEM pictures of original sample 1, 3 and 4, which didn't go through reflow. The lamellar eutectic structure can be observed. There are mainly two phases in the three alloys, and we call them light phase and dark phase. Distinctively, three phases, including a darker phase can be observed (circled in red) in Fig. 3 (d) of sample 3(Sn-40Bi-42.5In-4Pb). there are mainly two phases in these alloys. But a darker phase can be found (circled in red) in Fig. 4 (b) of sample 1(Sn-35Bi-47.5In-4Pb). There are mainly two phases, Bi-rich phase, and Sn-rich phase. In and Pb atoms are evenly distributed in the sample. Table 1 shows the EDX results of each phase in the samples. We also did XRD to further confirm the EDX results, as shown in Fig. 6 . The XRD result is complicated to analyze, and there are some weak peaks can't be indexed, which we tend to believe they are from the Sn-rich phase. We will discuss this point in Section 4.2. According to the EDX and XRD results, BiIn2 and Bi3In5 are the main components in Bi-rich phase, with some Sn and Pb atoms substituting In and Bi atoms. 12 The original sample 1 has the BiIn2 phase, because the content of In is the largest in sample 1, and that phase changes to two phases of BiIn2 and Bi3In5 after reflow. The original sample 3 has the Bi3In5 phase, and two Sn-rich phases before reflow. Sample 4 seems to have stable Bi3In5 and Sn-rich phases, which don't go through phase transformation during reflow. Tab. 1. The EDX results of sample 1, 3, and 4 before and after reflow. Fig. 6 . The XRD spectrum of (a) sample 1, (b) sample 3, and (c) sample 4 after reflow. We did TEM to investigate the Sn-rich phase, as shown in Fig.7 . Fig. 7(b) and (c) 15 show the electron diffraction pattern of sample 4, the corresponding areas are circled in Fig. 7(a) . The indexed result indicates that the phase has FCC structure, with the average lattice parameter as a = b = c = 0.478 nm. The TEM-EDX result show that the sample is mainly a mixing structure of Sn and In. Fig. 7 . The TEM diffraction pattern of sample 4 (Sn-37.5Bi-45In-4Pb). The tensile test results are shown in Fig. 8 , and the average elongations of samples 16 Fig. 9 . Different from the other four samples, a tip can be observed for the QEA sample after tensile fracturing. We also did the shear test. The shear strength of eutectic sample 4, Sn58Bi and pure 17 Sn solder joints is measured as 26.14  3.71 MPa, 43.65  2.52 MPa and 48.32  1.74 MPa, respectively. The measured shear strength of Sn58Bi and Sn solder joints is in the range of reported data [15] [16] [17] [18] [19] . Compared with Sn58Bi and Sn solder joint, the QEA solder has a lower shear strength. IV. Discussion With the DSC and the microstructure characterization results by SEM and XRD, we have a deeper understanding in the relationship between the microstructure and melting point. There are two endothermic peaks in sample 1 and two exothermic peaks in sample 3, which are due to the formation of a new phase. In sample 1, due to the excess of In atoms, the originally formed phase is BiIn2, and after reflow, it transforms to BiIn2 and Bi3In5. In sample 3, due to the excess of Bi atoms, the originally formed phase is Bi3In5. In addition to that phase, there is some high-In Sn-rich and low-In Sn-rich. After reflow, the high-In phase and low-In phase merges to a more stable Sn-rich phase. Compared to sample 1 and sample 3, sample 4 has one endothermic peak and one exothermic peak, with only 3 o C melting range, which can be regard as a eutectic alloy. The alloy has stable Bi3In5 and Sn-rich phase, both in melting and solidification processes. Thus, in the eutectic structure, the two phases are Bi3In5 and Sn-rich phase. The Sn-rich phase has been investigated and is believed to be a high entropy phase. 18 In the XRD result, in addition to the peaks from BiIn2 and Bi3In5, we have some peaks that are unable to be indexed. One interesting point is that, the undercooling of this eutectic phase in solidification is around 3 o C. Compared with SAC 305, Sn37Pb and Sn58Bi, the reported undercooling for the alloys is about 10 to 20 o C [20] [21] [22] , and we can conclude the undercooling of the QEA is quite small, implying the nucleation of the two phases in the eutectic alloy is easy. Bi3In5 is an IMC phase, with a melting point around 88 o C [23] . The high entropy phase of SnBiInPb, may reduce the nucleation barrier, and thus reduce the undercooling. Because of the small undercooling, the eutectic structure is quite different from the traditional eutectic lamellar structure. As shown in Fig. 3 (e) 19 and (f), the eutectic microstructure is "cloud-like", rather than "lamellar-like". This unique eutectic structure may contribute to the high plasticity that we are going to discuss below. 4.2 High plasticity of the QEA As shown in Fig. 8(a) However, we noted that sample 1 and sample 3 also have a low melting point, around 60 o C, but their elongation is much lower than QEA. One possible explanation is that, the QEA has a very small undercooling, leading to the "cloud-like" eutectic microstructure, rather than lamellar. For sample 1 and sample 3, because they are not perfect eutectic, and the undercooling is larger than 5 o C, they have lamellar microstructure. The lamellar microstructure should have less sliding mobility than "cloud-like" microstructure, due to anisotropic. Another possible reason is that the QEA sample has more weight of SnBiInPb high entropy phase, so that it can play a key role in the whole alloy. A detailed research on the high elongation of QEA will be carried out later. The small undercooling of QEA is meaningful in the design of solders. This is because a large undercooling means a large nucleation barrier, which is undesirable. 20 For a small undercooling, the difference is small on which phase nucleates first. Technically, the advantage of small undercooling is that we can have all the solder joints solidified almost at the same time. Another advantage is that the relationship between small undercooling, cloud-like microstructure, and plasticity may help us to improve the ductility of solder alloy in the application of mobile devices. In this work, our data base is not big enough so far, and we have only 300 data in the set of Sn, Bi, In and Pb. The result given by ML is not perfect, however, it's quite near our goal. By trials with limited times, we can find the perfect eutectic solder alloy. Although this eutectic SnBiInPb solder may be too soft in solder application, nevertheless, it's promising to be applied to the solid-liquid diffusion bonding process [24, 25] . We are now working on a much bigger data base, and the developed ML model can predict the eutectic composition from almost any Sn-based element combinations. We are confident that ML can overcome the challenge of searching for the future eutectic solders, which seems to be impossible by experimental search in the old days. We applied ML to the search for a new QEA. The composition given by the ML model is Sn-35Bi-47.5In-4Pb. Experimental data confirm that this alloy is near eutectic. By modifying the composition to Sn-37.5Bi-45In-4Pb, we obtained a QEA. This alloy has tensile strength of 22.55 MPa and elongation of 1.17. The high elongation is explained by the "cloud-like" microstructure. This soft and low melting point QEA may be applied to the solid-liquid diffusion bonding in packaging technology. This work has shown the potential to design the future solder alloys by ML. Warpage and thermal characterization of fan-out wafer-level packaging Interconnect quality and reliability of 3D packaging, 3D Microelectronic Packaging Scaling effect of interfacial reaction on intermetallic compound formation in Sn/Cu pillar down to 1 μm diameter Low melting point solders based on Sn, Bi, and In elements A review: On the development of low 22 melting temperature Pb-free solders A brief review on high-temperature, Pb-free dieattach materials A highentropy alloy as very low melting point solder for advanced electronic packaging Machine-learning phase prediction of highentropy alloys Machine learning assisted design of high entropy alloys with desired property Uncovering the eutectics design by machine learning in the Al-Co-Cr-Fe-Ni high entropy system Tensile Properties of Sn-Bi Lead-Free Solder Alloys, Solid State Phenomena Correlation between microstructure and mechanical properties of Sn-Bi-X solders Effects of cooling speed on microstructure and tensile properties of Sn-Ag-Cu alloys Tensile properties of Sn-Ag based lead-free solders and strain rate sensitivity Characterization of eutectic Sn-Bi solder joints The effect of aging on microstructure, room temperature deformation, and fracture of Sn-Bi/Cu solder joints Elevated temperature aging of solder joints based on Effects on joint microstructure and shear strength Influence of reflow and thermal aging on the shear strength and fracture behavior of Sn-3.5 Ag solder/Cu joints Mechanical properties of Pb-free SnAg solder joints Evolution of microstructure, thermal and creep properties of Nidoped Sn-0.5 Ag-0.7 Cu low-Ag solder alloys for electronic applications Determination of the melting and solidification characteristics of solders using differential scanning calorimetry The influence of Ni and Zn additions on microstructure and phase transformations in Sn-0 Cu/Cu solder joints Bi-In binary phase diagrams Sn-3.0 Ag-0.5 Cu/Sn-58Bi composite solder joint assembled using a low-temperature reflow process for PoP technology The application of Bi-based solders for low temperature reflow to reduce cost while improving SMT yields in client computing systems