Development of AMOEBA Force Field for 1,3-Dimethylimidazolium Based Ionic Liquids

Impact of an Ionic Liquid Solution on Horseradish Peroxidase Activity

Horseradish peroxidase (HRP) is an enzyme that oxidizes pollutants from wastewater. A previous report indicated that peroxidases can have an enhancement in initial enzymatic activity in an aqueous solution of 0.26 M 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) at neutral pH. However, the atomistic details remain elusive. In the enzymatic landscape of HRP, compound II (Cpd II) plays a key role and involves a histidine (H42) residue. Cpd II exists as oxoferryl (2a) or hydroxoferryl (2b(FeIV)) forms, where 2a is the predominantly observed form in experimental studies. Intriguingly, the ferric 2b(FeIII) form seen in synthetic complexes has not been observed in HRP. Here, we have investigated the structure and dynamics of HRP in pure water and aqueous [EMIm][EtSO4] (0.26 M), as well as the reaction mechanism of 2a to 2b conversion using polarizable molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations. When HRP is solvated in aq [EMIm][EtSO4], the catalytic water displaces, and H42 directly orients over the ferryl moiety, allowing a direct proton transfer (PT) with a significant energy barrier reduction. Conversely, in neat water, the reaction of 2a to 2b follows the previously reported mechanism. We further investigated the deprotonated form of H42. Analysis of the electric fields at the active site indicates that the aq [EMIm][EtSO4] medium facilitates the reaction by providing a more favorable environment compared with the system solvated in neat water. Overall, the atomic level supports the previous experimental observations and underscores the importance of favorable electric fields in the active site to promote catalysis.

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL-water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]⁺ cation and nitrate [NO₃]^- anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV^E, and liquid structures gr. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid-water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid-water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ∼0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. The ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

AMOEBA Polarizable Force Field for Molecular Dynamics Simulations of Glyme Solvents

Classical molecular dynamics (MD) simulations of electrolyte systems are important to gain insight into the atom-scale properties that determine the battery-relevant performance. The recent Tinker-HP software release enables efficient and accurate MD simulations with the AMOEBA polarizable force field. In this work, we developed a procedure to construct a universal AMOEBA model for the solvent family of glymes (glycol methyl ethers), which involves a refinement scheme for valence parameters by fitting the AMOEBA-derived atomic forces to those computed at the DFT level. The refined AMOEBA model provides a good description of both local and nonlocal properties in terms of the spectroscopic response of glyme molecules, as well as the liquid glyme density and dielectric constant. In addition, the complexation energies of alkali and alkaline-earth metal cations with tetraglyme molecules obtained from AMOEBA calculations are in good agreement with DFT results, demonstrating the suitability of the developed AMOEBA model for an accurate simulation of glyme-based battery electrolytes. We also expect the procedure to be transferable to the development of AMOEBA models for other battery electrolyte systems.

How Does Electronic Polarizability or Scaled-Charge Affect the Interfacial Properties of Room Temperature Ionic Liquids?

The room temperature ionic liquid (RTIL) air-liquid interface plays an important role in many applications. Herein, we present molecular dynamics simulation results for the air-liquid interface of a common RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) imide, [C4mim][NTf2]. To elucidate the effects of electronic polarizability and scaled-charge ions on the properties of the RTIL air-liquid interface, we employ three different kinds of force fields: a nonpolarizable force field (FF) with united ion charges (FixQ), a nonpolarizable FF with scaled-charge by 0.8 (ScaleQ), and a polarizable FF (Drude). To identify whether the ions reside at the interface or not, the method of identification of the truly interfacial molecules is used. The structural and dynamical properties in the interfacial, subinterfacial, and central layers are evaluated. In general for bulk liquids, the FixQ model predicts too-ordered structures and too-sluggish dynamics, while the ScaleQ model can serve as a simple cure. However, the ScaleQ model cannot reproduce the results of the Drude model at the interface, due to an inappropriate scaled-down charge near the interface.

Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions

