Thermodynamic Assessment of Molten Bix-Sn1−x (x = 0.1 to 0.9) Alloys and Microstructural Characterization of Some Bi-Sn Solder Alloys

Abstract

Properties such as lower melting temperature, good tensile strength, good reliability, and well creep resistance, together with low production cost, make the system Bi-Sn an ideal candidate for fine soldering in applications such as reballing or reflow. The first objective of the work was to determine the thermodynamic quantities of Bi and Sn using the electromotive force measurement method in an electrolytic cell (Gibbs’ enthalpies of the mixture, integral molar entropies, and the integral molar excess entropies were determined) at temperatures of 600 K and 903 K. The second objective addressed is the comprehensive characterization of three alloy compositions that were selected and elaborated, namely Bi25Sn75, Bi50Sn50, and Bi75Sn25, and morphological and structural investigations were carried out on them. Optical microscopy and SEM-EDS characterization revealed significant changes in the structure of the elaborated alloys, with all phases being uniformly distributed in the Bi50Sn50 and Bi75Sn25 alloys. These observations were confirmed by XRD and EDP-XRFS analyses. Diffractometric analysis reveals the prevalence of metallic Bi and traces of Sn, the formation of the Sn_0.3Bi_0.7, Sn_0.95Bi_0.05 compounds, and SnO and SnO₂ phases.

1. Introduction

The intensification of research in the field of lead-free soldering alloys started since 1 July 2006, when the new directives of the European Union (European Union Waste Electrical and Electronic Equipment Directive—WEEE and Restriction of Hazardous Substances Direction—RoHs) came into force and the European electronics industry had to be lead-free. Lead and lead compounds represent one of the 17 main categories of chemical substances and produce major problems that can threaten people’s lives or cause serious environmental imbalances [1,2].

Several conditions must be met for a soldering alloy to be appropriate: melting range that is suitable for the service temperature range; mechanical properties that are compatible with service conditions; metallurgical compatibility with surrounding metallization; low rate of intermetallic compound (IMC) formation at the service temperature; and acceptable wettability on surrounding metallization [3,4,5,6,7,8,9,10].

The demand for this category of alloys is constantly increasing due to the growth of the electronic equipment production sector, starting from applications for the common consumer (TV, PC, freezer, washing machine, mobile phone, notebook, etc.) to high-end applications such as in the field of satellite communications. Even if the replacement of traditional alloys with lead content is difficult, with their properties being difficult to match or even surpass, their drawbacks were further enhancing the properties of lead-free alloy systems by adding small amounts of alloying elements such as Bi, Cu, In, Ag, Al, Ga, Sb, Cr, Ni, and Ge to develop ternary and even quaternary Pb free systems [11,12,13,14,15].

Newly developed lead-free alloy systems have been discovered that can successfully replace them [16], such as those in the binary Sn-Cu, Sn-Ag, Sn-Bi, Sn-Zn, Sn-In, and Sn-Sb systems or the ternary Sn-Zn-Bi and Sn-Ag-Cu systems (SAC solders), with acceptable properties (melting temperature, mechanical properties, microstructure, wettability, solder ability, reliability, and cost), thus offering a positive trend in terms of finding new variants of accessible lead-free soldering alloys with properties similar to those prohibited [6,7,8,9,10,11,13]. Other investigated alloys are from the Bi-Ag, Bi-Sb, Bi-Sn, Bi-Zn, or Bi-Sn-Zn systems, where further investigations are necessary for better understanding. The addition of Bi decreases the melting point of the system and improves the wettability, and up to a certain extent, the mechanical properties could also be improved [7,17,18,19,20].

In both the conventional electronic industry (the industry that produces stand-alone equipment) and industrial electronics (equipment and components incorporated into products), the Bi-Sn alloy system is frequently used worldwide. The electronic industry is a major consumer of special materials such as semiconductors, insulators, and metals used in the construction of various electronic devices.

According to current research practice and the international research approach, the trend in materials science is a fundamental comprehensive approach to understanding material properties as accurately as possible, especially the physico-chemical characteristics essential for the intended applications. For this reason, we are increasingly moving toward a fundamental approach to material properties based on fine structure in correlation with other influencing factors that contribute to the measured values of those properties. In this regard, techniques and methods for investigating material structures at the atomic level promote exceptional scientific results, leading to the development of advanced materials for all practical areas of interest (e.g., metallurgy, medicine, electronics, etc.) [21,22,23,24,25,26,27,28,29,30]. The phenomenology associated with quantities such as Gibbs’ enthalpy, thermodynamic activity, entropy, free enthalpy, etc., must be seen from the perspective of material structure (phase stability analysis or phase transformations of an alloy system). These data are essential because they can provide us with valuable information, even if, in most cases, they are based on empirical calculation methods and on the problems of developing alloy systems or refining them, while knowing what problems face the development of the category of alloys for soldering.

The novelty of this research is the performance of detailed calculations regarding certain thermodynamic functions based on new research carried out on alloys from the Bi-Sn system by reviewing the data at an established temperature (600 K) as well as at a temperature at which no studies have been carried out until now (903 K). Another element of novelty is the development and study of new compositions of Bi-Sn alloys for solders (Bi25Sn75, Bi50Sn50, and Bi25Sn75) and the performance of detailed investigations on the micro-structural changes that occurred following an increase in the content of bismuth in the alloy mass (the appearance of phases and compounds that are not quite studied in the specialized literature regarding these types of alloys).

2. Materials and Methods

2.1. Elaboration of the Bi-Sn Solder Alloys

The elaboration of soldering alloys was carried out in an electric induction furnace, with a capacity of 2 kg, equipped with a PID temperature controller 220 V model REX-C100 (RKC Instruments Inc., Tokyo, Japan), which has an accuracy of ±0.5% FS ±1 °C (274.15 K). The maximum temperature reached by the furnace is 1100 °C (1373 K), and its schematic presentation is shown in Figure 1.

Bi, with a purity of 99.6%, and Sn, with a purity of 99.8%, were used to prepare the Sn-Bi alloys. Sn and Bi primary were introduced in the preheated crucible. The furnace (S Y Electric Melting Furnace, Italy) temperature was set to 903 K, and it reached the temperature of 903 K in approximately 25 min. The elaboration process proceeded with the preparation of materials and raw materials, encompassing both primary preparation and charging preparation stages.

