Polarizable Model Research Articles

Advances in Molecular Modeling of Ion-Protein Interaction Systems Towards Accurate Electrostatics: Methods and Applications

Metal ions are ubiquitous in complex with proteins and play key roles in protein structure and function. The ion-protein interactions are electrostatic delicate in nature, however, the description of electrostatic interactions could be problematic in conventional additive fixed-charge force fields. With empowered computational sources going beyond the common approximations, many efforts have been done to take account in more elaborate electrostatic description in molecular modeling using more sophisticated physical models and dilated algorithms/implementations. Here we review rencent progress in advanced polarizable models and new impletment approaches towards accurate electrostatics and highlight some successful application cases in ion-protein interaction systems in recent years.

Computational and Experimental Characterization of rDNA and rRNA G-Quadruplexes.

DNA G-quadruplexes in human telomeres and gene promoters are being extensively studied for their role in controlling the growth of cancer cells. G-quadruplexes have been unambiguously shown to exist both in vitro and in vivo, including in the guanine (G)-rich DNA genes encoding pre-ribosomal RNA (pre-rRNA), which is transcribed in the cell's nucleolus. Recent studies strongly suggest that these DNA sequences ("rDNA"), and the transcribed rRNA, are a potential anticancer target through the inhibition of RNA polymerase I (Pol I) in ribosome biogenesis, but the structures of ribosomal G-quadruplexes at atomic resolution are unknown and very little biophysical characterization has been performed on them to date. In the present study, circular dichroism (CD) spectroscopy is used to show that two putative rDNA G-quadruplex sequences, NUC 19P and NUC 23P and their counterpart rRNAs, predominantly adopt parallel topologies, reminiscent of the analogous telomeric quadruplex structures. Based on this information, we modeled parallel topology atomistic structures of the putative ribosomal G-quadruplexes. We then validated and refined the modeled ribosomal G-quadruplex structures using all-atom molecular dynamics (MD) simulations with the CHARMM36 force field in the presence and absence of stabilizing K+. Motivated by preliminary MD simulations of the telomeric parallel G-quadruplex (TEL 24P) in which the K+ ion is expelled, we used updated CHARMM36 force field K+ parameters that were optimized, targeting the data from quantum mechanical calculations and the polarizable Drude model force field. In subsequent MD simulations with optimized CHARMM36 parameters, the K+ ions are predominantly in the G-quadruplex channel and the rDNA G-quadruplexes have more well-defined, predominantly parallel-topology structures as compared to rRNA. In addition, NUC 19P is more structured than NUC 23P, which contains extended loops. Results from this study set the structural foundation for understanding G-quadruplex functions and the design of novel chemotherapeutics against these nucleolar targets and can be readily extended to other DNA and RNA G-quadruplexes.

OS100: A Benchmark Set of 100 Digitized UV-Visible Spectra and Derived Experimental Oscillator Strengths.

Excited-state quantum chemical calculations usually report excitation energies and oscillator strengths, f, for each electronic transition. On the other hand, UV-visible spectrophotometric experiments measure energy-dependent molar extinction/attenuation coefficients, ε(v), that give absorption band line shapes when plotted. ε(v) and f are related, but this relation is complicated by broadening and solvation effects. We fitted and integrated 100 experimental UV-visible spectra to obtain 164 fexp values for absorption bands appearing in these spectra. The 100 UV-visible spectra belong to solvated organic molecules ranging in size from 6-34 atoms. We estimated uncertainties in the fitting to indicate confidence level in the reported fexp values. The corresponding computed oscillator strengths (fcomp) were obtained with time-dependent density functional theory and a polarizable continuum solvent model. By expressing experimental and computed absorption strengths using a common quantity, we directly compared fcomp and fexp. Although fcomp and fexp are well correlated (linear regression R2 = 0.921), fcomp in most cases overestimated fexp (regression slope = 1.34). The agreement between absolute fcomp and fexp values was substantially improved by accounting for a solvent refractive index factor, as suggested in some derivations in the literature. The 100 digitized UV-visible spectra are included as plain text files in the Supporting Information to aid in benchmarking computational or machine learning methods that aim to simulate realistic UV-visible absorption spectra.