Polarizable molecular dynamics simulations are a fast progressing field in the scientific research of ionic liquids. The fundamentals of polarizable simulations, as well as their application to ionic liquids, were summarized in a review [Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, Jr., A. D.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem. Rev. 2019, 119, 7940–7995] in 2019. Since then, new methods to treat intermolecular interaction of induced dipoles in these highly charged systems were developed. This concerns the damping of these interactions and additional charge transfer as well as the prediction of ionic materials with ultrahigh refractive indices. In addition to the progress of the polarizable force fields, also thermostats and barostats for polarizable simulations evolved recently.

Fine-Tuning the Polarizable CL&Pol Force Field for the Deep Eutectic Solvent Ethaline

Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular, for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable, and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low q-vector values (2.4 nm^-1) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nanoscale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against ab initio data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL&Pol force field. The results were also compared with the nonpolarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs.

Development of a Polarizable Force Field for Molecular Dynamics Simulations of Lithium-Ion Battery Electrolytes: Sulfone-Based Solvents and Lithium Salts

A many-body polarizable force field (PFF) was developed for molecular dynamics (MD) simulations of sulfone-based solvents and lithium salts. Development of the polarizable force field included parameterization of atomic polarizabilities, electrostatic interactions, and van der Waals interactions of electrolyte components. 1λ⁶-thiolane-1,1-dione or sulfolane (SLF) compound was selected as one of the most appropriate solvents for high-voltage battery electrolytes. Atomic polarizabilities for the sulfolane solvent and lithium salts were obtained by means of a combination of quantum mechanics (QM) and molecular mechanics (MM) approaches using the isotropic atomic dipole polarizable (IADP) model. High-quality atomic polarizabilities were refined for 10 atomic types. Intermolecular interactions of Li⁺ ions with SLF were parameterized to reproduce the binding energies at the MP2/aug-cc-pvDZ level of theory in the gas phase. Intermolecular interactions of Li⁺ ions with polyatomic anions, such as nitrate [NO₃]^-, tetrafluoroborate [BF₄]^-, perchlorate [ClO₄]^-, hexafluorophosphate [PF₆]^-, bis(fluorosulfonyl)imide [FSI]^-, and bis(trifluoromethylsulfonyl)imide [TFSI]^-, were parameterized employing a similar methodology. A series of molecular dynamics simulations was performed for sulfolane-based electrolytes at several different lithium salt concentrations. Thermodynamic, structural, and transport properties were evaluated to validate the force field parameters against available simulation and experimental data. Transport properties of sulfolane were significantly improved as compared with those obtained from MD simulations using a nonpolarizable force field (NFF). A newly developed polarizable potential was shown to reproduce Li⁺ ion dynamics as a function of salt concentration. Faster diffusion of Li⁺ ions, among other electrolyte components, was obtained for high salt concentration electrolytes.

Does Explicit Polarizability Improve Simulations of Phase Behavior of Ionic Liquids?

Molecular dynamics simulations are performed for a test set of 20 aprotic ionic liquids to investigate whether including an explicit polarizability model in the force field leads to higher accuracy and reliability of the calculated phase behavior properties, especially the enthalpy of fusion. A classical nonpolarizable all-atom optimized potentials for liquid simulations (OPLS) force-field model developed by Canongia Lopes and Pádua (CL&P) serves as a reference level of theory. Polarizability is included either in the form of Drude oscillators, resulting in the CL&P-D models, or in the framework of the atomic multipole optimized energetics for biomolecular application (AMOEBA) force field with polarizable atomic sites. Benchmarking of the calculated fusion enthalpy values against the experimental data reveals that overall the nonpolarizable CL&P model and polarizable CL&P-D models perform similarly with average deviations of about 30%. However, fusion enthalpies from the CL&P-D models exhibit a stronger correlation with their experimental counterparts. The least successful predictions are interestingly obtained from AMOEBA (deviation ca. 60%), which may indicate that a reparametrization of this force-field model is needed to achieve improved predictions of the fusion enthalpy. In general, all FF models tend to underestimate the fusion enthalpies. In addition, quantum chemical calculations are used to compute the electronic cohesive energies of the crystalline phases of the ionic liquids and of the interaction energies within the ion pair. Significant positive correlations are found between the fusion enthalpy and the cohesive energies. The character of the present anions predetermines the magnitude of individual mechanistic components of the interaction energy and related enthalpic and cohesive properties.

