Jesse G. McDaniel

Physically Motivated, Robust, ab Initio Force Fields for CO2 and N2

We present a novel methodology for developing physically motivated, first-principles polarizable force fields and apply these techniques to the specific cases of CO(2) and N(2). Exchange, electrostatic, induction, and dispersion interaction parameters were fit to symmetry adapted perturbation theory (SAPT) dimer energy calculations, with explicit terms to account for each of the dominant fundamental interactions between molecules. Each term is represented by a physically appropriate functional form and fitted individually based on the results of the SAPT decomposition. The resulting CO(2) and N(2) force field was benchmarked against a diverse set of experimental data, including the second virial coefficient, density, scattering structure factor, heat capacity, enthalpy of vaporization, and vapor-liquid coexistence curves. In general, excellent agreement with experimental data is obtained with our model. Due to the physical nature of their construction, these force fields are robust and transferable to environments for which they were not specifically parametrized, including gas mixtures, and we anticipate applications in modeling CO(2)/N(2) adsorption in polar and/or heterogeneous media.

First-Principles, Physically Motivated Force Field for the Ionic Liquid [BMIM][BF4]

Molecular simulations play an important role in establishing structure-property relations in complex fluids such as room-temperature ionic liquids. Classical force fields are the starting point when large systems or long times are of interest. These force fields must be not only accurate but also transferable. In this work, we report a physically motivated force field for the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]) based on symmetry-adapted perturbation theory. The predictions (from molecular dynamics simulations) of the liquid density, enthalpy of vaporization, diffusion coefficients, viscosity, and conductivity are in excellent agreement with experiment, with no adjustable parameters. The explicit energy decomposition inherent in the force field enables a quantitative analysis of the important physical interactions in these systems. We find that polarization is crucial and there is little evidence of charge transfer. We also argue that the often used procedure of scaling down charges in molecular simulations of ionic liquids is unphysical for [BMIM][BF4]. Because all intermolecular interactions in the force field are parametrized from first-principles, we anticipate good transferability to other ionic liquid systems and physical conditions.

Next-Generation Force Fields from Symmetry-Adapted Perturbation Theory

Symmetry-adapted perturbation theory (SAPT) provides a unique set of advantages for parameterizing next-generation force fields from first principles. SAPT provides a direct, basis-set superposition error free estimate of molecular interaction energies, a physically intuitive energy decomposition, and a seamless transition to an asymptotic picture of intermolecular interactions. These properties have been exploited throughout the literature to develop next-generation force fields for a variety of applications, including classical molecular dynamics simulations, crystal structure prediction, and quantum dynamics/spectroscopy. This review provides a brief overview of the formalism and theory of SAPT, along with a practical discussion of the various methodologies utilized to parameterize force fields from SAPT calculations. It also highlights a number of applications of SAPT-based force fields for chemical systems of particular interest. Finally, the review ends with a brief outlook on the future opportunities and challenges that remain for next-generation force fields based on SAPT.

Ab Initio Force Fields for Imidazolium-Based Ionic Liquids

We develop ab initio force fields for alkylimidazolium-based ionic liquids (ILs) that predict the density, heats of vaporization, diffusion, and conductivity that are in semiquantitative agreement with experimental data. These predictions are useful in light of the scarcity of and sometimes inconsistency in experimental heats of vaporization and diffusion coefficients. We illuminate physical trends in the liquid cohesive energy with cation chain length and anion. These trends are different than those based on the experimental heats of vaporization. Molecular dynamics prediction of the room temperature dynamics of such ILs is more difficult than is generally realized in the literature due to large statistical uncertainties and sensitivity to subtle force field details. We believe that our developed force fields will be useful for correctly determining the physics responsible for the structure/property relationships in neat ILs.

An efficient multi-scale lattice model approach to screening nano-porous adsorbents

We present a multi-scale, hierarchical, approach for developing lattice models to estimate adsorption in nano-porous sorbents, derived on the basis of underlying atomistic potentials. This approach is a generalization of earlier work in zeolites (where the specific adsorption sites are easily definable) to encompass both specific as well as diffuse adsorption; the latter often dominates in the case of nano-porous metal-organic frameworks (MOFs). In conjunction with appropriately coarse grained guest-guest interactions, we demonstrate that our lattice approach offers semi-quantitative to quantitative agreement as compared to fully atomistic simulation from the low pressure regime through saturation. However, it also yields orders-of-magnitude acceleration versus the latter, thus enabling high-throughput screenings of both non-polar and polar adsorbates with high efficiency. We also show how our lattice model can be extended to facilitate rapid, qualitative screening of transport properties via appropriate calibration. Although our example applications focus on CO(2) adsorption in MOFs, this approach is readily generalizable to various nano-porous materials (MOFs, zeolites...) and guest adsorbates (CO(2), H(2), hydrocarbons).

Comment on “Isolating the non-polar contributions to the intermolecular potential for water-alkane interactions” [J. Chem. Phys. 141, 064905 (2014)]