A polarizable coarse-grained model for metal, metal oxide and composite metal/metal oxide nanoparticles and its applications.

We present a selected set of exemplifying applications of the novel polarizable coarse-grained model [see the first part] to various outstanding problems in the physics and chemistry of nanoparticles: electrostatic potential around silver and gold nanoparticles; spontaneous and external electric field-driven self-organization of gold and silver nanoparticle systems; and physisorption of carbon dioxide on titanium dioxide nanoparticles decorated with a gold catalyst. In the first application, the developed model has shown capabilities of predicting long-range potential with accuracy comparable to the tight-binding density functional theory methods. Furthermore, the electrostatic potential analysis in hot spot regions allowed us to identify twin defect lines in a silver nanostar as a promising candidate for an enhancer in surface-enhanced Raman spectroscopy. In the second application, the developed model has facilitated the elucidation of the microscopic mechanisms responsible for the self-organization of gold and silver nanoparticles. Analysis of Monte-Carlo simulations established that the self-organization process is driven by van der Waals interactions in the absence of an external electric field, and that it becomes gradually driven by electrostatic interactions in the presence of an external electric field with increasing strength of the external electric field. In the third application, the developed model combined with Monte-Carlo simulations has identified the dominant mechanism responsible for carbon dioxide transfer to the catalytic sites. Analysis of the obtained results indicates that surface diffusion is the dominant mechanism for the transport of carbon dioxide to the catalytic sites, and only in exceptional situations, direct physisorption becomes a competitive mechanism with the surface diffusion mechanism. These successful applications of the developed model indicate its wide range of applicability to various problems in the chemistry and physics of nanoparticles.

A polarizable coarse-grained model for metal, metal oxide and composite metal/metal oxide nanoparticles: development and implementation.

We present a polarizable coarse-grained model for metal, metal oxide, and composite metal/metal oxide nanoparticles with well-defined crystalline surfaces. The developed model uses a low-resolution polarizable "surface beads" representation of the nanoparticle's geometry and pairwise cross nanoparticle potential consisting of van der Waals and electrostatic interaction terms. The electrostatic interaction term of the cross nanoparticle potential incorporates a crucial physical aspect of electrostatic interaction into the metal and metal oxide systems, such as induced surface charges, making it possible to explore the nanoparticles' behavior in complex environments as well as investigate the interplay between electrostatic and van der Waals interactions in nanoparticle systems. The iterative stability, computational scaling, and performance of the presented model was tested on selected systems of gold, titanium dioxide, and composite gold/titanium dioxide nanoparticle systems. The model exhibits robust iterative stability and is able to converge the charge equilibration equation for fluctuating induced charges and dipoles within 10-60 "tug-tow" iterations in challenging situations, like crowded nanoparticle systems or nanoparticle systems in extreme external electric fields. The computation scaling of the presented model is semi-linear with respect to the number of nanoparticles in the system. It slightly varies depending on the size distribution of nanoparticles in a specific nanoparticle system. The computation cost of the model is significantly lower than that of conventional atomistic polarizable force field models and enables the treatment of large nanoparticle systems that are beyond the reach of currently existing atomistic force field models.

XXVI International Symposium Nanophysics and Nanoelectronics", Nizhny Novgorod, March 14-March 17, 2022 Effect of chemical passivation of GaAs(001) surface on anisotropy and orientation of gold nanoclusters formed on it and their plasmons