Determining the atomic charge of calcium ion requires the information of its coordination geometry in an EF-hand motif

It is challenging to parameterize the force field for calcium ions (Ca²⁺) in calcium-binding proteins because of their unique coordination chemistry that involves the surrounding atoms required for stability. In this work, we observed a wide variation in Ca²⁺ binding loop conformations of the Ca²⁺-binding protein calmodulin, which adopts the most populated ternary structures determined from the molecular dynamics simulations, followed by ab initio quantum mechanical (QM) calculations on all 12 amino acids in the loop that coordinate Ca²⁺ in aqueous solution. Ca²⁺ charges were derived by fitting to the electrostatic potential in the context of a classical or polarizable force field (PFF). We discovered that the atomic radius of Ca²⁺ in conventional force fields is too large for the QM calculation to capture the variation in the coordination geometry of Ca²⁺ in its ionic form, leading to unphysical charges. Specifically, we found that the fitted atomic charges of Ca²⁺ in the context of PFF depend on the coordinating geometry of electronegative atoms from the amino acids in the loop. Although nearby water molecules do not influence the atomic charge of Ca²⁺, they are crucial for compensating for the coordination of Ca²⁺ due to the conformational flexibility in the EF-hand loop. Our method advances the development of force fields for metal ions and protein binding sites in dynamic environments.

Systematic Parameterization and Simulation of Boronic Acid-β-Lactamase Aqueous Solution in Developing the ABEEMσπ Polarizable Force Field

Boronic acid, an inhibitor of β-lactamase, has begun to be applied to the treatment of biological infections and tumors. Scientists are working to develop new and more effective boronic acid. Molecular dynamics (MD) simulation provides a powerful auxiliary tool for drug design. However, the current force fields have no boron-related parameters. In this work, an atom-bond electronegativity equalization method at the σπ level (ABEEMσπ) polarizable force field (ABEEMσπ PFF) of boronic acid and β-lactamase has been developed to determine the potential functions and parameters. The interaction between boron and serine in β-lactamase is regarded as a bonded mode. The interaction between them is simulated by the Morse potential energy function, which is close to the experimental change of the stretching potential energy in a large range. The potential energy surfaces of the bond length, bond angle, and dihedral angle of boronic acid-β-lactamase have the same stability point and change trend as M06-2X/6-311G**. For 47 boronic acid-β-lactamase training molecules, the linear correlation coefficient (R) of the charge distribution between the ABEEMσπ PFF and HF/STO-3G is greater than 0.96. Attributed to the fact that the charge distribution of the ABEEMσπ PFF can fluctuate with the change of geometry and environment, the polarization effect and charge-transfer effect are well reflected. The binding ability of different boronic acids with the same β-lactamase is different. A total of 10 boronic acid-β-lactamase model molecules and 10 boronic acid-β-lactamase and water complexes are simulated. The order of binding energy of five large model molecules calculated by the ABEEMσπ PFF is consistent with that of the MP2 method. The binding energies of boronic acid-β-lactamase and water complexes are close to those of the MP2 method. The results of MD simulation of five aqueous boronic acid-β-lactamase complexes in the NVT ensemble verify the rationality of boron-related parameters of the ABEEMσπ PFF, which have a good application prospect. This study lays a solid theoretical foundation for further study of the inhibition of boronic acid on β-lactamase.

