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

J Phys Chem B. 2022 Jun 23;126(24):4611-4622. doi: 10.1021/acs.jpcb.2c02782. Epub 2022 Jun 13.

Abstract

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.