Experimental Measurements and Thermodynamic Optimization of the NaCl+RbCl Phase Diagram

In this study, the phase diagram of the NaCl+RbCl binary system was measured using differential scanning calorimetry, and the measured molar percentages of the binary mixture RbCl ranged from 12 to 97 mol%, updating and extending the experimental data for this binary system. The liquid phase of this system was thermodynamically modeled using a modified quasichemical model, and a computer optimization program for evaluation of the experimental data on phase equilibria and other thermodynamic data was developed based on this model. All thermodynamic models of the NaCl+RbCl binary system were constructed; calculations of the phase diagram, enthalpy of mixing, and activity coefficient of NaCl of the NaCl+RbCl binary system were completed using the obtained model parameters; and the calculated eutectic point positions were xRbCl = 0.567 and 550.2 °C. The calculated results were able to reproduce the various types of experimental data well, and the thermodynamic self-consistency between different types of experimental data was demonstrated.


Introduction
In recent years, four high-temperature molten salt materials, represented by nitrate, carbonate, chloride, and fluoride salts, have gained increased value in industrial applications with their excellent physicochemical properties. Compared with traditional heat transfer media (water and water vapor), molten salt materials have better characteristics, such as low vapor pressure, high boiling point, good thermal stability, relatively good thermal and electrical conductivities, higher specific heat, and high sensible and latent heat storage capacity. These four molten salts have not only led to the development of new technologies in energy fields such as nuclear power and solar energy but also play a key role in molten salt batteries, solid-state lasers, and material preparation.
The chemical stability of chloride salts is lower than that of fluoride salts but higher than that of nitrates and carbonates. Therefore, chlorinated salts are often used as heat transfer and heat storage media in advanced nuclear reactors, collector-type solar power, electrolytic production of rare earth metals, molten salt batteries, and other fields. Alkali metal chloride salts, such as LiCl+NaCl+MgCl 2 and LiCl+KCl+MgCl 2 , can be used as primary coolants in molten salt fast reactor (MSFR) [1]. Chlorinated mixed salts (e.g., NaCl+KCl+ZnCl 2 ) [2] are a high-temperature heat transfer fluid that can potentially be used in combination with advanced Brayton cycles for collector-type solar power systems. In the field of electrolytic reduction of rare earth metals, molten chloride LiCl+NaCl+KCl+YCl 3 [3] is used as an electrolyte to prepare pure yttrium, and molten chloride LiCl+KCl+MgCl 2 +YCl 3 is used as an electrolyte [4] to synthesize Mg+Li+Y alloys through electrolysis. The thermally activated cells currently being developed in the field of high-energy-density batteries are FeS 2 |KCl-LiCl|LiSi [5], CoS 2 |LiCl-LiBr-LiF|LiSi [6], Li|LiCl-KCl-RbCl-LiF|FeS 2 [7], etc. Alkali metal halides containing the cation Li + are suitable electrolyte materials because of their lower melting point (compared with alkali metal halides containing the cation Na + ), lower vapor pressure, and high ionic conductivity.
Ternary chloride salts containing the binary system NaCl+RbCl and Li+ ions are a potential molten salt cell. The addition of rare earth metal elements to the binary system NaCl+RbCl will be used for the electrolytic reduction of rare earth metals. The experimental and theoretical study of the phase equilibrium properties of the NaCl+RbCl binary system described in this paper will help engineers to understand the thermodynamic properties of chlorides containing Na + and Rb + cations, and will help with the design of new molten salt materials in the future.

