Properties of hydrogen bonds in the protic ionic liquid ethylammonium nitrate: DFT versus DFTB molecular dynamics

Storage and Diffusion of Carbon Dioxide in the Metal Organic Framework MOF-5─A Semi-empirical Molecular Dynamics Study

Metal-organic frameworks (MOFs) have attracted increasing attention due to their high porosity for exceptional gas storage applications. MOF-5 belongs to the family of isoreticular MOFs (IRMOFs) and consists of Zn4O6+ clusters linked by 1,4-benzenedicarboxylate. Due to the large number of atoms in the unit cell, molecular dynamics simulation based on density functional theory has proved to be too demanding, while force field models are often inadequate to model complex host-guest interactions. To overcome this limitation, an alternative semi-empirical approach using a set of approximations and extensive parametrization of interactions called density functional tight binding (DFTB) was applied in this work to study CO2 in the MOF-5 host. Calculations of pristine MOF-5 yield very good agreement with experimental data in terms of X-ray diffraction patterns as well as mechanical properties, such as the negative thermal expansion coefficient and the bulk modulus. In addition, different loadings of CO2 were introduced, and the associated self-diffusion coefficients and activation energies were investigated. The results show very good agreement with those of other experimental and theoretical investigations. This study provides detailed insights into the capability of semi-empirical DFTB-based molecular dynamics simulations of these challenging guest@host systems. Based on the comparison of the guest-guest pair distributions observed inside the MOF host and the corresponding gas-phase reference, a liquid-like structure of CO2 can be deduced upon storage in the host material.

Effects of Me-Solvent Interactions on the Structure and Infrared Spectra of MeTFSI (Me = Li, Na) Solutions in Carbonate Solvents-A Test of the GFN2-xTB Approach in Molecular Dynamics Simulations

We investigated the performance of the computationally effective GFN2-xTB approach in molecular dynamics (MD) simulations of liquid electrolytes for lithium/sodium batteries. The studied systems were LiTFSI and NaTFSI solutions in ethylene carbonate or fluoroethylene carbonate and the neat solvents. We focused on the structure of the electrolytes and on the manifestations of ion-solvent interactions in the vibrational spectra. The IR spectra were calculated from MD trajectories as Fourier transforms of the dipole moment. The results were compared to the data obtained from ab initio MD. The spectral shifts of the carbonyl stretching mode calculated from the GFN2-xTB simulations were in satisfactory agreement with the ab initio MD data and the experimental results for similar systems. The performance in the region of molecular ring vibrations was significantly worse. We also found some differences in structural data, suggesting that the GFN2-xTB overestimates interactions of Me ions with TFSI anions and Na+ binding to solvent molecules. We conclude that the GFN2-xTB method is an alternative worth considering for MD simulations of liquids, but it requires testing of its applicability for new systems.

State of anion in ethylammonium nitrate enclosed between micrometer-spaced glass plates as studied by 17O and 15N NMR

Some aprotic and protic ionic liquids (ILs) containing nitrate anion demonstrate unusual dynamic behavior of cations when these ILs are enclosed in micrometer-spaced layers between glass plates. We applied ¹⁷O and ¹⁵N NMR spectroscopy to discover the state and transformations of ¹⁷O and ¹⁵N isotopically enriched nitrate anion of ethylammonium nitrate (EAN) enclosed between glass plates. ¹⁵N NMR spectra demonstrated preferential orientation of the principal axes of the nitrate anions perpendicular to the normal of the glass surface. Therefore, isotropic ionic liquid EAN, when placed within a micrometer-spaced enclosure, forms an ordered phase, which is similar to a liquid crystal. The peculiarity of this phase is that the cations do not have a predominant orientation. Other features of this phase that are typical for liquid crystal phases are the changed local and translational dynamics in comparison with the isotropic state and slow transformation occurring under the action of an external magnetic field.

Hydrogen Bonding and Infrared Spectra of Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide/Water Mixtures: A View from Molecular Dynamics Simulations

