Electron transfer reaction of TEMPO-based organic radical batteries in different solvent environments: comparing quantum and classical approaches.

In this study, we delve into the complex electron transfer reactions associated with the redox-active (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO), a common component in organic radical batteries (ORBs). Our approach estimates quantum electron-transfer (ET) energies using Density Functional Theory (DFT) calculations by sampling from structures simulated classically. This work presents a comparative study of reorganization energies in ET reactions across different solvents. Furthermore, we investigate how changes in the electrolyte environment can modify the reorganization energy and, consequently, impact ET dynamics. We also explore the relationship between classical and quantum vertical energies using linear regression models. Importantly, this comparison between quantum and classical vertical energies underscores the role of quantum effects, like charge delocalization, in offering added stabilization post-redox reactions. These effects are not adequately represented by the classical vertical energy distribution. Our study shows that, although we find a significant correlation between the vertical energies computed by DFT and the classical force field, the regression parameters depend on the solvent, highlighting that classical methods should be benchmarked by DFT before applying them to novel electrolyte materials.

Molecular-Cling-Effect of Fluoroethylene Carbonate Characterized via Ethoxy(pentafluoro)cyclotriphosphazene on SiOx/C Anode Materials - A New Perspective for Formerly Sub-Sufficient SEI Forming Additive Compounds.

Effective electrolyte compositions are of primary importance in raising the performance of lithium-ion batteries (LIBs). Recently, fluorinated cyclic phosphazenes in combination with fluoroethylene carbonate (FEC) have been introduced as promising electrolyte additives, which can decompose to form an effective dense, uniform, and thin protective layer on the surface of electrodes. Although the basic electrochemical aspects of cyclic fluorinated phosphazenes combined with FEC were introduced, it is still unclear how these two compounds interact constructively during operation. This study investigates the complementary effect of FEC and ethoxy(pentafluoro)cyclotriphosphazene (EtPFPN) in aprotic organic electrolyte in LiNi0.5 Co0.2 Mn0.3 O ∥ SiOx /C full cells. The formation mechanism of lithium ethyl methyl carbonate (LEMC)-EtPFPN interphasial intermediate products and the reaction mechanism of lithium alkoxide with EtPFPN are proposed and supported by Density Functional Theory calculations. A novel property of FEC is also discussed here, called molecular-cling-effect (MCE). To the best knowledge, the MCE has not been reported in the literature, although FEC belongs to one of the most investigated electrolyte additives. The beneficial MCE of FEC toward the sub-sufficient solid-electrolyte interphase forming additive compound EtPFPN is investigated via gas chromatography-mass spectrometry, gas chromatography high resolution-accurate mass spectrometry, in situ shell-isolated nanoparticle-enhanced Raman spectroscopy, and scanning electron microscopy.

Polypropylene carbonate-based electrolytes as model for a different approach towards improved ion transport properties for novel electrolytes.

Linear poly(alkylene carbonates) such as polyethylene carbonate (PEC) and polypropylene carbonate (PPC) have gained increasing interest due to their remarkable ion transport properties such as high Li+ transference numbers. The cause of these properties is not yet fully understood which makes it challenging to replicate them in other polymer electrolytes. Therefore, it is critical to understand the underlying mechanisms in polycarbonate electrolytes such as PPC. In this work we present insights from impedance spectroscopy, transference number measurements, PFG-NMR, IR and Raman spectroscopy as well as molecular dynamics simulations to address this issue. We find that in addition to plasticization, the lithium ion coordination by the carbonate groups of the polymer is weakened upon gelation, leading to a rapid exhange of the lithium ion solvation shell and consequently a strong increase of the conductivity. Moreover, we study the impact of the anions by employing different conducting salts. Interestingly, while the total conductivity decreases with increasing anion size, the reverse trend can be observed for the lithium ion transference numbers. Via our holistic approach, we demonstrate that this behavior can be attributed to differences in the collective ion dynamics.

Mechanistic understanding of the correlation between structure and dynamics of liquid carbonate electrolytes: impact of polarization.

