Molecular dynamics simulations of concentrated aqueous electrolyte solutions

Toward a molecular understanding of the conductivity of lithium-ion conducting polyanion polymer electrolytes by molecular dynamics simulation

With the improved lithium-ion transference number near unity, the low conductivity of single lithium-ion conducting solid polymer electrolytes (SLIC-SPEs) still hinders their application in high-rate batteries. Though some empirical conclusions on the conducting mechanism of SLIC-SPEs have been obtained, a more comprehensive study on the quantitative relationship between the molecular structure factors and ionic conduction performance is expected. In this study, a model structure that contains adjustable main chain and anion groups in the polyethylene oxide (PEO) matrix was used to clarify the influence of molecular structural factors on ionic conductivity and electrochemical stability of SLIC-SPEs. The anionic group was further disassembled into the intermediate group and end group while the main chain structure was distinguished into different degrees of polymerization and various lengths of the spacers between anions. Therefore, a well-defined molecular structure was employed to describe its relationship with ionic conductivity. In addition, the dissociation degree of salts and mobility of ions changing with the molecular structure were also discussed to explore the fundamental causes of conductivity. It can be concluded that the anion group affects the conductivity mainly via the dissociation degree, while the main chain structure impacts the conductivity by both dissociation degree and mobility.

Surface charge-dependent hydrodynamic properties of an electroosmotic slip flow

The electroosmosis effects at the interface of an aqueous NaCl solution and a charged silicon surface are studied using a molecular dynamics (MD) method. Considering a plug-like electroosmotic flow, we identified a thin interfacial layer in the immediate vicinity of the charged surface, where the flow velocity experiences almost linear spatial variations. The thickness of this interfacial layer is found to be linearly dependent on the surface charge density, with a negative slope which slightly depends on the surface hydrophobicity while being independent of the salt concentration, electric field strength, and orientation of the surface lattice. It is also found that upon increasing the surface charge density, the effective slip length first increases up to a maximum amount and then follows an almost linear reduction. We found that increasing the salt concentration drastically reduces the surface charge at which the effective slip length reaches its maximum amount. For highly concentrated solutions, therefore, the effective slip length could be assumed to change linearly in the whole range of the surface charge density, with a slope which is proportional to the square root of the electric field strength divided by the depth of the potential well assigned to the surface atoms εwall. Also, in a wide range of the surface charge density, the slip velocity is found to be a constant fraction of the electroosmotic velocity, which could be measured experimentally. Finally, by comparing the electroosmotic velocities calculated from the Stokes equation (considering both the slip and no-slip boundary conditions) with our MD results, we found that the no-slip boundary condition, which is normally used in analytical calculations, leads to a very inaccurate result for the studied system.

Role of Electroosmosis in the Permeation of Neutral Molecules: CymA and Cyclodextrin as an Example

To quantify the flow of small uncharged molecules into and across nanopores, one often uses ion currents. The respective ion-current fluctuations caused by the presence of the analyte make it possible to draw some conclusions about the direction and magnitude of the analyte flow. However, often this flow appears to be asymmetric with respect to the applied voltage. As a possible reason for this asymmetry, we identified the electroosmotic flow (EOF), which is the water transport associated with ions driven by the external transmembrane voltage. As an example, we quantify the contribution of the EOF through a nanopore by investigating the permeation of α-cyclodextrin through CymA, a cyclodextrin-specific channel from Klebsiella oxytoca. To understand the results from electrophysiology on a molecular level, all-atom molecular dynamics simulations are used to detail the effect of the EOF on substrate entry to and exit from a CymA channel in which the N-terminus has been deleted. The combined experimental and computational results strongly suggest that one needs to account for the significant contribution of the EOF when analyzing the penetration of cyclodextrins through the CymA pore. This example study at the same time points to the more general finding that the EOF needs to be considered in translocation studies of neutral molecules and, at least in many cases, should be able to help in discriminating between translocation and binding events.

Multiscale modeling of lithium ion batteries: thermal aspects

The thermal behavior of lithium ion batteries has a huge impact on their lifetime and the initiation of degradation processes. The development of hot spots or large local overpotentials leading, e.g., to lithium metal deposition depends on material properties as well as on the nano- und microstructure of the electrodes. In recent years a theoretical structure emerges, which opens the possibility to establish a systematic modeling strategy from atomistic to continuum scale to capture and couple the relevant phenomena on each scale. We outline the building blocks for such a systematic approach and discuss in detail a rigorous approach for the continuum scale based on rational thermodynamics and homogenization theories. Our focus is on the development of a systematic thermodynamically consistent theory for thermal phenomena in batteries at the microstructure scale and at the cell scale. We discuss the importance of carefully defining the continuum fields for being able to compare seemingly different phenomenological theories and for obtaining rules to determine unknown parameters of the theory by experiments or lower-scale theories. The resulting continuum models for the microscopic and the cell scale are numerically solved in full 3D resolution. The complex very localized distributions of heat sources in a microstructure of a battery and the problems of mapping these localized sources on an averaged porous electrode model are discussed by comparing the detailed 3D microstructure-resolved simulations of the heat distribution with the result of the upscaled porous electrode model. It is shown, that not all heat sources that exist on the microstructure scale are represented in the averaged theory due to subtle cancellation effects of interface and bulk heat sources. Nevertheless, we find that in special cases the averaged thermal behavior can be captured very well by porous electrode theory.

First-passage-time analysis of atomic-resolution simulations of the ionic transport in a bacterial porin

We have studied the dynamics of chloride and potassium ions in the interior of the Outer membrane porin F (OmpF) under the influence of an external electric field. From the results of extensive all-atom molecular dynamics (MD) simulations of the system, we computed several first-passage-time (FPT) quantities to characterize the dynamics of the ions in the interior of the channel. Such FPT quantities obtained from MD simulations demonstrate that it is not possible to describe the dynamics of chloride and potassium ions inside the whole channel with a single constant diffusion coefficient. However, we showed that a valid, statistically rigorous description in terms of a constant diffusion coefficient D and an effective deterministic force F(eff) can be obtained after appropriate subdivision of the channel in different regions suggested by the x-ray structure. These results have important implications for popular simplified descriptions of channels based on the one-dimensional Poisson-Nernst-Planck equations. Also, the effect of entropic barriers on the diffusion of the ions is identified and briefly discussed.