Simulations of ab initio molecular dynamics have been performed for mixtures of ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) ionic liquid and water. Statistics of donors and acceptors of hydrogen bonds has revealed that with increasing water content, hydrogen bonds between EMIM cations and TFSI anions are replaced by bonds to water molecules. In the mixture of liquids, the total number of bonds (from EMIM cations or water molecules) formed by TFSI acceptors increases. IR spectra obtained from ab initio molecular dynamics trajectories are in good agreement with literature data for ionic liquid/water systems. Analysis of oscillations of individual C-H and O-H bonds has shown correlations between vibrational frequencies and hydrogen bonds formed by an EMIM cation or water molecule and has indicated that the changes in the IR spectrum result from the decreased number of water-water hydrogen bonds in the mixture. The tests of DFTB methodology with tailored parameterizations have yielded reasonably good description of the IR spectrum of bulk water, whereas available parameterizations have failed in satisfactory reproduction of the IR spectrum of EMIM-TFSI/water mixtures in the region above 3000 cm^-1.

A density functional theory based tight-binding study on the water effect on nanostructuring of choline chloride + ethylene glycol deep eutectic solvent

The effect of water on the properties of an archetypical type III deep eutectic solvent [choline chloride : ethyleneglycol (1:2)] is analyzed using ab initio molecular dynamics simulations in the 0 to 60 wt. % water content range. The properties of the mixed fluids are studied considering nanostructuring, intermolecular forces (hydrogen bonding), the energy of interactions, dynamic properties, and domain analysis. The reported results confirm that the change in the properties of the studied deep eutectic solvent is largely dependent on the amount of water. The competing effect of water molecules for the available hydrogen bonding sites determines the evolution of the properties upon water sorption. The main structural features of the considered deep eutectic were maintained even for large water contents; thus, its hydrophilicity could be used for tuning fluid physicochemical properties.

Assessing the Structure of Protic Ionic Liquids Based on Triethylammonium and Organic Acid Anions

We present a computational analysis of the short-range structure of three protic ionic liquids based on strong organic acids: trifluoracetate, methanesulfonate, and triflate of triethylammonium. Accurate ab initio computations carried out on the gas-phase dimers show that the protonation of triethylamine is spontaneous. We have identified the anion-cation binding motif that is due to the presence of a strong hydrogen bond and to electrostatic interactions. The strength of the hydrogen bond and the magnitude of the binding energy decrease in the order trifluoroacetate ≳ methanesulfonate > triflate. The corresponding simulations of the bulk phases, obtained using a semiempirical evaluation of the interatomic forces, reveal that on short timescales, the state of the three liquids remains highly ionized and that the gas-phase cation-/anion-binding motif is preserved while no other peculiar structural features seem to emerge.

Extension of the CL&Pol Polarizable Force Field to Electrolytes, Protic Ionic Liquids, and Deep Eutectic Solvents