During primary preparation, a thorough assessment of the raw materials and auxiliary materials was conducted for qualitative and quantitative verification. Additionally, within this stage, the raw materials were dosed, sized, and weighed using an electronic balance. The obtained samples (Bi25Sn75, Bi50Sn50, and Bi75Sn25) underwent a process of melting and casting within graphite crucibles and molds.

The melting temperature was maintained constant for 10 min and subjected to a refining operation with 0.5% flux (3 portions), immersed in the molten bath with a graphite rod, and subsequently homogenized before pouring into the graphite mold. The alloy was cleaned of the slag formed and poured into the mold at a temperature of 623 K using a crucible handling tong. Refined ammonium chloride (NH₄Cl) served as a protective flux to shield the molten metal from oxidizing gases within the furnace environment.

2.2. Determination of Thermodynamic Functions of the Bi-Sn System

The aim of this research was to obtain extensive knowledge about the thermodynamic behavior of Bi-Sn alloys in the molten state over a large temperature; the selection of the values of 600 K and 903 K has a certain motivation.

First, tests were carried out at a temperature of 600 K to verify and correlate the thermodynamic data with those obtained by Magnus and Mannheimer [31] and, at the same time, to verify the correctness of the measurement process using the electromotive force method (EMF) which ensures greater precision and accuracy of the results.

Secondly, tests were carried out to determine the thermodynamic functions at a temperature of 903 K since there are no specialized studies carried out by other researchers so far, thus representing a challenge for the team.

The thermodynamic properties of the Bi-Sn system were also investigated at temperatures of 673 K, 723 K, and 793 K using EMF measurements [23,32]. Another study on the Bi-Sn binary system concerns the measurement of the integral enthalpy of the mixture carried out at temperatures of 767 K and 855 K by calorimetric methods [33] or vapor pressure at a temperature of 1373 K [34].

Bi-Sn solder alloys (Sn-5Bi, Sn-15Bi, Sn-30Bi, Sn-45Bi, and Sn-58Bi) were elaborated in a ceramic crucible at a temperature of 873 K, cooled to 573 K and then chilled and cast in a graphite mold and was observed that the Bi addition modify the microstructure, interfacial behaviors, and joint strength [35].

Melting of Bi-Sn alloys to 903 K can offer several advantages in various applications: improved thermal resistance (advantageous in applications where parts must withstand high temperatures), dimensional stability (not to deform or contract excessively during cooling and use), improved mechanical performance (corrosion resistance or durability), ease of processing (easier to machine and form into various shapes and configurations), and versatility in demanding applications (where a reliable and durable material is required). For electronics applications, Bi-Sn alloys can provide good adhesion between electronic components and the circuit board without affecting their operation.

The tests to determine the thermodynamic functions were carried out using the same furnace mentioned previously, which was used for the elaboration of the alloys, in order not to influence the results. At the same time, the use of this furnace allowed the introduction and fixing of an H-shaped electrochemical cell used for thermodynamic determinations (vessel height—15 cm; cylinder inner diameter—1.8 cm; the thickness of quartz walls –0.2 cm; the length of the W wires—20 cm; connection channel diameter—1.8 cm; and the height of the connecting channel—2.5 cm).

Following the cell preparation, i.e., introducing the pure Bi electrode, Bi-Sn electrode, and the electrolyte (a mixture of 17% moles KCl and 83% moles SnCl₂) into the cell compartments, the cell was placed into the furnace and slowly heated until the furnace had reached the envisaged temperatures, namely 600 K and 903K. Once the working temperature was stabilized, the W conductors were introduced into the working cell, put into contact with the melted electrodes, and coupled to the digital multimeter [7,22].

The electrochemical system was maintained at a constant temperature until the cell’s electromotive voltage became constant, which is considered the equilibrium electromotive voltage (E) of the cell. The equilibrium value of the electromotive voltage was reached after approximately 20 min, and five measurements were carried out on each molar fraction in the concentration range [0.1 to 1.0] [7].

Based on the measured values, we were able to calculate the thermodynamic activity as well as the thermodynamic activity coefficients for Bi [7,22]. Using the calculated values for Bi depending on concentration, the activity of the Sn component (as_n) could be determined at constant temperature based on the Gibbs-Duhem relationship [36,37].

