The “ionic” to “molecular” transitions in AlCl3 and FeCl3 as predicted by an ionic interaction model

Tracing mechanistic pathways and reaction kinetics toward equilibrium in reactive molten salts

In the dynamic environment of multi-component reactive molten salts, speciation unfolds as a complex process, involving multiple competing reaction pathways that are likely to face free energy barriers before reaching the reaction equilibria. Herein, we unravel intricate speciation in the AlCl3-KCl melt compositions with rate theory and ab initio molecular dynamics simulations. We find that the compositions with 100 and 50 mol% AlCl3 exclusively comprise neutral Al2Cl6 dimers and charged AlCl4- monomers, respectively. In intermediate AlCl3-KCl compositions, the chemical speciation proves to be a very complex process, requiring over 0.5 nanosecond to reach an equilibrium distribution of multiple species. It is a consequence of the competitive formation and dissociation of additional species, including charged Al dimers, trimers, and tetramers. Here, the species formation occurs through ion exchange events, which we explain by computing free energy landscapes and employing a Marcus-like rate theory. We show that both interspecies and intraspecies ion exchanges are probable and are dictated by the local structural reorganization reflected in the change of local coulombic fields. The species distributions are validated by comparing computed Raman spectra and neutron structure factors with the available experimental data. We find an excellent simulation-experiment agreement in both cases. Nevertheless, Raman spectroscopy turns out to be particularly advantageous for distinguishing between unique species distributions because of the distinct vibrational signatures of different species. The mechanistic insight into reaction dynamics gained in this study will be essential for the advancement of molten salts as reactive media in high-temperature energy applications.

Multicomponent diffusion in molten LiCl-KCl: dynamical correlations and divergent Maxwell-Stefan diffusivities

Multicomponent diffusional mechanisms in the ternary LiCl-KCl system are elucidated using the Green-Kubo formalism and equilibrium molecular dynamics simulations. The Maxwell-Stefan (MS) diffusion matrix is evaluated from the Onsager dynamical matrix that contains the diffusion flux correlation functions. From the temporal behavior of the correlation functions, we observe that the Li-Li and Li-Cl ion pairs have a pronounced cage dynamics that remains noticeably strong even at high temperatures. Even though the Onsager coefficients, which are the time integrals of the diffusion flux correlation functions, portray a relatively smooth variation across various compositions and temperatures, we observe a sign change and a divergent-like behavior for the MS diffusivity of the K-Li ion pair at a temperature of ~1100 K for the eutectic composition, and at a KCl mole fraction of ~0.49 at 1043 K. Negative MS diffusivities, while unusual, are however shown to satisfy the nonnegative entropic constraints.

Molecular dynamics simulation of LiTFSI-acetamide electrolytes: structural properties

The liquid structures of nonaqueous electrolytes composed of lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) and acetamide, with LiTFSI/acetamide molar ratios of 1:2, 1:4, and 1:6, were studied by molecular dynamics simulations. The simulations indicate that the Li+ cations prefer to be six-coordinate by the sulfonyl oxygen atoms of the TFSI- anions and the carbonyl oxygen atoms of the acetamide molecules, rather than by the most electronegative nitrogen atom of the TFSI- anion. Therefore, close Li+-TFSI- contact pairs exist in the system. The TFSI- anion prefers to provide only one of four possible oxygen atoms to coordinate to the same Li+ cation. Three conformations (cis, trans, and gauche) of the TFSI- anions were found to coexist in the liquid electrolyte. At high salt concentrations, the TFSI- anions mainly adopt the gauche conformation in order to provide more oxygen atoms to coordinate to different Li+ cations, while simultaneously reducing the repulsion among the Li+ cations. On the other hand, the fraction of TFSI- anions adopting the cis conformation is largest for the system with the molar ratio of 1:6, in which many clusters, mainly composed of the Li+ cations and the TFSI- anions, are immersed in the acetamide molecules. The size and charge distribution of clusters were also investigated. In the system with the molar ratio of 1:2, nearly all of the ions in the PBC (periodic boundary conditions) box aggregate into a bulky cluster that gradually disassembles into small clusters with decreasing salt concentration. The addition of acetamide molecules was found to effectively relax the liquid electrolyte structure, and the system with the molar ratio of 1:4 was found to exhibit a more homogeneous liquid structure than the other two electrolyte systems with molar ratios of 1:2 and 1:6.

Molecular-network-ionic structure transitions in liquid AlCl(3) and ZnCl(2) halogenides under pressure

We present the in situ high-pressure-high-temperature x-ray diffraction study of the liquid AlCl(3) and ZnCl(2) halogenides having a quasi-molecular network structure in liquid state at normal pressure. These liquids are intermediate between pure covalent and ionic melts. Structural study of these liquid halogenides is indicative of a rapid and strong breakdown of an intermediate-range order in a tetrahedral network of melts for the initial pressure range, 0-2.5 GPa for AlCl(3) and 0-1.8 GPa for ZnCl(2), and points to rather sharp transitions in liquids with the formation of a short-range order structure similar to ionic melt structures around 4 GPa for AlCl(3) and 3 GPa for ZnCl(2). Thus, pseudo-covalent liquid halogenides like AlCl(3) and ZnCl(2) provide testimony to two phenomena under high pressures, namely, a gradual decay of structural correlations in the tetrahedral network of the melt and a sharp transition from molecular-network to ionic structure in liquid on further compression. Such a two-stage structural transformation under pressure is the general feature for a wide class of simple melts, including most of the pseudo-covalent halogenides.