The principal role of chemical passivation of GaAs surface in the formation on it of oriented anisotropic nanoclusters of gold is discussed. The nanoclusters are fabricated by thermal annealing of a gold film deposited onto GaAs(001) surface passivated as a preliminary by a monolayer of nitrogen or sulfur atoms. These atoms, bonded chemically to gallium atoms of the crystal surface, form a crystal lattice and prevent the chemical interaction of Au with GaAs. As a result of annealing, the arrays of anisotropic (elongated) nanoclusters of chemically pure Au oriented preferably in crystal [110] direction are formed on passivated GaAs(001) surface. The presence of strong anisotropy and orientation of Au clusters on passivated GaAs surfaces is established by the methods of probe diagnostics and of optical reflectance anisotropy spectroscopy and polarized reflection spectroscopy. Using an optical model of plasmonic polarizability of elongated Au spheroids, it is shown that the spectral features observed in polarized reflection originate from anisotropic plasmons of Au nanoclusters polarized mainly in direction [110] of crystal. Keywords: semiconductor surface, nitride passivation, gold nanoclusters, anisotropic plasmons, polarized reflectance.

XXVI Международный симпозиум Нанофизика и наноэлектроника", Нижний Новгород, 14-17 марта 2022 г. Влияние химической пассивации поверхности GaAs(001) на анизотропию и ориентацию образующихся на ней нанокластеров золота и их плазмонов

A principal role of chemical passivation of GaAs surface in the formation of oriented anisotropic on-surface nanoclusters of gold is discussed. The gold nanoclusters are fabricated by thermal annealing of the gold film deposited onto GaAs(001) surface passivated preliminary by an atomic monolayer of nitrogen or sulphur. These atoms forming a crystal lattice are bonded chemically to gallium atoms of the crystal surface to prevent chemical interaction of Au with GaAs. Annealing Au on passivated GaAs(001) surface results in formation of arrays of anisotropic (elongated) nanolusters of chemically pure Au oriented preferably in crystal direction. The presence of strong anisotropy and orientation of Au clusters on passivated GaAs surfaces is established by the methods of probe diagnostics and the optical techniques of reflectance anisotropy spectroscopy and polarized reflection spectroscopy. Within an optical model of plasmon polarizability of elongated Au spheroids, it is shown that the spectral features observed in polarized reflection originate from anisotropic plasmons polarized in Au nanoclusters predominantly in direction on GaAs(001) surface.

Predicting the structure configuration and Raman analysis of caffeine molecules encapsulated into single-walled carbon nanotubes: Evidence for charge transfer

A new hybrid nano-system constituted of a single-walled carbon nanotube (SWNT) filled by caffeine (Caff) molecule(s) is proposed as a potential candidate for organic solar cells. The stability of this hybrid system with either a single or two Caff molecule(s) encapsulated into SWNTs has been investigated. In particular, the optimal SWNT diameter is discussed for each configuration. Raman spectra have been calculated using an approach combining the density functional theory, molecular mechanics, bond polarizability model and the spectral moment’s method. We have analyzed the influence of the nanoconfinement on the Raman-active modes of Caff molecule and those of SWNT. The analysis of the nanoconfinement effect on the Raman active modes of Caff molecules and those of SWNTs (radial breathing mode (RBM) and tangential mode (TM)) confirm the structural stability and indicate an evidence for a charge transfer between the Caff molecules and the nanotubes.

On the Cartesian Representation of the Molecular Polarizability Tensor Surface by Polynomial Fitting to Ab Initio Data.

We describe an approach to constructing an analytic Cartesian representation of the molecular dipole polarizability tensor surface in terms of polynomials in interatomic distances with a training set of ab initio data points obtained from a molecular dynamics (MD) simulation or by any other available means. The proposed formulation is based on a perturbation treatment of the unmodified point dipole polarizability model of Applequist [ J. Am. Chem. Soc. 1972, 94, 2952] and is shown here to be, by construction (i) free of short-range or other singularities or discontinuities, (ii) symmetric and translationally invariant, and (iii) nonreliant on a body-fixed coordinate system. Permutational invariance of like nuclei is demonstrated to be readily applicable, making this approach useful for highly fluxional and reactive systems. Derivation of the method is described in detail, adding brief didactic numerical examples of H2 and H2O and concluding with an MD simulation of the Raman spectrum of H5O2+ at 300 K with the polarizability tensor fitted to CCSD(T)/aug-cc-pVTZ data obtained using the HBB-4B potential [ J. Chem. Phys. 2005, 122, 044308].

Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method