Liquid electrolyte design and modelling is an essential part of the development of improved lithium ion batteries. For mixed organic carbonates (ethylene carbonate (EC) and ethyl-methyl carbonate (EMC) mixtures)-based electrolytes with LiPF6 as salt, we have compared a polarizable force field with the standard non-polarizable force field with and without charge rescaling to model the structural and dynamic properties. The result of our molecular dynamics simulations shows that both polarizable and non-polarizable force fields have similar structural factors, which are also in agreement with X-ray diffraction experimental results. In contrast, structural differences are observed for the lithium neighborhood, while the lithium-anion neighbourhood is much more pronounced for the polarizable force field. Comparison of EC/EMC coordination statistics with Fourier transformed infrared spectroscopy (FTIR) shows the best agreement for the polarizable force field. Also for transport quantities such as ionic conductivities, transference numbers, and viscosities, the agreement with the polarizable force field is by far better for a large range of salt concentrations and EC : EMC ratios. In contrast, for the non-polarizable variants, the dynamics are largely underestimated. The excellent performance of the polarizable force field is explored in different ways to pave the way to a realistic description of the structure-dynamics relationships for a wide range of salt and solvent compositions for this standard electrolyte. In particular, we can characterize the distinct correlation terms between like and unlike ions, relate them to structural properties, and explore to which degree the transport in this electrolyte is mass or charge limited.

Hydrodynamic interactions in ion transport-Theory and simulation.

We present a hydrodynamic theory describing pair diffusion in systems with periodic boundary conditions, thereby generalizing earlier work on self-diffusion [B. Dünweg and K. Kremer, J. Chem. Phys. 99, 6983-6997 (1993) and I.-C. Yeh and G. Hummer, J. Phys. Chem. B 108, 15873-15879 (2004)]. Its predictions are compared with Molecular Dynamics simulations for a liquid carbonate electrolyte and two ionic liquids, for which we characterize the correlated motion between distinct ions. Overall, we observe good agreement between theory and simulation data, highlighting that hydrodynamic interactions universally dictate ion correlations. However, when summing over all ion pairs in the system to obtain the cross-contributions to the total cationic or anionic conductivity, the hydrodynamic interactions between ions with like and unlike charges largely cancel. Consequently, significant conductivity contributions only arise from deviations from a hydrodynamic flow field of an ideal fluid, which is from the local electrolyte structure as well as the relaxation processes in the subdiffusive regime. In the case of ionic liquids, the momentum-conservation constraint additionally is vital, which we study by employing different ionic masses in the simulations. Our formalism will likely also be helpful to estimate finite-size effects of the conductivity or of Maxwell-Stefan diffusivities in simulations.

Controlling Li+ transport in ionic liquid electrolytes through salt content and anion asymmetry: a mechanistic understanding gained from molecular dynamics simulations.

In this work, we report the results from molecular dynamics simulations of lithium salt-ionic liquid electrolytes (ILEs) based either on the symmetric bis[(trifluoromethyl)sulfonyl]imide (TFSI-) anion or its asymmetric analogue 2,2,2-(trifluoromethyl)sulfonyl-N-cyanoamide (TFSAM-). Relating lithium's coordination environment to anion mean residence times and diffusion constants confirms the remarkable transport behaviour of the TFSAM--based ILEs that has been observed in recent experiments: for increased salt doping, the lithium ions must compete for the more attractive cyano over oxygen coordination and a fragmented landscape of solvation geometries emerges, in which lithium appears to be less strongly bound. We present a novel, yet statistically straightforward methodology to quantify the extent to which lithium and its solvation shell are dynamically coupled. By means of a Lithium Coupling Factor (LCF) we demonstrate that the shell anions do not constitute a stable lithium vehicle, which suggests for this electrolyte material the commonly termed "vehicular" lithium transport mechanism could be more aptly pictured as a correlated, flow-like motion of lithium and its neighbourhood. Our analysis elucidates two separate causes why lithium and shell dynamics progressively decouple with higher salt content: on the one hand, an increased sharing of anions between lithium limits the achievable LCF of individual lithium-anion pairs. On the other hand, weaker binding configurations naturally entail a lower dynamic stability of the lithium-anion complex, which is particularly relevant for the TFSAM--containing ILEs.

