Charge-Scaling Effect in Ionic Liquids from the Charge-Density Analysis ofN,N′-Dimethylimidazolium Methylsulfate

Diffusion and structure of propylene carbonate-metal salt electrolyte solutions for post-lithium-ion batteries: From experiment to simulation

The diffusion of cations in organic solvent solutions is important for the performance of metal-ion batteries. In this article, pulsed field gradient nuclear magnetic resonance experiments and fully atomistic molecular dynamic simulations were employed to study the temperature-dependent diffusive behavior of various liquid electrolytes representing 1M propylene carbonate solutions of metal salts with bis(trifluoromethylsulfonyl)imide (TFSI-) or hexafluorophosphate (PF6-) anions commonly used in lithium-ion batteries and beyond. The experimental studies revealed the temperature dependence of the diffusion coefficients for the propylene carbonate (PC) solvent and for the anions following an Arrhenius type of behavior. It was observed that the PC molecules are the faster species. For the monovalent cations (Li+, Na+, K+), the PC solvent diffusion was enhanced as the cation size increased, while for the divalent cations (Mg2+, Ca2+, Sr2+, Ba2+), the opposite trend was observed, i.e., the diffusion coefficients decreased as the cation size increased. The anion diffusion in LiTFSI and NaTFSI solutions was found to be similar, while in electrolytes with divalent cations, a decrease in anion diffusion with increasing cation size was observed. It was shown that non-polarizable charge-scaled force fields could correspond perfectly to the experimental values of the anion and PC solvent diffusion coefficients in salt solutions of both monovalent (Li+, Na+, K+) and divalent (Mg2+, Ca2+, Sr2+, Ba2+) cations at a range of operational temperatures. Finally, after calculating the radial distribution functions between cations, anions, and solvent molecules, the increase in the PC diffusion coefficient established with the increase in cation size for monovalent cations was clearly explained by the large hydration shell of small Li+ cations, due to their strong interaction with the PC solvent. In solutions with larger monovalent cations, such as Na+, and with a smaller solvation shell of PC, the PC diffusion is faster due to more liberated solvent molecules. In the salt solutions with divalent cations, both the anion and the PC diffusion coefficients decreased as the cation size increased due to an enhanced cation-anion coordination, which was accompanied by an increase in the amount of PC in the cation solvation shell due to the presence of anions.

Alkylcyanoborate Anions: Building Blocks for Fluorine-Free Low-Viscosity, Electrochemically and Thermally Stable Ionic Liquids

A set of mixed-substituted potassium alkylcyano- and alkylcyanofluoroborates has been synthesized using easily accessible starting compounds and characterized by elemental analysis, NMR and vibrational spectroscopy, and mass spectrometry. In addition, single-crystal structures of salts of the cyanoborate anions have been derived from X-ray diffraction experiments. The 1-ethyl-3-methylimidazolium room temperature ionic liquids ([EMIm]+ -RTILs) with the new borate anions have been synthesized and their physicochemical properties, that is, high thermal and electrochemical stability, low viscosity, and high conductivity, have been compared to the properties of related [EMIm]+ -RTILs. The influence of the different alkyl substituents at boron has been assessed. The exemplary study on the properties with the [EMIm]+ -ILs with the mixed water stable alkylcyanoborate anions points towards the potential of these fluorine-free borate anions, in general.

Molecular Modeling of Water-in-Salt Electrolytes: A Comprehensive Analysis of Polarization Effects and Force Field Parameters in Molecular Dynamics Simulations

Accurate modeling of highly concentrated aqueous solutions, such as water-in-salt (WiS) electrolytes in battery applications, requires proper consideration of polarization contributions to atomic interactions. Within the force field molecular dynamics (MD) simulations, the atomic polarization can be accounted for at various levels. Nonpolarizable force fields implicitly account for polarization effects by incorporating them into their van der Waals interaction parameters. They can additionally mimic electron polarization within a mean-field approximation through ionic charge scaling. Alternatively, explicit polarization description methods, such as the Drude oscillator model, can be selectively applied to either a subset of polarizable atoms or all polarizable atoms to enhance simulation accuracy. The trade-off between simulation accuracy and computational efficiency highlights the importance of determining an optimal level of accounting for atomic polarization. In this study, we analyze different approaches to include polarization effects in MD simulations of WiS electrolytes, with an example of a Na-OTF solution. These approaches range from a nonpolarizable to a fully polarizable force field. After careful examination of computational costs, simulation stability, and feasibility of controlling the electrolyte properties, we identify an efficient combination of force fields: the Drude polarizable force field for salt ions and non-polarizable models for water. This cost-effective combination is sufficiently flexible to reproduce a broad range of electrolyte properties, while ensuring simulation stability over a relatively wide range of force field parameters. Furthermore, we conduct a thorough evaluation of the influence of various force field parameters on both the simulation results and technical requirements, with the aim of establishing a general framework for force field optimization and facilitating parametrization of similar systems.