In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. We show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.

Energetics of quantum vacuum friction: Field fluctuations

Quantum fluctuations can induce a friction on a neutral but polarizable particle and cause it to radiate energy even if the particle is moving in free space filled with blackbody radiation, and is not in contact with or close to any surface or other object. We explore the energetics of such a particle moving uniformly in vacuum, continuing our previous investigations of quantum friction. The intrinsic polarizability of the particle is considered to be purely real before it is dressed by radiation. The particle is then guaranteed to be in the nonequilibrium steady state (NESS), where it absorbs and emits energy at the same rate. We first calculate the quantum frictional power and force on the particle in the rest frame of the blackbody radiation from first principles, namely the Maxwell-Heaviside equations and the Lorentz force law. Then we provide a simpler method of obtaining the same quantities in the rest frame of the particle by using the principle of virtual work. The equivalence of the two approaches is illustrated. The formulas we derive for quantum vacuum frictional power and force are fully relativistic and applicable to finite temperature. In NESS, the quantum vacuum frictional force on the particle is shown to be a true drag, independent of the model for polarizability and the polarization state of the particle. Finally, we give an estimate of the quantum vacuum friction on a gold atom and comment on the feasibility of detecting such quantum vacuum frictional effects.

Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential.

Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C4 term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C4 term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C4 term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.

Water Nanoconfined in a Hydrophobic Pore: Molecular Dynamics Simulations of Transmembrane Protein 175 and the Influence of Water Models.

Water molecules within biological ion channels are in a nano-confined environment and therefore exhibit behaviours which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously de-wet to form a ‘vapour lock’ if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/de-wetting behaviour in the TMEM175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/de-wetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/de-wetting behaviours, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nano-confined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/de-wetting processes. These findings are of importance for the design of functionalised biomimetic nanopores (for e.g. sensing or desalination), as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.

Comparative study of quantum confinements effect present in Silicon Nanowires using absorption and Raman spectroscopy

Silver-induced chemical etching (SIE) has been widely explored for developing silicon (Si) based photovoltaic devices with its benefits for low-cost and large-area fabrication of silicon nanowires (SiNWs) of high aspect ratios. Structures and optical properties of SiNWs fabricated through chemical etching are significantly affected by experimental as well as environmental conditions of etching. A comparative study of quantum confinement (QC) effect due to variation in nano crystallite size (NCS) of SiNWs samples, fabricated by SIE technique has been investigated here using ultraviolet visible (UV-VIS) spectroscopy and Raman spectroscopy. An enhanced optical band gap has been observed from the samples under UV excitation at room temperature due to QC effect. To improve the understanding of basic mechanism of QC effect present in samples has been analyzed through characteristic features of Raman line-shape like red shift, asymmetry ratio and broadening as compared to its bulk counterpart. The NCS of SiNWs present in the samples has been calculated by QC based Brus model, bond polarizability (BP) model and Raman line-shape fitting model. This study prove itself to explore a new direction towards a systematic building block for fabricating optoelectronic devices based on light - matter interactions present at the nano-scale. • Silicon nanowires (Si-NWs) onn type-Si(100) by silver-induced chemical etching. • Results obtained, indicate that the fabrication of size controlled Si-NWs. • Quantum confinements effectobserved in Si-NWs using Brus, bond polarizability and line shape model. • Estimated size from Brus and bond polarizability model using PL and Raman data were around 1.7–3.2 nm. • Theoretical estimation of size from line shape model were 2.66–3.16 nm.

Global Optimization of the Lennard-Jones Parameters for the Drude Polarizable Force Field.