After obtaining these thermodynamic functions, were calculated further the partial molar enthalpy of the mixture ($∆{\overline{H}}_i$), partial molar entropies of the mixture (($∆{\overline{S}}_i$), integral molar free energies ($∆G_i$) and the integral molar excess energies ($∆G_i^E$), integral molar enthalpies ($∆H_i$) and integral molar excess enthalpies ($∆H_i^E$) and integral molar entropies ($∆S_i$) and the integral molar excess entropies ($∆S_i^E$). The obtained results were processed using the Mathcad 15 analysis software and transposed in Excel for a better graphic display.

2.3. Sample Characterization

Conventional cast samples of binary Bi-Sn alloys underwent characterization through optical microscopy using the Optika model B383 MET optical microscope, (Optika, Ponteranica, Italy) which was equipped with a digital camera and image processing software. Prior to microscopic examination, the samples underwent metallographic preparation to reveal their structural properties. The resulting micrographs depicted variations in the microstructure corresponding to the percentage composition of bismuth within each sample.

The microstructural characterization was conducted through scanning electron microscopy (SEM) FEI QUANTA 250 (FEI Company, Hillsboro, OR, USA), in high vacuum mode, backscattered electron techniques (BSE) using an Angular Backscattered Detector. Point analysis, elemental mapping, and semi-quantitative analysis were performed through energy dispersive X-ray spectroscopy (EDS) using the EDAX Element EDS Analysis System consisting of ELEMENT Silicon Drift Detector Fixed and ELEMENT EDS Analysis Software Suite (EDAX APEX™, DigitalMicrograph 3.6). For microstructural investigations, metallographic specimens were prepared by sectioning, mounting, grinding, fine grinding, polishing, and etching with a H₂O₂ (30%):H₂O mixture.

For comprehensive elemental and structural characterization of Bi-Sn alloys, X-ray fluorescence spectrometry (XRFS) and X-ray diffractometry (XRD) were identified as the most effective analytical approaches. These instrumental techniques require minimal sample preparation while providing extensive information on elemental and phase compositions. The Energy Dispersive Polarization X-ray Fluorescence (EDP-XRFS) method utilized for the characterization of Bi-Sn alloys employed the XEPOS spectrometer (SPECTRO Analytical GmbH). The obtained spectra were processed using the analytical program TurboQuant-Alloys.

X-ray diffraction investigations were conducted using the Bruker D8 Advance X-ray diffractometer in a Bragg-Brentano configuration for the interval 2θ = 10…90° (Δ_2θ = 0.02°, λ_Cu = 1.5415 Å). The structural data required for comparison were taken from the ICDD (International Center for Diffraction Data), formally JCPDS—(Joint Committee on Powder Diffraction Standards) files. The identification of the nature of the phase was performed by considering the positioning of the diffraction lines, respectively, the interplanar distances d_{h k l} associated with the diffraction maxima through the Bragg relation. Based on the intensities of the diffraction lines, the volume concentration of the phases in the analyzed sample is established in a semi-quantitative way, conducted using the Bruker-EVA program version 4.2.1.10, and was indexed based on the ICDD-PDF4+ database.

3. Results and Discussion

3.1. Thermodynamical Assessment

After calculating the activity coefficients, the partial molar enthalpy of the mixture was calculated based on the following relationship:

(1)$∆{\overline{H}}_i = ∆{\overline{G}}_i+T∆{\overline{S}}_i,$

where

$∆{\overline{H}}_i$—the partial molar enthalpy of the mixture for element i [J/mol];

$∆{\overline{G}}_i$—the partial molar free enthalpy of the mixture for element i [J/mol];

$∆{\overline{S}}_i$—the partial-molar entropy of the mixture for element i [J/mol].

T—temperature, [K].

• -. The values of the partial molar enthalpies of the mixture ($∆{\overline{H}}_i$) of Bi and Sn at 600 K and 903 K, calculated based on Equation (1), are presented in Table 1.

• -. The calculation of the values of partial molar entropies for the mixture was performed based on the following relationship (2):

(2)$∆\overline{S_i}=-\frac{\partial ∆\overline{G_i}}{\partial T}=-Rln{x}_i-\frac{\partial }{\partial T}\left(RTln{y}_i\right)=-Rln{x}_i-\frac{\partial ∆{\overline{G}}_i^E}{\partial T}$

where

$∆\overline{S_i}$—the partial molar entropy of the mixture for element i [J/mol];

$∆{\overline{G}}_i$—the partial molar free enthalpy of the mixture for element i [J/mol];

$∆{\overline{G}}_i^E$—the partial molar excess free enthalpy of the mixture for element i [J/mol];

R—the universal gas constant [J/mol·K];

T—temperature, [K].

The partial molar entropies of the mixture ($∆{\overline{S}}_i$) for Bi and Sn, calculated using Equation (2), are presented in Table 2.

• -. The integral molar free energies ($∆G_i$) and the integral molar excess energies ($∆G_i^E$) at 600K and 903K have been calculated using Equations (3) and (4), and they are presented in Table 3.

(3) $\mathrm{\Delta }{G}_i=\sum x_i\mathrm{\Delta }{\overline{G}}_i=x_{Bi}\mathrm{\Delta }{\overline{G}}_{Bi}+x_{Sn}\mathrm{\Delta }{\overline{G}}_{Sn}$

(4)$\mathrm{\Delta }{G}^E=\mathrm{\Delta }G-\mathrm{\Delta }{G}^{id}=RT(x_{Bi}ln{\gamma }_{Bi}+x_{Sn}ln{\gamma }_{Sn})$

where

$\mathrm{\Delta }{G}_i$—the molar-free enthalpy of the mixture for element i [J/mol];

$∆{\overline{G}}_i$—the partial molar free enthalpy of the mixture for element i [J/mol];

$\mathrm{\Delta }{G}^E$—the integral molar excess free enthalpy of the mixture for element i [J/mol];

$\mathrm{\Delta }G$—free energy variation [J/mol];

$\mathrm{\Delta }{G}^{id}$—ideal free energy variation [J/mol];

R—the universal gas constant [J/mol·K];

T—temperature, [K].

Table 3

The integral molar free and the integral molar excess energies at T = 600 K and T = 903 K.

	Temperature 600 K		Temperature 903 K
X Bi	$∆\mathit{G}$[J/mol]	$∆{\mathit{G}}^{\mathit{E}}$[J/mol]	$∆\mathit{G}$[J/mol]	$∆{\mathit{G}}^{\mathit{E}}$[J/mol]
0.0	0	0	0	0
0.1	−6425	−35,620	−602	−9820
0.2	−5802	−32,060	−1177	−11,850
0.3	−5292	−28,590	−1739	−14,450
0.4	−4906	−21,990	−2249	−17,400
0.5	−4679	−18,920	−2853	−20,550
0.6	−4672	−16,060	−3593	−23,830
0.7	−4946	−13,490	−4697	−27,180
0.8	−5674	−11,370	−6548	−30,570
0.9	−7157	−10,000	−10,530	−33,970
1.0	0	0	0	0

• -. The integral molar enthalpies ($∆H_i$) and integral molar excess enthalpies ($∆H_i^E$) at temperatures of 600 K and 903 K were calculated using Equations (5) and (6) and are presented in Table 4.

(5) $\mathrm{\Delta }{H}_i=\sum x_i\mathrm{\Delta }{\overline{H}}_i=x_{Bi}\mathrm{\Delta }{\overline{H}}_{Bi}+x_{Sn}\mathrm{\Delta }{\overline{H}}_{Sn}$

(6)$∆H^E=∆H^M-∆H^{id}=∆H^M$

where

$\mathrm{\Delta }{H}_i$—the integral molar enthalpies of the mixture [J/mol];

$\mathrm{\Delta }{\overline{H}}_i$—the partial molar enthalpies of the mixture [J/mol];

$∆H^E$—the variation excess enthalpy [J/mol];

$∆H^M$—the variation of the integral molar enthalpy [J/mol];

$∆H^{id}$—ideal free enthalpy variation [J/mol].

Table 4

The integral molar enthalpies and integral molar excess enthalpies of Bi and Sn at T = 600 K and T = 903 K.

	Temperature 600 K		Temperature 903 K
X Bi	$∆\mathit{H}$[J/mol]	$∆{\mathit{H}}^{\mathit{E}}$[J/mol]	$∆\mathit{H}$[J/mol]	$∆{\mathit{H}}^{\mathit{E}}$[J/mol]
0.0	0	0	0	0
0.1	−19,960	−18,140	−13,050	−17,130
0.2	−11,360	−14,020	−8486	−11,370
0.3	−8794	−11,410	−5593	−8019
0.4	−6838	−9188	−3628	−5280
0.5	−5248	−7031	−2278	−2672
0.6	−3908	−4865	−1342	−95
0.7	−2769	−2791	−698	2391
0.8	−1756	−870	−295	4616
0.9	−817	736	−68	6381
1.0	0	0	0	0

• -. The integral molar entropies ($∆S_i$) and the integral molar excess entropies ($∆S_i^E$) of the alloy Bi–Sn at temperatures of 600 K and 903 K, they were calculated using Equations (7) and (8) and are presented in Table 5.

(7) $\mathsf{\Delta }{S}_i=\sum x_i\mathsf{\Delta }{\overline{S}}_i=x_{Bi}\mathsf{\Delta }{\overline{S}}_{Bi}+x_{Sn}\mathsf{\Delta }{\overline{S}}_{Sn}$

(8)$\mathsf{\Delta }{S}^E=\mathsf{\Delta }S-\mathsf{\Delta }{S}^{id}=\frac{\mathsf{\Delta }{H}^E-\mathsf{\Delta }{G}^E}{T}+R(x_{Bi}ln{x}_{Bi}+x_{Sn}ln{x}_{Sn})=\frac{\mathsf{\Delta }{H}^E-\mathsf{\Delta }{G}^E}{T}$

where

$\mathsf{\Delta }{S}_i$—integral molar entropies of the mixture for element i [J/mol];

$∆\overline{S_i}$—the partial molar entropy of the mixture for element i [J/mol];

$\mathsf{\Delta }{S}_i^E$—integral excess molar entropy of the mixture for element i [J/mol];

$∆H^E$—the variation excess enthalpy [J/mol];

$\mathsf{\Delta }{G}^E$—the integral molar excess free enthalpy of the mixture for element i [J/mol];

R—the universal gas constant [J/mol·K];

T—temperature, [K].

Table 5

The integral molar entropies and the integral molar excess entropies at T = 600 K and T = 903 K.

	Temperature 600 K		Temperature 903 K
X Bi	$∆\mathit{S}$[J/mol]	$∆{\mathit{S}}^{\mathit{E}}$[J/mol]	$∆\mathit{S}$[J/mol]	$∆{\mathit{S}}^{\mathit{E}}$[J/mol]
0.0	0	0	0	0
0.1	−1.189	1.514	−1.175	1.484
0.2	0.922	5.082	0.823	4.784
0.3	3.850	8.929	2.954	7.726
0.4	6.186	11.781	5.881	10.690
0.5	7.469	13.232	6.841	12.467
0.6	7.632	13.228	6.326	12.563
0.7	6.726	11.805	5.765	10.679
0.8	4.890	9.051	3.909	8.611
0.9	2.606	5.309	1.806	4.629
1.0	0	0	0	0

3.2. Characterization of Bi-Sn Solder Alloys

3.2.1. Optical Microscopy

Due to the slow cooling of the samples, micro-cracks have been formed, and the grain size is large. This aspect can be improved through rapid or forced cooling of the samples and, of course, achieving a finer structure if the alloy is intended for use in “chip-off/chip-on” or “re-balling” techniques for repairing processors or memory devices in electronics. However, these techniques represent a small percentage of the electronics market.

According to the resulting microstructures, significant changes in the structure of the elaborated alloys are observed (Figure 2a–c). From a structure primarily based on the Sn matrix with few areas where rare zones with intermetallic constituents have been formed, we have transitioned to a much more balanced structure in this regard. This new structure exhibits a greater quantity of intermetallic constituents, metallic phases, and even the formation of compounds with polyhedral or tetragonal shapes.

In the microstructures, we can observe the α phase, which is rich in bismuth (gray areas), and the β phase, which is rich in tin (white areas). No voids or discontinuities are visible in the analyzed section. This is because, unlike other metals, Bi expands upon solidification. The presence of these two phases indicates significant diffusion during the cooling of the melt.

According to the phase diagram of Bi-Sn [19], the solubility limit of Bi in Sn is approximately 6 at. % for a temperature of 138 °C. Considering that the percentages of Bi used in the alloy compositions during fabrication greatly exceed this value, bismuth is observed to precipitate from the tin-rich matrix, resulting in the observed structures.

All phases are uniformly distributed in the Bi75Sn25 alloy structure without the formation of conglomerates and/or clusters that would cause heterogeneity in chemical composition. Instead, a fine and uniformly distributed multiphase structure is evident. Fine grains of Bi, some of them with rounded shapes, can be observed to have precipitated within the Sn matrix.

These observations are confirmed by XRD, ED-XRFS, and SEM-EDS investigations.

3.2.2. SEM-EDS Investigations

To highlight the effect of bismuth on the microstructure of Bi-Sn alloys (Bi25Sn75, Bi50Sn50, and Bi75Sn25), scanning electron microscopy (SEM) investigations coupled with energy dispersive spectrometry (EDS) micro-compositional analyzes were carried out.

Analyzing the SEM-EDS images in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, it is found that Bi has an obvious effect on the formation of the alloy structure and implicitly on the formation of phases and metallographic constituents. Bismuth distribution is homogeneous in all analyzed microstructures. In all studied alloys, Sn-Bi or Bi-Sn solid solutions were formed, composed of β-Sn phases in which fine Bi particles are embedded.

Several EDS analyses were carried out (Figure 4, Figure 6 and Figure 8) at the points marked on the microstructures to determine the micro-composition of the analyzed regions. They are presented together with the SEM microstructures.

Following the SEM-EDS analysis, the microstructural transformations in the studied alloys are quite visible. The darker areas represent b-Sn primary phases, and the lighter areas are eutectic structures between Sn-Bi and Bi-Sn. By increasing the bismuth content, the eutectic structure increases and the b-Sn phase decreases. The Bi25Sn75 alloy shows a typical hypoeutectic structure, with areas rich in bismuth; by equalizing the metal content, the Bi50Sn50 alloy shows a eutectic structure with bismuth spread homogeneously throughout the mass of the alloy, and in the last alloy, Bi75Sn25, by increasing the bismuth content, it is created a typical hypereutectic structure, with large bismuth grains between which different Bi-Sn phases are formed.

Also, from the microstructures in Figure 5 (Bi50Sn50 alloy) and Figure 7 (Bi75Sn25 alloy), the formation of fine needle-shaped beta phases embedded in the alpha phase dendritic matrix can be observed.

3.2.3. Structural Analysis through XRD and ED-XRFS

Three representative alloys were investigated, namely Bi25Sn75, Bi50Sn50, and Bi75Sn25. The Bi-Sn samples provided three XRF spectra (Figure 9, Figures 11 and 13) obtained through the secondary scattering of Rh X-ray radiation on Al₂O₃ targets—the blue curve; Mo targets—the red curve; and HOPG targets—the magenta curve.

It can be observed that the fluorescence lines of Sn and Bi are dominant. Additionally, the XRF spectra also exhibit accompanying lines due to Compton scattering on secondary targets. The elemental concentrations provided by the TurboQuant-Alloys analytical program are presented in Table 6, Table 7 and Table 8, in which the concentrations of the main elements—Bi and Sn can be observed with values very close to those calculated for loading in the furnace. The other elements that appear come from the materials used, which do not have a very high degree of purity. Most of the values represent negligible amounts.

The diffractogram from Figure 10 shows the diffraction peaks of the Bi metal, Sn_0.3Bi_0.7, and Sn_0.95Bi_0.05 compounds in the obtained Bi25Sn75 alloy. This demonstrates a strong interaction between metals, leading to a precipitation of intermetallic compounds in the β Sn matrix (according to EDP-XRF spectra from Figure 9), creating a rather solid solution of Bi in Sn.

The Sn_0.95Bi_0.05 compound is mentioned in a paper of [38] where Bi-added PbSn solders were used as joint samples fabricated from NbTi and Bi2223 wires using an in-situ sheath-dissolution method. In the solder alloy, a major phase (Pb_0.63Sn_0.1Bi_0.27) and a minor phase (Sn_0.95Bi_0.05) were shown. According to the authors, this phase is identified as small islands that do not show superconductivity at liquid He temperature.

The semi-quantitative phase analysis reveals that the sample consists of two phases rich in Sn, accounting for approximately 90% of Sn and 9% of segregated Bi. The 25% of Bi is found in two states—diffused and segregated—within Sn, indicating the potential formation of a Bi-Sn eutectic.

The EDP-XRFS test conducted on the Bi50Sn50 alloy sample provided fluorescence spectra with more intense lines in the low-energy region of the spectra due to the increased proportion of Bi in the alloy (Figure 11).

The diffractometric investigations carried out on the Bi50Sn50 alloy sample reveal an increase in the intensity of the diffraction lines associated with the metallic Bi phase (Figure 12) and a decrease in the intensity of the lines of the Sn_0.95Bi_0.05 compound. In the diffractogram, it is evident that the Sn_0.3Bi_0.7 compound no longer appears, but the orthorhombic SnO₂ phase (cassiterite) is observed.

The diffractogram in Figure 12 does not show lines from phases other than those mentioned earlier, indicating that the alloy does not contain crystallized impurities.

The XRF spectra for the Bi75-Sn25 alloy sample (Figure 13), obtained in the same manner as the previous spectra, exhibit intense characteristic spectral lines, which facilitated a highly accurate elemental analysis, as shown in Table 8.

The diffractometric analysis confirms the presence of metallic Bi in the Bi75Sn25 alloy sample (Figure 14), the formation of the Sn_0.95Bi_0.05 compound, along with traces of metallic Sn. The orthorhombic SnO₂ (cassiterite) phase no longer appears; it was replaced by the tetragonal SnO phase. This phase transformation occurs at temperatures above 573 K. These tin oxides are part of a class of materials with high electrical conductivity, being of interest in many fields, but the most important ones are electronics or optoelectronics.

The elemental composition of the Bi75-Sn25 alloy sample closely aligns with the projected composition of the standard. Residual elements identified in the sample exhibit values near the detection limit of the spectrometer. The origin of characteristic spectral lines attributed to these residual elements may stem from impurities present in the precursors or artifacts like the Compton scattering of specific Sn and Bi fluorescence photons.

4. Conclusions

According to the first objective of the paper, we successfully obtained temperature-dependent thermodynamic properties—enthalpy, entropy, and Gibbs free energy—for Sn and Bi at temperatures of 600 K and 903 K. The obtained experimental values revealed consistent trends, demonstrating agreement between the obtained and expected values. This confirms the accuracy of the experimental procedure. The derived thermodynamic functions enable precise calculations at any temperature–concentration point of interest with high accuracy.

According to the second objective of the work, three types of alloys, namely Bi25-Sn75, Bi50-Sn50, and Bi75-Sn25 were cast and investigated using optical microscopy, SEM-EDS, XRD, and EDP-XRFS methods.

• -. Microstructural variations were observed depending on the bismuth content, with structures exhibiting Sn-based matrices, Bi-based matrices, and equilibrated structures with numerous intermetallic constituents. Bismuth has an obvious effect on the formation of the alloy structure. The absence of conglomerates or clusters causing heterogeneity in chemical composition was noted. Metallic phases were uniformly distributed, and the formation of polyhedral, tetragonal, or orthorhombic compounds was observed.

• -. The SEM-EDS characterization revealed significant changes in the structure of the elaborated alloys. It starts with a hypoeutectic structure (Bi25Sn75), and as the bismuth content is increased, a hypereutectic structure (Bi75Sn25) is reached. The formation of fine needle-shaped beta phases embedded in the alpha phase dendritic matrix is highlighted for the Bi75Sn25 alloy.

• -. The EDP-XRFS analyses revealed that the concentrations of the main elements, Bi and Sn, can be observed with values that are very close to those calculated, as well as the presence of negligible amounts of other elements. Sample one exhibited dominant fluorescence lines of bismuth and tin, with semi-quantitative phase analysis indicating the presence of two Sn-rich phases and segregated bismuth. Sample two showed more intense fluorescence lines in the low-energy region due to increased bismuth content, with diffractometric investigations confirming the intended alloy state. Sample three closely matched the projected standard composition, with identified residual elements attributed to impurities or artifacts.

• -. Diffractometric analysis revealed the prevalence of metallic bismuth and the formation of the Sn_0.3Bi_0.7 and Sn_0.95Bi_0.05 compounds, along with traces of metallic tin, thus demonstrating strong interactions between Bi and Sn leading to the formation of different solid solutions. The presence of SnO is undoubted, as its diffraction lines overlap with the lines of other indexed phases (Figure 14). Also, the lines of SnO₂ overlap except for a line of very small intensity (Figure 13).

Further research will be carried out to establish the optimal Bi-Sn alloy for soldering purposes, including mechanical tests, melting temperature variation, and wetting properties depending on circuit board characteristics.

Author Contributions

Conceptualization, F.N., E.-L.N. and G.I.; methodology, F.N., I.P., G.I. and M.G.; software, F.N., I.P., M.-M.S., R.-M.Ş., N.Ş., G.I. and M.G.; validation, F.N. and G.I.; formal analysis, M.-M.S., M.-I.P., M.B., C.-D.S., A.-A.I., N.Ş., F.N. and M.G.; investigation, I.P., M.B., C.-D.S., A.-A.I., N.Ş., E.-L.N. and M.G.; resources, M.-I.P., M.B., G.I. and F.N.; data curation, F.N., I.P., R.-M.Ş. and G.I.; writing—original draft preparation, F.N., R.-M.Ş. and I.P.; writing—review and editing, F.N., G.I., I.P. and M.G.; visualization, F.N. and E.-L.N.; supervision, G.I. and F.N.; project administration, G.I. and F.N.; funding acquisition, A.-A.I. and M.-I.P. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Footnotes

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.

Figures and Tables

Figure 1. The drawing of the installation used in the experiments [7].

Figure 2. The optical microstructures for the following alloy samples (magnification 50×): (a) Bi25Sn75; (b) Bi50Sn50; (c) Bi75Sn25.

Figure 3. SEM images on the microstructure of Bi25Sn75 alloy sample.

Figure 4. EDS analysis in different spots of Bi25Sn75 alloy sample.

Figure 5. SEM images on the microstructure of Bi50Sn50 alloy sample.

Figure 6. EDS analysis in different spots of Bi50Sn50 alloy sample.

Figure 7. SEM images on the microstructure of Bi75Sn25 alloy sample.

Figure 8. EDS analysis in different spots of Bi75Sn25 alloy sample.

Figure 9. The EDP-XRF spectra correspond to the Bi25Sn75 alloy sample.

Figure 10. The diffractogram of the Bi25Sn75 alloy.

Figure 11. The EDP-XRF spectra correspond to the Bi50Sn50 alloy sample.

Figure 12. The diffractogram of the Bi50Sn50 alloy.

Figure 13. The EDP-XRF spectra correspond to the Bi75Sn25 alloy sample.

Figure 14. The diffractogram of the Bi75Sn25 alloy.

Table 1

Partial molar enthalpies of Bi and Sn la T = 600 K and T = 903 K.

	Temperature 600 K		Temperature 903 K
X Bi	$∆{\overline{\mathit{H}}}_{\mathit{B}\mathit{i}}$[J/mol]	$∆{\overline{\mathit{H}}}_{\mathit{S}\mathit{n}}$[J/mol]	$∆{\overline{\mathit{H}}}_{\mathit{B}\mathit{i}}$[J/mol]	$∆{\overline{\mathit{H}}}_{\mathit{S}\mathit{n}}$[J/mol]
0.0	−6.093 × 103	−5.812 × 103	−9.166 × 103	−4.054 × 103
0.1	−3.389 × 103	−2.593 × 103	−5.100 × 103	0.302 × 103
0.2	−1.188 × 103	−1.188 × 103	−1.788 × 103	3.674 × 103
0.3	0.498 × 103	0.117 × 103	0.751 × 103	6.533 × 103
0.4	1.978 × 103	2.270 × 103	2.977 × 103	9.077 × 103
0.5	3.469 × 103	4.365 × 103	5.217 × 103	1.131 × 103
0.6	5.265 × 103	6.244 × 103	7.911 × 103	1.341 × 103
0.7	8.023 × 103	8.045 × 103	1.207 × 104	1.528 × 103
0.8	1.398 × 104	9.724 × 103	2.104 × 104	1.682 × 103
0.9	1.634 × 104	1.118 × 104	2.459 × 104	1.871 × 104
1.0	0	0	0	0

Table 2

The partial molar entropies of Bi and Sn at 903 K.

	Temperature 903 K
X Bi	$∆{\overline{\mathit{S}}}_{\mathit{B}\mathit{i}}$ [J/mol]	$∆{\overline{\mathit{S}}}_{\mathit{S}\mathit{n}}$ [J/mol]
0.0	27.227	20.721
0.1	23.304	18.631
0.2	14.194	16.964
0.3	10.540	15.047
0.4	8.671	13.000
0.5	7.354	10.975
0.6	6.335	8.831
0.7	5.154	6.696
0.8	3.498	4.524
0.9	1.743	2.245
1.0	0	0

Table 6

XRF-EDP Test Report for Bi25Sn75 alloy sample.

Sample Name 1—Bi25Sn75Description Method—TurboQuant-Alloys
Z	Symbol	Element	Norm. Int.	Concentration	Abs. Error
13	Al	Aluminum	61.6353	0.0775%	0.0014%
15	P	Phosphorus	285.3130	0.08956%	0.00050%
20	Ca	Calcium	3071.7985	0.08837%	0.0010%
28	Ni	Nickel	100.9788	0.0527%	0.0032%
29	Cu	Copper	154.5501	0.0829%	0.0032%
50	Sn	Tin	30,350.5791	75.0850%	0.06%
72	Hf	Hafnium	105.3269	0.0157%	0.0084%
73	Ta	Tantalum	35.1758	0.0153%	0.0024%
74	W	Tungsten	63.0604	0.0623%	0.0031%
79	Au	Gold	534.9385	0.02587%	0.0054%
83	Bi	Bismuth	41,692.0316	24.317%	0.003%
		Other		0.0178%
		Sum		99.93%

Table 7

XRF-EDP Test Report for Bi50Sn50 alloy sample.

Sample Name 2—Bi50Sn50Description Method—TurboQuant-Alloys
Z	Symbol	Element	Norm. Int.	Concentration	Abs. Error
13	Al	Aluminum	88.0285	0.01747%	0.0032%
15	P	Phosphorus	583.8644	0.02200%	0.0010%
20	Ca	Calcium	1977.1571	0.08571%	0.0013%
28	Ni	Nickel	84.0540	0.0540%	0.0074%
29	Cu	Copper	185.2806	0.0354%	0.0076%
30	Zn	Zinc	24.0341	0.01265%	0.00094%
40	Zr	Zirconium	358.2148	0.05063%	0.0018%
47	Ag	Silver	50.5271	0.0137%	0.031%
50	Sn	Tin	18,162.4573	50.062%	0.05%
51	Sb	Antimony	26.0649	0.1180%	0.0027%
72	Hf	Hafnium	89.0883	0.0183%	0.021%
74	W	Tungsten	76.8009	0.1045%	0.0080%
78	Pt	Platinum	2567.5260	0.034%	0.073%
81	Tl	Thallium	845.4254	0.0402%	0.0062%
82	Pb	Lead	506.6176	0.0273%	0.0062%
83	Bi	Bismuth	97,237.1198	49.035%	0.01%
		Other		0.1768%
		Sum		99.91%

Table 8

XRF-EDP Test Report for Bi75Sn25 alloy sample.

Sample Name 3—Bi75Sn25Description Method—TurboQuant-Alloys
Z	Symbol	Element	Norm. Int.	Concentration	Abs. Error
13	Al	Aluminum	145.4822	0.0531%	0.0094%
20	Ca	Calcium	762.9267	0.06833%	0.0019%
28	Ni	Nickel	98.3900	0.0110%	0.019%
29	Cu	Copper	101.9181	0.0122%	0.018%
30	Zn	Zinc	39.4357	0.0354%	0.0025%
48	Cd	Cadmium	4.1658	0.0472%	0.0062%
50	Sn	Tin	5420.5119	24.9102%	0.05%
51	Sb	Antimony	13.5239	0.0119%	0.0035%
52	Te	Tellurium	8.0419	0.01203%	0.00058%
83	Bi	Bismuth	163,456.8712	74.7006%	0.02%
		Other		0.0918%
		Sum		99.95%

References

1. European Commission. Restriction of Hazardous Substances in Electrical and Electronic Equipment (RoHS), Energy, Climate Change, Environment. Available online: https://environment.ec.europa.eu/topics/waste-and-recycling/rohs-directive_en (accessed on 5 January 2024).

2. European Commission. Waste from Electrical and Electronic Equipment (WEEE), Energy, Climate Change, Environment. Available online: https://environment.ec.europa.eu/topics/waste-and-recycling/waste-electrical-and-electronic-equipment-weee_en (accessed on 5 January 2024).

3. Sujong, K.; Wonsik, H.; Hyunbin, N.; Namhyun, K. Growth Behavior of Intermetallic Compounds in Various Solder Joints Induced by Electromigration. J. Weld. Join.; 2021; 39, pp. 89-102.

4. Ko, Y.K.; Ko, Y.H.; Lee, C.W. Micro-bump Joining Technology for 3-Dimensional Chip Stacking. J. Korean Soc. Precis. Eng.; 2014; 31, pp. 865-871. [DOI: https://dx.doi.org/10.7736/KSPE.2014.31.10.865]

5. Arkhipov, P.A.; Grishenkova, O.V.; Kholkina, A.S. Thermodynamic characteristics of liquid metallic alloys containing lead. antimony and bismuth. J. Mol. Liq.; 2021; 335, 116071. [DOI: https://dx.doi.org/10.1016/j.molliq.2021.116071]

6. Park, J.M.; Kim, S.H.; Jeong, M.H.; Hand, M.; Park, Y.B. Effect of Cu-Sn intermetallic compound reactions on the Kirkendall void growth characteristics in Cu/Sn/Cu microbumps. Jpn. J. Appl. Phys.; 2014; 53, 05HA06. [DOI: https://dx.doi.org/10.7567/JJAP.53.05HA06]

7. Niculescu, E.; Iacob, G.; Niculescu, F.; Pencea, I.; Buzatu, M.; Petrescu, M.I.; Marcu, D.M.; Turcu, R.N.; Geanta, V.; Butu, M. Experimental Determination of the Activities of Liquid Bi-Sn Alloys. J. Phase Equilib. Diff.; 2021; 42, pp. 278-289. [DOI: https://dx.doi.org/10.1007/s11669-021-00880-2]

8. Ren, G.; Wilding, I.J.; Collins, M.N. Alloying Influences on Low Melt Temperature SnZn and SnBi Solder Alloys for Electronic Interconnections. J. Alloys Compd.; 2016; 665, pp. 251-260. [DOI: https://dx.doi.org/10.1016/j.jallcom.2016.01.006]

9. Miao, H.W.; Duh, J.G. Microstructure evolution in Sn–Bi and Sn–Bi–Cu solder joints under thermal aging. Mater. Chem. Phys.; 2001; 71, pp. 255-271. [DOI: https://dx.doi.org/10.1016/S0254-0584(01)00298-X]

10. Suksongkarm, P.; Rojananan, S.; Rojananan, S. Bismuth Formation in Lead-Free Cu–Zn–Si Yellow Brass with Various Bismuth–Tin Alloy Additions. Mater. Trans.; 2017; 58, pp. 1754-1760. [DOI: https://dx.doi.org/10.2320/matertrans.M2017227]

11. Gao, Y.; Bian, X.; Qiu, X.; Jia, Y.; Yi, J.; Wang, G. Investigation of Microstructure and Mechanical Properties of SAC105 Solders with Sb, In, Ni, and Bi Additions. Materials; 2023; 16, 4059. [DOI: https://dx.doi.org/10.3390/ma16114059] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37297193]

12. Chen, S.; Huang, C.M.; Pecht, M. A review of lead-free solders for electronics applications. Microelectron. Reliab.; 2017; 75, pp. 77-95. [DOI: https://dx.doi.org/10.1016/j.microrel.2017.06.016]

13. Huang, M.; Zhao, N.; Liu, S.; He, Y. Drop failure modes of Sn–3.0Ag–0.5Cu solder joints in wafer level chip scale package. Trans. Nonferrous Met. Soc. China; 2016; 26, pp. 1663-1669. [DOI: https://dx.doi.org/10.1016/S1003-6326(16)64272-3]

14. Chong, D.Y.R.; Che, F.X.; Pang, J.H.L.; Ng, K.; Tan, J.Y.N.; Low, P.T.H. Drop impact reliability testing for lead-free and lead-based soldered IC packages. Microelectron. Reliab.; 2006; 46, pp. 1160-1171. [DOI: https://dx.doi.org/10.1016/j.microrel.2005.10.011]

15. Suraski, D.; Seelig, K. The current status of lead-free solder alloys. IEEE Trans. Adv. Packag.; 2001; 24, pp. 244-248. [DOI: https://dx.doi.org/10.1109/6104.980031]

16. Ng, W.; Weng, C. Evolution of Pb-Free Solders. Recent Progress in Soldering Material; Mohamad, A.A. IntechOpen Limited: London, UK, 2017.

17. Braga, M.H.; Vizdal, J.; Kroupa, A.; Ferreira, J.A. The experimental study of the Bi–Sn, Bi–Zn and Bi–Sn–Zn systems. Calphad; 2007; 31, pp. 468-478. [DOI: https://dx.doi.org/10.1016/j.calphad.2007.04.004]

18. Manasijević, D.; Vřešťál, J.; Minić, D.; Kroupa, A.; Źivković, D.; Źivković, Ź. Phase Equilibria and Thermodynamics of the Bi-Sb-Sn Ternary System. J. Alloys Compd.; 2007; 438, pp. 150-157. [DOI: https://dx.doi.org/10.1016/j.jallcom.2006.08.021]

19. Vizdal, J.; Braga, M.H.; Kroupa, A.; Richter, K.W.; Soares, D.; Malheiros, L.F.; Ferreira, J. Thermodynamic Assessment of the Bi-Sn-Zn System. Calphad; 2007; 31, pp. 438-448. [DOI: https://dx.doi.org/10.1016/j.calphad.2007.05.002]

20. El-Ashram, T. Formation of Supersaturated Solid Solutions Bi-Ag and Bi-Zn by Rapid Solidification Using Melt Spinning Technique. Mater. Sci. Appl.; 2015; 6, pp. 183-188. [DOI: https://dx.doi.org/10.4236/msa.2015.62021]

21. Gaskell, D.R. Introduction to the Thermodynamics of Materials; CRC Press: Boca Raton, CA, USA, 2008.

22. Niculescu, E.L.; Buzatu, M.; Petrescu, M.I.; Butu, M.; Iacob, G.; Niculescu, F.; Florea, B.; Marcu, D.F.; Stancel, C.D.; Stanescu, M.M. Thermodynamic calculation of the binary systems Bi-Sn by implementing a JAVA interface. UPB Sci. Bull. B Chem. Mater. Sci.; 2021; 83, pp. 219-228.

23. Asryan, N.; Mikula, A. Thermodynamic properties of Bi–Sn melts. Inorg. Mater.; 2004; 40, pp. 457-461. [DOI: https://dx.doi.org/10.1023/B:INMA.0000023961.34926.67]

24. Yang, C.; Chen, F.; Gierlotka, W.; Chen, S.W.; Hsieh, K.C.; Huang, L.L. Thermodynamic properties and phase equilibria of Sn–Bi–Zn ternary alloys. Mater. Chem. Phys.; 2008; 112, pp. 94-103. [DOI: https://dx.doi.org/10.1016/j.matchemphys.2008.05.034]

25. Kirshon, Y.; Ben Shalom, S.; Emuna, M.; Greenberg, Y.; Lee, J.; Makov, G.; Yahel, E. Thermophysical Measurements in Liquid Alloys and Phase Diagram Studies. Materials; 2019; 12, 3999. [DOI: https://dx.doi.org/10.3390/ma12233999]

26. Stlen, S.; Grande, T. Chemical Thermodynamics of Materials: Macroscopic and Microscopic Aspects; John Wiley & Sons: Hoboken, NJ, USA, 2004.

27. Lukas, H.L.; Fries, S.G.; Sundman, B. Computational Thermodynamics; Cambridge University Press: Cambridge, UK, 2007.

28. Niculescu, F.; Constantin, I.; Buzatu, M.; Marcu, D.F.; Csaki, I. The thermodynamic proprieties of the Pb-Sb system. UPB Sci. Bull. B Chem. Mater. Sci.; 2016; 78, pp. 167-174.

29. Wojciech, G.; Cho-yu, L.; Prasert, C.; Arifur, R.; Ting-Nan, K. Thermodynamic Re-optimization of the Binary Pb-Sb System. J. Phase Equilib. Diff.; 2013; 34, pp. 421-424.

30. Elliott, J.F.; Chipman, J. The Thermodynamic Properties of Liquid Ternary Cadmium Solutions. J. Am. Chem. Soc.; 1951; 73, pp. 2682-2693. [DOI: https://dx.doi.org/10.1021/ja01150a075]

31. Magnus, A.; Mannheimer, M. The Measurements of the Enthalpy of Mixing of Liquid Metals. Z. Physik. Chem.; 1926; 121, pp. 267-278. [DOI: https://dx.doi.org/10.1515/zpch-1926-12121]

32. Asryan, N.; Mikula, A. Thermodynamic Properties of Liquid Bi-Sn Alloys. Int. J. Mater. Res.; 2004; 95, pp. 132-135. [DOI: https://dx.doi.org/10.3139/146.017925]

33. Kumar, M.R.; Mohan, S.; Behera, C.K. Thermodynamic Accessment Experimentally on Bi-Sn System by Calorimeter. Mater. Today Proc.; 2018; 5, pp. 27777-27785. [DOI: https://dx.doi.org/10.1016/j.matpr.2018.10.013]

34. Yazawa, A.; Koike, K. Activity Measurement in Sn-Pb, Sn-Bi, Sn-Tl and Sn-Pb-Bi Alloys at 110 °C. J. Min. Met. Inst. Japan; 1968; 184, pp. 1593-1596.

35. Wang, F.; Huang, Y.; Zhang, Z.; Yan, C. Interfacial Reaction and Mechanical Properties of Sn-Bi Solder joints. Materials; 2017; 10, 920. [DOI: https://dx.doi.org/10.3390/ma10080920]

36. Hultgren, R. Selected Values of the Thermodynamic Properties of the Elements; ASM: Metals Park, OH, USA, 1973.

37. Hultgren, R. Selected Values of the Thermodynamic Properties of Binary Alloys; ASM: Metals Park, OH, USA, 1973.

38. Ryo, M.; Hirotsugu, I.; Aichi, Y.; Hiroshi, H.; Gen, N.; Masashi, T.; Hiroyuki, T.; Yoshihiko, T. Superconducting joints using Bi-added PbSn solders. Appl. Phys. Express; 2017; 10, pp. 93-102.