The manuscript by Ballal et al.(Ref 1) presents an interesting study demonstrating the inability of popular force fields with standard combination rules to accurately describe water/alkane interactions. The authors find that the Lorentz-Berthelot combination rules on the SPC/E water and TraPPE alkane potentials give a cross interaction that fails to predict the (low-water content) water solubility in various alkanes. Realizing that both explicit polarization as well as the static octupole moment of methane are missing in these potentials, the authors examine the effect of these terms, but are still unable to resolve the discrepancy. They conclude with the statement that “the research community lacks a complete picture of water-alkane interactions at the molecular level.

Conformational and Dynamic Properties of Poly(ethylene oxide) in an Ionic Liquid: Development and Implementation of a First-Principles Force Field.

The conformational properties of polymers in ionic liquids are of fundamental interest but not well understood. Atomistic and coarse-grained molecular models predict qualitatively different results for the scaling of chain size with molecular weight, and experiments on dilute solutions are not available. In this work, we develop a first-principles force field for poly(ethylene oxide) (PEO) in the ionic liquid 1-butyl 3-methylimidazolium tetrafluoroborate ([BMIM][BF4]) using symmetry adapted perturbation theory (SAPT). At temperatures above 400 K, simulations employing both the SAPT and OPLS-AA force fields predict that PEO displays ideal chain behavior, in contrast to previous simulations at lower temperature. We therefore argue that the system shows a transition from extended to more compact configurations as the temperature is increased from room temperature to the experimental lower critical solution temperature. Although polarization is shown to be important, its implicit inclusion in the OPLS-AA force is sufficient to describe the structure and energetics of the mixture. The simulations emphasize the difference between ionic liquids from typical solvents for polymers.

Dynamics of Water in Gemini Surfactant-Based Lyotropic Liquid Crystals

The dynamics of water confined to nanometer-sized domains is important in a variety of applications ranging from proton exchange membranes to crowding effects in biophysics. In this work, we study the dynamics of water in gemini surfactant-based lyotropic liquid crystals (LLCs) using molecular dynamics simulations. These systems have well characterized morphologies, for example, hexagonal, gyroid, and lamellar, and the surfaces of the confining regions can be controlled by modifying the headgroup of the surfactants. This allows one to study the effect of topology, functionalization, and interfacial curvature on the dynamics of confined water. Through analysis of the translational diffusion and rotational relaxation, we conclude that the hydration level and resulting confinement length scale is the predominate determiner of the rates of water dynamics, and other effects, namely, surface functionality and curvature, are largely secondary. This novel analysis of the water dynamics in these LLC systems provides an important comparison for previous studies of water dynamics in lipid bilayers and reverse micelles.

Importance of hydrophobic traps for proton diffusion in lyotropic liquid crystals

The diffusion of protons in self-assembled systems is potentially important for the design of efficient proton exchange membranes. In this work, we study proton dynamics in a low-water content, lamellar phase of a sodium-carboxylate gemini surfactant/water system using computer simulations. The hopping of protons via the Grotthuss mechanism is explicitly allowed through the multi-state empirical valence bond method. We find that the hydronium ion is trapped on the hydrophobic side of the surfactant-water interface, and proton diffusion then proceeds by hopping between surface sites. The importance of hydrophobic traps is surprising because one would expect the hydronium ions to be trapped at the charged headgroups. The physics illustrated in this system should be relevant to the proton dynamics in other amphiphilic membrane systems, whenever there exist exposed hydrophobic surface regions.

Ion Correlation and Collective Dynamics in BMIM/BF4 Based Organic Electrolytes: From Dilute Solutions to the Ionic Liquid Limit

Quantifying ion association and collective dynamical processes in organic electrolytes is essential for fundamental property interpretation and optimization for electrochemical applications. The extent of ion correlation depends on both the ion concentration and dielectric strength of the solvent; ions may be largely uncorrelated in sufficiently high-dielectric solvents at low concentration, but properties of concentrated electrolytes are dictated by correlated and collective ion processes. In this work, we utilize molecular dynamics simulations to characterize ion association and collective ion dynamics in organic electrolytes composed of binary mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate [BMIM+][BF4-] and 1,2-dichloroethane, acetone, acetonitrile, and water solvents. We illustrate different physical regimes of characteristically distinct ion correlations for the systematic range of electrolyte concentrations and solvent dielectric strengths. Dilute electrolytes composed of low-dielectric solvents exhibit significant counterion correlation in the form of ion pairing and clustering driven by both weak screening and relatively low solvation energies. This regime is characterized by enhanced ion coordination numbers and near equality of cation and anion diffusion coefficients, despite the significantly different ion sizes. In contrast, ion correlation in highly concentrated electrolytes is dominated by the anti-correlated motion of both like-charge and opposite-charge ions, approaching neat ionic liquid behavior. We show that the cross-over of these correlation regimes is clearly illuminated by quantifying the fractional self and distinct contributions to the net ionic conductivity. For organic electrolytes composed of low-dielectric solvents, we conclude that significant ion correlation exists at all concentrations but the nature of the correlation changes markedly from the dilute electrolyte to the pure ionic liquid limit.