Differential Scanning Calorimetry
Differential scanning calorimetry (DSC), as a thermal analysis technique, achieves the equality of sample and reference temperatures at any moment and follows the procedure to raise and lower the temperature. Because the sample and reference require different heating amounts, there is a heating amount difference. After a DSC measurement, a DSC curve is plotted. It represents heat flow as a function of temperature or time. There are two types of DSC, the power compensated DSC and the heat flow DSC, and they use different measurement methods.
Power-compensated DSC adopts the "zero equilibrium" principle, which is obtained when the temperature difference ∆T between sample and reference is equal to 0. This principle should be permanently maintained during measurements. Power-compensated DSC consists of two separate heating control systems. Through precise control of the heating power, one system ensures that the reference and sample are heated or cooled at a uniform heating or cooling rate, and the other system ensures compliance with the zero equilibrium principle must, permanently providing ∆T = 0. The sample and reference furnaces are individually heated or cooled according to the same temperature-time program. During the cooling process, the compensation heating wire rapidly heats the reference, and during the heating process, the compensation heating wire starts to heat the specimen, always keeping the temperature difference ∆T equal to 0.During the experimental measurement, the difference between the power applied to the investigated sample and to the reference is measured as a function of temperature or time. The working principle of heat-flow DSC is similar to that of the power-compensated DSC. The difference between them is that the sample and reference are simultaneously heated on an electrically heating platform by a single heater. The heating platform is made of a material with good thermal conductivity, which allows for rapid heat transfer between the sample and reference. This high thermal conductivity material ensures that the temperature difference between the reference and sample is close to 0 when the thermal effect occurs. Figure 1 shows a schematic diagram of the heat-flow-type and power-compensationtype DSC. Power-compensated DSC has separate heaters for the sample and reference. They are separately heated at the command of the power compensation system. In a heatflow DSC, the sample and reference are in the same heating chamber and are placed on a heating platform with good thermal conductivity. The heater under the heating platform simultaneously heats the sample and the reference. Figure 2 depicts a schematic diagram of a typical heat-flow DSC. The sample and reference crucibles are placed on a copper conduction heating platform; two nickel-chromium alloy disks are placed underneath the platform, directly beneath the sample and reference. The thermocouple connection points for measuring temperature and temperature difference are connected under the nickel-chromium alloy disk. OR PEER REVIEW 3 of 12 mocouple connection points for measuring temperature and temperature difference are connected under the nickel-chromium alloy disk.

DSC Curve Analysis
Currently, differential scanning calorimetry (DSC) is the most used method, and a large amount of experimental data on the phase diagrams of binary and multiple chloride systems have been accumulated through experiments. However, the main disagreement lies in how to analyze and determine the characteristic points in the phase diagram based on the obtained DSC curves. This has become an important part of experimental accuracy. The analysis of differential thermal analysis (DTA) curves of a large number of alloys by Boettinger et al. [8], which we also adopted in this study, is now commonly used internationally.
Essentially, power-compensated DSC is the real DSC. It can comply with the zero equilibrium principle and guarantee that the temperature difference between the reference and sample is equal to 0 at all times. The basic configurations of DTA and heat-flow-type DSC are the same. Their heating sources for controlling the temperatures of the reference and sample are different. The heating sources of DTA and heat-flow-type DSC are the furnace wall and the separate heater respectively. Both the DTA and the heat flow type DSC can quantitatively measure the heat of phase transition. Their differences mainly exist in sensitivity and accuracy. Therefore, the analytical methods of DTA and DSC curves are basically the same. mocouple connection points for measuring temperature and temperature differ connected under the nickel-chromium alloy disk.

DSC Curve Analysis
Currently, differential scanning calorimetry (DSC) is the most used metho large amount of experimental data on the phase diagrams of binary and multip ride systems have been accumulated through experiments. However, the ma greement lies in how to analyze and determine the characteristic points in the p agram based on the obtained DSC curves. This has become an important part o mental accuracy. The analysis of differential thermal analysis (DTA) curves o number of alloys by Boettinger et al. [8], which we also adopted in this study commonly used internationally.
Essentially, power-compensated DSC is the real DSC. It can comply with equilibrium principle and guarantee that the temperature difference between t ence and sample is equal to 0 at all times. The basic configurations of D heat-flow-type DSC are the same. Their heating sources for controlling the temp of the reference and sample are different. The heating sources of DTA and heat-fl DSC are the furnace wall and the separate heater respectively. Both the DTA and flow type DSC can quantitatively measure the heat of phase transition. Their dif mainly exist in sensitivity and accuracy. Therefore, the analytical methods of D DSC curves are basically the same.

