Self-Consistent Scheme Combining MD and Order-N DFT Methods: An Improved Set of Nonpolarizable Force Fields for Ionic Liquids

Amide cyclodextrin that recognises monophosphate anions in harmony with water molecules

Anion recognition in water by synthetic host molecules is a popular and challenging topic. It has been considered difficult because the water molecules compete for the recognition units. In this study, we have successfully created a novel macrocycle that achieves precise recognition through multipoint hydrogen bonding in harmony with water molecules. Specifically, an N-methylpyridinium amide β-cyclodextrin (β-CD) derivative 1(OTf)7 was synthesized, whose amide groups are directly attached to each pyranose ring. The pyridinium amide CD encapsulated a monophosphate anion in water, but it did not show interactions with sulfonates or carboxylates, thus a remarkable selectivity was demonstrated. Two monophosphates with different substituents, phenyl phosphate (PhOPO3 2-) and adamantyl phosphate (AdOPO3 2-), exhibited interesting contrasting pictures in the inclusion process, which were revealed by a combination of NOESY experiments, ITC measurements, and MD simulations. PhOPO3 2- was positioned slightly "upper" (closer to the pyridinium amide side) in 17+ with the oxygen atom of the phosphate ester R-O-P involved in the hydrogen bonds with the amide N-H, and configurational entropy plays a key role in the inclusion. Meanwhile, AdOPO3 2- was positioned "lower" (closer to the methoxy rim of CD) with the terminal -PO3 2- forming hydrogen bonds with the amides, and the hydrophobic effect is a major contributing driving force of the inclusion. The molecular design presented herein to achieve the precise recognition in water and clarification of the detailed mechanisms including the hydration phenomenon greatly contribute to the development of functional molecules that work in aqueous environments.

Aquatic Functional Liquid Crystals: Design, Functionalization, and Molecular Simulation

Aquatic functional liquid crystals, which are ordered molecular assemblies that work in water environment, are described in this review. Aquatic functional liquid crystals are liquid-crystalline (LC) materials interacting water molecules or aquatic environment. They include aquatic lyotropic liquid crystals and LC based materials that have aquatic interfaces, for example, nanoporous water treatment membranes that are solids preserving LC order. They can remove ions and viruses with nano- and subnano-porous structures. Columnar, smectic, bicontinuous LC structures are used for fabrication of these 1D, 2D, 3D materials. Design and functionalization of aquatic LC sensors based on aqueous/LC interfaces are also described. The ordering transitions of liquid crystals induced by molecular recognition at the aqueous interfaces provide distinct optical responses. Molecular orientation and dynamic behavior of these aquatic functional LC materials are studied by molecular dynamics simulations. The molecular interactions of LC materials and water are key of these investigations. New insights into aquatic functional LC materials contribute to the fields of environment, healthcare, and biotechnology.

Treating Polarization Effects in Charged and Polar Bio-Molecules Through Variable Electrostatic Parameters

Polarization plays important roles in charged and hydrogen bonding containing systems. Much effort ranging from the construction of physics-based models to quantum mechanism (QM)-based and machine learning (ML)-assisted models have been devoted to incorporating the polarization effect into the conventional force fields at different levels, such as atomic and coarse grained (CG). The application of polarizable force fields or polarization models was limited by two aspects, namely, computational cost and transferability. Different from physics-based models, no predetermining parameters were required in the QM-based approaches. Taking advantage of both the accuracy of QM calculations and efficiency of molecular mechanism (MM) and ML, polarization effects could be treated more efficiently while maintaining the QM accuracy. The computational cost could be reduced with variable electrostatic parameters, such as the charge, dipole, and electronic dielectric constant with the help of linear scaling fragmentation-based QM calculations and ML models. Polarization and entropy effects on the prediction of partition coefficient of druglike molecules are demonstrated by using both explicit or implicit all-atom molecular dynamics simulations and machine learning-assisted models. Directions and challenges for future development are also envisioned.

Nonpolarizable Force Fields through the Self-Consistent Modeling Scheme with MD and DFT Methods: From Ionic Liquids to Self-Assembled Ionic Liquid Crystals

