Engineering topics
Kim, Hojong
Publications and source records attributed to Kim, Hojong.
Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI
The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.
Determination of Kinetic Properties of Ni(II) Ions in Molten LiF-NaF-KF via Voltammetry
Kinetic properties of Ni(II) in eutectic LiF-NaF-KF (FLiNaK) molten salt were determined at T = 748–823 K using cyclic voltammetry (CV), square wave voltammetry (SWV), and chronoamperometry (CA) measurements using a glassy C working electrode, Ni(II)/Ni reference electrode, and Ni counter electrode. Reduction of Ni(II) to Ni(s) was determined to be a single step, two-electron transfer process. Diffusivity values were calculated using the Berzins and Delahay equation and semi-integral electroanalysis from the CV measurements as well as using the Cottrell equation from the CA measurements. Diffusivity of Ni(II) in molten FLiNaK at T = 748–823 K was determined to be 3.27 × 10 –7 –3.04 × 10 –6 cm 2 s –1 with an activation energy of 62–104 kJ mol –1 . The estimated kinetic properties varied appreciably among methodologies possibly due to inherent assumptions in theory regarding reversibility and unit activity during Ni metal deposition on the working electrode.
Electrochemical properties of Gd(III) ions in LiCl-KCl-GdCl 3 at 723–1023 K
Electrochemical behavior of Gd(III) ions in molten LiCl-KCl-GdCl 3 was investigated at 723–1023 K via cyclic voltammetry using tungsten as a working electrode, Gd-Bi (mole fraction, x Gd = 0.16) as a reference electrode, and Gd-Bi (x Gd = 0.02) as a counter electrode. A single reduction–oxidation wave was observed, confirming a single-step, 3-electron transfer Gd(III)/Gd transition. The cathodic peak potential exhibited minimal change (<13 mV) over a wide range of scan rates (0.05–0.30 V s –1 ), indicating facile charge transfer kinetics (i.e., a reversible electrode process). A nucleation overpotential associated with solid Gd deposition was observable at low temperatures (T < 823 K). The mass transport properties of Gd(III) ions were estimated using the Berzins and Delahay relation based on diffusion-limiting peak current. The diffusivity values were determined to be D Gd(III) = 0.5–2.7 × 10 –5 cm 2 s –1 at 723–1023 K with an associated activation energy of E a = 33.9 (±1.0) kJ mol –1 . Finally, the two-phase [liquid + GdBi] Gd-Bi alloy reference electrode experienced less than 0.5 mV of drift over 5 days of repeated electrochemical measurements, indicating high stability.
Electrochemical recovery of Nd using liquid metals (Bi and Sn) in LiCl-KCl-NdCl 3
Highly efficient recovery of Nd into liquid metals of Bi and Sn was achieved in molten LiCl-KCl-NdCl 3 electrolyte at 773–973 K by leveraging the strong interactions of Nd with liquid metals. Based on the emf measurements of Nd-Sn and Nd-Bi alloys, the activity values of Nd were determined as low as 1.1–5.8×10 –13 in both liquid metals at 973 K while the solubility of Nd was found to be 1.46 mol% in Sn and 5.65 mol% in Bi. Both liquid metals demonstrated high round-trip coulombic efficiencies (>99.3%) during deposition-removal cycles of 10–50 mA cm –2 and high recovery capacity up to approximately 20 mol% Nd beyond the solubility limit. In addition, a high Nd recovery yield (84–90%) with respect to the applied charge was confirmed based on chemical analysis of electrolysis products in Bi after constant current electrolysis (–50 mA cm –2 ) at 873–973 K. Overpotentials during the Nd deposition process were attributed to charge-transfer and mass-transport resistances based on the current-potential curve and electrochemical impedance spectroscopy. The charge-transfer kinetics of Nd deposition into liquid metals was facile with high exchange current densities at ~220 mA cm –2 . Finally, the exceptionally high recovery efficiency for Nd in the molten chloride is thought to result from strong chemical interactions (i.e., low activity) of Nd in liquid metals that encourage one-step reduction, i.e., Nd 3+ + 3e → Nd(in Bi or Sn) by effectively suppressing side reaction pathways from multivalent states (Nd 2+ and Nd 3+ ).
Thermodynamic properties of the Nd-Bi system via emf measurements, $\mathrm{DFT}$ calculations, machine learning, and $\mathrm{CALPHAD}$ modeling
Thermodynamic properties of the Nd-Bi system were investigated using a combination of experimental measurements, first-principles calculations based on density functional theory (DFT), data mining and machine learning (DM + ML) predictions, and calculation of phase diagrams (CALPHAD) modeling. The electromotive force (emf) of Nd-Bi alloys in molten LiCl-KCl-NdCl 3 at 773–973 K was measured via coulometric titration of Nd into Bi for the determination of thermochemical properties such as activity coefficients and solubilities of Nd in Bi. A new peritectic reaction of [liquid + NdBi 2 = Nd 3 Bi 7 ] at 774 K was confirmed using differential scanning calorimetry, structural (X-ray diffraction), and microstructural (scanning electron microscopy) analyses. The unknown crystal structure of NdBi2 was suggested to be a mixture of the anti-La 2 Sb configuration and the La 2 Te-type configuration based on ML predictions for over 26,000 data-mined AB 2 -type configurations together with DFT-based verifications. Using the newly acquired experimental data and DFT-based calculations, the thermodynamic description of the Nd-Bi system was remodeled, and a more complete Nd-Bi phase diagram was calculated, including the Nd 3 Bi 7 compound, invariant transition reactions, and liquidus temperatures.
DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations
In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .
Thermodynamic properties of Gd-Bi alloys determined by emf measurements in LiCl-KCl-GdCl 3 electrolyte
Thermodynamic properties of binary Gd-Bi alloys (mole fraction, x Gd = 0.02–0.35) were determined using electromotive force (emf) measurements in molten LiCl-KCl-GdCl 3 electrolyte in complement with structural, microstructural, and thermal characterization. X-ray diffractometry (XRD) and scanning electron microscopy (SEM) identified the phase constituents of Gd-Bi alloys (x Gd = 0.02–0.35) as Bi matrix and GdBi compound. According to differential scanning calorimetry (DSC) measurements up to 1250 K, a eutectic transition [L = Bi + GdBi] was observed at 543 K but the reported peritectic transition [L + GdBi = GdBi 2 ] was not detected, suggesting the absence of the GdBi 2 compound. Based on structural, microstructural, and thermal characterization of Gd-Bi alloys (x Gd = 0.02–0.35), the GdBi 2 compound was not detected, and its thermodynamic stability is questionable. The emf of Gd-Bi alloy (x Gd = 0.16) at 700–1048 K relative to Gd(s) was measured by electrodepositing pure Gd metal at 25 K increments. Further, the emf values of Gd-Bi alloys were determined via coulometric titration of Gd into liquid Bi using a two-phase [L + GdBi] alloy as a reference electrode at 773–973 K. The solubility of Gd in liquid Bi was also estimated to be 0.50 mol% (773 K), 0.99 mol% (873 K), and 2.21 mol% (973 K) and the excess partial molar Gibbs energy (Δ$\mathrm{\bar{G}^{ex}_{Gd}}$) of liquid Gd-Bi alloys was as low as 65 kJ mol –1 , indicating strong chemical interactions between Gd and Bi.
Multi-element liquid metal battery
An electrochemical cell including: a negative electrode including calcium and an alkali metal; a positive electrode including one or more elements selected from the group consisting of Al, Si, Zn, Ga, Ge, Cd, In, Sn, Sb, Hg, Tl, Pb, Bi, Te, Bi, Pb, Sb, Zn, Sn and Mg; and an electrolyte including a salt of calcium and a salt of the alkali metal. The electrolyte is configured to allow the cations of the calcium and alkali metal to be transferred from the negative electrode to the positive electrode during discharging and to be transferred from the positive electrode to the negative electrode during charging. The electrolyte exists as a liquid phase and one or both of the negative electrode and the positive electrode exists as liquid or partially liquid phases at operating temperatures of the electrochemical cell.