DSC Curve Analysis
Currently, differential scanning calorimetry (DSC) is the most used method, and a large amount of experimental data on the phase diagrams of binary and multiple chloride systems have been accumulated through experiments. However, the main disagreement lies in how to analyze and determine the characteristic points in the phase diagram based on the obtained DSC curves. This has become an important part of experimental accuracy. The analysis of differential thermal analysis (DTA) curves of a large number of alloys by Boettinger et al. [8], which we also adopted in this study, is now commonly used internationally.
Essentially, power-compensated DSC is the real DSC. It can comply with the zero equilibrium principle and guarantee that the temperature difference between the reference and sample is equal to 0 at all times. The basic configurations of DTA and heat-flow-type DSC are the same. Their heating sources for controlling the temperatures of the reference and sample are different. The heating sources of DTA and heat-flow-type DSC are the furnace wall and the separate heater respectively. Both the DTA and the heat flow type DSC can quantitatively measure the heat of phase transition. Their differences mainly exist in sensitivity and accuracy. Therefore, the analytical methods of DTA and DSC curves are basically the same.
The ideal case of the DTA curve for the melting and solidification process of a pure substance is shown in Figure 3. The melting point of a pure substance can be determined by the extrapolation onset point of the heating peak and by the extrapolation termination point of the cooling peak. For binary molten salts or binary alloy materials, take the phase diagram of the A-B binary system in Figure 4 as an example. The phase diagram of the A-B binary system is a simple phase diagram with one eutectic point. The eutectic temperatures and liquidus temperatures can be determined by the onset point of the first endothermic peak and the peak point of the second endothermic peak in the endothermic process. They can also be determined using the termination points of the first and second exothermic peaks of the DTA curve during the exothermic process. This analytical method was proposed by Boettinger et al. [8] after analyzing the DTA curves of a large number of alloys. In this study, this method was used to analyze the measured DSC curves.

Materials 2022, 15, x FOR PEER REVIEW
The ideal case of the DTA curve for the melting and solidification process o substance is shown in Figure 3. The melting point of a pure substance can be dete by the extrapolation onset point of the heating peak and by the extrapolation term point of the cooling peak. For binary molten salts or binary alloy materials, take th diagram of the A-B binary system in Figure 4as an example. The phase diagram A-B binary system is a simple phase diagram with one eutectic point. The eutec peratures and liquidus temperatures can be determined by the onset point of endothermic peak and the peak point of the second endothermic peak in the endo process. They can also be determined using the termination points of the first and exothermic peaks of the DTA curve during the exothermic process. This an method was proposed by Boettinger et al. [8] after analyzing the DTA curves o number of alloys. In this study, this method was used to analyze the measur curves.

Experimental Instruments and Materials
The thermal analysis measurements completed in this study were performed synchronous thermal analyzer (TGA-DSC) (NETZSCH, instrument model: STA Jupiter, Selb, Germany). It allows the measurement of thermal effects at temp ranging from room temperature to 1450°C. Pure NaCl and RbCl, with a purity o and 99.95%, respectively, were purchased from Aladdin. The pure RbCl was hygroscopic. In order to dry the pure RbCl, it was dried in a vacuum-drying o temperature of 200 °C for 2 h. After completing drying, the pure RbCl was stored glove box. We prepared about 40 mg of the binary mixture and placed it in num-rhodium crucible. Nitrogen was used as a protective gas, and the tempera point of the cooling peak. For binary molten salts or binary alloy mater diagram of the A-B binary system in Figure 4as an example. The ph A-B binary system is a simple phase diagram with one eutectic point peratures and liquidus temperatures can be determined by the onse endothermic peak and the peak point of the second endothermic peak process. They can also be determined using the termination points of t exothermic peaks of the DTA curve during the exothermic proce method was proposed by Boettinger et al. [8] after analyzing the DTA number of alloys. In this study, this method was used to analyze t curves.

