1. Introduction
Electrochemical processing utilizing a high temperature molten salt electrolyte, namely Pyroprocessing, is a widely sought technology for used nuclear fuel processing to recover actinides and minimize radioactive waste [1-7]. To be prepared for a large-scale facility deployment, sensor and accounting technologies have received significant attention to establish material accountability consistent with traditional international safeguard methods [8-13]. Within this system, the electrorefiner (ER) is an electrochemical cell that dissolves used metallic fuel anodically and collects metallic actinides cathodically. The molten salt electrolyte in the ER contains the majority of special nuclear materials in the pyroprocessing system [14-20]. Molten salt electrolyte composition sensors are of interest for safeguards due to the typical accumulation of plutonium and other actinides in this salt.
A variety of methods have been pursued to assess the inventory of elements in the molten electrolyte of the pyroprocessing system [21-40]. Among those, electrochemical methods can form the basis of such sensors, since the various species have different electrochemical properties, and there are electrode materials compatible with the harsh environmental conditions found in the ER. Such electrochemistry- based techniques as cyclic voltammetry (CV) and chronopotentiometry (CP), provide simple yet powerful means to quantify species concentrations in molten salts let alone their traditional usages for identifying electrochemical reactions, and measuring electrochemical properties such as diffusion and activity coefficients. One can find numerous applications of these electrochemical probing techniques to molten salts containing uranium [21-25] as well as species that contribute to the composition of the used fuel; plutonium [22,26], americium [22], lanthanum [22], neodymium [22], samarium [27-30], thorium [31], gadolinium [32], cerium [33-35], erbium [36], and neptunium [22,37].
While the aforementioned techniques rely on direct current (DC), electrochemical impedance spectroscopy (EIS) adopts alternating current (AC). Impedance techniques are useful to a variety of industrial applications, e.g. fuel cells [41], solar salt systems [42] and anti-corrosion coatings [43]. As an analytical tool, via frequency and amplitude control, the EIS measurements provide the frequency domain characterization and a unique electrochemical perspective. For example, EIS methods are valuable for elucidating reactions and associated electrochemical mechanisms at electrode/ electrolyte interfaces [44]. More specifically, such parameters as reaction rate constant [45-47], corrosion rate [48-50], electrolyte conductivity [51-54], double-layer capacitance [55,56], and exchange current density [57-59] can be determined by EIS techniques.
A selection of reported applications of EIS for the characterization of electrochemical reactions in LiCl-KCl molten salts used in pyroprocessing applications are briefly described herein. Reddy et al. [60] reported the EIS investigation of uranium chloride in LiCl-KCl eutectic at 435℃. In this application EIS spectra in the range of 0.1 to 65 kHz were obtained at potentials ranging from -0.15 to -1.55 V relative to a Ag/AgCl reference electrode. Four equivalent circuit models were used across the potential range to analyze the EIS data. The analysis focused on qualitatively explaining the diffusion, absorption and charge transfer processes associated with the electrochemical reactions occurring at the working electrode. Ghosh et al. [61] used EIS data to calculate the heterogeneous rate constant associated with the reduction of CeO+ at a glassy carbon electrode at 550℃ in MgCl2-NaCl-KCl eutectic. An equivalent circuit was used to determine the charge transfer resistance from an EIS spectrum which had an observable Warburg behavior at frequencies down to 100 Hz. Similarly, Vandarkuzhali et al. [62] used EIS data to calculate the heterogeneous rate constant associated with the reduction of lanthanum at a tungsten electrode at 525℃ in LiCl-KCl using data down to 10 Hz that exemplified a Warburg behavior at -0.25 V. Ferry and Picard [63] investigated the Ti4+/Ti3+ reduction at 470℃ by analyzing EIS data with an equivalent circuit and obtaining the rate constant for this reaction.
Although there have been efforts to characterize the electrochemical properties with EIS, to the best of the authors’ knowledge, there has been no reported attempt to correlate impedance response to the concentration of a given species in a LiCl-KCl molten salt. The intent of this paper is to introduce the use of EIS for the prescribed purpose. Also provided in this paper is the application of CV technique to estimate the species concentration with the aid of the Randles- Sevcik equation.
To this end, SmCl3 in the range of 0.5 to 3 wt% in LiCl- KCl eutectic was characterized using CV and EIS at 500℃. CV was used to calculate an average diffusion coefficient for Sm3+ in the concentration range, while EIS was used to calculate a relationship between polarization resistance and bulk SmCl3 concentration. Comparison of calculated bulk SmCl3 concentrations based on CV and EIS with experimental bulk concentrations demonstrated the applicability of these methods for monitoring concentration changes in LiCl-KCl melts.
In the effort reported here, the Sm3+/Sm2+ redox reaction was characterized using EIS. An equivalent circuit composed of a series of resistor-capacitor pairs was used to model spectra. The model was also used to estimate the polarization resistance by extrapolation. The influences of applied potential and SmCl3 concentration were determined by comparing the resistor-time constant pairs and the polarization resistance.
2. Experimental
The experiments were conducted in an MBraun glovebox containing an argon atmosphere with H2O and O2 concentrations less than 5 ppm. A Kerr Lab Electro-Melt furnace was located in the glovebox and used to heat the molten salts to 500℃. Eutectic LiCl-KCl salt (99.99%) was obtained from AAPL and anhydrous SmCl3 (99.99%) was obtained from Alfa Aesar. The use of anhydrous salts and under low H2O and O2 conditions is meant to minimize the effects of oxide and oxychloride formation in the salt. A tapered cylindrical glassy carbon crucible obtained from HTW Germany was used for containing the salts, which was 16.8 cm in height with a wall thickness of 3 mm. The tapered cylinder had outside diameters of 45 mm and 41 mm at the top and bottom, respectively. The glassy carbon crucible was placed inside the graphite crucible of the furnace.
The electrode assembly consisted of working, counter, and reference electrodes together with a thermocouple for temperature measurement. The working electrode was a 2.0 mm diameter tungsten rod (99.95%, 30 cm in length) obtained from Alfa Aesar. The counter electrode was a 3.0 mm diameter glassy carbon rod (30 cm in length) obtained from HTW Germany. The reference electrode was constructed using a Pyrex tube (10 mm outside diameter and 1.0 mm wall thickness) unto which a Pyrex rod was attached. One end of a 1.0 mm diameter silver wire (99.9%), obtained from Acros Organics, was made into a spiral sufficiently small in diameter to be inserted into the tube. 0.135 g of AgCl beads and 1.00 g of the LiCl-KCl eutectic were placed in the tube. Upon melting, a 5.0 mol% AgCl mixture was obtained that served as the reference electrode. An alumina tube (2 mm and 3 mm inner and outer diameters, respectively) closed at one end was used to sheath the K-type thermocouple.
The electrode assembly was made by placing the working and counter electrodes and the silver wire connect of the reference electrode in alumina tubes of 2 mm and 3 mm inner and outer diameters, respectively, as shown in Fig. 1. The tungsten working electrode was secured in the alumina tube such that 14 mm of it was exposed at one end, the end that would be the working electrode. An alumina paste was used to seal the alumina tube. For the purpose of calculations, the geometric surface area (0.911 cm2) of the working electrode is used. The glassy carbon counter electrode was not secured in the alumina tube such that it could contact the bottom of the crucible, effectively making the entire crucible the counter electrode.
The CV and EIS experiments were conducted using a Princeton Applied Research VersaSTAT 4-400 potentiostat with VersaStudio software. The electrical leads from the potentiostat were connected to electrical feed-through fittings in the glovebox. The leads from these fittings were connected to the electrodes using alligator clips.
The experimental procedure initiated with mixture preparation. For the first experiment, 80 g of LiCl-KCl eutectic and 0.41 g of SmCl3 were introduced into the crucible, with the mixture comprised of 0.5 wt% (3.23 × 10-5 mol·cm-3) SmCl3. The furnace set-point was adjusted to 500℃ with a ramp-rate of 4℃ min-1. Starting from room temperature, this resulted with a heating time of 115 min. The electrode assembly was lowered into the crucible to a position where the working and reference electrodes and thermocouple were approximately 1 cm above the surface of the salt mixture when the furnace temperature was 350℃. At a furnace temperature of 450℃, the salt was assumed to be molten and the electrode assembly was lowered slowly into the melt. Upon contact of the crucible by the glassy carbon counter electrode rod, the assembly was secured in position. The assembly was made such that the exposed end of the working electrode would be totally immersed; approximately 60 % of the Pyrex tube of the reference electrode was also immersed. The thermocouple was immersed to the same extent as the working electrode. The molten salt with electrode assembly immersed was left for 6 h to allow the mixture to become homogeneous before beginning the measurements. This consisted of three scans at a given scan rate (0.1 to 1.5 V·s-1).
The EIS experiments were conducted afterwards using the same three-electrode set up with the frequency range of 10 kHz to 0.01 Hz and perturbation amplitude of 10 mV. EIS measurements were acquired at bias potential of -0.84, -0.93 and -1.03 V. Six EIS spectra were acquired at each bias potential. A time period of 1 min was allowed between EIS measurements at the same bias potential while a period of 3 min was allowed when the bias potential was changed. It was found that electrode cleaning steps between measurements was not necessary because the reduced species is soluble. Upon completion of the electrochemical experiments, the furnace was switched off and the electrode assembly was retracted from the molten salt to a position approximately 2 cm above the surface of the melt. The salt was allowed to freeze completely. The experimental procedure was repeated except that only SmCl3 was added to the existing salt mixture to increment its concentration to 1.0 wt% (6.31 × 10-5 mol·cm-3), and subsequently to 1.5 (9.47 × 10-5 mol·cm-3), 2.0 (1.26 × 10-4 mol·cm-3), 2.5 (1.58 × 10-4 mol·cm-3) and 3.0 wt% (1.89 × 10-4 mol·cm-3). It is noted the concentration of SmCl3 will be presented as wt % in the remainder of this manuscript, however the molar concentration has been provided above for convenience.
3. Results
The CV results associated with the 1.0, 2.0 and 3.0 wt% concentrations of SmCl3 in the eutectic are shown in Fig. 2 where each scan was the third of three scans for a given scan rate.
The potential range of 0 V to -1.6 V was used for the scans in order to observe the reduction of Sm3+ to Sm2+ and the oxidation of Sm2+ to Sm3+. The reduction of Sm2+ to Sm occurs at potential more negative than the reduction of Li+, a constituent of the eutectic. The voltammograms consisted of an anodic wave occurring at approximately -0.84 V and a cathodic wave at approximately -1.03 V. There were neither pre-waves nor post-waves associated with either peak. The anodic peak potential, , corresponded to the electrooxidation of Sm2+,
and the cathodic peak potential, , the electro-reduction of Sm3+,
Although the experimental conditions were tightly controlled, it is possible that oxide or oxychloride formation occurred, albeit to a very small extent. However, as shown by Fig. 2., there are no unidentified current peaks in the CV plot. Further, it has been reported that oxide and oxychloride of samarium are electrochemically inactive in the potential window of LiCl-KCl eutectic [64].
The values of and as functions of concentration are shown in Fig. 3 where each value was the value associated with the third CV scans for a given concentration. It was observed there was no significant influence on and potentials by the concentration. The electrode potential E1/2 was obtained by
and is also shown in Fig. 3 as a function of SmCl3 concentration. The values of , and E1/2 for the first EIS experiment were based on the 0.5 wt% CV data and were -0.84, -1.03 and -0.93 V, respectively. These values were used for the remaining EIS experiments as there was consistency in the and values as shown in Fig. 3. The value of E1/2 at -0.93 V was close to the standard potential of the redox Sm3+/Sm2+ system in LiCl-KCl that was reported as -0.97 ± 0.02 V relative to Ag/AgCl in literature [27].
The Nyquist and Bode phase angle plots of the EIS data associated with the 1.0 wt% SmCl3 system obtained at a
bias potential of -0.93 V, E1/2, is shown in Fig. 4 for the frequency range of 10 kHz to 0.01 Hz. There was inductive behavior in the high frequency region of the Nyquist plot not easily observed due to scale of the plot, however it is more clearly indicated by the positive phase angle in the Bode phase diagram. There was a small amount of scatter at the low frequency region observed in the Nyquist plot. The repeatability of the EIS spectra is demonstrated by the overlaying of the experimental data as highlighted in the Bode phase angle plot.
Nyquist plots of the EIS data for the 1.0 wt% SmCl3 system are shown in Fig. 5 where the bias potentials were -0.84 V and -1.03 V. These potentials corresponded to the and , respectively. Inductive loops were present in both data sets although not observable in this figure. The repeatability of the data at the potential was demonstrated by overlaying the EIS spectra. There was repeatability of the EIS data down to 0.1 Hz for the potential with observable scatter at the low frequency region for frequencies lower than 0.1 Hz. Comparison of the data for the 1.0 wt% SmCl3 system at the , and E1/2 potentials indicated that there was more scatter at low frequencies for the most negative potential as compared to the and E1/2 potentials. This trend was consistent with the different concentrations of SmCl3 investigated.
The Nyquist and Bode phase angle plots of the EIS data associated with the 1.0 wt% SmCl3 system obtained at a bias potential of -0.93 V, E1/2, is shown in Fig. 4 for the frequency range of 10 kHz to 0.01 Hz. There was inductive behavior in the high frequency region of the Nyquist plot not easily observed due to scale of the plot, however it is more clearly indicated by the positive phase angle in the Bode phase diagram. There was a small amount of scatter at the low frequency region observed in the Nyquist plot. The repeatability of the EIS spectra is demonstrated by the overlaying of the experimental data as highlighted in the Bode phase angle plot.
Nyquist plots of the EIS data for the 1.0 wt% SmCl3 system are shown in Fig. 5 where the bias potentials were -0.84 V and -1.03 V. These potentials corresponded to the and , respectively. Inductive loops were present in both data sets although not observable in this figure. The repeatability of the data at the potential was demonstrated by overlaying the EIS spectra. There was repeatability of the EIS data down to 0.1 Hz for the potential with observable scatter at the low frequency region for frequencies lower than 0.1 Hz. Comparison of the data for the 1.0 wt% SmCl3 system at the , and E1/2 potentials indicated that there was more scatter at low frequencies for the most negative potential as compared to the and E1/2 potentials. This trend was consistent with the different concentrations of SmCl3 investigated.
The Nyquist and Bode phase angle representations of the EIS data associated with spectra obtained at the different concentrations of SmCl3 at a bias potential of -0.93 V are shown in Fig. 6. A smaller value in the real axis was observed at the low frequency region and was associated with a larger concentration. There was overlay of the data in the Bode phase angle plot at frequencies less than 1 Hz and greater than 1000 Hz. The influence of the concentrations was observable in the frequency range from 1 Hz to 1000 Hz where the response was predominantly capacitive.
4. Discussion
The analysis of the CV and EIS data detailed in this section is focused at demonstrating the applicability of these methods for monitoring the concentration of SmCl3 in the LiCl-KCl eutectic. The CV approach required the calculation of the diffusion coefficient of Sm3+, while the EIS approach required obtaining a relationship between the polarization resistance and the concentration. The diffusion coefficient of the Sm3+ ion, was calculated using the Randles-Sevcik equation
where D is the diffusion coefficient (cm2·s-1), Ip is the cathodic current (C·s-1), c0 is the bulk concentration (mol· cm-3), S is the surface area of the electrode (cm2), and n is the number of electrons transferred, in this case n=1. In this equation F is the Faraday constant (C·mol-1), T is the temperature of the electrolyte, which was 500℃ (773 K), and R is the universal gas constant (J·K·mol-1). The plot of Ip as a function of ν0.5 is shown in Fig. 7 with the concentration of SmCl3 as a parameter.
As expected from Equation 4, there was a linear relationship between Ip and ν0.5 for a given concentration. Linear fits by regression of the data associated with a given concentration yielded the overlaid line on a data set. The slope of a given line was used to calculate the value of , shown in Fig. 8. as a function of concentration. The linearity of the curves in Fig. 7 along with the consistency of peak potentials shown in Fig. 2 indicates the system is reversible and supports the use of Equation 4. The average and standard deviation of the diffusion coefficient, , for the range of 0.5 to 3.0 wt% (0.03 mmol·cm-3 to 0.20 mmol·cm-3) was (0.647 ± 0.076) × 10-5 cm2·s-1. There appears to be no systematic trend with respect to concentration.
The calculated values of the bulk concentration are shown as a function of the experimental c0 in Fig. 9 for a scan rate of 1 V·s-1. The calculated c0 values were obtained using Equation 4 together with the value and the experimental Ip value associated with a given concentration. There was agreement between the calculated and experimental c0 values as demonstrated by the proximity of the data points to the diagonal dashed line. The largest error relative to the experimental c0 values was associated with the lowest concentration at approximately 9% while the relative errors for the other concentrations were less than 2%.
A reported value for the diffusion coefficient of Sm3+ ions, , is 1.30 × 10-5 cm2·s-1 that was obtained at a concentration of 0.04 mol·kg-1, approximately 0.065 mmol·cm-3, using a similar CV approach [27]. The asterisk superscript denotes the value of diffusivity was taken from literature. was different than the value of , which is reported here as 0.647 × 10-5 cm2·s-1. It was also significantly different than the value of 0.759 × 10-5 cm2·s-1 (reported here) that was associated with a concentration of 0.0637 mmol·cm-3, which is close to that of 0.065 mmol·cm-3. The influence of the literature diffusion coefficient value was investigated by using this value together with the experimental Ip values to calculate c0 values. These calculated c0 values are shown as a function of the experimental c0 values in Fig. 9. It can be seen there was less agreement between the experimental and calculated concentrations with an increase in SmCl3 concentration. The value which was calculated using the Ip values associated with concentrations ranging from 0.0317 to 0.195 mmol·cm-3 yielded very good agreement between the experimental and calculated bulk concentration values.
It is interesting that the other reported values for at 500℃ and 0.065 mmol·cm-3 were 1.41 × 10-5 and 0.98 × 10-5 cm2·s-1 obtained using semi-integral and CP methods, respectively. This demonstrates the calculation of the diffusion coefficient is sensitive to the experimental method. Also, it is likely that variation in reported diffusion coefficients is a result of the difficulty associated with accurately measuring the true surface area of the working electrode and the concentration of the electrochemically active species at the electrode/electrolyte interface. Table 1 lists values for D found in literature. Although there is no significant influence of the concentration on the calculated value of in this study, the average value over the range of concentrations was better suited for use in monitoring concentration changes. These findings suggest it may be best practice to determine and utilize a diffusion coefficient that is calculated for a given system rather than relying solely on values from literature.
For the purpose of correlating polarization resistance to concentration of SmCl3, analysis of the EIS data was performed using an equivalent model consisting of a solution resistance in series with several Voigt elements (a measurement model) as shown in Fig. 10. The measurement model method (curve fitting) is a valid approach to implicitly analyze impedance spectra, whereas a process model is phenomena-based [65,66]. In this study, the measurement model was sufficient to estimate the polarization resistance. The model was fitted to data in the range from 0.01 to 1000 Hz where the impedance response is predominantly capacitive.
There was no Warburg element included, as a diffusion- limited process was not observed at the low frequency end of the EIS data. Diffusion effects are usually indicated by a 45-degree line in the low frequency region of a Nyquist plot. The impedance of such a model can be expressed as
where Ri represents the resistance (Ω·cm-2) for Voigt element i [65]. The characteristic time constant associated with a Voigt element is τi = RiCi, where Ci is the capacitance (farad) of Voigt element i. The characteristic frequency fi (Hz) and characteristic time constant are related by 2πfiτi = 1. The model as shown in Fig. 10 has been used to distinguish between bias and stochastic errors associated with impedance data [67]. The model was used to analyze the EIS data using a complex non-linear least squares algorithm where modulus weighting was used. This involved minimizing the function
where n is the number of data points, m is the number of parameters, and p is the vector of parameters. The variables Zr(ωi;p) and Zj(ωi;p) are the calculated real and imaginary values, respectively, which are functions of the angular frequency ωi (rad·s-1) and p. The measured real and imaginary values were Zr(ωi) and Zj(ωi), respectively. A down-hill simplex method was used in the regression procedure. Analytical software was written using FORTRAN to analyze the data. The Nyquist and Bode phase angle representations of the second EIS spectrum associated with the 1.0 wt% SmCl3 system at a bias potential of -0.93 V is shown in Fig. 11. for the frequency range of 1000 Hz to 0.01 Hz.
The optimal fit of the model to the data is included in Fig. 11, where the optimal fit was a six Voigt element model. This was determined by sequentially fitting one, two, three, four, five, six and then seven elements. The agreement between the model and the data improved with each addition of an element. However, the inclusion of a seventh element resulted in parameter errors that included zero. There was good agreement between the six-element model and the data observable in the Nyquist and Bode representations. Regression of the model to the individual measured spectra resulted in chi-squared values that were all less than three (χ2 < 3). These values of χ2 compare favorably with the critical value of 35, which is associated with 50 degrees of freedom and 95% probability.
The values of Ri and τi listed in Table 2 were obtained from fitting the second EIS spectra for the 1.0 wt% SmCl3 system to the six-element Voigt model. The parameters associated with this model were used to calculate the impedance in the frequency range 0.01 to 0.0001 Hz, thereby extrapolating to the polarization resistance as shown in Fig. 11 by the dotted line. The six spectra associated with the 1.0 wt% SmCl3 system were each fit to the six-element Voigt model with the Ri values shown in Fig. 12 as a function of the associated values of τi where the error bars represent the errors obtained from the regression procedure.
Repeatability was observed of the overlay of Ri values for a given τi for the six spectra analyzed. Characterization of the EIS response for a given concentration was given by the average of Ri and τi values from the six spectra. The average values Ri and τi calculated from the six parameter sets associated with the six spectra of the 1.0 wt% SmCl3 system are listed in Table 2 where the standard deviation of a given average value is calculated from the six regressed parameter values.
The Nyquist and Bode phase representations of the EIS data associated with the 1.0 wt% system are shown in Fig. 11 with the applied potential as a parameter. Included is the fit of the six-element Voigt model to each data set and extrapolation to the polarization resistance. The low frequency data for the -1.03 V spectrum was more scattered than the -0.84 V spectrum based on the Nyquist representation. The data for these two spectra overlaid in the Bode phase angle representation. The Bode phase angle representation for the -0.93 V spectrum was different than that of the other applied potentials for frequencies greater than 1 Hz. The Ri and τi model parameter values obtained from the six spectra associated with the bias potentials of -0.84, -0.93 and -1.03 V are shown in Fig. 13.
The similarity in the EIS response among the three bias potentials was demonstrated in the time constants being approximately the same for a given element. This indicated that the EIS responses of each bias potential was composed of similar line shapes suggesting that the same mechanisms were involved at the three potentials. The influence of the bias potential was observed in the resistor value for elements with a τi value greater than 0.5 s where the Ri value for the applied E1/2 bias potential was less than the values for and potentials which were similar. This suggested the resistance to any reactions being characterized by the Voigt elements was less at E1/2 as compared to and . This would have been expected as the reversible reaction would be most favorable.
The trend of a smaller Ri value for the E1/2 bias potential as compared with the potential was also observed for the 1.5, 2.0, 2.5 and 3.0 wt% SmCl3 systems. A smaller Ri value was associated with a given τi value for the EIS data at the potential of -1.03 V for the 1 wt% and 1.5 wt% systems. There was scatter in the frequencies less than 0.1 Hz for the systems with higher concentrations and these were not analyzed. The plot of Rp as a function of the applied potential with the concentration of SmCl3 as a parameter is shown in Fig. 14. Each data point represents the average of the six Rp parameters obtained from the six spectra associated with a given condition and the error bar represents one standard deviation.
The EIS data associated with the -1.03 V potential and greater than 1.5 wt% SmCl3 were not analyzed. The Rp value was the lowest for the E1/2 potential of -0.93 V at 1.0 and 1.5 wt% SmCl3 as compared to the -1.03 V and -0.84 V potentials, and , respectively. The Rp for the E1/2 value was lower than the -0.84 V potential for concentrations greater than 1.5 wt % SmCl3. This was consistent where the least resistance was associated with the E1/2 bias potential was also observed with the Ri parameters.
The Nyquist and Bode phase representations of the EIS data associated with -0.93 V applied potential are shown in Fig. 15 with the concentration of SmCl3 as a parameter. The fit of the six-element Voigt model to each data set is included together with the extrapolation to approximate the polarization resistance.
The EIS data for the 0.5 wt% SmCl3 system was analyzed for the 1000 Hz to 0.1 Hz. The influence of the concentration was seen on the Nyquist plot, where a smaller semicircular loop and Rp value was associated with a larger concentration. The influence of the concentration was observable in the Bode phase angle plot, where for frequencies greater than 0.1 Hz a more negative phase angle was associated with a larger concentration for a given frequency. The Bode phase angle plot also indicated that the phase angle data at the low frequency end overlaid for frequencies less than 0.1 Hz. The agreement between the six-element Voigt model and the EIS data was demonstrated by the overlay of the model on the data observable both in the Nyquist and Bode phase angle plots.
The value of Ri as a function of τi is shown in Fig. 16 with concentration as a parameter for the -0.93 V applied potential condition. The τi value was approximately equal for a given element except for the values for the 0.5 wt% system, which was slightly less than the value for the other systems for a given element. This indicated that the mechanisms involved in the EIS response were the same for each concentration, as would be expected. Concentration influenced the Ri value where a smaller Ri value was associated with a larger concentration, for a given element. This indicated that the resistance of a mechanism characterized by an element was less when there was an increase in the concentration, as expected.
The values of Rp and Rsol as functions of concentration are shown in Fig. 17(a) and 17(b), respectively. The Rp and Rsol values at a given concentration associated with the six spectra obtained for that concentration overlay on each other and are not distinguishable. The larger scatter in the values associated with the 0.5 wt% system was attributed to the lowest frequency being measured for these spectra was 0.1 Hz, while for the other concentrations it was 0.01 Hz. Overlaid on these plots are the average values, Rp and Rsol, as functions of concentration. It is observed from Fig. 17(a) that there was a monotonic decrease in Rp with concentration as highlighted by the trend associated with the Rp values.
A monotonic trend was also observed for the Rsol values however, this excluded the largest concentration. The lower value of Rp and Rsol associated with a higher concentration was expected. A higher concentration of SmCl3 was more favorable for the reversible reaction with there being a lower resistance exhibited and the solution resistance was lowered due to an increase in ionic conductivity.
The experimental c0 values were cast as a function of the Rp values as shown in Fig. 18 from which an equation of the form
was derived by fitting a power function to the data. The Rp values associated with a given experimental c0 values were used to calculate a c0 values using this equation. The calculated c0 values were cast as a function of the experimental c0 value in Fig. 9. There was agreement between the experimental and calculated c0 values with the relative errors ranging from 3% to 13% for the concentration range investigated.
5. Conclusions
This paper reports both cyclic voltammetry (CV) and electrical impedance spectroscopy (EIS) were effective for the measurement of SmCl3 concentration (0.5 to 3.0 wt%) in molten eutectic LiCl-KCl at 500℃. The CV approach with the Randles-Sevcik equation requires knowledge of the diffusion coefficient, which can either be obtained from the literature or via CV of samples with known concentrations. In this study, an average diffusion coefficient for Sm3+ in the salt () was determined via the latter approach. For the EIS approach, a polarization resistance/concentration relationship was determined. This demonstrated the consistency of using electrochemical based methods for monitoring the concentration changes in a molten salt such as the LiCl- KCl system. We recommend that an averaged diffusion coefficient of a species in the concentration range should be used for the application of the developed CV approach to concentration estimation. Also, it is recommended for better performance of the developed CV method that diffusion coefficients should be measured experimentally for the given salt matrix rather than relying upon values published in the literature for similar salt matrices. The EIS approach has benefits for implementation, since it does not require an independent measurement of the diffusion coefficient. Note that this study involves a soluble-soluble redox transition. Further study is needed to evaluate the application of EIS to redox transitions with an insoluble reduced species.