Current Status of AMOEBA-IL: A Multipolar/Polarizable Force Field for Ionic Liquids

Computational simulations of ionic liquid solutions have become a useful tool to investigate various physical, chemical and catalytic properties of systems involving these solvents. Classical molecular dynamics and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of IL systems have provided significant insights at the atomic level. Here, we present a review of the development and application of the multipolar and polarizable force field AMOEBA for ionic liquid systems, termed AMOEBA–IL. The parametrization approach for AMOEBA–IL relies on the reproduction of total quantum mechanical (QM) intermolecular interaction energies and QM energy decomposition analysis. This approach has been used to develop parameters for imidazolium– and pyrrolidinium–based ILs coupled with various inorganic anions. AMOEBA–IL has been used to investigate and predict the properties of a variety of systems including neat ILs and IL mixtures, water exchange reactions on lanthanide ions in IL mixtures, IL–based liquid–liquid extraction, and effects of ILs on an aniline protection reaction.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Polarizable ab initio QM/MM Study of the Reaction Mechanism of N-tert-Butyloxycarbonylation of Aniline in [EMIm][BF₄]

N-t e r t-butoxycarbonylation of amines in solution (water, organic solvents, or ionic liquids) is a common reaction for the preparation of drug molecules. To understand the reaction mechanism and the role of the solvent, quantum mechanical/molecular mechanical simulations using a polarizable multipolar force field with long⁻range electrostatic corrections were used to optimize the minimum energy paths (MEPs) associated with various possible reaction mechanisms employing the nudged elastic band (NEB) and the quadratic string method (QSM). The calculated reaction energies and energy barriers were compared with the corresponding gas-phase and dichloromethane results. Complementary Electron Localization Function (ELF)/NCI analyses provide insights on the critical structures along the MEP. The calculated results suggest the most likely path involves a sequential mechanism with the rate⁻limiting step corresponding to the nucleophilic attack of the aniline, followed by proton transfer and the release of CO 2 without the direct involvement of imidazolium cations as catalysts.

Molecular dynamics investigation of water-exchange reactions on lanthanide ions in water/1-ethyl-3-methylimidazolium trifluoromethylsufate ([EMIm][OTf])

We report a kinetic study of the water exchange on lanthanide ions in water/[1-ethyl-3-methylimidazolium][trifluoromethylsufate] (water/[EMIm][OTf]). The results from ¹⁷O-NMR measurements show that the water-exchange rates in water/[EMIm][OTf] increase with decreasing size of the lanthanide ions. This trend for water-exchange is similar to the previously reported trend in water/1-ethyl-3-methylimidazolium ethyl sulfate (water/[EMIm][EtSO₄]) but opposite to that in water. To gain atomic-level insight into these water-exchange reactions, molecular dynamics simulations for lanthanide ions in water/[EMIm][OTf] have been performed using the atomic-multipole-optimized-energetics-for-biomolecular-application polarizable force field. Our molecular dynamics simulations reproduce the experimental water-exchange rates in terms of the trend and provide possible explanations for the observed experimental behavior. The smaller lanthanide ions in water/[EMIm][OTf] undergo faster water exchange because the smaller lanthanide ions coordinate to the first shell [OTf]^- anions more tightly, resulting in a stronger screening effect for the second-shell water. The screening effect weakens the interaction of the lanthanide ions with the second-shell water molecules, facilitating the dissociation of water from the second-shell and subsequent association of water molecules from the outer solvation shells.

Tinker-HP: a massively parallel molecular dynamics package for multiscale simulations of large complex systems with advanced point dipole polarizable force fields

