Molecular dynamics study of response of liquid N,N-dimethylformamide to externally applied electric field using a polarizable force field

The behavior of Liquid N,N-dimethylformamide subjected to a wide range of externally applied electric fields (from 0.001 V/nm to 1 V/nm) has been investigated through molecular dynamics simulation. To approach the objective the AMOEBA polarizable force field was extended to include the interaction of the external electric field with atomic partial charges and the contribution to the atomic polarization. The simulation results were evaluated with quantum mechanical calculations. The results from the present force field for the liquid at normal conditions were compared with the experimental and molecular dynamics results with non-polarizable and other polarizable force fields. The uniform external electric fields of higher than 0.01 V/nm have a significant effect on the structure of the liquid, which exhibits a variation in numerous properties, including molecular polarization, local cluster structure, rotation, alignment, energetics, and bulk thermodynamic and structural properties.

Ion clustering in molecular dynamics simulations of sodium iodide solutions

Model systems of sodium iodide dissolved in dimethyl ether or 1,2-dimethoxyethane (glyme) were studied in order to investigate the structural and dynamic properties of ionic solutions in small and polymeric ethers. Full molecular dynamics simulations were performed at a range of different salt concentrations. An algorithm was designed which assigns ions to clusters and then calculates all the terms which contribute to ionic conductivity. In dilute solutions, free ions are the most common ionic species, followed by ion pairs. As the concentration increases, pairs become the most common species, with significant concentrations of clusters with 3 through 6 ions. Changing the solvent from dimethyl ether to glyme significantly decreases the ion clustering due to the chelate effect in which the two oxygens on a solvent stabilize an associated cation. The conductivity in stable systems is shown to be primarily the result of the movement of free ions and the relative movement of ions within neutral pairs. The Nernst-Einstein relation, commonly used in the discussion of polymer electrolytes, is shown to be inadequate to quantitatively describe conductivity in the model systems.

Insight into local structure and molecular dynamics in organic solid-state ionic conductors

Elucidating the rate and geometry of molecular dynamics is particularly important for unravelling ion-conduction mechanisms in electrochemical materials. The local molecular motions in the plastic crystal 1-ethyl-1-methylpyrrolidinium tetrafluoroborate ([C2 mpyr][BF4 ]) are studied by a combination of quantum chemical calculations and advanced solid-state nuclear magnetic resonance spectroscopy. For the first time, a restricted puckering motion with a small fluctuation angle of 25° in the pyrrolidinium ring has been observed, even in the low-temperature phase (-45 °C). This local molecular motion is deemed to be particularly important for the material to maintain its plasticity, and hence, its ion mobility at low temperatures.

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 approach to the structural characterization and transport properties of poly(acrylonitrile)/N,N-dimethylformamide solutions

Poly(acrylonitrile) (PAN) in N,N-dimethylformamide (DMF) is a popular solution for producing large variety of polymer products. To precisely describe the behaviours of PAN and DMF in the synthesis processes, it is significant to call for more details about the structure, some thermodynamic and dynamical properties of PAN-DMF solutions. A PAN-DMF solution was simulated via molecular dynamics with an all-atom OPLS type potential in both the NPT and NVT ensembles. The simulation results were evaluated with quantum mechanical calculations (MP2/6-311 ++G(d,p) and counterpoise procedure) and were compared with available experimental results. The liquid structure was illustrated with pair correlation functions and transport and dynamics properties were calculated with the mean-square displacements MSD and the velocity autocorrelation functions. The strong H-bonds of C≡N « H-C=O, CH » O=C-H and CH2 O=C-H, with distances of 2.55 Å, 2.55 Å and 2.65 Å, respectively, were found. The largest interaction energy of - 7.157 kcal/mol between DMF molecules and PAN molecules was found at 4.9 Å center-of-mass distance. A potential profile of intermolecular interaction of DMF with PAN along the interaction distance was presented, clearly showing an increase of DMF vaporisation heat when it getting close to PAN molecules. This provided very useful information to analyse the vaporisation behaviours of DMF at the microscopic level, which is essential to comprehensively understand molecular rearrangements towards the design of synthetic processes. The impact of the presence of the PAN on the DMF solution properties were also benchmarked with pure DMF solution.

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.

Molecular dynamics simulations of highly concentrated salt solutions : structural and transport effects in polymer electrolytes

Structural, thermodynamic and transport properties have been calculated in concentrated non-aqueous NaI solutions using molecular dynamics simulations. Although the solvent has been represented by a simplistic Stockmayer fluid (spherical particles with point dipoles), the general trends observed are still a useful indication of the behavior of real non-aqueous electrolyte systems. Results indicate that in low dielectric media, significant ion pairing and clustering occurs. Contact ion pairs become more prominent at higher temperatures, independent of the dielectric strength of the solvent. Thermodynamic analysis shows that this temperature behavior is predominantly entropically driven. Calculation of ionic diffusivities and conductivities in the NaI/ether system confirms the clustered nature of the salt, with the conductivities significantly lower than those predicted from the Nernst-Einstein relation. In systems where the solvent-ion interactions increase relative to ion-ion interactions (lower charge or higher solvent dipole moment), less clustering is observed and the transport properties indicate independent motion of the ions, with higher calculated conductivities. The solvent in this system is the most mobile species, in comparison with the polymer electrolytes where the solvent is practically immobile.

Highly concentrated salt solutions : molecular dynamics simulations of structure and transport

Molecular dynamics (MD) simulations in NaI solutions, where the solvent has been represented by the Stockmayer fluid, were performed as a function of temperature, salt concentration, and solvent dipole strength. At higher temperatures contact ion pairs become more prevalent, regardless of solvent strength. An examination of the temperature dependence of the potential of mean force demonstrates the entropic nature of this effect. The transport properties calculated in the simulations are dependent on the balance between solvent dielectric constant and ion charge. In systems with a large solvent dipole moment, the ions appear to be independently mobile, and deviations from Nernst–Einstein behavior are small. In systems of smaller solvent dipole moment or greater ion charge, the ions form clusters, and large deviations from Nernst–Einstein behavior are observed.