RESUMO
MgCl2 hydrates are considered as high-potential candidates for seasonal heat storage materials. These materials have high storage capacity and fast dehydration kinetics. However, as a side reaction to dehydration, hydrolysis may occur. Hydrolysis is an irreversible reaction, which produces HCl gas thus affecting the durability of heat storage systems. In this study, we present the parameterization of a reactive force field (ReaxFF) for MgCl2 hydrates to study the dehydration and hydrolysis kinetics of MgCl2·H2O and MgCl2·2H2O. The ReaxFF parameters have been derived by training against quantum mechanics data obtained from Density Functional Theory (DFT) calculations consisting of bond dissociation curves, angle bending curves, reaction enthalpies, and equation of state. A single-parameter search algorithm in combination with a Metropolis Monte Carlo algorithm is successfully used for this ReaxFF parameterization. The newly developed force field is validated by examining the elastic properties of MgCl2 hydrates and the proton transfer reaction barrier, which is important for the hydrolysis reaction. The bulk moduli of MgCl2·H2O and MgCl2·2H2O obtained from ReaxFF are in close agreement with the bulk moduli obtained from DFT. A barrier of 20.24 kcal mol(-1) for the proton transfer in MgCl2·2H2O is obtained, which is in good agreement with the barrier (19.55 kcal mol(-1)) obtained from DFT. Molecular dynamics simulations using the newly developed ReaxFF on 2D-periodic slabs of MgCl2·H2O and MgCl2·2H2O show that the dehydration rate increases more rapidly with temperature in MgCl2·H2O than in MgCl2·2H2O, in the temperature range 300-500 K. The onset temperature of HCl formation, a crucial design parameter in seasonal heat storage systems, is observed at 340 K for MgCl2·H2O, which is in agreement with experiments. The HCl formation is not observed for MgCl2·2H2O. The diffusion coefficient of H2O through MgCl2·H2O is lower than through MgCl2·2H2O, and can become a rate-limiting step. The diffusion coefficient increases with temperature and follows the Arrhenius law both for MgCl2·H2O and MgCl2·2H2O. These results indicate the validity of the ReaxFF approach for studying MgCl2 hydrates and provide important atomistic-scale insight of reaction kinetics and H2O transport in these materials.
RESUMO
Salt hydrates store solar energy in chemical form via a reversible dehydration-hydration reaction. However, as a side reaction to dehydration, hydrolysis (HCl formation) may occur in chloride based salt hydrates (specially in MgCl2 hydrates), affecting the durability of the storage system. The mixture of CaCl2 and MgCl2 hydrates has been shown experimentally to have exceptional cycle stability and improved kinetics. However, the optimal operating conditions for the mixture are unknown. To understand the appropriate balance between dehydration and hydrolysis kinetics in the mixtures, it is essential to gain in-depth insight into the mixture components. We present a GGA-DFT level study to investigate the various gaseous structures of CaCl2 hydrates and to understand the relative stability of their conformers. The hydration strength and relative stability of conformers are dominated by electrostatic interactions. A wide network of intramolecular homonuclear and heteronuclear hydrogen bonds is observed in CaCl2 hydrates. Equilibrium product concentrations are obtained during dehydration and hydrolysis reactions under various temperature and pressure conditions. The trend of the dehydration curve with temperature in CaCl2 hydrates is similar to the experiments. Comparing these results to those of MgCl2 hydrates, we find that CaCl2 hydrates are more resistant towards hydrolysis in the temperature range of 273-800 K. Specifically, the present study reveals that the onset temperatures of HCl formation, a crucial design parameter for MgCl2 hydrates, are lower than for CaCl2 hydrates except for the mono-hydrate.