Connecting the quantum and classical mechanics simulation world: Applications of reactive step molecular dynamics simulations.

This article presents the application of the reactive step molecular dynamics simulation method [M. Biedermann, D. Diddens, and A. Heuer, J. Chem. Theory Comput. 17, 1074 (2021)] toward two different atomistic, chemically reactive systems. During reactive steps, transitions from reactant to product molecules are modeled according to physically correct transition probabilities based on quantum chemical information about the reactions such as molecular reaction rates via instant exchange of the employed force field and a subsequent, short relaxation of the structure. In the first application, we study the follow-up reactions of singly reduced ethylene carbonate (EC) radicals in EC solution, first, via extensive ab initio molecular dynamics simulations and, second, with the reactive step algorithm. A direct comparison of both simulation methods shows excellent agreement. Then, we employ the reactive step algorithm to simulate the enolate formation of 2-methylcyclopropanone with the base lithium diisopropylamine. Thereby, we can demonstrate that the reactive step algorithm is also capable of capturing effects from kinetic vs thermodynamic control of chemical reactions during simulation.

Cation-Assisted Lithium-Ion Transport for High-Performance PEO-based Ternary Solid Polymer Electrolytes.

N-alkyl-N-alkyl pyrrolidinium-based ionic liquids (ILs) are promising candidates as non-flammable plasticizers for lowering the operation temperature of poly(ethylene oxide) (PEO)-based solid polymer electrolytes (SPEs), but they present limitations in terms of lithium-ion transport, such as a much lower lithium transference number. Thus, a pyrrolidinium cation was prepared with an oligo(ethylene oxide) substituent with seven repeating units. We show, by a combination of experimental characterizations and simulations, that the cation's solvating properties allow faster lithium-ion transport than alkyl-substituted analogues when incorporated in SPEs. This proceeds not only by accelerating the conduction modes of PEO, but also by enabling new conduction modes linked to the solvation of lithium by a single IL cation. This, combined with favorable interfacial properties versus lithium metal, leads to significantly improved performance on lithium-metal polymer batteries.

3D structure of the electric double layer of ionic liquid-alcohol mixtures at the electrochemical interface.

Mixtures of the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate with amphiphilic cosolvents, such as methanol and ethanol, nanoconfined between graphene walls are studied by means of molecular dynamics simulations and the results are compared with those of the pure ionic liquid and its mixtures with water confined in the same conditions. We investigate the adsorption of cosolvent molecules at the graphene walls as well as their distribution across the system. The results show that, due to a higher affinity of the polar groups to be close to the anions in combination with the electrostatic and excluded volume interactions, there exists a high tendency of the OH groups to lie close to the anode, inducing small changes in the first cation layer. The orientation of cosolvent molecules is found to be closely related to the alignment of the molecular dipole moment. We also investigate the lateral ionic distribution in the layers close to the electrodes, which shows a structural transition from liquid-like lamellar ordering to solid-like hexagonal patterns as the size of the cosolvent molecules increases leading to smaller position fluctuations of the ions. The dependence of the specific patterns on the nature of the electrodes is also studied. This study strongly suggests that the ionic patterns formed in the first ionic layers next to the charged interfaces are universal since their existence does not crucially depend on the atomic composition of the interfacial material, but only on the net charge density of the considered ionic layer, which significantly changes the ionic mobility in this region.

ForConX: A forcefield conversion tool based on XML.

The force field conversion from one MD program to another one is exhausting and error-prone. Although single conversion tools from one MD program to another exist not every combination and both directions of conversion are available for the favorite MD programs Amber, Charmm, Dl-Poly, Gromacs, and Lammps. We present here a general tool for the force field conversion on the basis of an XML document. The force field is converted to and from this XML structure facilitating the implementation of new MD programs for the conversion. Furthermore, the XML structure is human readable and can be manipulated before continuing the conversion. We report, as testcases, the conversions of topologies for acetonitrile, dimethylformamide, and 1-ethyl-3-methylimidazolium trifluoromethanesulfonate comprising also Urey-Bradley and Ryckaert-Bellemans potentials. © 2017 Wiley Periodicals, Inc.