We present Tinker-HP, a massively MPI parallel package dedicated to classical molecular dynamics (MD) and to multiscale simulations, using advanced polarizable force fields (PFF) encompassing distributed multipoles electrostatics. Tinker-HP is an evolution of the popular Tinker package code that conserves its simplicity of use and its reference double precision implementation for CPUs. Grounded on interdisciplinary efforts with applied mathematics, Tinker-HP allows for long polarizable MD simulations on large systems up to millions of atoms. We detail in the paper the newly developed extension of massively parallel 3D spatial decomposition to point dipole polarizable models as well as their coupling to efficient Krylov iterative and non-iterative polarization solvers. The design of the code allows the use of various computer systems ranging from laboratory workstations to modern petascale supercomputers with thousands of cores. Tinker-HP proposes therefore the first high-performance scalable CPU computing environment for the development of next generation point dipole PFFs and for production simulations. Strategies linking Tinker-HP to Quantum Mechanics (QM) in the framework of multiscale polarizable self-consistent QM/MD simulations are also provided. The possibilities, performances and scalability of the software are demonstrated via benchmarks calculations using the polarizable AMOEBA force field on systems ranging from large water boxes of increasing size and ionic liquids to (very) large biosystems encompassing several proteins as well as the complete satellite tobacco mosaic virus and ribosome structures. For small systems, Tinker-HP appears to be competitive with the Tinker-OpenMM GPU implementation of Tinker. As the system size grows, Tinker-HP remains operational thanks to its access to distributed memory and takes advantage of its new algorithmic enabling for stable long timescale polarizable simulations. Overall, a several thousand-fold acceleration over a single-core computation is observed for the largest systems. The extension of the present CPU implementation of Tinker-HP to other computational platforms is discussed.

Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

Computational investigation of O2 diffusion through an intra-molecular tunnel in AlkB; influence of polarization on O2 transport

E. Coli AlkB catalyzes the direct dealkylation of various alkylated bases in damaged DNA. The diffusion of molecular oxygen to the active site in AlkB is an essential step for the oxidative dealkylation activity. Despite detailed studies on the stepwise oxidation mechanism of AlkB, there is no conclusive picture of how O2 molecules reach the active site of the protein. Yu et al. (Nature, 439, 879) proposed the existence of an intra-molecular tunnel based on their initial crystal structures of AlkB. We have employed computational simulations to investigate possible migration pathways inside AlkB for O2 molecules. Extensive molecular dynamics (MD) simulations, including explicit ligand sampling and potential of mean force (PMF) calculations, have been performed to provide a microscopic description of the O2 delivery pathway in AlkB. Analysis of intra-molecular tunnels using the CAVER software indicates two possible pathways for O2 to diffuse into the AlkB active site. Explicit ligand sampling simulations suggests that only one of these tunnels provides a viable route. The free energy path for an oxygen molecule to travel along each of these tunnels has been determined with AMBER and AMOEBA. Both PMFs indicate passive transport of O2 from the surface of the protein. However, the inclusion of explicit polarization shows a very large barrier for diffusion of the co-substrate out of the active site, compared with the non-polarizable potential. In addition, our results suggest that the mutation of a conserved residue along the tunnel, Y178, has dramatic effects on the dynamics of AlkB and on the transport of O2 along the tunnel.

Toward individually tunable compound eyes with transparent graphene electrode