Experimental Instruments and Materials
The thermal analysis measurements completed in this study were synchronous thermal analyzer (TGA-DSC) (NETZSCH, instrument m Jupiter, Selb, Germany). It allows the measurement of thermal effec ranging from room temperature to 1450°C. Pure NaCl and RbCl, with and 99.95%, respectively, were purchased from Aladdin. The pure hygroscopic. In order to dry the pure RbCl, it was dried in a vacuum temperature of 200 °C for 2 h. After completing drying, the pure RbCl w

Experimental Instruments and Materials
The thermal analysis measurements completed in this study were performed using a synchronous thermal analyzer (TGA-DSC) (NETZSCH, instrument model: STA 449 F3 Jupiter, Selb, Germany). It allows the measurement of thermal effects at temperatures ranging from room temperature to 1450 • C. Pure NaCl and RbCl, with a purity of 99.99% and 99.95%, respectively, were purchased from Aladdin. The pure RbCl was slightly hygroscopic. In order to dry the pure RbCl, it was dried in a vacuum-drying oven at a temperature of 200 • C for 2 h. After completing drying, the pure RbCl was stored in a dry glove box. We prepared about 40 mg of the binary mixture and placed it in a platinumrhodium crucible. Nitrogen was used as a protective gas, and the temperature rise and fall rate were set to 10 K·min −1 . Two heating and cooling cycles were completed. The first cycle was performed to obtain a homogenized binary salt mixture. The DSC curves obtained from the first cycle were not used for analysis. The DSC curves obtained from the second heating and cooling cycle were mainly analyzed.

Experimental Results and Discussion
3.1. Experimental Results for Pure Substances Figure 5 presents the DSC curves of pure NaCl and RbCl. As shown in Figure 5a, for pure NaCl, the melting temperatures obtained from the curve analysis are 803.4 • C (heating curve) and 800.1 • C (cooling curve). As shown in Figure 5b, for pure RbCl, the obtained melting temperatures are 721.7 • C (heating curve) and 719.2 • C (cooling curve). These values are listed in Table 1. It can be seen that the effect of supercooling on the melting point is negligible. rials 2022, 15, x FOR PEER REVIEW obtained from the first cycle were not used for analysis. The DSC curves obtained the second heating and cooling cycle were mainly analyzed.  Table 1. It can be seen that the effect of supercooling melting point is negligible.      Figure 6 illustrates the experimentally obtained DSC curves for the NaCl+RbCl mixtures with different mole fractions x RbCl . Table 1 presents the analytically determined eutectic and liquidus temperatures. As shown in Figure 6, basically all of the cooling DSC curves contain two exothermic peaks, and some of the heating curves have only one heat endothermic peak. The eutectic temperature of the NaCl+RbCl binary system can be clearly reflected by the endothermic and exothermic peaks. The liquidus line can also be determined by the exothermic peaks. However, some endothermic peaks cannot reflect the liquidus temperature. The reason may be that the compositions of some NaCl+RbCl mixtures are very close to the eutectic composition, and the endothermic peaks of the eutectic reaction are too strong to cover up the endothermic peak for reflecting the liquidus lines in the heating process, such as the heating DSC curves with an x RbCl equal to 52.5 and 57.0 mol%. It is also possible that anthropogenic or nonanthropogenic interference affected the formation of peaks, and a garbled signal appeared. The heating DSC curves for x RbCl equal to 40.7 and 67.6 mol% reflect this situation. , x FOR PEER REVIEW 6 of 12 Figure 6 illustrates the experimentally obtained DSC curves for the NaCl+RbCl mixtures with different mole fractions xRbCl. Table 1 presents the analytically determined eutectic and liquidus temperatures. As shown in Figure 6, basically all of the cooling DSC curves contain two exothermic peaks, and some of the heating curves have only one heat endothermic peak. The eutectic temperature of the NaCl+RbCl binary system can be clearly reflected by the endothermic and exothermic peaks. The liquidus line can also be determined by the exothermic peaks. However, some endothermic peaks cannot reflect the liquidus temperature. The reason may be that the compositions of some NaCl+RbCl mixtures are very close to the eutectic composition, and the endothermic peaks of the eutectic reaction are too strong to cover up the endothermic peak for reflecting the liquidus lines in the heating process, such as the heating DSC curves with an xRbCl equal to 52.5 and 57.0 mol%. It is also possible that anthropogenic or nonanthropogenic interference affected the formation of peaks, and a garbled signal appeared. The heating DSC curves for xRbCl equal to 40.7 and 67.6 mol% reflect this situation.  As can be seen in Table 1, due to the effect of supercooling, in most ca transition temperature determined by the heating curve is higher than th by the cooling curve. Only when xRbCl= 30.6 mol%, its liquidus temperatu by the heating curve is lower than that determined by the cooling curve. As can be seen in Table 1, due to the effect of supercooling, in most cases, the phase transition temperature determined by the heating curve is higher than that determined by the cooling curve. Only when x RbCl = 30.6 mol%, its liquidus temperature determined by the heating curve is lower than that determined by the cooling curve. Therefore, for this mole fraction, only the liquidus temperature determined by the cooling curve was used in this study. For the two DSC curves that are very close to the pure RbCl, such as x RbCl = 96 and 97 mol%, the endothermic and exothermic peaks on the DSC curve reflecting the eutectic reaction become very small, especially on the cooling DSC curve, where the exothermic peak reflecting the eutectic reaction is almost invisible. This indicates that the eutectic reaction is no longer occurring here, so for the thermodynamic optimization in this study, we used only the liquidus temperatures.

