Conformational isomerisms and nano-aggregation in substituted alkylammonium nitrates ionic liquids: An x-ray and computational study of 2-methoxyethylammonium nitrate

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

New insights into chloromethyl-oxirane and chloromethyl-thiirane in liquid and solid phase from low-temperature infrared spectroscopy and ab initio modeling

A detailed study of the conformational landscape of chloromethyl-oxirane and chloromethyl-thiirane is here reported. The equilibrium of the three different conformers of the two molecules was assessed, using a joint approach of experimental and theoretical methods. High quality infrared spectroscopy measurements of the liquid and of the crystalline phases were interpreted with the aid of ab initio Molecular Dynamics (AIMD) simulations, anharmonic frequencies and free energy calculations, obtaining a very good reproduction of the experimental data. The modulation of the conformational equilibrium upon the addition of polar and non-polar solvents was computationally evaluated and results found a confirmation in experimental measures.

New Insights into the Structure of Glycols and Derivatives: A Comparative X-Ray Diffraction, Raman and Molecular Dynamics Study of Ethane-1,2-Diol, 2-Methoxyethan-1-ol and 1,2-Dimethoxy Ethane

In this study, we report a detailed experimental and theoretical investigation of three glycol derivatives, namely ethane-1,2-diol, 2-methoxyethan-1-ol and 1,2-dimethoxy ethane. For the first time, the X-ray spectra of the latter two liquids was measured at room temperature, and they were compared with the newly measured spectrum of ethane-1,2-diol. The experimental diffraction patterns were interpreted very satisfactorily with molecular dynamics calculations, and suggest that in liquid ethane-1,2-diol most molecules are found in gauche conformation, with intramolecular hydrogen bonds between the two hydroxyl groups. Intramolecular H-bonds are established in the mono-alkylated diol, but the interaction is weaker. The EDXD study also evidences strong intermolecular hydrogen-bond interactions, with short O···O correlations in both systems, while longer methyl-methyl interactions are found in 1,2-dimethoxy ethane. X-ray studies are complemented by micro Raman investigations at room temperature and at 80 °C, that confirm the conformational analysis predicted by X-ray experiments and simulations.

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.

Charge Environment and Hydrogen Bond Dynamics in Binary Ionic Liquid Mixtures: A Computational Study

The extent of charge transfer between the cation and the anion in a room-temperature ionic liquid depends on the basicity of the anion. Ion charges determined in the condensed state via density functional theory calculations capture this effect rather well, and charges derived in such a manner have been employed in force field-based molecular dynamics simulations to quantitatively reproduce several physical properties of the liquids. However, the issue of transferability of cation charges in mixtures of ionic liquids, say with one type of cation and two different anion types needs to be addressed. Herein, we demonstrate that the cation charge in such a mixture varies linearly with anion composition, a result that ties in rather well with X-ray photoelectron spectroscopic experiments. The variation in cation charge with bulk anion composition is shown to be a result of changes in its coordination environment. Cations surrounded by a higher proportion of more basic anions possess lower charges than those surrounded by less basic anions. Time scales for the exchange of anion types for the occupation of hydrogen bonding sites around the cation have been determined and are seen to be constituted by three processes-breakage of existing hydrogen bond, diffusion to the hydrogen bonding site and displacement of the incumbent anion from its site in the cation coordination shell. These time scales explain the differences observed between infrared and NMR spectroscopic experiments in ionic liquid mixtures rather well.

Structure and dynamics of propylammonium nitrate-acetonitrile mixtures: An intricate multi-scale system probed with experimental and theoretical techniques

In this article, we report the study of structural and dynamical properties for a series of acetonitrile/propylammonium nitrate mixtures as a function of their composition. These systems display an unusual increase in intensity in their X-ray diffraction patterns in the low-q regime, and their ¹H-NMR diffusion-ordered NMR spectroscopy (DOSY) spectra display unusual diffusivities. However, the magnitude of both phenomena for mixtures of propylammonium nitrate is smaller than those observed for ethylammonium nitrate mixtures with the same cosolvent, suggesting that the cation alkyl tail plays an important role in these observations. The experimental X-ray scattering data are compared with the results of molecular dynamics simulations, including both ab initio studies used to interpret short-range interactions and classical simulations to describe longer range interactions. The higher level calculations highlight the presence of a strong hydrogen bond network within the ionic liquid, only slightly perturbed even at high acetonitrile concentration. These strong interactions lead to the symmetry breaking of the NO₃^- vibrations, with a splitting of about 88 cm^-1 in the ν₃ antisymmetric stretch. The classical force field simulations use a greater number of ion pairs, but are not capable of fully describing the longest range interactions, although they do successfully account for the observed concentration trend, and the analysis of the models confirms the nano-inhomogeneity of these kinds of samples.

Bio ionic liquids and water mixtures: a structural study

In this study we have explored, by means of ab initio molecular dynamics, a subset of three different water/cho⁺–phe⁻ mixtures.

X-ray and molecular dynamics studies of butylammonium butanoate–water binary mixtures

The nanostructural organisation of mixtures of the ionic liquid (butylammonium butanoate) and water at several mole fractions of water has been investigated using small and wide angle X-ray scattering (S-WAXS) and molecular dynamics (MD) simulations.

Direct Comparison of Atomistic Molecular Dynamics Simulations and X-ray Scattering of Polymerized Ionic Liquids

The design of solid-state electrolytes for electrochemical applications that utilize polymerized ionic liquids (polyILs) would greatly benefit from a molecular-level understanding of structure-function relationships. We herein use atomistic molecular dynamics simulations to investigate the structural properties of a homologous series of poly(n-alkyl-vinylimidzolium bistrifluoromethylsulfonylimide) poly(nVim Tf2N) and present the first direct comparison of the structure factors obtained from X-ray scattering and simulations. Excellent agreement is found in terms of peak position and shape. The backbone-to-backbone correlation length increases at a rate of 1 Å/CH₂. The longer alkyl chains lead to the longer backbone-to-backbone separation and the larger nonpolar nanodomains. This quantitative comparison of atomistic simulations to X-ray scattering will lead to a fundamental understanding in structure and morphology of polyILs and pave a path forward toward the rational design of future polyILs for electrochemical devices.

Direct calculation of the X-ray structure factor of ionic liquids

A conceptually simple and computationally efficient direct method to calculate the total X-ray structure factor of ionic liquids from molecular simulations is advocated to be complementary to the popular Fourier transform (FT) method.

Is a medium-range order pre-peak possible for ionic liquids without an aliphatic chain?

The pre-peak in the X-ray scattering patterns of ionic liquids is not always due to alkyl chains. In choline–proline (a bio-compatible ionic liquid), it is due to second-shell cation–cation distances.

Cholinium-amino acid based ionic liquids: a new method of synthesis and physico-chemical characterization

Fourteen cholinium-amino acid based room temperature ionic liquids were prepared using a cleaner synthetic method. Chemicophysical properties were well correlated with the wide range of amino acid chemical structures.