Charge transfer and polarisability in ionic liquids: a case study

The practical use of ionic liquids (ILs) is benefiting from a growing understanding of the underpinning structural and dynamic properties, facilitated through classical molecular dynamics (MD) simulations. The predictive and explanatory power of a classical MD simulation is inextricably linked to the underlying force field. A key aspect of the forcefield for ILs is the ability to recover charge based interactions. Our focus in this paper is on the description and recovery of charge transfer and polarisability effects, demonstrated through MD simulations of the widely used 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C₄C₁im][NTf₂] IL. We study the charge distributions generated by a range of ab initio methods, and present an interpolation method for determining atom-wise scaled partial charges. Two novel methods for determining the mean field (total) charge transfer from anion to cation are presented. The impact of using different charge models and different partial charge scaling (unscaled, uniformly scaled, atom-wise scaled) are compared to fully polarisable simulations. We study a range of Drude particle explicitly polarisable potentials and shed light on the performance of current approaches to counter known problems. It is demonstrated that small changes in the charge description and MD methodology can have a significant impact; biasing some properties, while leaving others unaffected within the structural and dynamic domains.

Coupled hydroxyl and ether functionalisation in EAN derivatives: the effect of hydrogen bond donor/acceptor groups on the structural heterogeneity studied with X-ray diffractions and fixed charge/polarizable simulations

We present a study by energy-dispersive X-ray diffraction of liquid 2-(2-hydroxyethoxy)ethan-1-ammonium nitrate, NH₃CH₂CH₂(OCH₂CH₂OH)⁺NO₃^- (22HHEAN). This ionic liquid is derived from the parent ethylammonium nitrate (EAN) with an ether link in the chain and a hydroxyl group in the terminal position. The absence of peaks at low-q values in the experimental diffraction curve indicates that the added polar groups and the high conformational isomerism of the cations alter strongly the nanosegregation of the parent EAN liquid. Aggregation between ionic species may involve hydrogen bonding between cations and anions and a variety of intermolecular hydrogen bonds between cations. Diffraction patterns are compared with the results of molecular dynamics simulations with two different force fields: the fixed point charge force field (GAFF) with different charge scaling protocols and the polarizable AMOEBA force field. Most point charge models lead to the appearance of a quite evident low q-peak which decreases gradually, when the percentage and type of the scaling (uniform vs. non-uniform) are increased. In the polarisable model and in the model where only anion charges are scaled to 20%, instead, the pre-peak is absent in agreement with our experiments.

Doxycycline hydrate and doxycycline hydrochloride dihydrate – crystal structure and charge density analysis

High-resolution low-temperature X-ray diffraction experiments for doxycycline monohydrate and hydrochloride dihydrate have been performed. Translation-Libration-Screw (TLS) analysis for both crystal forms as well as the data from neutron diffraction experiment for hydrochloride combined with the Hansen-Coppens formalism resulted in precise charge density distribution models for both the zwitterionic monohydrate and a protonated hydrochloride crystal forms. Their detailed topological analysis suggested that the electron structure of doxycycline’s amide moiety undergoes significant changes during protonation due to formation of a very strong resonance-assisted hydrogen bond. A notably increased participation of amide nitrogen atom and hydrogen-accepting oxygen atom in the resonance upon doxycycline protonation was observed. A comparison of TLS- and neutron data-derived hydrogen parameters confirmed the experimental neutron data to be vital for proper description of intra- and inter-molecular interactions in this compound. Finally, calculated lattice and interaction energies quantified repulsive Dox-Dox interactions in the protonated crystal form of the antibiotic, relating with a good solubility of doxycycline hydrochloride relative to its hydrate.

