Predicting the conditions in which a compound adopts a metastable structure when it crystallizes out of solution is an unsolved and fundamental problem in materials synthesis, and one which, if understood and harnessed, could enable the rational design of synthesis pathways toward or away from metastable structures. Crystallization of metastable phases is particularly accessible via low-temperature solution-based routes, such as chimie douce and hydrothermal synthesis, but although the chemistry of the solution plays a crucial role in governing which polymorph forms, how it does so is poorly understood. Here, we demonstrate an ab initio technique to quantify thermodynamic parameters of surfaces and bulks in equilibrium with an aqueous environment, enabling the calculation of nucleation barriers of competing polymorphs as a function of solution chemistry, thereby predicting the solution conditions governing polymorph selection. We apply this approach to resolve the long-standing “calcite–aragonite problem”––the observation that calcium carbonate precipitates as the metastable aragonite polymorph in marine environments, rather than the stable phase calcite––which is of tremendous relevance to biomineralization, carbon sequestration, paleogeochemistry, and the vulnerability of marine life to ocean acidification. We identify a direct relationship between the calcite surface energy and solution Mg–Ca ion concentrations, showing that the calcite nucleation barrier surpasses that of metastable aragonite in solutions with Mg:Ca ratios consistent with modern seawater, allowing aragonite to dominate the kinetics of nucleation. Our ability to quantify how solution parameters distinguish between polymorphs marks an important step toward the ab initio prediction of materials synthesis pathways in solution.

Broader context: Electrical energy storage is a key technology for a clean energy economy, but currently requires significant improvement in energy density beyond the capabilities of traditional Li-ion batteries. Mg-ion batteries offer an exciting alternative in terms of the amount of energy that can be delivered, safety, manufacturing and disposal costs, with limited environmental impact. The electrochemical functions of the Mg-ion battery ultimately depend on the choice of the electrolyte, which is limited by the peculiar chemistry of Mg. To date, very few electrolytes can reversibly plate and strip Mg. The Magnesium Aluminum Chloride Complex (MACC) electrolyte can reversibly plate and strip Mg with significantly higher voltage (3.1 V) as compared to other electrolytes, but there is a pressing need to address critical questions about the structural evolution of this electrolyte during electrochemical cycling. Abstract: Non-aqueous Mg-ion batteries offer a promising way to overcome safety, costs, and energy density limitations of state-of-the-art Li-ion battery technology. We present a rigorous analysis of the Magnesium Aluminum Chloro Complex (MACC) in tetrahydrofuran (THF), one of the few electrolytes that can reversibly plate and strip Mg. We use ab initio calculations and classical molecular dynamics simulations to interrogate the MACC electrolyte composition with the goal of addressing two urgent questions that have puzzled battery researchers: i) the functional species of the electrolyte, and ii) the complex equilibria regulating the MACC speciation after prolonged electrochemical cycling, a process termed as conditioning, and after prolonged inactivity, a process called aging. A general computational strategy to untangle the complex structure of electrolytes, ionic liquids and other liquid media is presented. The analysis of formation energies and grand-potential phase diagrams of Mg-Al-Cl-THF suggests that the MACC electrolyte bears a simple chemical structure with few simple constituents, namely the electro-active species MgCl + and AlCl − 4 in equilibrium with MgCl2 and AlCl3. Knowledge of the stable species of the MACC electrolyte allows us to determine the most important equilibria occurring during electrochemical cycling. We observe that Al deposition is always preferred to Mg deposition, explaining why freshly synthesized MACC cannot operate and needs to undergo preparatory conditioning. Similarly, we suggest that aluminum displacement and depletion from the solution upon electrolyte resting (along with continuous MgCl2 regeneration) represents one of the causes of electrolyte aging. Finally, we compute the NMR shifts from shielding tensors of selected molecules and ions providing fingerprints to guide future experimental investigations. arXiv:1511.02504v1 [cond-mat.mtrl-sci] 8 Nov 2015