Molecular dynamics (MD) simulations based on atomic models play an important role in the drug-discovery process to screen molecules, estimate binding free energies, and optimize lead compounds in chemical space. Accurate computations of thermodynamic and kinetic properties using MD simulations are highly dependent on the accuracy of the underlying atomic force field. In this context, going beyond the nonpolarizable fixed-charge model by accounting explicitly for induced polarization is highly desirable. The CHARMM polarizable force field based on classical Drude oscillators, in which an auxiliary charged particle is attached via a harmonic spring to its parent nucleus, offers both a computationally convenient and rigorous framework to model explicitly induced electronic polarization in MD simulations. For any molecule of interest, electrostatic partial charges, atomic polarizabilities, and Thole shielding factors, as well as bonded parameters can either be determined from ab initio calculations or ascribed from the knowledge-based library of the CHARMM Generalized force field. While this approach is fairly reliable in general, it is well understood that the overall accuracy of the models with respect to thermodynamic properties such as bulk density, enthalpies, and solvation free energies is particularly sensitive to the nonbonded Lennard-Jones (LJ) parameters. In the present study, we systematically refined the set of LJ parameters for the atom types available in the Drude force field to best match the experimental thermodynamic properties for 416 small druglike organic molecules. To further test the transferability of the optimized parameters, the hydration free energy of 372 molecules was computed. The calculations resulted in a small average error of 0.46 kcal/mol and a Pearson R of 0.9, representing a significant improvement over the additive GAFF force field in our previous study, where an average error of ∼2 kcal/mol was obtained. Such an improvement is consistent with the ability of the polarizable Drude model to more accurately model interactions in different environments. The effort provides a roadmap for the global optimization of force field parameters using experimental data. It is hoped that the present effort will further the application of the Drude polarizable force field in molecular simulations including drug design and discovery.

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λ6-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 [NO3]-, tetrafluoroborate [BF4]-, perchlorate [ClO4]-, hexafluorophosphate [PF6]-, 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.

Computing the Liquidus of Binary Monatomic Salt Mixtures with Direct Simulation and Alchemical Free Energy Methods.

We describe and validate a free-energy-based method for computing the liquidus for binary solid-liquid phase diagrams in molecular simulations of monatomic salts. The method is demonstrated by calculating the liquidus for LiCl-KCl and MgCl2-KCl salt mixtures with the polarizable ion model (PIM). The free-energy-based method is cross-validated with direct coexistence simulations. Both techniques show excellent agreement with one another. Though the predictions of the PIM disagree with experiments, we use our free-energy-based approach to decouple the contributions of liquid mixture nonidealities and pure component solid-liquid equilibrium to the phase diagram. In both mixtures, the PIM accurately reproduces the liquid phase nonidealities but fails to predict the liquidus because it does not accurately predict the pure component melting temperature of LiCl or MgCl2.

Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study.

Picosecond fast motions and their involvement in the biochemical processes such as protein–ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein–ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme–triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues’ spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein–ligand binding chemistry and dynamics.

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.

Investigating the Accuracy of Water Models through the Van Hove Correlation Function.

We present molecular-simulation-based calculations of the Van Hove correlation function (VHF) of water using multiple modeling approaches: classical molecular dynamics with simple three-site nonpolarizable models, with a polarizable model, and with a reactive force field; density functional tight-binding molecular dynamics; and ab initio molecular dynamics. Due to the many orders of magnitude difference in the computational cost of these approaches, we investigate how small and short the simulations can be while still yielding sufficiently accurate and interpretable results for the VHF. We investigate the accuracy of the different models by comparing them to recently published inelastic X-ray scattering measurements of the VHF. We find that all of the models exhibit qualitative agreement with the experiments, and in some models and for some properties, the agreement is quantitative. This work lays the foundation for future simulation approaches to calculating the VHF for aqueous solutions in bulk and under nanoconfinement.