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.

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.

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.

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.

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.

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.

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.

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.

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.

Insights into the Solubility of Poly(vinylphenothiazine) in Carbonate-Based Battery Electrolytes.

Organic materials are promising candidates for next-generation battery systems. However, many organic battery materials suffer from high solubility in common battery electrolytes. Such solubility can be overcome by introducing tailored high-molecular-weight polymer structures, for example, by cross-linking, requiring enhanced synthetic efforts. We herein propose a different strategy by optimizing the battery electrolyte to obtain insolubility of non-cross-linked poly(3-vinyl-N-methylphenothiazine) (PVMPT). Successive investigation and theoretical insights into carbonate-based electrolytes and their interplay with PVMPT led to a strong decrease in the solubility of the redox polymer in ethylene carbonate/ethyl methyl carbonate (3:7) with 1 M LiPF6. This allowed accessing its full theoretical specific capacity by changing the charge/discharge mechanism compared to previous reports. Through electrochemical, spectroscopic, and theoretical investigations, we show that changing the constituents of the electrolyte significantly influences the interactions between the electrolyte molecules and the redox polymer PVMPT. Our study demonstrates that choosing the ideal electrolyte composition without chemical modification of the active material is a successful strategy to enhance the performance of organic polymer-based batteries.

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.

rs@md: Introducing Reactive Steps at the Molecular Dynamics Simulation Level.

A concept is presented to extend molecular dynamics simulations by the so-called reactive steps, during which transitions from reactant to product molecules are performed with physically correct transition probabilities. This goes along with an instant exchange of the employed force field. We provide a detailed mathematical derivation for how the acceptance probability for such reactive steps can be computed from molecular reaction rates and introduce a simulation program that performs such reactive step molecular dynamics simulations. Our program is designed in a modular fashion and can thus be extended to any conventional molecular dynamics program. Furthermore, the working principle of these reaction rate-based reactive step simulations is demonstrated by applying them to a reactive model system based on associating and dissociating Lennard-Jones particles and compared to a similar approach from Nagaoka et al. which uses the Metropolis Monte Carlo scheme for the reactive steps. Overall, we find that our approach not only recovers the correct thermodynamics but also ensures proper kinetics, that is, the correct time evolution of the system.

Dioxolanone-Anchored Poly(allyl ether)-Based Cross-Linked Dual-Salt Polymer Electrolytes for High-Voltage Lithium Metal Batteries.

Novel cross-linked polymer electrolytes (XPEs) are synthesized by free-radical copolymerization induced by ultraviolet (UV)-light irradiation of a reactive solution, which is composed of a difunctional poly(ethylene glycol) diallyl ether oligomer (PEGDAE), a monofunctional reactive diluent 4-vinyl-1,3-dioxolan-2-one (VEC), and a stock solution containing lithium salt (lithium bis(trifluoromethanesulfonyl)imide, LiTFSI) in a carbonate-free nonvolatile plasticizer, poly(ethylene glycol) dimethyl ether (PEGDME). The resulting polymer matrix can be represented as a linear polyethylene chain functionalized with cyclic carbonate (dioxolanone) moieties and cross-linked by ethylene oxide units. A series of XPEs are prepared by varying the [O]/[Li] ratio (24 to 3) of the stock solution and thoroughly characterized using physicochemical (thermogravimetric analysis-mass spectrometry, differential scanning calorimetry, NMR, etc.) and electrochemical techniques. In addition, quantum chemical calculations are performed to elucidate the correlation between the electrochemical oxidation potential and the lithium ion-ethylene oxide coordination in the stock solution. Later, lithium bis(fluorosulfonyl)imide (LiFSI) salt is incorporated into the electrolyte system to produce a dual-salt XPE that exhibits improved electrochemical performance, a stable interface against lithium metal, and enhanced physical and chemical characteristics to be employed against high-voltage cathodes. The XPE membranes demonstrated excellent resistance against lithium dendrite growth even after reversibly plating and stripping lithium ions for more than 1000 h with a total capacity of 0.5 mAh cm-2. Finally, the XPE films are assembled in a lab-scale lithium metal battery configuration by using carbon-coated LiFePO4 (LFP) or LiNi0.8Co0.15Al0.05O2 (NCA) as a cathode and galvanostatically cycled at 20, 40, and 60 °C. Remarkably, at 20 °C, the NCA-based lithium metal cells displayed excellent cycling stability and good capacity retention (>50%) even after 1000 cycles.