Local Chain Segregation and Entanglements in a Confined Polymer Melt.

The reptation mechanism, introduced by de Gennes and Edwards, where a polymer diffuses along a fluffy tube, defined by the constraints imposed by its surroundings, convincingly describes the relaxation of long polymers in concentrated solutions and melts. We propose that the scale for the tube diameter is set by local chain segregation, which we study analytically. We show that the concept of local segregation is especially operational for confined geometries, where segregation extends over mesoscopic domains, drastically reducing binary contacts, and provide an estimate of the entanglement length. Our predictions are quantitatively supported by extensive molecular dynamics simulations on systems consisting of long, entangled chains.

Aqueous ionic liquids and their influence on peptide conformations: denaturation and dehydration mechanisms.

Low concentrated aqueous ionic liquids (ILs) and their influence on protein structures have attracted a lot of interest over the last few years. This can be mostly attributed to the fact that aqueous ILs, depending on the ion species involved, can be used as protein protectants or protein denaturants. Atomistic molecular dynamics (MD) simulations are performed in order to study the influence of different aprotic ILs on the properties of a short hairpin peptide. Our results reveal distinct binding and denaturation effects for 1-ethyl-3-methylimidazolium (EMIM) in combination with different anions, namely, chloride (CL), tetrafluoroborate (BF4) and acetate (ACE). The simulation outcomes demonstrate that the studied ILs with larger anions reveal a more pronounced accumulation behavior of the individual ion species around the peptide, which is accomplished by a stronger dehydration effect. We can relate these findings to the implications of the Kirkwood-Buff theory, which provides a thermodynamic explanation for the denaturation strength in terms of the IL accumulation behavior. The results for the spatial distribution functions, the binding energies and the local/bulk partition coefficients are in good agreement with metadynamics simulations in order to determine the energetically most stable peptide conformations. The free energy landscapes indicate a decrease of the denaturation strength in the order EMIM/ACE, EMIM/BF4 and EMIM/CL, which coincides with a decreasing size of the anion species. An analysis of the potential binding energies reveals that this effect is mainly of enthalpic nature.

Counterintuitive trends of the wetting behavior of ionic liquid-based electrolytes on modified lithium electrodes.

The demand for high energy densities has brought rechargeable lithium metal batteries back into the research focus. Ionic liquids (ILs) are considered as suitable electrolyte components for these systems. In this work, the wetting behavior of 1-ethyl-3-methylimidazolium bis((trifluoromethyl)sulfonyl)imide ([C2MIm]TFSI), 1-butyl-3-methylimidazolium bis-((trifluoromethyl)sulfonyl)imide ([C4MIm]TFSI), 1-hexyl-3-methylimidazolium bis((trifluoromethyl)sulfonyl)imide ([C6MIm]TFSI), and N-butyl-N-methylpyrrolidinium bis((trifluoromethyl)sulfonyl)imide (Pyr14TFSI) on mechanically modified lithium electrodes, with and without lithium bis((trifluoromethyl)sulfonyl)imide (LiTFSI) conducting salt, is investigated and is compared to an organic carbonate-based electrolyte. Three different patterns were chosen for the lithium modification, enabling a surface area increase of 12%, 20%, and 56% for the modified lithium electrodes. Especially for pure ILs, the contact angle on lithium was significantly larger with higher surface areas of the lithium electrodes. Since the addition of LiTFSI remarkably decreased the contact angles of the ILs on the modified lithium surfaces, it could be shown that the effect of LiTFSI can be attributed to a decreased surface tension. This observation could be explained by an interruption of the ordering of ionic liquid cations and anions, which is supported by Raman spectroscopy and molecular dynamics (MD) simulations.

Chain end mobilities in polymer melts--a computational study.

