Equilibration strategies in molecular dynamics simulation of polyethers

A comparison of the NVT and NVE ensemble simulations of tetraglyme showed that, in terms of energy, temperature and most of the structural features the results were very similar. However, major differences were observed in dynamic properties, ie in the mean square displacement and in the O--O distances. A fast equilibration method suitable for amorphous polymer systems is also detailed. This was accomplished by the reassignment of the velocity distribution to the particles, after the interruption of an NVE simulation.

Lithium-doped plastic crystal electrolytes exhibiting fast ion conduction for secondary batteries

Rechargeable lithium batteries have long been considered an attractive alternative power source for a wide variety of applications. Safety and stability1 concerns associated with solvent-based electrolytes has necessitated the use of lithium intercalation materials (rather than lithium metal) as anodes, which decreases the energy storage capacity per unit mass. The use of solid lithium ion conductors - based on glasses, ceramics or polymers - as the electrolyte would potentially improve the stability of a lithium metal anode while alleviating the safety concerns. Glasses and ceramics conduct via a fast ion mechanism, in which the lithium ions move within an essentially static framework. In contrast, the motion of ions in polymer systems is similar to that in solvent-based electrolytes - motion is mediated by the dynamics of the host polymer, thereby restricting the conductivity to relatively low values. Moreover, in the polymer systems, the motion of the lithium ions provides only a small fraction of the overall conductivity2, which results in severe concentration gradients during cell operation, causing premature failure3. Here we describe a class of materials, prepared by doping lithium ions into a plastic crystalline matrix, that exhibit fast lithium ion motion due to rotational disorder and the existence of vacancies in the lattice. The combination of possible structural variations of the plastic crystal matrix and conductivities as high as 2 3 1024 S cm21 at 60 8C make these materials very attractive for secondary battery applications.

Solid-state NMR as a probe of anion binding: molecular dynamics and associations in a [5]polynorbornane bisurea host complexed with terephthalate

A range of solid-state NMR techniques is used to characterise a molecular host:guest complex consisting of a [5]polynorbornane bisurea host binding a terephthalate dianion guest. Detailed information is obtained on the molecular dynamics and associations from the point of view of both the host and guest molecules. The formation of the complex in the solid state is confirmed using (1)H 2D exchange NMR, and the 180° flipping of the (2)H-labelled terephthalate guest and its eventual expulsion from the complex at elevated temperatures are quantified using variable-temperature (2)H spin-echo experiments. Two-dimensional (1)H-(13)C HETCOR spectra obtained under fast magic angle spinning conditions (60 kHz) show a high resolution despite the poor crystallinity of the solid complex, and clearly reveal changes in the rigidity of the host molecule when complexed. Short-range intra- and intermolecular (1)H-(1)H proximities are also detected using 2D SQ-DQ correlation methods, providing insight into the molecular packing in the solid phase.

CISR-ODE, A C plus plus framework with ODE solver for code based system dynamics simulation

Ordinary differential equations are used for modelling a wide range of dynamic systems. Even though there are many graphical software applications for this purpose, a fully customised solution for all problems is code-level programming of the model and solver. In this project, a free and open source C++ framework is designed to facilitate modelling in native code environment and fulfill the common simulation needs of control and many other engineering and science applications. The solvers of this project are obtained from ODEINT and specialised for Armadillo matrix library to provide an easy syntax and a fast execution. The solver code is minimised and its modification for users have become easier. There are several features added to the solvers such as controlling maximum step size, informing the solver about sudden input change and forcing custom times into the results and calling a custom method at these points. The comfort of the model designer, code readability, extendibility and model isolation have been considered in the structure of this framework. The application manages the output results, exporting and plotting them. Modifying the model has become more practical and a portion of corresponding codes are updated automatically. A set of libraries is provided for generation of output figures, matrix hashing, control system functions, profiling, etc. In this paper, an example of using this framework for a classical washout filter model is explained.

Polymer architecture effect on sodium ion transport in PSTFSI-based ionomers: a molecular dynamics study

In this study, we investigated the effect of polymer architecture on the ion dynamics and local structure to understand the factors that might lead to the design of highly conductive and mechanically robust polyelectrolytes. Molecular dynamic simulations were undertaken on the sodium poly[(4-styrenesulfonyl) (trifluoromethanesulfonyl) imide] P(STFSINa) homopolymer and its copolymers with either ether or styrene spacer groups to investigate the spacer length and polarity dependence of Na-ion transport. Using a scaled charge model, we observed a continuous ion aggregate network in the homopolymer, which facilitates the fast ion dynamics despite the rigid polymer matrix. The longest spacer groups disrupt this percolating ionic network differently, with the ether group being more disruptive than the styrene group, and leading to more discrete ionic aggregates. The copolymer with the ether spacer was also found to result in an alternative Na-ion diffusion mechanism.

Molecular dynamics study of a dual-cation ionomer electrolyte

The poly(N1222)xLi1-x[AMPS] ionomer system with dual cations has previously shown decoupled Li ion dynamics from polymer segmental motions, characterized by the glass transition temperature, which can result in a conductive electrolyte material whilst retaining an appropriate modulus (ie. stiffness) so that it can suppress dendrite formation, thereby improving safety when used in lithium metal batteries. To understand this ion dynamics behavior, molecular dynamics techniques have been used in this work to simulate structure and dynamics in these materials. These simulations confirm that the Li ion transport is decoupled from the polymer particularly at intermediate N1222+ concentrations. At 50 mol% N1222+ concentration the polymer backbone is more rigid than for higher N1222+ concentrations, but with increasing temperature Li ion transport is more significant than polymer or quaternary ammonium cation motions. Here we suggest an ion hopping mechanism for Li+, arising from structural rearrangement of ionic clusters that could explain its decoupled behavior. Higher temperatures favor an aggregated ionic structure as well as enhancing these hopping motions. The simulations discussed here provide an atomic-level understanding of ion dynamics that could contribute to designing an improved ionomer with fast ion transport and mechanical robustness.