Search

Multi-criteria optimization is used for developing molecular models for ethylene carbonate (EC) and propylene carbonate (PC), organic solvents commonly used in Li-ion batteries. The molecular geometry and partial charges of the solvents are obtained from quantum mechanical calculations. Using a novel optimization strategy that combines systematic variations of the Lennard-Jones parameters with a reduced units approach, the models are fitted to experimental data on the liquid density, vapor pressure, relative permittivity, and self-diffusion coefficient. Since no experimental data for the self-diffusion coefficient of pure EC were available in the literature, they are measured in this work using a gradient-based nuclear magnetic resonance technique. For all pure component properties, excellent agreement between experiment and simulation is obtained. Moreover, the predictive capabilities of the new solvent models are assessed by comparison to experimental data for the liquid density and relative permittivity of mixtures of EC and PC. In addition, molecular models for the anions P F 6 − , B F 4 − , and C l O 4 − in solutions of their lithium electrolytes in PC are developed using experimental data on the solution densities. Finally, the self-diffusion coefficients of LiPF6 in PC and in aqueous solution are predicted and compared, showing that diffusion is much slower in the organic solution due to the formation of larger solvent shells around the ions. Furthermore, an analysis of the radial distribution functions in these solutions suggests that the ions have much less impact on the structure of the solvent PC than on water. [ABSTRACT FROM AUTHOR]

A set of molecular models for the alkali nitrates (LiNO3, NaNO3, KNO3, RbNO3, and CsNO3) in aqueous solutions is presented and used for predicting the thermophysical properties of these solutions with molecular dynamics simulations. The set of models is obtained from a combination of a model for the nitrate anion from the literature with a set of models for the alkali cations developed in previous works of our group. The water model is SPC/E and the Lorentz–Berthelot combining rules are used for describing the unlike interactions. This combination is shown to yield fair predictions of thermophysical and structural properties of the studied aqueous solutions, namely the density, the water activity and the mean ionic activity coefficient, the self-diffusion coefficients of the ions, and radial distribution functions, which were studied at 298 K and 1 bar; except for the density of the solutions of all five nitrates and the activity properties of solutions of NaNO3, which were also studied at 333 K. For calculating the water the activity and the mean ionic activity coefficient, the OPAS (osmotic pressure for the activity of selvents) method was applied. The new models extend an ion model family for the alkali halides developed in previous works of our group in a consistent way. [ABSTRACT FROM AUTHOR]

The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation. [ABSTRACT FROM AUTHOR]

A method for determining activity coefficients by molecular dynamics simulations is presented. It is an extension of the OPAS (osmotic pressure for the activity of the solvent) method in previous work for studying the solvent activity in electrolyte solutions. That method is extended here to study activities of all components in mixtures of molecular species. As an example, activity coefficients in liquid mixtures of water and methanol are calculated for 298.15 K and 323.15 K at 1 bar using molecular models from the literature. These dense and strongly interacting mixtures pose a significant challenge to existing methods for determining activity coefficients by molecular simulation. It is shown that the new method yields accurate results for the activity coefficients which are in agreement with results obtained with a thermodynamic integration technique. As the partial molar volumes are needed in the proposed method, the molar excess volume of the system water + methanol is also investigated. [ABSTRACT FROM AUTHOR]

A method for determining the activity of the solvent in electrolyte solutions by molecular dynamics simulations is presented. The electrolyte solution is simulated in contact with the pure solvent. Between the two phases, there is a virtual membrane, which is permeable only for the solvent. In the simulation, this is realized by an external field which acts only on the solutes and confines them to a part of the simulation volume. The osmotic pressure, i.e., the pressure difference between both phases, is obtained with high accuracy from the force on the membrane, so that reliable data on the solvent activity can be determined. The acronym of the new method is therefore OPAS (osmotic pressure for activity of solvents). The OPAS method is verified using tests of varying complexity. This includes a comparison of results from the OPAS method for aqueous NaCl solutions to results from the literature which were obtained with other molecular simulation methods. Favorable agreement is observed not only for the solvent activity but also for the activity coefficient of NaCl, which is obtained by application of the Gibbs-Duhem equation. [ABSTRACT FROM AUTHOR]