Error Analysis of DSC Measurements
Errors in DSC measurement mainly come from the systematic error of DSC instruments and personal errors during the operating procedure (including sampling and data processing). The systematic error of DSC instruments includes thermal analysis accuracy and DSC curve analysis error. NETZSCH's software can be used to mark the start point, peak, inflection point, and endpoint temperatures of peaks, and to perform automatic peak searches. The sources of personal errors mainly come from sample storage, dehumidification, weighing, and mixing. Table 2 presents all thermodynamic data (∆H 0 298.15 K ,S 0 298.15 K , and C p (T)) for the condensed pure compounds of the NaCl+RbCl systems. The molar Gibbs free energy of a stoichiometric phase can be written as: The modified quasichemical model [10,11] was used to describe the liquid phase of the NaCl+RbCl system. The nonconfigurational Gibbs energy change ∆g AB containing the parameters of this model is an important term in the expression for the Gibbs free energy of the solution. ∆g AB can be expanded as a polynomial as:

Thermodynamic Model
where ∆g 0 AB and g ij AB are the parameters of the model that can be functions of temperature. These parameters can be determined by optimizing the phase equilibrium and 8 of 11 thermodynamic data, which were measured by the experimental method.χ AB and χ BA are component variables, and they can be defined as: where X AA , X AB , and X BB represent the pair fractions of cation-cation pairs. There are two terminal solid solutions in the NaCl+RbCl system: the terminal solid solution α (near the pure NaCl) and the terminal solid solution β (near the pure RbCl). For the solid solution α, the molar fraction x RbCl is very small, and its excess Gibbs energy expression can be written as: where L α is a temperature-independent coefficient. Solid solution β, it can be considered that the Na+ cations enter into the interspaces of the RbCl crystal structure. Using the classical polynomial, its excess Gibbs energy expression can be written as: where a, b, and c are temperature-independent coefficients. The values of L α , a, b, and c all need to be determined by thermodynamic optimization.

Calphad Method
The essence of the Calphad method is the computer coupling of phase diagrams and thermochemistry. A classic book describing the generation and application of the Calphad techniques was authored by Saunders and Miodownik [12]. The experience and skill for the thermodynamic optimization of phase diagrams were provided by Hari Kumar and Wollants [13]. After learning a great deal about the Calphad method from the above references, we developed a computer program to perform all the optimization and calculations in this work. The calculating flow of our computer program was based on the flowchart of the Calphad method in [13]. This computer program has been used to perform thermodynamic evaluation and optimization for some binary and ternary chloride molten salt systems [14][15][16]. In this program, the minimum Gibbs energy of the system was found through the algorithms of golden section search and parabolic interpolation, and the optimal model parameters were obtained by Nelder-Mead simplex algorithm [17].