The Rouse model can be regarded as the standard model to describe the dynamics of a short polymer chain under melt conditions. In this contribution, we explicitly check one of the fundamental assumptions of this model, namely, that of a uniform friction coefficient for all monomers, on the basis of MD simulation data of a poly(ethylene oxide) (PEO) melt. This question immediately arises from the fact that in a real polymer melt, the terminal monomers have on average more intermolecular neighbors than the central monomers, and one would expect that exactly these details affect the precise value of the friction coefficient. The mobilities are determined by our recently developed statistical method, which provides detailed insights into the local polymer dynamics. Moreover, it yields complementary information to that obtained from the mean square displacement (MSD) or the Rouse mode analysis. It turns out that the Rouse assumption of a uniform mobility is fulfilled to a good approximation for the PEO melt. However, a more detailed analysis reveals that the underlying microscopic dynamics are highly affected by different contributions from intra- and intermolecular excluded volume interactions, which cannot be taken into account by a modified friction coefficient. Minor deviations occur only for the terminal monomers on larger time scales, which can be attributed to the presence of two different escape mechanisms from their first coordination sphere. These effects remain elusive when studying the dynamics with the MSD only.

Effects of ionic liquids on cation dynamics in amorphous polyethylene oxide electrolytes.

We perform extensive molecular dynamics simulations of a poly(ethylene oxide)-based polymer electrolyte material containing lithium bis(trifluoromethanesulfonyl)imide salt for a wide temperature regime above and below the experimental crystallization temperature with and without N-methyl-N-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide ionic liquid (IL). The impact of the IL-concentration on the cation dynamics is studied. The increase of the cation mobility upon addition of IL is significant but temperature-independent. This can be related to distinct variations of the underlying transport properties as expressed within the previously introduced transport model of polymer electrolytes. Even for the largest IL concentration the transport model perfectly predicts the non-trivial time-dependence of the cationic mean square displacement for all temperatures. Finally, we compare our numerical and theoretical findings with the results of recent nuclear magnetic resonance experiments. In this way we can exclusively relate the strong experimentally observed dependence of the low-temperature Li-diffusivity on the IL concentration to the impact of IL on crystallization.

Machine Learning Isotropic g Values of Radical Polymers.

Methods for electronic structure computations, such as density functional theory (DFT), are routinely used for the calculation of spectroscopic parameters to establish and validate structure-parameter correlations. DFT calculations, however, are computationally expensive for large systems such as polymers. This work explores the machine learning (ML) of isotropic g values, giso, obtained from electron paramagnetic resonance (EPR) experiments of an organic radical polymer. An ML model based on regression trees is trained on DFT-calculated g values of poly(2,2,6,6-tetramethylpiperidinyloxy-4-yl methacrylate) (PTMA) polymer structures extracted from different time frames of a molecular dynamics trajectory. The DFT-derived g values, gisocalc, for different radical densities of PTMA, are compared against experimentally derived g values obtained from in operando EPR measurements of a PTMA-based organic radical battery. The ML-predicted giso values, gisopred, were compared with gisocalc to evaluate the performance of the model. Mean deviations of gisopred from gisocalc were found to be on the order of 0.0001. Furthermore, a performance evaluation on test structures from a separate MD trajectory indicated that the model is sensitive to the radical density and efficiently learns to predict giso values even for radical densities that were not part of the training data set. Since our trained model can reproduce the changes in giso along the MD trajectory and is sensitive to the extent of equilibration of the polymer structure, it is a promising alternative to computationally more expensive DFT methods, particularly for large systems that cannot be easily represented by a smaller model system.

Ternary Solid Polymer Electrolytes at the Electrochemical Interface: A Computational Study.

Polymer-based solid-like gel electrolytes have emerged as a promising alternative to improve battery performance. However, there is a scarcity of studies on the behavior of these media at the electrochemical interface. In this work, we report classical MD simulations of ternary polymer electrolytes composed of poly(ethylene oxide), a lithium salt [lithium bis(trifluoromethanesulfonyl)imide], and different ionic liquids [1-butyl-1-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide and 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide] confined between two charged and uncharged graphene-like surfaces. The molecular solvation of Li+ ions and their diffusion as well as the polymer conformational picture were characterized in terms of the radial distribution functions, coordination numbers, number density profiles, orientations, displacement variance, polymer radius of gyration, and polymer end-to-end distance. Our results show that the layering behavior of the ternary electrolyte in the interfacial region leads to a decrease of Li+ mobility in the direction perpendicular to the electrodes and high energy barriers that hinder lithium cations from coming into direct contact with the graphene-like surface. The nature of the ionic liquid and its concentration were found to influence the structural and dynamic properties at the electrode/electrolyte interface, the electrolyte with low amounts of the pyrrolidinium-based ionic liquid being that with the best performance since it favors the migration of Li+ cations toward the negative electrode when compared to the imidazolium-based one.