Short-Range Order in Liquid Aluminum Chloride: Ab Initio Molecular Dynamics Simulations and Quantum-Chemical Calculations

We use ab initio molecular dynamics simulations based on density-functional theory and quantum-chemistry calculations on molecular clusters to examine the structure of liquid AlCl3. In the past, conflicting descriptions of the short-range-order in molten AlCl3, based on either edge-sharing dimers or corner-sharing oligomers, have been proposed. This liquid also poses a simulation challenge, due to the possibility of ring-like trimers which can be metastable on the order of >10 ps. Simulations which begin with monomers, either random or ordered, appear to be able to produce proper ratios of ring-trimer to dimer-plus-tail molecular structures without the need to achieve long-time scale chemical equilibrium. Single-molecule calculations lend further support to the conclusion that the liquid is composed largely of edge-sharing dimers.

Ionic Liquids from Car−Parrinello Simulations, Part I: Liquid AlCl3

The properties of isolated AlCl3 clusters and the bulk system are investigated by means of static and dynamic electronic structure methods. We find important structural motifs with the edge connectivity dominant in a dimer and the corner connectivity dominant in a trimer. Furthermore, the trimer cluster exhibits an interesting ring structure with large cooperative effects relative to the dimer. Comparing the found structural motifs in isolated molecule calculations with the structure of the liquid allows us to determine the dominance of edge connectivity in the liquid. The size of the clusters present in the liquid indicates indeed that the dimer is the most abundant species, but there are also trimers, tetramers, and pentamers present. From the local dipole analysis both for the isolated clusters as well as for the liquid, further proof for the edge connectivity is given. However, all results point to the fact that there is also some small percentage of corner connectivity present that might be attributed to the most stable corner-connected cluster, namely the trimer. Importantly, we find that energetic considerations of isolated (static) clusters only do not represent the findings in liquid phase. Instead, a quantum cluster equilibrium approach or simulations are needed.

Raman spectra of ionic liquids: Interpretation via computer simulation

Theoretical Raman spectra of the complex-forming ionic liquids LaCl3 and ScCl3, derived from molecular dynamics computer simulations, are presented. These simulations, which use polarizable ion interaction models, have previously been shown to predict structural properties in excellent agreement with diffraction experiments. The dependence of the polarizability of the melt on the ionic positions, which determines the Raman spectrum through the time dependence of the polarizability correlation function, is modeled on the basis of ab initio electronic structure calculations carried out on alkali chlorides. New simulation techniques are introduced in order to allow the spectrum to be calculated with acceptable statistics. The calculated spectra are in semiquantitative agreement with experimental data. The distinctive bands which appear in the spectra of such complex melts are linked to the vibrations of the transient coordination complexes which form in these melts and new interpretations for the origin of several well-known features are proposed. The simulations thus enable a link between the structure of a melt as perceived through Raman spectroscopy and through diffraction experiments to be made.

Gauge symmetry, chirality and parity effects in four-particle systems: Coulomb's law as a universal function for diatomic molecules

Following recent work in search for a universal function (Van Hooydonk, Eur. J. Inorg. Chem., (1999), 1617), we test four symmetric +/- a(n)Rn potentials for reproducing molecular potential energy curves (PECs). Classical gauge symmetry for 1/R-potentials results in generic left right asymmetric PECs. A pair of symmetric perturbed Coulomb potentials is quantitatively in accordance with observed PECs. For a bond, a four-particle system, charge inversion (a parity effect, atom chirality) is the key to explain this shape generically. A parity adapted Hamiltonian reduces from ten to two terms and to a soluble Bohr-like formula, a Kratzer (1 - Re/R)2 potential. The result is similar to the combined action of spin and wave function symmetry upon the Hamiltonian in Heitler-London theory. Analytical perturbed Coulomb functions varying with (1 - Re/R) scale attractive and repulsive branches of PECs for 13 bonds H2, HF, LiH, KH, AuH, Li2, LiF, KLi, NaCs, Rb2, RbCs, Cs2 and I2 in a single straight line. The 400 turning points for 13 bonds are reproduced with a deviation of 0.007 A at both branches. For 230 points at the repulsive side, the deviation is 0.003 A. The perturbed electrostatic Coulomb law is a universal molecular function. Ab initio zero molecular parameter functions give PECs of acceptable quality, just using atomic ionisation energies. The function can be used as a model potential for inverting levels and gives a first principle's comparison of short- and long-range interactions, important for the study of cold atoms. Wave-packet dynamics, femto-chemistry applied to the crossing of covalent and ionic curves, can provide evidence for this theory. We anticipate this scale/shape invariant scheme applies to smaller scales in nuclear and high-energy particle physics. For larger gravitational scales (Newton 1/R potentials), problems with super-unification are discussed. Reactions between hydrogen and antihydrogen, feasible in the near future, will probably produce normal H2.