Results and Discussion
In this study, all available phase equilibrium property data and thermodynamic data for the binary NaCl+RbCl system were carefully evaluated. The phase diagram data of the NaCl+RbCl system measured in this study were also simultaneously evaluated. These experimental data used in the thermodynamic optimization process were the activity of the NaCl of the NaCl+RbCl system [18], enthalpy of mixing [19], and phase diagram [20,21]. The optimization of the phase diagram and thermodynamic properties of the NaCl+RbCl system was performed using a computer program we developed. Firstly, the data of mixing enthalpy, activity, and eutectic point were used to determine the initial values of the model parameters. Secondly, all experimental data were used to optimize the model parameters.
Pelton and Flengas [20] measured the phase diagram of the NaCl+RbCl system using the step-cooling curve method in 1970. The molar percentages of RbCl of the compositions for binary mixtures they reported ranged from 10 to 90 mol%. The one reported eutectic point was x RbCl = 0.56 at 550 • C. In this study, the phase diagram of the NaCl+RbCl system was also measured using the DSC method, and the molar percentages of RbCl ranged from 12 to 97 mol%. More data for the eutectic temperature are also given. Short and Roy [21] measured the solubility limits of the two terminal solid solutions in the NaCl+RbCl system by using the interdiffusion method. They found that, at 515 • C, the solubility limit of RbCl in solid solution α (where NaCl is the main body) is 0.8 mol%, and the solubility limit of NaCl in solid solution β (where RbCl is the main body) is 6 mol%. Figure 7 shows a comparison of the calculated phase diagram of the NaCl+RbCl system with the corresponding experimental data. The calculated values basically reproduced the experimental data. The calculated solubility lines for the two solid solutions below the eutectic temperature also reproduced the experimental data given by Short and Roy. The calculated eutectic points are x RbCl = 0.567 and 550.2 • C, which are also in better agreement with the results recommended by Pelton and Flengas.
Materials 2022, 15, x FOR PEER REVIEW 10 of 12 NaCl+RbCl system was also measured using the DSC method, and the molar percentages of RbCl ranged from 12 to 97 mol%. More data for the eutectic temperature are also given. Short and Roy [21] measured the solubility limits of the two terminal solid solutions in the NaCl+RbCl system by using the interdiffusion method. They found that, at 515 °C, the solubility limit of RbCl in solid solution α (where NaCl is the main body) is 0.8 mol%, and the solubility limit of NaCl in solid solution β (where RbCl is the main body) is 6 mol%. Figure 7 shows a comparison of the calculated phase diagram of the NaCl+RbCl system with the corresponding experimental data. The calculated values basically reproduced the experimental data. The calculated solubility lines for the two solid solutions below the eutectic temperature also reproduced the experimental data given by Short and Roy. The calculated eutectic points are xRbCl= 0.567 and 550.2 °C, which are also in better agreement with the results recommended by Pelton and Flengas. Figure 7. The calculated phase diagram of the NaCl+RbCl system at atmospheric pressure. □,Heating process; ■cooling process phase diagram points measured in this study; (•)Pelton and Flengas [20]; (★) Short and Roy [21]. Figures 8 and 9 show the calculated enthalpy of mixing and activity of the liquid phase of the NaCl+RbCl system, respectively, in which we compared the corresponding experimental data. It can be seen that the calculated enthalpy of mixing of the liquid phase at 810 °C is in good agreement with the experimental data [19]. Only the calculated activity of NaCl somewhat deviates from the experimental data [18], and the average relative error between the calculated and experimental values is 2.86%. Figure 7. The calculated phase diagram of the NaCl+RbCl system at atmospheric pressure. , Heating process; cooling process phase diagram points measured in this study; (•) Pelton and Flengas [20]; ( ) Short and Roy [21]. Figures 8 and 9 show the calculated enthalpy of mixing and activity of the liquid phase of the NaCl+RbCl system, respectively, in which we compared the corresponding experimental data. It can be seen that the calculated enthalpy of mixing of the liquid phase at 810 • C is in good agreement with the experimental data [19]. Only the calculated activity of NaCl somewhat deviates from the experimental data [18], and the average relative error between the calculated and experimental values is 2.86%.
Materials 2022, 15, x FOR PEER REVIEW 10 of 12 NaCl+RbCl system was also measured using the DSC method, and the molar percentages of RbCl ranged from 12 to 97 mol%. More data for the eutectic temperature are also given. Short and Roy [21] measured the solubility limits of the two terminal solid solutions in the NaCl+RbCl system by using the interdiffusion method. They found that, at 515 °C, the solubility limit of RbCl in solid solution α (where NaCl is the main body) is 0.8 mol%, and the solubility limit of NaCl in solid solution β (where RbCl is the main body) is 6 mol%. Figure 7 shows a comparison of the calculated phase diagram of the NaCl+RbCl system with the corresponding experimental data. The calculated values basically reproduced the experimental data. The calculated solubility lines for the two solid solutions below the eutectic temperature also reproduced the experimental data given by Short and Roy. The calculated eutectic points are xRbCl= 0.567 and 550.2 °C, which are also in better agreement with the results recommended by Pelton and Flengas. Figure 7. The calculated phase diagram of the NaCl+RbCl system at atmospheric pressure. □,Heating process; ■cooling process phase diagram points measured in this study; (•)Pelton and Flengas [20]; (★) Short and Roy [21]. Figures 8 and 9 show the calculated enthalpy of mixing and activity of the liquid phase of the NaCl+RbCl system, respectively, in which we compared the corresponding experimental data. It can be seen that the calculated enthalpy of mixing of the liquid phase at 810 °C is in good agreement with the experimental data [19]. Only the calculated activity of NaCl somewhat deviates from the experimental data [18], and the average relative error between the calculated and experimental values is 2.86%.