Multi-valent (MV) battery architectures based on pairing a Mg metal anode with a high-voltage (∼ 3 V) intercalation cathode offer a realistic design pathway toward significantly surpassing the energy storage performance of traditional Li-ion based batteries, but there are currently only few electrolyte systems that support reversible Mg deposition. Using both static first-principles calculations and ab initio molecular dynamics, we perform a comprehensive adsorption study of several salt and solvent species at the interface of Mg metal with an electrolyte of Mg 2+ and Cl − dissolved in liquid tetrahydrofuran (THF). Our findings not only provide a picture of the stable species at the interface, but also explain how this system can support reversible Mg deposition and as such we provide insights in how to design other electrolytes for Mg plating and stripping. The active depositing species are identified to be (MgCl) + monomers coordinated by THF, which exhibit preferential adsorption on Mg compared to possible passivating species (such as THF solvent or neutral MgCl2 complexes). Upon deposition, the energy to desolvate these adsorbed complexes and facilitate charge-transfer is shown to be small (∼ 61 -46.2 kJ mol −1 to remove 3 THF from the strongest adsorbing complex), and the stable orientations of the adsorbed but desolvated (MgCl) + complexes appear favorable for charge-transfer. Finally, observations of Mg-Cl dissociation at the Mg surface at very low THF coordinations (0 and 1) suggest that deleterious Cl incorporation in the anode may occur upon plating. In the stripping process, this is beneficial by further facilitating the Mg removal reaction.

The melting point, enthalpy of fusion, and thermodynamic stability of two crystal polymorphs of the ionic liquid 1-n-butyl-3-methylimidazolium chloride are calculated using a thermodynamic integration-based atomistic simulation method. The computed melting point of the orthorhombic phase ranges from 365 to 369 K, depending on the classical force field used. This compares reasonably well with the experimental values, which range from 337 to 339 K. The computed enthalpy of fusion ranges from 19 to 29 kJ/mol, compared to the experimental values of 18.5-21.5 kJ/mol. Only one of the two force fields evaluated in this work yielded a stable monoclinic phase, despite the fact that both give accurate liquid state densities. The computed melting point of the monoclinic polymorph was found to be 373 K, which is somewhat higher than the experimental range of 318-340 K. The computed enthalpy of fusion was 23 kJ/mol, which is also higher than the experimental value of 9.3-14.5 kJ/mol. The simulations predict that the monoclinic form is more stable than the orthorhombic form at low temperature, in agreement with one set of experiments but in conflict with another. The difference in free energy between the two polymorphs is very small, due to the fact that a single trans-gauche conformational difference in an alkyl sidechain distinguishes the two structures. As a result, it is very difficult to construct simple classical force fields that are accurate enough to definitively predict which polymorph is most stable. A liquid phase analysis of the probability distribution of the dihedral angles in the alkyl chain indicates that less than half of the dihedral angles are in the gauche-trans configuration that is adopted in the orthorhombic crystal. The low melting point and glass forming tendency of this ionic liquid is likely due to the energy barrier for conversion of the remaining dihedral angles into the gauche-trans state. The simulation procedure used to perform the melting point calculations is an extension of the so-called pseudo-supercritical path sampling procedure. This study demonstrates that the method can be effectively applied to quite complex systems such as ionic liquids and that the appropriate choice of tethering potentials for a key step in the thermodynamic path can enable first order phase transitions to be avoided.

We present an adaptable method to compute the solubility limit of solids by molecular simulation, which avoids the difficulty of reference state calculations. In this way, the method is highly adaptable to molecules of complex topology. Results are shown for solubility calculations of sodium chloride in water and light alcohols at atmospheric conditions. The pseudosupercritical path integration method is used to calculate the free energy of the solid and gives results that are in good agreement with previous studies that reference the Einstein crystal. For the solution phase calculations, the self-adaptive Wang-Landau transition-matrix Monte Carlo method is used within the context of an expanded isothermal-isobaric ensemble. The method shows rapid convergence properties and the uncertainty in the calculated chemical potential was 1% or less for all cases. The present study underpredicts the solubility limit of sodium chloride in water, suggesting a shortcoming of the molecular models. Importantly, the proper trend for the chemical potential in various solvents was captured, suggesting that relative solubilities can be computed by the method.