A key to achieve the accuracy of molecular dynamics (MD) simulation is the set of force fields used to express the atomistic interactions. In particular, the electrostatic interaction remains the main issue for the precise simulation of various ionic soft materials from ionic liquids to their supramolecular compounds. In this study, we test the nonpolarizable force fields of ionic liquids (ILs) and self-assembled ionic liquid crystals (ILCs) for which the intermolecular charge transfer and intramolecular polarization are significant. The self-consistent modeling scheme is adopted to refine the atomic charges of ionic species in a condensed state through the use of density functional theory (DFT) under the periodic boundary condition. The atomic charges of the generalized amber force field (GAFF) are effectively updated to express the electrostatic properties of ionic molecules obtained by the DFT calculation in condensed phase, which improves the prediction accuracy of ionic conductivity with the obtained force field (GAFF-DFT). The derived DFT charges then suggest that the substitution of a hydrophobic liquid-crystalline moiety into IL-based cations enhances the charge localization of ionic groups in the amphiphilic molecules, leading to the amplification of the electrostatic interactions among the hydrophilic/ionic groups in the presence of hydrophobic moieties. In addition, we focus on an ion-conductive pathway hidden in the self-assembled nanostructure. The MD results indicate that the ionic groups of cation and anion interact strongly for keeping the bicontinuous nanosegregation of ionic nanochannel. The partial fractions of hydrophilic/ionic and hydrophobic nanodomains are then quantified with the volume difference from referenced IL systems, while the calculated ionic conductivity decreases in the self-assembled ILCs more than the occupied volume of ionic nanodomains. These analyses suggest that the mobility of ions in the self-assembled ILCs remains quite restricted even with small tetrafluoroborate anions because of strong attractive interaction among ionic moieties.

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.

Establishing the accuracy of density functional approaches for the description of noncovalent interactions in ionic liquids

We test a number of dispersion corrected versatile Generalized Gradient Approximation (GGA) and meta-GGA functionals for their ability to predict the interactions of ionic liquids, and show that most can achieve energies within 1 kcal mol^-1 of benchmarks. This compares favorably with an accurate dispersion corrected hybrid, ωB97X-V. Our tests also reveal that PBE (Perdew-Burke-Ernzerhof GGA) calculations using the plane-wave projector augmented wave method and Gaussian Type Orbitals (GTOs) differ by less than 0.6 kJ mol^-1 for ionic liquids, despite ions being difficult to evaluate in periodic cells - thus revealing that GTO benchmarks may be used also for plane-wave codes. Finally, the relatively high success of explicit van der Waals density functionals, compared to elemental and ionic dispersion models, suggests that improvements are required for low-cost dispersion correction models of ions.

Molecular insights on confined water in the nanochannels of self-assembled ionic liquid crystal

Self-assembled ionic liquid crystals can transport water and ions via the periodic nanochannels, and these materials are promising candidates as water treatment membranes. Molecular insights on the water transport process are, however, less investigated because of computational difficulties of ionic soft matters and the self-assembly. Here we report specific behavior of water molecules in the nanochannels by using the self-consistent modeling combining density functional theory and molecular dynamics and the large-scale molecular dynamics calculation. The simulations clearly provide the one-dimensional (1D) and 3D-interconnected nanochannels of self-assembled columnar and bicontinuous structures, respectively, with the precise mesoscale order observed by x-ray diffraction measurement. Water molecules are then confined inside the nanochannels with the formation of hydrogen bonding network. The quantitative analyses of free energetics and anisotropic diffusivity reveal that, the mesoscale geometry of 1D nanodomain profits the nature of water transport via advantages of dissolution and diffusion mechanisms inside the ionic nanochannels.

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.

Spatial-Decomposition Analysis of Electrical Conductivity in Mixtures of Ionic Liquid and Sodium Salt for Sodium-Ion Battery Electrolytes

Ionic liquid (IL)-based electrolytes are a promising material for the development of sodium-ion batteries, and their performance can be quantified by electrical conductivity. In this highly concentrated ionic system, the correlated motions of ion pairs are influential on the ionic transport properties. Herein, all-atom analyses are conducted through molecular dynamics simulations to bridge the macroscopically observable electrical conductivity with the molecular pictures of correlated motion of ion pairs. The analysis is applied to three mixtures of IL with sodium salt that are relevant to the electrolyte for a sodium-ion battery: [1-ethyl-3-methylimidazolium, Na][bis(fluorosulfonyl)amide] ([C₂C₁im, Na][FSA]), [N-methyl-N-propylpyrrolidinium, Na][FSA] ([C₃C₁pyrr, Na][FSA]), and [K, Na][FSA]. The computational results on electrical conductivities are in agreement with the experimental reports, and their dependency on temperature and sodium-ion composition is reproduced well. The overall contributions from cross-correlated motions are found to be negative in all the IL mixtures; thus, the total conductivities are less than their Nernst-Einstein estimates. The spatial view of cross-correlated motions is further obtained by decomposing the time correlation functions of velocities according to the distances between ion pairs. It is observed that ion pairs are moving in the same direction for ∼0.3 ps when they were initially within the first coordination shell, followed by motions toward opposite directions. The cross-correlation terms are also dissected into local and nonlocal components, and it is commonly seen for all the ion pairs that the local component is negative for cation-anion pairs and is positive for cation-cation and anion-anion pairs. The motions of ion pairs are accompanied by a "backflow" that manifests in the form of the nonlocal component whose sign is opposite to the corresponding local component. In fact, the contributions of the correlated motions of ions to the electrical conductivity are not localized to contact pairs and extend spatially beyond the first coordination shell of the cation-anion pairs.