Conclusions
In this study, differential scanning calorimetry analysis of the various compositions in the NaCl+RbCl system was performed. A thermodynamic data set for the NaCl+RbCl system was generated to reproduce all available experimental data. The following conclusions were obtained: (1) The phase diagram of the NaCl+RbCl system was measured using DSC, and the measured molar percentages of RbCl of the binary mixture ranged from 12 to 97 mol%, which enriches the experimental data of the phase diagram of the NaCl+RbCl system.
(2)With the help of a computational program that we developed, all available experimental data on phase equilibria and other thermodynamic data, also including the measurements in this work, were evaluated. The obtained model parameters were used to complete the calculations of the phase diagram, enthalpy of mixing, and activity for the NaCl+RbCl system. The model parameters ∆g AB of the modified quasichemical model in Equation (2) obtained by optimization are: ∆g NaRb/Cl = −1137.35 + 1.11T/K + (−34.22 − 0.75T/K)χ NaRb +(55.08 + 0.42T/K)χ RbNa J/mol (7) For solid solution α (NaCl is the main body), the optimized excess Gibbs energy parameter ∆ xs G α in Equation (5) is obtained as: For solid solution β (RbCl is the main body), the optimized excess Gibbs energy parameter ∆ xs G β in Equation (6) is obtained as: ∆ xs G β NaRb/Cl = x NaCl x RbCl (9587.93 + 92749.24x NaCl + 2014.85x RbCl ) J/mol (9)

Conclusions
In this study, differential scanning calorimetry analysis of the various compositions in the NaCl+RbCl system was performed. A thermodynamic data set for the NaCl+RbCl system was generated to reproduce all available experimental data. The following conclusions were obtained: (1) The phase diagram of the NaCl+RbCl system was measured using DSC, and the measured molar percentages of RbCl of the binary mixture ranged from 12 to 97 mol%, which enriches the experimental data of the phase diagram of the NaCl+RbCl system.
(2) With the help of a computational program that we developed, all available experimental data on phase equilibria and other thermodynamic data, also including the measurements in this work, were evaluated. The obtained model parameters were used to complete the calculations of the phase diagram, enthalpy of mixing, and activity for the NaCl+RbCl system.
(3) The phase diagram measurements and thermodynamic optimization obtained in this study provide a foundation for the subsequent completion of the phase diagram measurements and thermodynamic optimization of the multicomponent system containing NaCl+RbCl.