The polarizable CL&Pol force field presented in our previous study, Transferable, Polarizable Force Field for Ionic Liquids (J. Chem. Theory Comput. 2019, 15, 5858, DOI: http://doi.org/10.1021/acs.jctc.9b0068910.1021/acs.jctc.9b00689), is extended to electrolytes, protic ionic liquids (PIL), deep eutectic solvents (DES), and glycols. These systems are problematic in polarizable simulations because they contain either small, highly charged ions or strong hydrogen bonds, which cause trajectory instabilities due to the pull exerted on the induced dipoles. We use a Tang-Toennies (TT) function to dampen, or smear, the interactions between charges and induced dipole at a short range involving small, highly charged atoms (such as hydrogen or lithium), thus preventing the "polarization catastrophe". The new force field gives stable trajectories and is validated through comparison with experimental data on density, viscosity, and ion diffusion coefficients of liquid systems of the above-mentioned classes. The results also shed light on the hydrogen-bonding pattern in ethylammonium nitrate, a PIL, for which the literature contains conflicting views. We describe the implementation of the TT damping function, of the temperature-grouped Nosé-Hoover thermostat for polarizable molecular dynamics (MD) and of the periodic perturbation method for viscosity evaluation from non-equilibrium trajectories in the LAMMPS MD code. The main result of this work is the wider applicability of the CL&Pol polarizable force field to new, important classes of fluids, achieving robust trajectories and a good description of equilibrium and transport properties in challenging systems. The fragment-based approach of CL&Pol will allow ready extension to a wide variety of PILs, DES, and electrolytes.

Multi-Interactions in Ionic Liquids for Natural Product Extraction

Natural products with a variety of pharmacological effects are important sources for commercial drugs, and it is very crucial to develop effective techniques to selectively extract and isolate bioactive natural components from the plants against the background of sustainable development. Ionic liquids (ILs) are a kind of designable material with unique physicochemical properties, including good thermal stability, negligible vapor pressure, good solvation ability, etc. ILs have already been used in pharmaceuticals for extraction, purification, drug delivery, etc. It has been reported that multi-interactions, like hydrogen bonding, hydrophobic interactions, play important roles in the extraction of bioactive components from the plants. In this review, recent progress in the understanding of scientific essence of hydrogen bonding, the special interaction, in ILs was summarized. The extraction of various natural products, one important area in pharmaceutical, by conventional and functional ILs as well as the specific roles of multi-interactions in this process were also reviewed. Moreover, problems existing in bioactive compound extraction by ILs and the future developing trends of this area are given, which might be helpful for scientists, especially beginners, in this field.

Theoretical Insights into the Structure of the Aminotris(Methylenephosphonic Acid) (ATMP) Anion: A Possible Partner for Conducting Ionic Media

We present a computational characterisation of Aminotris(methylenephosphonic acid) (ATMP) and its potential use as an anionic partner for conductive ionic liquids (ILs). We argue that for an IL to be a good candidate for a conducting medium, two conditions must be fulfilled: (i) the charge must be transported by light carriers; and (ii) the system must maintain a high degree of ionisation. The result trends presented herein show that there are molecular ion combinations that do comply with these two criteria, regardless of the specific system used. ATMP is a symmetric molecule with a total of six protons. In the bulk phase, breaking the symmetry of the fully protonated state and creating singly and doubly charged anions induces proton transfer mechanisms. To demonstrate this, we used molecular dynamics (MD) simulations employing a variable topology approach based on the reasonably reliable semiempirical density functional tight binding (DFTB) evaluation of the atomic forces. We show that, by choosing common and economical starting compounds, we can devise a viable prototype for a highly conductive medium where charge transfer is achieved by proton motion.

Ion Pair Structures and Hydrogen Bonding in RnNH4-n Alkylammonium Ionic Liquids with Hydrogen Sulfate and Mesylate Anions by DFT Computations

Density function theory calculations are employed to study the interaction of amines bearing different numbers of alkyl substituents of different sizes on the nitrogen atom with sulfuric and methanesulfonic acids. The proton affinities of the studied amines are calculated, and it is shown that the higher the value is, the more probable is its protonation. The most stable structures of the ion pairs resulting from the acid-base interaction are obtained and characterized. The geometric parameters of the ion pairs and the characteristics derived from the NBO and QTAIM analysis show that there are hydrogen bonding interactions between the cation and the anion. The hydrogen bonding character of the ion pairs and the strength of the interaction between the ions strongly depend on the nature of the cation itself. The interaction between the ions in the ion pairs weakens with the increase in the cation size. The trend of change in the structural parameters of the H-bonds and energetic characteristics in the cation series for the studied ion pairs is not dependent on the nature of the anion.

QM/MM simulations of organic phosphorus adsorption at the diaspore-water interface

Phosphorus (P) immobilization and thus its availability for plants are mainly affected by the strong interaction of phosphates with soil components especially soil mineral surfaces. The related reactions have been studied extensively via sorption experiments especially by carrying out adsorption of ortho-phosphates onto Fe-oxide surfaces. But a molecular-level understanding of the P-binding mechanisms at the mineral-water interface is still lacking, especially for forest eco-systems. Therefore, the current contribution provides an investigation of the molecular binding mechanisms for two abundant phosphates in forest soils, inositol hexaphosphate (IHP) and glycerolphosphate (GP), at the diaspore mineral surface. Here a hybrid electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) based molecular dynamics simulation has been applied to explore the diaspore-IHP/GP-water interactions. The results provide evidence for the formation of different P-diaspore binding motifs involving monodentate (M) and bidentate (B) for GP and two (2M) as well as three (3M) monodentates for IHP. The interaction energy results indicated the abundance of the GP B motif compared to the M one. The IHP 3M motif has a higher total interaction energy compared to its 2M motif, but exhibits a lower interaction energy per bond. Compared to GP, IHP exhibited stronger interaction with the surface as well as with water. Water was found to play an important role in controlling these diaspore-IHP/GP-water interactions. The interfacial water molecules form moderately strong H-bonds (HBs) with GP and IHP as well as with the diaspore surface. For all the diaspore-IHP/GP-water complexes, the interaction of water with the diaspore exceeds that with the studied phosphates. Furthermore, some water molecules form covalent bonds with diaspore Al atoms while others dissociate at the surface to protons and hydroxyl groups leading to proton transfer processes. Finally, the current results confirm the previous experimental conclusions indicating the importance of the number of phosphate groups, HBs, and proton transfers in controlling the P-binding at soil mineral surfaces.