We present tunable compound eyes made of ionic liquid lenses, of which both curvatures (R ₁ and R ₂ in the lensmaker's equation) can be individually changed using electrowetting on dielectric (EWOD) and applied pressure. Flexible graphene is used as a transparent electrode and is integrated on a flexible polydimethylsiloxane (PDMS)/parylene hybrid substrate. Graphene electrodes allow a large lens aperture diameter of between 2.4 mm and 2.74 mm. Spherical aberration analysis is performed using COMSOL to investigate the optical property of the lens under applied voltage and pressure. The final lens system shows a resolution of 645.1 line pair per millimeter. A prototype of a tunable lens array is proposed for the application of a compound eye.

Induced polarization restricts the conformational distribution of a light-harvesting molecular triad in the ground state

Free energy surface of the light-harvesting triad employing a non-polarizable force field (NFF) and a polarizable force field (PFF) shows that induced polarization limits the motion of rotation about chemical bonds as well as bending at the porphyrin, which are prominent using the NFF, thus limiting the conformational space of the triad.

Truncated Conjugate Gradient: An Optimal Strategy for the Analytical Evaluation of the Many-Body Polarization Energy and Forces in Molecular Simulations

We introduce a new class of methods, denoted as Truncated Conjugate Gradient(TCG), to solve the many-body polarization energy and its associated forces in molecular simulations (i.e. molecular dynamics (MD) and Monte Carlo). The method consists in a fixed number of Conjugate Gradient (CG) iterations. TCG approaches provide a scalable solution to the polarization problem at a user-chosen cost and a corresponding optimal accuracy. The optimality of the CG-method guarantees that the number of the required matrix-vector products are reduced to a minimum compared to other iterative methods. This family of methods is non-empirical, fully adaptive, and provides analytical gradients, avoiding therefore any energy drift in MD as compared to popular iterative solvers. Besides speed, one great advantage of this class of approximate methods is that their accuracy is systematically improvable. Indeed, as the CG-method is a Krylov subspace method, the associated error is monotonically reduced at each iteration. On top of that, two improvements can be proposed at virtually no cost: (i) the use of preconditioners can be employed, which leads to the Truncated Preconditioned Conjugate Gradient (TPCG); (ii) since the residual of the final step of the CG-method is available, one additional Picard fixed point iteration (“peek”), equivalent to one step of Jacobi Over Relaxation (JOR) with relaxation parameter ω, can be made at almost no cost. This method is denoted by TCG-n(ω). Black-box adaptive methods to find good choices of ω are provided and discussed. Results show that TPCG-3(ω) is converged to high accuracy (a few kcal/mol) for various types of systems including proteins and highly charged systems at the fixed cost of four matrix-vector products: three CG iterations plus the initial CG descent direction. Alternatively, T(P)CG-2(ω) provides robust results at a reduced cost (three matrix-vector products) and offers new perspectives for long polarizable MD as a production algorithm. The T(P)CG-1(ω) level provides less accurate solutions for inhomogeneous systems, but its applicability to well-conditioned problems such as water is remarkable, with only two matrix-vector product evaluations.

General van der Waals potential for common organic molecules

This work presents a systematic development of a new van der Waals potential (vdW2016) for common organic molecules based on symmetry-adapted perturbation theory (SAPT) energy decomposition. The Buf-14-7 function, as well as Cubic-mean and Waldman-Hagler mixing rules were chosen given their best performance among other popular potentials. A database containing 39 organic molecules and 108 dimers was utilized to derive a general set of vdW parameters, which were further validated on nucleobase stacking systems and testing organic dimers. The vdW2016 potential is anticipated to significantly improve the accuracy and transferability of new generations of force fields for organic molecules.

Simulations of the water exchange dynamics of lanthanide ions in 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) and water

Molecular dynamics simulations have been carried out to explain the water-exchange rates of lanthanide ions in water and water/[EMIm][EtSO₄] observed from ¹⁷O-NMR experiments. Our simulations are in agreement with experimental results with respect to water-exchange trends.

Computational approaches to understanding reaction outcomes of organic processes in ionic liquids

The utility of using a combined experimental and computational approach for understanding ionic liquid media, and their effect on reaction outcome, is highlighted through a number of case studies.

Observation of aliphatic C–H⋯X hydrogen bonds in imidazolium ionic liquids

Measurement of H/D isotope effects on the¹⁹F resonances of sequentially deuterated 1-n-butyl-3-methylimidazolium hexafluorophosphate and tetrafluoroborate isotopologues reveals the formation of aliphatic C–H⋯X hydrogen bonds between the fluorinated anions and protons along theN-alkyl sidechains of the cations.

Simulations of room temperature ionic liquids: from polarizable to coarse-grained force fields

This perspective article summarizes the recent advances in the classical molecular modelling of room temperature ionic liquids.