Li+ transport in poly(ethylene oxide) based electrolytes: neutron scattering, dielectric spectroscopy, and molecular dynamics simulations.

The dynamics of Li(+) transport in polyethylene oxide (PEO) and lithium bis(trifluoromethanesulfonyl)imde mixtures are investigated by combining neutron spin-echo (NSE) and dielectric spectroscopy with molecular dynamics (MD) simulations. The results are summarized in a relaxation time map covering wide ranges of temperature and time. The temperature dependence of the dc conductivity and the dielectric α relaxation time is found to be identical, indicating a strong coupling between both. The relaxation times obtained from the NSE measurements at 0.05 Å(-1)

Study of a High-Voltage NMC Interphase in the Presence of a Thiophene Additive Realized by Operando SHINERS.

Improving the electrochemical properties and cycle life of high-voltage cathodes in lithium-ion batteries requires a deep understanding of the structural properties and failure mechanisms at the cathode electrolyte interphase (CEI). We present a study implementing an advanced Raman spectroscopy technique to specifically address the compositional features of interphase during cell operation. Our operando technique, shell-isolated nanoparticle-enhanced Raman spectroscopy (SHINERS), provides a reliable platform to investigate the dynamics of the interphase structure and elucidate the compositional changes near the cathode surface. To improve the CEI properties, thiophene was introduced and investigated as an effective, high-voltage film-forming additive by largely diminishing the capacity fading triggered at high potentials in LiNi1/3Co1/3Mn1/3O2 cathodes. While the cells without thiophene show severe capacity fading, cells with an optimized concentration of thiophene exhibit a significant performance improvement. Operando SHINERS detects the presence of a stable CEI. The results suggest that the composition of the CEI is dominated by polythiophene and copolymerization products of ethylene carbonate with thiophene, which protects the electrolyte components from further decomposition. The formation mechanism of the polymeric film was modeled using quantum chemistry calculations, which shows good agreement with the experimental data.

Multimodal investigation of electronic transport in PTMA and its impact on organic radical battery performance.

Organic radical batteries (ORBs) represent a viable pathway to a more sustainable energy storage technology compared to conventional Li-ion batteries. For further materials and cell development towards competitive energy and power densities, a deeper understanding of electron transport and conductivity in organic radical polymer cathodes is required. Such electron transport is characterised by electron hopping processes, which depend on the presence of closely spaced hopping sites. Using a combination of electrochemical, electron paramagnetic resonance (EPR) spectroscopic, and theoretical molecular dynamics as well as density functional theory modelling techniques, we explored how compositional characteristics of cross-linked poly(2,2,6,6-tetramethyl-1-piperidinyloxy-4-yl methacrylate) (PTMA) polymers govern electron hopping and rationalise their impact on ORB performance. Electrochemistry and EPR spectroscopy not only show a correlation between capacity and the total number of radicals in an ORB using a PTMA cathode, but also indicates that the state-of-health degrades about twice as fast if the amount of radical is reduced by 15%. The presence of up to 3% free monomer radicals did not improve fast charging capabilities. Pulsed EPR indicated that these radicals readily dissolve into the electrolyte but a direct effect on battery degradation could not be shown. However, a qualitative impact cannot be excluded either. The work further illustrates that nitroxide units have a high affinity to the carbon black conductive additive, indicating the possibility of its participation in electron hopping. At the same time, the polymers attempt to adopt a compact conformation to increase radical-radical contact. Hence, a kinetic competition exists, which might gradually be altered towards a thermodynamically more stable configuration by repeated cycling, yet further investigations are required for its characterisation.