Accurate prediction of the structure and vibrational spectra of ionic liquid clusters with the generalized energy-based fragmentation approach: critical role of ion-pair-based fragmentation

The GEBF method with the ion-pair-based fragmentation has been developed to facilitate ab initio calculations of general ionic liquid clusters.

Poly[n]catenanes: Synthesis of molecular interlocked chains

As the macromolecular version of mechanically interlocked molecules, mechanically interlocked polymers are promising candidates for the creation of sophisticated molecular machines and smart soft materials. Poly[n]catenanes, where the molecular chains consist solely of interlocked macrocycles, contain one of the highest concentrations of topological bonds. We report, herein, a synthetic approach toward this distinctive polymer architecture in high yield (~75%) via efficient ring closing of rationally designed metallosupramolecular polymers. Light-scattering, mass spectrometric, and nuclear magnetic resonance characterization of fractionated samples support assignment of the high-molar mass product (number-average molar mass ~21.4 kilograms per mole) to a mixture of linear poly[7-26]catenanes, branched poly[13-130]catenanes, and cyclic poly[4-7]catenanes. Increased hydrodynamic radius (in solution) and glass transition temperature (in bulk materials) were observed upon metallation with Zn².

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.

Photoswitching of Azobenzene-Based Reverse Micelles above and at Subzero Temperatures As Studied by NMR and Molecular Dynamics Simulations

We designed and studied the structure, dynamics, and photochemistry of photoswitchable reverse micelles (RMs) composed of azobenzene-containing ammonium amphiphile 1 and water in chloroform at room and subzero temperatures by NMR spectroscopy and molecular dynamics simulations. The NMR and diffusion coefficient analyses showed that micelles containing either the E or Z configuration of 1 are stable at room temperature. Depending on the water-to-surfactant molar ratio, the size of the RMs remains unchanged or is slightly reduced because of the partial loss of water from the micellar cores upon extensive E → Z or Z → E photoisomerization of the azobenzene group in 1. Upon freezing at 253 or 233 K, E-1 RMs partially precipitate from the solution but are redissolved upon warming whereas Z-1 RMs remain fully dissolved at all temperatures. Light-induced isomerization of 1 at low temperatures does not lead to the disintegration of RMs remaining in the solution; however, its scope is influenced by a precipitation process. To obtain a deeper molecular view of RMs, their structure was characterized by MD simulations. It is shown that RMs allow for amphiphile isomerization without causing any immediate significant structural changes in the micelles.

Experimental and theoretical studies on solvation in aqueous solutions of ionic liquids carrying different side chains: the n-butyl-group versus the methoxyethyl group

We used solvatochromic compounds to probe solvation in mixtures of water, W, and four ionic liquids (ILs), 1-R-3-methylimidazoliumX, where R =n-butyl or methoxyethyl and X = acetate and chloride.

A New Volume-Based Approach for Predicting Thermophysical Behavior of Ionic Liquids and Ionic Liquid Crystals

Volume-based prediction of melting points and other properties of ionic liquids (ILs) relies on empirical relations with volumes of ions in these low-melting organic salts. Here we report an accurate way to ionic volumes by Bader's partitioning of electron densities from X-ray diffraction obtained via a simple database approach. For a series of 1-tetradecyl-3-methylimidazolium salts, the volumes of different anions are found to correlate linearly with melting points; larger anions giving lower-melting ILs. The volume-based concept is transferred to ionic liquid crystals (ILs that adopt liquid crystalline mesophases, ILCs) for predicting the domain of their existence from the knowledge of their constituents. For 1-alkyl-3-methylimidazolium ILCs, linear correlations of ionic volumes with the occurrence of LC mesophase and its stability are revealed, thus paving the way to rational design of ILCs by combining suitably sized ions.

Understanding the hydrogen bonds in ionic liquids and their roles in properties and reactions