Fluorinated Cyclic Phosphorus(III)-Based Electrolyte Additives for High Voltage Application in Lithium-Ion Batteries: Impact of Structure-Reactivity Relationships on CEI Formation and Cell Performance.

Two selected and designed fluorinated cyclic phosphorus(III)-based compounds, namely 2-(2,2,3,3,3-pentafluoropropoxy)-1,3,2-dioxaphospholane (PFPOEPi) and 2-(2,2,3,3,3-pentafluoro-propoxy)-4-(trifluormethyl)-1,3,2-dioxaphospholane (PFPOEPi-1CF3), were synthesized and comprehensively characterized for high voltage application in lithium-ion batteries (LIBs). Cyclic voltammetry (CV) and constant current cycling were conducted, followed by post mortem analysis of the NMC111 electrode surface via scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS). To support and complement obtained experimental results, density functional theory (DFT) calculations and molecular dynamics (MD) simulations were performed. Theoretical and experimental findings show that the considered phospholane molecule class enables high voltage LIB application by sacrificial decomposition on the cathode surface and involvement in the formation of a cathode electrode interphase (CEI) via polymerization reaction. In addition, obtained results point out that the introduction of the CF3 group has a significant influence on the formation and dynamics of the CEI as well as on the overall cell performance, as the cell impedance as well as the thickness of the CEI is increased compared to the cells containing PFPOEPi, which results in a decreased cycling performance. This systematic approach allows researchers to understand the structure-reactivity relationship of the newly synthesized compounds and helps to further tailor the vital physicochemical properties of functional electrolyte additives relevant for high voltage LIB application.

Improved lithium ion dynamics in crosslinked PMMA gel polymer electrolyte.

Since PMMA-based gel polymer electrolytes could substitute PVDF-HFP based gels currently used in Li-ion batteries at lower financial and environmental costs, we investigate here the solvation and transport properties of the lithium ions in a crosslinked PMMA-based gel polymer electrolyte by a combination of thermal and electrochemical methods, Raman spectroscopy, pulse field gradient (PFG) and electrophoretic NMR (eNMR) techniques, as well as ab initio calculations. The conductivity of the gel containing 10 wt% polymer is only reduced by 14% relative to the liquid electrolyte. In addition, the co-solvation by polymer functional groups, a priori expected to slow lithium transport relatively to the anion, has instead a positive effect on lithium transport. Indeed, the ester groups not only participate in lithium solvation and increase ionic dissociation, but since this interaction is rather weak, rather than lowering the lithium diffusion relatively to other species, it mainly decorrelates lithium transport from anionic mobility. Compared to its liquid fraction, the gels show, at the same time, better dissociation and a higher lithium transference number, which results in a higher cationic conductivity, despite the overall conductivity loss.

Microscopic Structure of Compacted Polyelectrolyte Complexes: Insights from Molecular Dynamics Simulations.

We utilize atomistic molecular dynamics (MD) simulations to study local structural changes inside a polyelectrolyte complex consisting of poly(styrenesulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) upon densification, in analogy to ultracentrifugation in experiments. In particular, we focus on the water content and on the reinforcement of the PSS-PDADMA network for various external accelerations. We demonstrate that apart from the formation of mesoscopic pores observed experimentally also the microscopic structure and the local relaxation processes likely affect the unique rheological properties of compacted polyelectrolyte complexes, as densification increases both the number of PSS-PDADMA coordinations and the intermixing of PSS and PDADMA. These processes slow down local rearrangements, thus further stabilizing the compacted state. We find that the concept of binary PSS-PDADMA salt bonds-relevant for theoretical models-is not strictly valid in the dense limit.

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.

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.