Ionic liquids (ILs) have many potential applications in the chemical industry. In order to understand ILs, their molecular details have been extensively investigated. Intuitively, electrostatic forces are solely important in ILs. However, experiments and calculations have provided strong evidence for the existence of H-bonds in ILs and their roles in the properties and applications of ILs. As a structure-directing force, H-bonds are responsible for ionic pairing, stacking and self-assembling. Their geometric structure, interaction energy and electronic configuration in the ion-pairs of imidazolium-based ILs and protic ionic liquids (PILs) show a great number of differences compared to conventional H-bonds. In particular, their cooperation with electrostatic, dispersion and π interactions embodies the physical nature of H-bonds in ILs, which anomalously influences their properties, leading to a decrease in their melting points and viscosities and thus fluidizing them. Using ILs as catalysts and solvents, many reactions can be activated by the presence of H-bonds, which reduce the reaction barriers and stabilize the transition states. In the dissolution of lignocellulosic biomass by ILs, H-bonds exhibit a most important role in disrupting the H-bonding network of cellulose and controlling microscopic ordering into domains. In this article, a critical review is presented regarding the structural features of H-bonds in ILs and PILs, the correlation between H-bonds and the properties of ILs, and the roles of H-bonds in typical reactions.

Probing the effects of the ester functional group, alkyl side chain length and anions on the bulk nanostructure of ionic liquids: a computational study

The effects of ester addition on nanostructural properties of biodegradable ILs.

Quantification of Conventional and Nonconventional Charge-Assisted Hydrogen Bonds in the Condensed and Gas Phases

Charge-assisted hydrogen bonds (CAHBs) play critical roles in many systems from biology through to materials. In none of these areas has the role and function of CAHBs been explored satisfactorily because of the lack of data on the energy of CAHBs in the condensed phases. We have, for the first time, quantified three types of CAHBs in both the condensed and gas phases for 1-(2'-hydroxylethyl)-3-methylimidazolium acetate ([C2OHmim][OAc]). The energy of conventional OH···[OAc](-) CAHBs is ∼10 kcal·mol(-1), whereas nonconventional C(sp2)H···[OAc](-) and C(sp3)H···[OAc](-) CAHBs are weaker by ∼5-7 kcal·mol(-1). In the gas phase, the strength of the nonconventional CAHBs is doubled, whereas the conventional CAHBs are strengthened by <20%. The influence of cooperativity effects on the ability of the [OAc](-) anion to deprotonate the imidazolium cation is evaluated. The ability to quantify CAHBs in the condensed phase on the basis of easier accessible gas-phase estimates is highlighted.

Hydrogen bonding in ionic liquids

Ionic liquids (IL) and hydrogen bonding (H-bonding) are two diverse fields for which there is a developing recognition of significant overlap. Doubly ionic H-bonds occur when a H-bond forms between a cation and anion, and are a key feature of ILs. Doubly ionic H-bonds represent a wide area of H-bonding which has yet to be fully recognised, characterised or explored. H-bonds in ILs (both protic and aprotic) are bifurcated and chelating, and unlike many molecular liquids a significant variety of distinct H-bonds are formed between different types and numbers of donor and acceptor sites within a given IL. Traditional more neutral H-bonds can also be formed in functionalised ILs, adding a further level of complexity. Ab initio computed parameters; association energies, partial charges, density descriptors as encompassed by the QTAIM methodology (ρBCP), qualitative molecular orbital theory and NBO analysis provide established and robust mechanisms for understanding and interpreting traditional neutral and ionic H-bonds. In this review the applicability and extension of these parameters to describe and quantify the doubly ionic H-bond has been explored. Estimating the H-bonding energy is difficult because at a fundamental level the H-bond and ionic interaction are coupled. The NBO and QTAIM methodologies, unlike the total energy, are local descriptors and therefore can be used to directly compare neutral, ionic and doubly ionic H-bonds. The charged nature of the ions influences the ionic characteristics of the H-bond and vice versa, in addition the close association of the ions leads to enhanced orbital overlap and covalent contributions. The charge on the ions raises the energy of the Ylp and lowers the energy of the X-H σ* NBOs resulting in greater charge transfer, strengthening the H-bond. Using this range of parameters and comparing doubly ionic H-bonds to more traditional neutral and ionic H-bonds it is clear that doubly ionic H-bonds cover the full range of weak through to very strong H-bonds.

Invariom approach as a new tool in electron density studies of ionic liquids: a model case of 1-butyl-2,3-dimethylimidazolium chloride BDMIM[Cl]

In this comparative study of an ionic liquid (IL), BDMIM[Cl], the invariom approximation emerges as a useful tool in search for ‘structure–property’ relations in ILs, as entangling complex interplay of interionic interactions that operate in them.