1
|
Li Y, Li Z, Misra RP, Liang C, Gillen AJ, Zhao S, Abdullah J, Laurence T, Fagan JA, Aluru N, Blankschtein D, Noy A. Molecular transport enhancement in pure metallic carbon nanotube porins. NATURE MATERIALS 2024; 23:1123-1130. [PMID: 38937586 DOI: 10.1038/s41563-024-01925-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2023] [Accepted: 05/16/2024] [Indexed: 06/29/2024]
Abstract
Nanofluidic channels impose extreme confinement on water and ions, giving rise to unusual transport phenomena strongly dependent on the interactions at the channel-wall interface. Yet how the electronic properties of the nanofluidic channels influence transport efficiency remains largely unexplored. Here we measure transport through the inner pores of sub-1 nm metallic and semiconducting carbon nanotube porins. We find that water and proton transport are enhanced in metallic nanotubes over semiconducting nanotubes, whereas ion transport is largely insensitive to the nanotube bandgap value. Molecular simulations using polarizable force fields highlight the contributions of the anisotropic polarizability tensor of the carbon nanotubes to the ion-nanotube interactions and the water friction coefficient. We also describe the origin of the proton transport enhancement in metallic nanotubes using deep neural network molecular dynamics simulations. These results emphasize the complex role of the electronic properties of nanofluidic channels in modulating transport under extreme nanoscale confinement.
Collapse
Affiliation(s)
- Yuhao Li
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Zhongwu Li
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Rahul Prasanna Misra
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Chenxing Liang
- Walker Department of Mechanical Engineering, Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, USA
| | - Alice J Gillen
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
- Vivani Medical Inc., Emeryville, CA, USA
| | - Sidi Zhao
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
- School of Engineering, University of California Merced, Merced, CA, USA
| | - Jobaer Abdullah
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
- School of Natural Sciences, University of California Merced, Merced, CA, USA
| | - Ted Laurence
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Jeffrey A Fagan
- Materials Science and Engineering Division, National Institute of Standards and Technology (NIST), Gaithersburg, MD, USA
| | - Narayana Aluru
- Walker Department of Mechanical Engineering, Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, USA.
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA.
| | - Aleksandr Noy
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA.
- School of Natural Sciences, University of California Merced, Merced, CA, USA.
| |
Collapse
|
2
|
Rezaei M, Sakong S, Groß A. Sodium Triflate Water-in-Salt Electrolytes in Advanced Battery Applications: A First-Principles-Based Molecular Dynamics Study. ACS APPLIED MATERIALS & INTERFACES 2024; 16:32169-32188. [PMID: 38862108 PMCID: PMC11212028 DOI: 10.1021/acsami.4c01449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 04/03/2024] [Accepted: 05/29/2024] [Indexed: 06/13/2024]
Abstract
Offering a compelling combination of safety and cost-effectiveness, water-in-salt (WiS) electrolytes have emerged as promising frontiers in energy storage technology. Still, there is a strong demand for research and development efforts to make these electrolytes ripe for commercialization. Here, we present a first-principles-based molecular dynamics (MD) study addressing in detail the properties of a sodium triflate WiS electrolyte for Na-ion batteries. We have developed a workflow based on a machine learning (ML) potential derived from ab initio MD simulations. As ML potentials are typically restricted to the interpolation of the data points of the training set and have hardly any predictive properties, we subsequently optimize a classical force field based on physics principles to ensure broad applicability and high performance. Performing and analyzing detailed MD simulations, we identify several very promising properties of the sodium triflate as a WiS electrolyte but also indicate some potential stability challenges associated with its use as a battery electrolyte.
Collapse
Affiliation(s)
- Majid Rezaei
- Institute
of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
| | - Sung Sakong
- Institute
of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
| | - Axel Groß
- Institute
of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
- Helmholtz
Institute Ulm (HIU) for Electrochemical Energy Storage, Helmholtzstraße 11, 89069 Ulm, Germany
| |
Collapse
|
3
|
Mohandas N, Bawari S, Shibuya JJT, Ghosh S, Mondal J, Narayanan TN, Cuesta A. Understanding electrochemical interfaces through comparing experimental and computational charge density-potential curves. Chem Sci 2024; 15:6643-6660. [PMID: 38725490 PMCID: PMC11077530 DOI: 10.1039/d4sc00746h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Accepted: 04/21/2024] [Indexed: 05/12/2024] Open
Abstract
Electrode-electrolyte interfaces play a decisive role in electrochemical charge accumulation and transfer processes. Theoretical modelling of these interfaces is critical to decipher the microscopic details of such phenomena. Different force field-based molecular dynamics protocols are compared here in a view to connect calculated and experimental charge density-potential relationships. Platinum-aqueous electrolyte interfaces are taken as a model. The potential of using experimental charge density-potential curves to transform cell voltage into electrode potential in force-field molecular dynamics simulations, and the need for that purpose of developing simulation protocols that can accurately calculate the double-layer capacitance, are discussed.
Collapse
Affiliation(s)
- Nandita Mohandas
- Tata Institute of Fundamental Research-Hyderabad Hyderabad 500046 India
- Advanced Centre for Energy and Sustainability (ACES), School of Natural and Computing Sciences, University of Aberdeen AB24 3UE Aberdeen Scotland UK
| | - Sumit Bawari
- Tata Institute of Fundamental Research-Hyderabad Hyderabad 500046 India
| | - Jani J T Shibuya
- Advanced Centre for Energy and Sustainability (ACES), School of Natural and Computing Sciences, University of Aberdeen AB24 3UE Aberdeen Scotland UK
| | - Soumya Ghosh
- Tata Institute of Fundamental Research-Hyderabad Hyderabad 500046 India
| | - Jagannath Mondal
- Tata Institute of Fundamental Research-Hyderabad Hyderabad 500046 India
| | | | - Angel Cuesta
- Advanced Centre for Energy and Sustainability (ACES), School of Natural and Computing Sciences, University of Aberdeen AB24 3UE Aberdeen Scotland UK
- Centre for Energy Transition, University of Aberdeen AB24 3FX Aberdeen Scotland UK
| |
Collapse
|
4
|
Corsini C, M Correa C, Scaglione N, Costa Gomes M, Padua A. How Do Deep Eutectic Solvents Form Porous Liquids? The Example of Methyltriphenylphosphonium Bromide: Glycerol and ZIF-8. J Phys Chem B 2024. [PMID: 38433612 DOI: 10.1021/acs.jpcb.3c08490] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2024]
Abstract
Porous liquids are new materials that provide permanent porosity in the liquid phase through the dispersion of nanoporous solid particles in a bulky solvent. Herein, we aim at understanding how new sustainable solvents such as deep eutectic solvent (DES) can be used to form porous stable suspensions for the capture of gases of interest for sustainable chemistry. The properties of an ionic DES, methyltriphenylphosphonium bromide/glycerol in a 1:3 molar composition, and its behavior at the interface with a metal-organic framework (MOF), ZIF-8, are here investigated by polarizable molecular dynamics simulations. The structural and dynamic properties of the DES are analyzed in the bulk liquid and in the interfacial regions with the MOF, namely, in the accessible cavities exposed at the surface. The porosity of the suspension is maintained, and it is caused not only by the Coulomb cohesive energy between cations and anions, which prevents the small anions from entering the pores, but also by the glycerol molecules not penetrating the small aperture of the ZIF-8 structure. This was further verified by simulating a system composed of glycerol and ZIF-8. Simulations with CO2 show its partition between the DES and the MOF, with higher concentrations registered in the MOF cavities.
Collapse
Affiliation(s)
- Chiara Corsini
- Laboratoire de Chimie de l'ENS Lyon, CNRS and Université de Lyon, 46 allée d'Italie, 69364 Lyon, France
| | - Cintia M Correa
- Laboratoire de Chimie de l'ENS Lyon, CNRS and Université de Lyon, 46 allée d'Italie, 69364 Lyon, France
| | - Nicolas Scaglione
- Laboratoire de Chimie de l'ENS Lyon, CNRS and Université de Lyon, 46 allée d'Italie, 69364 Lyon, France
| | - Margarida Costa Gomes
- Laboratoire de Chimie de l'ENS Lyon, CNRS and Université de Lyon, 46 allée d'Italie, 69364 Lyon, France
| | - Agilio Padua
- Laboratoire de Chimie de l'ENS Lyon, CNRS and Université de Lyon, 46 allée d'Italie, 69364 Lyon, France
| |
Collapse
|
5
|
Luo S, Misra RP, Blankschtein D. Water Electric Field Induced Modulation of the Wetting of Hexagonal Boron Nitride: Insights from Multiscale Modeling of Many-Body Polarization. ACS NANO 2024; 18:1629-1646. [PMID: 38169482 DOI: 10.1021/acsnano.3c09811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Understanding the behavior of water contacting two-dimensional materials, such as hexagonal boron nitride (hBN), is important in practical applications, including seawater desalination and energy harvesting. Water, being a polar solvent, can strongly polarize the hBN surface via the electric fields that it generates. However, there is a lack of molecular-level understanding about the role of polarization effects at the hBN/water interface, including its effect on the wetting properties of water. In this study, we develop a theoretical framework that introduces an all-atomistic polarizable force field to accurately model the interactions of water molecules with hBN surfaces. The force field is then utilized to self-consistently describe the water-induced polarization of hBN using the classical Drude oscillator model, including predicting the hBN-water binding energies which are found to be in excellent agreement with diffusion Monte Carlo (DMC) predictions. By carrying out molecular dynamics (MD) simulations, we demonstrate that the polarizable force field yields a water contact angle on multilayered hBN which is in close agreement with the recent experimentally reported values. Conversely, an implicit modeling of the hBN-water polarization energy utilizing a Lennard-Jones (LJ) potential, a commonly utilized approximation in previous MD simulation studies, leads to a considerably lower water contact angle. This difference in the predicted contact angles is attributed to the significant energy-entropy compensation resulting from the incorporation of polarization effects at the hBN-water interface. Our work highlights the importance of self-consistently modeling the hBN-water polarization energy and offers insights into the wetting-related interfacial phenomena of water on polarizable materials.
Collapse
Affiliation(s)
- Shuang Luo
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Rahul Prasanna Misra
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
6
|
Jedlinska ZM, Riggleman RA. The effect of monomer polarizability on the stability and salt partitioning in model coacervates. SOFT MATTER 2023; 19:7000-7010. [PMID: 37668019 DOI: 10.1039/d3sm00706e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
Coacervation of charged polymer chains has been a topic of major interest in both polymer and biological sciences, as it is a subset of a phenomenon called liquid-liquid phase separation (LLPS). In this process the polymer-rich phase separates from the polymer-lean supernatant while still maintaining its liquid-like properties. LLPS has been shown to play a crucial role in cellular homeostasis by driving the formation of membraneless organelles. It also has the potential to be harnessed to aid in novel therapeutical applications. Recent studies have demonstrated that there is no one simple mechanism which drives LLPS, which is instead a result of the combined effect of electrostatic, dipolar, hydrophobic, and other weak interactions. Using coarse-grained polymer simulations we investigate the relatively unexplored effects of monomer polarizability and spatially varying dielectric constant on LLPS propensity, and these factors affect the properties of the resulting condensates. In order to produce spatial variations in the dielectric constant, all our simulations include explicit solvent and counterions. We demonstrate that polarizability has only a minor effect on the bulk behaviour of the condensates but plays a major role when ion partitioning and microstructure are considered. We observe that the major contribution comes from the nature of the neutral blocks as endowing them with an induced dipole changes their character from hydrophobic to hydrophilic. We hypothesize that the results of this work can aid in guiding future studies concerned with LLPS by providing a general framework and by highlighting important factors which influence LLPS.
Collapse
Affiliation(s)
- Zuzanna M Jedlinska
- Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, USA
| | - Robert A Riggleman
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, USA.
| |
Collapse
|
7
|
Rezaei M, Sakong S, Groß A. Molecular Modeling of Water-in-Salt Electrolytes: A Comprehensive Analysis of Polarization Effects and Force Field Parameters in Molecular Dynamics Simulations. J Chem Theory Comput 2023; 19:5712-5730. [PMID: 37528639 DOI: 10.1021/acs.jctc.3c00171] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/03/2023]
Abstract
Accurate modeling of highly concentrated aqueous solutions, such as water-in-salt (WiS) electrolytes in battery applications, requires proper consideration of polarization contributions to atomic interactions. Within the force field molecular dynamics (MD) simulations, the atomic polarization can be accounted for at various levels. Nonpolarizable force fields implicitly account for polarization effects by incorporating them into their van der Waals interaction parameters. They can additionally mimic electron polarization within a mean-field approximation through ionic charge scaling. Alternatively, explicit polarization description methods, such as the Drude oscillator model, can be selectively applied to either a subset of polarizable atoms or all polarizable atoms to enhance simulation accuracy. The trade-off between simulation accuracy and computational efficiency highlights the importance of determining an optimal level of accounting for atomic polarization. In this study, we analyze different approaches to include polarization effects in MD simulations of WiS electrolytes, with an example of a Na-OTF solution. These approaches range from a nonpolarizable to a fully polarizable force field. After careful examination of computational costs, simulation stability, and feasibility of controlling the electrolyte properties, we identify an efficient combination of force fields: the Drude polarizable force field for salt ions and non-polarizable models for water. This cost-effective combination is sufficiently flexible to reproduce a broad range of electrolyte properties, while ensuring simulation stability over a relatively wide range of force field parameters. Furthermore, we conduct a thorough evaluation of the influence of various force field parameters on both the simulation results and technical requirements, with the aim of establishing a general framework for force field optimization and facilitating parametrization of similar systems.
Collapse
Affiliation(s)
- Majid Rezaei
- Institute of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, Oberberghof 7, 89081 Ulm, Germany
- Helmholtz Institute Ulm (HIU) for Electrochemical Energy Storage, Helmholtzstraße 11, 89069 Ulm, Germany
| |
Collapse
|
8
|
Wang Y, Wang K, Zhao X, Xu X, Sun T. Influence of pH on the self-assembly of diphenylalanine peptides: molecular insights from coarse-grained simulations. SOFT MATTER 2023; 19:5749-5757. [PMID: 37462931 DOI: 10.1039/d3sm00739a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/03/2023]
Abstract
Nanostructures fabricated from peptide self-assemblies are attracting increasing attention owing to their possible applications in biology and nanotechnology. A known example is an aromatic dipeptide (diphenylalanine, FF) which is extracted from Alzheimer's β-amyloid polypeptide as the core recognition motif for molecular self-assembly. Many studies have been carried out to organize FF peptides into various functional ordered nanostructures. For potential applications of self-assembled FF-based nanomaterials, it becomes important to consider some influencing factors (e.g., solvents, peptide concentrations, pH, temperature, etc.) on the self-assembly process. Among these factors, the effect of pH on the self-assembly process of FF peptides into assembled nanostructures through simulation studies is the main focus of the present work. In the current study, we have investigated the assembly pathway of 1000 FF peptides and qualitatively evaluated the morphological changes of FF-based nanostructures at different pH values by performing extensive coarse-grained molecular dynamics (CG-MD) simulations. Structural analyses suggest that FF peptides can spontaneously assemble into nanotubes with different shapes under acidic, neutral and basic conditions. Based on the analysis of FF nanostructure formation pathways in different pH solutions, the self-assembly of the nanotube involves the aggregation of molecules to form a bilayer, the curling of a bilayer to form a vesicle and the transformation of a vesicle into a tubular structure. It is noted that a flat hollow columnar structure is observed as a special intermediate state during the transformation process of a vesicle-like to a tube-like structure. Energetic analysis suggests that the aggregation of FF peptides is driven by the vdW interactions but the aggregation shape is mainly affected by the electrostatic interactions. Overall, this study provides further understanding of the self-assembly behavior of aromatic short peptide derivatives in different pH solutions.
Collapse
Affiliation(s)
- Yan Wang
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| | - Kang Wang
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| | - Xinyi Zhao
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| | - Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou, Jiangsu 213001, China
| | - Tingting Sun
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| |
Collapse
|
9
|
Pieroni M, Madeddu F, Di Martino J, Arcieri M, Parisi V, Bottoni P, Castrignanò T. MD-Ligand-Receptor: A High-Performance Computing Tool for Characterizing Ligand-Receptor Binding Interactions in Molecular Dynamics Trajectories. Int J Mol Sci 2023; 24:11671. [PMID: 37511429 PMCID: PMC10380688 DOI: 10.3390/ijms241411671] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Revised: 07/15/2023] [Accepted: 07/18/2023] [Indexed: 07/30/2023] Open
Abstract
Molecular dynamics simulation is a widely employed computational technique for studying the dynamic behavior of molecular systems over time. By simulating macromolecular biological systems consisting of a drug, a receptor and a solvated environment with thousands of water molecules, MD allows for realistic ligand-receptor binding interactions (lrbi) to be studied. In this study, we present MD-ligand-receptor (MDLR), a state-of-the-art software designed to explore the intricate interactions between ligands and receptors over time using molecular dynamics trajectories. Unlike traditional static analysis tools, MDLR goes beyond simply taking a snapshot of ligand-receptor binding interactions (lrbi), uncovering long-lasting molecular interactions and predicting the time-dependent inhibitory activity of specific drugs. With MDLR, researchers can gain insights into the dynamic behavior of complex ligand-receptor systems. Our pipeline is optimized for high-performance computing, capable of efficiently processing vast molecular dynamics trajectories on multicore Linux servers or even multinode HPC clusters. In the latter case, MDLR allows the user to analyze large trajectories in a very short time. To facilitate the exploration and visualization of lrbi, we provide an intuitive Python notebook (Jupyter), which allows users to examine and interpret the results through various graphical representations.
Collapse
Affiliation(s)
- Michele Pieroni
- Department of Computer Science, "Sapienza" University of Rome, V. le Regina Elena 295, 00161 Rome, Italy
| | - Francesco Madeddu
- Department of Computer Science, "Sapienza" University of Rome, V. le Regina Elena 295, 00161 Rome, Italy
| | - Jessica Di Martino
- Department of Ecological and Biological Sciences, Tuscia University, Viale dell'Università s.n.c., 01100 Viterbo, Italy
| | - Manuel Arcieri
- Department of Health Technology, Technical University of Denmark, Anker Engelunds Vej 101, 2800 Kongens Lyngby, Denmark
| | - Valerio Parisi
- Department of Physics, "Sapienza" University of Rome, P. le Aldo Moro, 5, 00185 Rome, Italy
| | - Paolo Bottoni
- Department of Computer Science, "Sapienza" University of Rome, V. le Regina Elena 295, 00161 Rome, Italy
| | - Tiziana Castrignanò
- Department of Ecological and Biological Sciences, Tuscia University, Viale dell'Università s.n.c., 01100 Viterbo, Italy
| |
Collapse
|
10
|
Buyukdagli S. Impact of the inner solute structure on the electrostatic mean-field and strong-coupling regimes of macromolecular interactions. Phys Rev E 2023; 107:064604. [PMID: 37464605 DOI: 10.1103/physreve.107.064604] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Accepted: 06/04/2023] [Indexed: 07/20/2023]
Abstract
The structural diversity of the solute molecules involved in biomolecular processes necessitates the characterization of the forces between charged macromolecules beyond the point-ion description. From the field-theoretic partition function of an electrolyte confined between two anionic membranes, we derive a contact-value identity valid for general intramolecular solute structure and electrostatic coupling strength. In the electrostatic mean-field regime, the inner charge spread of the solute particles is shown to induce the twofold enhancement of the short-range Poisson-Boltzmann level membrane repulsion and a longer-range depletion attraction. Our contact theorem indicates that the twofold repulsion enhancement by solute size is equally present in the opposite strong-coupling regime of linear and spherical solute molecules. Upon the inclusion of the dielectric contrast between the electrolyte and the interacting membranes, the emerging polarization forces substantially amplify the solute specificity of the macromolecular interactions. Namely, the finite size of the dumbbell-like solute particles composed of similar terminal charges weakens the intermembrane repulsion. However, the extended structure of the solute molecules carrying opposite elementary charges such as ionized atoms and zwitterionic molecules enhances the membrane repulsion by several factors. We also show that these polarization forces can extend the range of the solute structure effects up to intermembrane distances exceeding the solute size by an order of magnitude. This radical alteration of the intermembrane interactions by the salt structure identifies the solute specificity as a key ingredient of the thermodynamic stability in colloidal systems.
Collapse
|
11
|
Dodin A, Geissler PL. Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics. J Chem Theory Comput 2023; 19:2906-2917. [PMID: 37130215 DOI: 10.1021/acs.jctc.3c00278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Drude oscillator potentials are a popular and computationally efficient class of polarizable models that represent each polarizable atom as a positively charged Drude core harmonically bound to a negatively charged Drude shell. We show that existing force fields that place all non-Coulomb forces on the Drude core and none on the shell inadvertently couple the dipole to non-Coulombic forces. This introduces errors where interactions with neutral particles can erroneously induce atomic polarization, leading to spurious polarizations in the absence of an electric field, exacerbating violations of equipartition in the employed Carr-Parinello scheme. A suitable symmetrization of the interaction potential that correctly splits the force between the Drude core and shell can correct this shortcoming, improving the stability and numerical performance of Drude oscillator-based simulations. The symmetrization procedure is straightforward and only requires the rescaling of a few force field parameters.
Collapse
Affiliation(s)
- Amro Dodin
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Phillip L Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| |
Collapse
|
12
|
Gullbrekken Ø, Røe IT, Selbach SM, Schnell SK. Charge Transport in Water-NaCl Electrolytes with Molecular Dynamics Simulations. J Phys Chem B 2023; 127:2729-2738. [PMID: 36921121 PMCID: PMC10068734 DOI: 10.1021/acs.jpcb.2c08047] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023]
Abstract
A systematic description of microscopic mechanisms is necessary to understand mass transport in solid and liquid electrolytes. From Molecular Dynamics (MD) simulations, transport properties can be computed and provide a detailed view of the molecular and ionic motions. In this work, ionic conductivity and transport numbers in electrolyte systems are computed from equilibrium and nonequilibrium MD simulations. Results from the two methods are compared with experimental results, and we discuss the significance of the frame of reference when determining and comparing transport numbers. Two ways of computing ionic conductivity from equilibrium simulations are presented: the Nernst-Einstein approximation or the Onsager coefficients. The Onsager coefficients take ionic correlations into account and are found to be more suitable for concentrated electrolytes. Main features and differences between equilibrium and nonequilibrium simulations are discussed, and some potential anomalies and critical pitfalls of using nonequilibrium molecular dynamics to determine transport properties are highlighted.
Collapse
Affiliation(s)
- Øystein Gullbrekken
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Ingeborg Treu Røe
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Sverre Magnus Selbach
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| | - Sondre Kvalvåg Schnell
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, Trondheim NO-7491, Norway
| |
Collapse
|
13
|
Singh A, Mason TG, Lu Z, Hill AJ, Pas SJ, Teo BM, Freeman BD, Izgorodina EI. Structural elucidation of polydopamine facilitated by ionic liquid solvation. Phys Chem Chem Phys 2023; 25:14700-14710. [PMID: 36806848 DOI: 10.1039/d2cp05439f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/13/2023]
Abstract
Minimal understanding of the formation mechanism and structure of polydopamine (pDA) and its natural analogue, eumelanin, impedes the practical application of these versatile polymers and limits our knowledge of the origin of melanoma. The lack of conclusive structural evidence stems from the insolubility of these materials, which has spawned significantly diverse suggestions of pDA's structure in the literature. We discovered that pDA is soluble in certain ionic liquids. Using these ionic liquids (ILs) as solvents, we present an experimental methodology to solvate pDA, enabling us to identify pDA's chemical structure. The resolved pDA structure consists of self-assembled supramolecular aggregates that contribute to the increasing complexity of the polymer. The underlying molecular energetics of pDA solvation and a macroscopic picture of the disruption of the aggregates using IL solvents have been investigated, along with studies of the aggregation mechanism in water.
Collapse
Affiliation(s)
- Abhishek Singh
- School of Chemistry, Monash University, Clayton, Melbourne, VIC 3800, Australia. .,IITB-Monash Research Academy, Bombay 400076, India
| | - Thomas G Mason
- School of Chemistry, Monash University, Clayton, Melbourne, VIC 3800, Australia.
| | - Zhenzhen Lu
- School of Chemistry, Monash University, Clayton, Melbourne, VIC 3800, Australia.
| | - Anita J Hill
- Manufacturing, CSIRO, Clayton, VIC 3168, Australia
| | - Steven J Pas
- Maritime Division, Defence Science and Technology Group, Department of Defence, 506 Lorimer St Fisherman's Bend, VIC 3207, Australia
| | - Boon Mia Teo
- School of Chemistry, Monash University, Clayton, Melbourne, VIC 3800, Australia.
| | - Benny D Freeman
- Department of Chemical Engineering, The University of Texas at Austin, Austin, TX 78712, USA
| | | |
Collapse
|
14
|
Chen S, Voth GA. How Does Electronic Polarizability or Scaled-Charge Affect the Interfacial Properties of Room Temperature Ionic Liquids? J Phys Chem B 2023; 127:1264-1275. [PMID: 36701801 PMCID: PMC9924258 DOI: 10.1021/acs.jpcb.2c07981] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2022] [Revised: 01/01/2023] [Indexed: 01/27/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sijia Chen
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, The James Franck Institute,
and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, The James Franck Institute,
and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois60637, United States
| |
Collapse
|
15
|
Li Z, Misra RP, Li Y, Yao YC, Zhao S, Zhang Y, Chen Y, Blankschtein D, Noy A. Breakdown of the Nernst-Einstein relation in carbon nanotube porins. NATURE NANOTECHNOLOGY 2023; 18:177-183. [PMID: 36585518 DOI: 10.1038/s41565-022-01276-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Accepted: 10/17/2022] [Indexed: 06/17/2023]
Abstract
For over 100 years, the Nernst-Einstein relation has linked a charged particle's electrophoretic mobility and diffusion coefficient. Here we report experimental measurements of diffusion and electromigration of K+ ions in narrow 0.8-nm-diameter single-walled carbon nanotube porins (CNTPs) and demonstrate that the Nernst-Einstein relation in these channels breaks down by more than three orders of magnitude. Molecular dynamics simulations using polarizable force fields show that K+ ion diffusion in CNTPs in the presence of a single-file water chain is three orders of magnitude slower than bulk diffusion. Intriguingly, the simulations also reveal a disintegration of the water chain upon application of electric fields, resulting in the formation of distinct K+-water clusters, which then traverse the CNTP at high velocity. Finally, we show that although individual ion-water clusters still obey the Nernst-Einstein relation, the overall relation breaks down because of two distinct mechanisms for ion diffusion and electromigration.
Collapse
Affiliation(s)
- Zhongwu Li
- Materials Science Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
- Jiangsu Key Laboratory for Design and Manufacture of Micro-Nano Biomedical Instruments, School of Mechanical Engineering, Southeast University, Nanjing, China
- School of Mechatronic Engineering, China University of Mining and Technology, Xuzhou, China
| | - Rahul Prasanna Misra
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Yuhao Li
- Materials Science Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Yun-Chiao Yao
- Materials Science Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
- School of Natural Sciences, University of California Merced, Merced, CA, USA
| | - Sidi Zhao
- Materials Science Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
- School of Engineering, University of California Merced, Merced, CA, USA
| | - Yuliang Zhang
- Materials Science Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Yunfei Chen
- Jiangsu Key Laboratory for Design and Manufacture of Micro-Nano Biomedical Instruments, School of Mechanical Engineering, Southeast University, Nanjing, China
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA.
| | - Aleksandr Noy
- Materials Science Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA.
- School of Natural Sciences, University of California Merced, Merced, CA, USA.
| |
Collapse
|
16
|
Balzer C, Frischknecht AL. Explicit Polarization in Coarse-Grained Simulations of Ionomer Melts. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c01608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Affiliation(s)
- Christopher Balzer
- Division of Chemistry and Chemical Engineering, California Institute of Technology, 1200 E California Blvd, Pasadena, California91125, United States
| | - Amalie L. Frischknecht
- Center for Integrated Nanotechnologies, Sandia National Laboratories, P.O. Box 5800
MS 1303, Albuquerque, New Mexico87185-1303, United States
| |
Collapse
|
17
|
Philippi F, Rauber D, Palumbo O, Goloviznina K, McDaniel J, Pugh D, Suarez S, Fraenza CC, Padua A, Kay CWM, Welton T. Flexibility is the key to tuning the transport properties of fluorinated imide-based ionic liquids. Chem Sci 2022; 13:9176-9190. [PMID: 36093026 PMCID: PMC9384794 DOI: 10.1039/d2sc03074h] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2022] [Accepted: 07/02/2022] [Indexed: 11/21/2022] Open
Abstract
Ionic liquids are becoming increasingly popular for practical applications such as biomass processing and lithium-ion batteries. However, identifying ionic liquids with optimal properties for specific applications by trial and error is extremely inefficient since there are a vast number of potential candidate ions. Here we combine experimental and computational techniques to determine how the interplay of fluorination, flexibility and mass affects the transport properties of ionic liquids with the popular imide anion. We observe that fluorination and flexibility have a large impact on properties such as viscosity, whereas the influence of mass is negligible. Using targeted modifications, we show that conformational flexibility provides a significant contribution to the success of fluorination as a design element. Contrary to conventional wisdom, fluorination by itself is thus not a guarantor for beneficial properties such as low viscosity.
Collapse
Affiliation(s)
- Frederik Philippi
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London White City Campus London W12 0BZ UK
| | - Daniel Rauber
- Department of Chemistry, Saarland University Campus B2.2 Saarbrücken Germany
| | - Oriele Palumbo
- Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi Piazzale Aldo Moro 5 00185 Rome Italy
| | | | - Jesse McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology Atlanta Georgia 30332-0400 USA
| | - David Pugh
- Department of Chemistry, King's College London 7 Trinity Street London SE1 1DB UK
| | - Sophia Suarez
- Department of Physics, Brooklyn College of CUNY Brooklyn New York 11210 USA
| | - Carla C Fraenza
- Department of Physics and Astronomy, Hunter College of CUNY New York 10065 USA
| | - Agilio Padua
- Laboratoire de Chimie, École Normale Supérieure de Lyon, CNRS 69364 Lyon France
| | - Christopher W M Kay
- Department of Chemistry, Saarland University Campus B2.2 Saarbrücken Germany
- London Centre for Nanotechnology, University College London 17-19 Gordon Street London WC1H 0AH UK
| | - Tom Welton
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London White City Campus London W12 0BZ UK
| |
Collapse
|
18
|
Reis GSA, de Souza RM, Ribeiro MCC. Molecular Dynamics Simulation Study of the Far-Infrared Spectrum of a Deep Eutectic Solvent. J Phys Chem B 2022; 126:5695-5705. [PMID: 35858287 DOI: 10.1021/acs.jpcb.2c03277] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Deep eutectic solvents (DESs) are similar to ionic liquids (IL) in terms of physicochemical properties and technical uses. In ILs, far-infrared (FIR) spectroscopy has been utilized to reveal ionic interactions and even to produce a signature of the strengthening of the cation-anion hydrogen bond. However, for the situation of the DES, where the mixing of a salt and a molecular species makes the interplay between multiple intermolecular interactions even more complex, a full investigation of FIR spectra is still absent. In this work, the FIR spectrum of the DES, often referred to as ethaline, which is a 1:2 mixture of choline chloride and ethylene glycol, is calculated using classical molecular dynamics (MD) simulations and compared to experimental data. To explore the induced dipole effect on the computed FIR spectrum, MD simulations were run with both nonpolarizable and polarizable models. The calculation satisfactorily reproduces the position of the peak at ∼110 cm-1 and the bandwidth seen in the experimental FIR spectrum of ethaline. The MD simulations show that the charge current is the most important contributor to the FIR spectrum, but the cross-correlation between the charge current and dipole reorientation also plays a role in the polarizable model. The dynamics of the chloride-ethylene glycol correlation span a wide frequency range, with a maximum at ∼150 cm-1, but it participates as a direct mechanism only in the charge current-dipole reorientation cross-term. Anion correlations, whose dynamics are regulated via correlation with both ethylene glycol and choline, make the most significant contribution to the charge current mechanism. The MD simulations were also utilized to investigate the effect on the FIR spectrum of adding water to the DES and switching to a 1:1 composition.
Collapse
Affiliation(s)
- Gabriela S A Reis
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05513-970 São Paulo, São Paulo, Brazil
| | - Rafael M de Souza
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05513-970 São Paulo, São Paulo, Brazil
| | - Mauro C C Ribeiro
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05513-970 São Paulo, São Paulo, Brazil
| |
Collapse
|
19
|
Szabadi A, Schröder C. Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2022. [DOI: 10.1142/s2737416521420035] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
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.
Collapse
Affiliation(s)
- András Szabadi
- University of Vienna, Faculty of Chemistry, Department of Computational Biological Chemistry, A-1090 Vienna, Austria
| | - Christian Schröder
- University of Vienna, Faculty of Chemistry, Department of Computational Biological Chemistry, A-1090 Vienna, Austria
| |
Collapse
|
20
|
Guan Y, Clark R, Philippi F, Zhang X, Welton T. How do external forces related to mass and charge affect the structures and dynamics of an ionic liquid?. J Chem Phys 2022; 156:204312. [DOI: 10.1063/5.0091322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ionic liquids (ILs) are novel promising materials widely used in various fields. Their structures and properties can be tuned by means of external perturbations, thus further broadening their applications. Herein, forces proportional to atomic mass (mass-related field) and atomic charge (electric field) are applied in molecular dynamics simulations to the IL 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide to investigate the origin of the resulting changes in structures and dynamics. The results show that both electric and mass-related fields cause the ion cages to expand and deform, eventually leading to their breakdown to produce a transformation of IL from cage structure to channel-like structure, which results in faster self-diffusion of ions in the directions of the applied force and to a lesser extent other directions. Further comparison of electric and mass-related fields demonstrates that only the electric fields reorientate cations to produce a hydrodynamically favoured conformation in the force direction which shows faster diffusion. The cis isomer of the anion is preferred in the presence of the electric fields, whereas applying the forces proportional to mass does not change the anion conformer equilibrium significantly. The results presented in this work aid in the understanding of how ions adjust their structures to adapt to external perturbations and facilitates the application of ILs as electrolytes.
Collapse
Affiliation(s)
- Yongji Guan
- School of Information Science and Engineering, Lanzhou University, China
| | - Ryan Clark
- Imperial College Department of Chemistry, United Kingdom
| | | | - Xiaoping Zhang
- School of Information Science and Engineering, Lanzhou University, China
| | - Thomas Welton
- Department of Chemistry, Imperial College London, United Kingdom
| |
Collapse
|
21
|
Zhang Z, Zofchak E, Krajniak J, Ganesan V. Influence of Polarizability on the Structure, Dynamic Characteristics, and Ion-Transport Mechanisms in Polymeric Ionic Liquids. J Phys Chem B 2022; 126:2583-2592. [DOI: 10.1021/acs.jpcb.1c10662] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zidan Zhang
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Everett Zofchak
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Jakub Krajniak
- Independent Researcher, os. Kosmonautow 13/56, 61-631 Poznan, Poland
| | - Venkat Ganesan
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
22
|
Bernardes CES. DLPGEN: Preparing Molecular Dynamics Simulations with Support for Polarizable Force Fields. J Chem Inf Model 2022; 62:1471-1478. [PMID: 35239343 DOI: 10.1021/acs.jcim.1c01431] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Preparing input files for molecular dynamics (MD) simulations can be a tedious task, particularly if different MD programs need to be used. Most simulation packages are accompanied by applications that handle this task, and, in some cases, software to perform interconversion between MD programs can be found. However, the conversion between different types of files is not always foolproof or the force field may not be fully supported, as quite often observed with polarization models. This work describes the program DLPGEN, which produces input files for the MD programs GROMACS, CHARMM, DL_POLY, and LAMMPS. The program can prepare polarizable force fields using a self-consistent field approach or with a dual thermostat-extended Lagrangian model. For GROMACS, a new polarization scheme suitable for the simulation of molecules containing more than one virtual particle is described. Results obtained for ethanol in the liquid state revealed that the system configurational energy, liquid density, and self-diffusion coefficient obtained with GROMACS are in good agreement with the data found with LAMMPS and CHARMM. In the case of DL_POLY, a problem with the Shells temperature control was found, suggesting that this program may not be suitable for the simulation of molecules containing multiple Drude particles.
Collapse
Affiliation(s)
- Carlos E S Bernardes
- Centro de Química Estrutural, Institute of Molecular Sciences, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Lisboa 1749-016, Portugal
| |
Collapse
|
23
|
Klajmon M, Červinka C. Does Explicit Polarizability Improve Molecular Dynamics Predictions of Glass Transition Temperatures of Ionic Liquids? J Phys Chem B 2022; 126:2005-2013. [PMID: 35195429 DOI: 10.1021/acs.jpcb.1c10809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations are used for predictions of the glass transition temperatures for a test set of five aprotic ionic liquids. Glass transitions are localized with the trend-shift method, analyzing volumetric and transport properties of bulk amorphous phases. A classical nonpolarizable all-atom OPLS force-field model developed by Canongia Lopes and Pádua (CL&P) is employed as a starting level of theory for all calculations. Alternative approaches of charge scaling and the Drude oscillator model, accounting for atomic polarizability either implicitly or explicitly, respectively, are used to investigate the sensitivity of the glass transition temperatures to induction effects. The former nonpolarizable model overestimates the glass transition temperature by tens of Kelvins (37 K on average). The charge-scaling technique yields a significant improvement, and the best estimations were achieved using polarizable simulations with the Drude model, which yielded an average deviation of 11 K. Although the volumetric data usually exhibit a lesser trend shift upon vitrification, their lower statistical uncertainty enables to predict the glass transition temperature with lower uncertainty than the ionic self-diffusivities, the temperature dependence of which is usually more scattered. Additional analyses of the simulated data were also performed, revealing that the Drude model predicts lower densities for most subcooled liquids but higher densities for the glasses than the original CL&P, and that the Drude model also invokes some longer-range organization of the subcooled liquid, greatly impacting the temperature trend of ionic self-diffusivities in the low-temperature region.
Collapse
Affiliation(s)
- Martin Klajmon
- Department of Physical Chemistry, Faculty of Chemical Engineering, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| | - Ctirad Červinka
- Department of Physical Chemistry, Faculty of Chemical Engineering, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| |
Collapse
|
24
|
Gäding J, Tocci G, Busch M, Huber P, Meißner RH. Impact of confinement and polarizability on dynamics of ionic liquids. J Chem Phys 2022; 156:064703. [DOI: 10.1063/5.0077408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Johannes Gäding
- Institute of Polymers and Composites, Hamburg University of Technology, 21073 Hamburg, Germany
| | - Gabriele Tocci
- Department of Chemistry, University of Zurich, 8057 Zürich, Switzerland
| | - Mark Busch
- Institute for Materials and X-Ray Physics, Hamburg University of Technology, 21073 Hamburg, Germany
- Deutsches Elektronen-Synchrotron DESY, Centre for X-Ray and Nano Science CXNS, 22607 Hamburg, Germany
- Centre for Hybrid Nanostructures CHyN, Hamburg University, 22761 Hamburg, Germany
| | - Patrick Huber
- Institute for Materials and X-Ray Physics, Hamburg University of Technology, 21073 Hamburg, Germany
- Deutsches Elektronen-Synchrotron DESY, Centre for X-Ray and Nano Science CXNS, 22607 Hamburg, Germany
- Centre for Hybrid Nanostructures CHyN, Hamburg University, 22761 Hamburg, Germany
| | - Robert H. Meißner
- Institute of Polymers and Composites, Hamburg University of Technology, 21073 Hamburg, Germany
- Helmholtz-Zentrum Hereon, Institute of Surface Science, 21502 Geesthacht, Germany
| |
Collapse
|
25
|
Philippi F, Goloviznina K, Gong Z, Gehrke S, Kirchner B, Pádua AAH, Hunt PA. Charge transfer and polarisability in ionic liquids: a case study. Phys Chem Chem Phys 2022; 24:3144-3162. [PMID: 35040843 DOI: 10.1039/d1cp04592j] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The practical use of ionic liquids (ILs) is benefiting from a growing understanding of the underpinning structural and dynamic properties, facilitated through classical molecular dynamics (MD) simulations. The predictive and explanatory power of a classical MD simulation is inextricably linked to the underlying force field. A key aspect of the forcefield for ILs is the ability to recover charge based interactions. Our focus in this paper is on the description and recovery of charge transfer and polarisability effects, demonstrated through MD simulations of the widely used 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C4C1im][NTf2] IL. We study the charge distributions generated by a range of ab initio methods, and present an interpolation method for determining atom-wise scaled partial charges. Two novel methods for determining the mean field (total) charge transfer from anion to cation are presented. The impact of using different charge models and different partial charge scaling (unscaled, uniformly scaled, atom-wise scaled) are compared to fully polarisable simulations. We study a range of Drude particle explicitly polarisable potentials and shed light on the performance of current approaches to counter known problems. It is demonstrated that small changes in the charge description and MD methodology can have a significant impact; biasing some properties, while leaving others unaffected within the structural and dynamic domains.
Collapse
Affiliation(s)
- Frederik Philippi
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, London W12 0BZ, UK
| | - Kateryna Goloviznina
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| | - Zheng Gong
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| | - Sascha Gehrke
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4+6, D-53115 Bonn, Germany.,Department of Physics and Astronomy, University College London, London, WC1E 6BT, UK
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| | - Agílio A H Pádua
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| | - Patricia A Hunt
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, London W12 0BZ, UK.,School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington, New Zealand.
| |
Collapse
|
26
|
Atwi R, Chen Y, Han KS, Mueller KT, Murugesan V, Rajput NN. An automated framework for high-throughput predictions of NMR chemical shifts within liquid solutions. NATURE COMPUTATIONAL SCIENCE 2022; 2:112-122. [PMID: 38177518 DOI: 10.1038/s43588-022-00200-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 01/20/2022] [Indexed: 01/06/2024]
Abstract
Identifying stable speciation in multi-component liquid solutions is fundamentally important to areas from electrochemistry to organic chemistry and biomolecular systems. Here we introduce a fully automated, high-throughput computational framework for the accurate prediction of stable species in liquid solutions by computing the nuclear magnetic resonance (NMR) chemical shifts. The framework automatically extracts and categorizes hundreds of thousands of atomic clusters from classical molecular dynamics simulations, identifies the most stable species in solution and calculates their NMR chemical shifts via density functional theory calculations. Additionally, the framework creates a database of computed chemical shifts for liquid solutions across a wide chemical and parameter space. We compare our computational results to experimental measurements for magnesium bis(trifluoromethanesulfonyl)imide Mg(TFSI)2 salt in dimethoxyethane solvent. Our analysis of the Mg2+ solvation structural evolutions reveals key factors that influence the accuracy of NMR chemical shift predictions in liquid solutions. Furthermore, we show how the framework reduces the performance of over 300 13C and 600 1H density functional theory chemical shift predictions to a single submission procedure.
Collapse
Affiliation(s)
- Rasha Atwi
- Department of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - Ying Chen
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Kee Sung Han
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Karl T Mueller
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Vijayakumar Murugesan
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Nav Nidhi Rajput
- Department of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, NY, USA.
| |
Collapse
|
27
|
Maglia de Souza R, Karttunen M, Ribeiro MCC. Fine-Tuning the Polarizable CL&Pol Force Field for the Deep Eutectic Solvent Ethaline. J Chem Inf Model 2021; 61:5938-5947. [PMID: 34797679 DOI: 10.1021/acs.jcim.1c01181] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
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.
Collapse
Affiliation(s)
- Rafael Maglia de Souza
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Avenida Professor Lineu Prestes 748, São Paulo 05508-070, Brazil
| | - Mikko Karttunen
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 3K7, Canada.,Centre for Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoi Prospect V.O. 31, St. Petersburg 199004, Russia
| | - Mauro Carlos Costa Ribeiro
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Avenida Professor Lineu Prestes 748, São Paulo 05508-070, Brazil
| |
Collapse
|
28
|
Kognole AA, Lee J, Park SJ, Jo S, Chatterjee P, Lemkul JA, Huang J, MacKerell AD, Im W. CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field. J Comput Chem 2021; 43:359-375. [PMID: 34874077 DOI: 10.1002/jcc.26795] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Revised: 11/10/2021] [Accepted: 11/25/2021] [Indexed: 12/18/2022]
Abstract
Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. In addition, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. The capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.
Collapse
Affiliation(s)
- Abhishek A Kognole
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Jumin Lee
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Sang-Jun Park
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Sunhwan Jo
- Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois, USA
| | - Payal Chatterjee
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia, USA
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Zhejiang, Hangzhou, China
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| |
Collapse
|
29
|
Tiribocchi A, Montessori A, Lauricella M, Bonaccorso F, Brown KA, Succi S. Microscale modelling of dielectrophoresis assembly processes. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200407. [PMID: 34455845 DOI: 10.1098/rsta.2020.0407] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/04/2021] [Indexed: 06/13/2023]
Abstract
This work presents a microscale approach for simulating the dielectrophoresis assembly of polarizable particles under an external electric field. The model is shown to capture interesting dynamical and topological features, such as the formation of chains of particles and their incipient aggregation into hierarchical structures. A quantitative characterization in terms of the number and size of these structures is also discussed. This computational model could represent a viable numerical tool to study the mechanical properties of particle-based hierarchical materials and suggest new strategies for enhancing their design and manufacture. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.
Collapse
Affiliation(s)
- A Tiribocchi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
| | - A Montessori
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
| | - M Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
| | - F Bonaccorso
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
- Department of Physics and INFN, University of Rome 'Tor Vergata', Via della Ricerca Scientifica, 1 00133 Rome, Italy
| | - K A Brown
- Department of Mechanical Engineering, Boston University, Boston, MA 02215, USA
| | - S Succi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
- Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| |
Collapse
|
30
|
Klajmon M, Červinka C. Does Explicit Polarizability Improve Simulations of Phase Behavior of Ionic Liquids? J Chem Theory Comput 2021; 17:6225-6239. [PMID: 34520200 DOI: 10.1021/acs.jctc.1c00518] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Martin Klajmon
- Department of Physical Chemistry, Faculty of Chemical Engineering, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| | - Ctirad Červinka
- Department of Physical Chemistry, Faculty of Chemical Engineering, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| |
Collapse
|
31
|
Kadaoluwa Pathirannahalage SP, Meftahi N, Elbourne A, Weiss ACG, McConville CF, Padua A, Winkler DA, Costa Gomes M, Greaves TL, Le TC, Besford QA, Christofferson AJ. Systematic Comparison of the Structural and Dynamic Properties of Commonly Used Water Models for Molecular Dynamics Simulations. J Chem Inf Model 2021; 61:4521-4536. [PMID: 34406000 DOI: 10.1021/acs.jcim.1c00794] [Citation(s) in RCA: 64] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Water is a unique solvent that is ubiquitous in biology and present in a variety of solutions, mixtures, and materials settings. It therefore forms the basis for all molecular dynamics simulations of biological phenomena, as well as for many chemical, industrial, and materials investigations. Over the years, many water models have been developed, and it remains a challenge to find a single water model that accurately reproduces all experimental properties of water simultaneously. Here, we report a comprehensive comparison of structural and dynamic properties of 30 commonly used 3-point, 4-point, 5-point, and polarizable water models simulated using consistent settings and analysis methods. For the properties of density, coordination number, surface tension, dielectric constant, self-diffusion coefficient, and solvation free energy of methane, models published within the past two decades consistently show better agreement with experimental values compared to models published earlier, albeit with some notable exceptions. However, no single model reproduced all experimental values exactly, highlighting the need to carefully choose a water model for a particular study, depending on the phenomena of interest. Finally, machine learning algorithms quantified the relationship between the water model force field parameters and the resulting bulk properties, providing insight into the parameter-property relationship and illustrating the challenges of developing a water model that can accurately reproduce all properties of water simultaneously.
Collapse
Affiliation(s)
- Sachini P Kadaoluwa Pathirannahalage
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia.,Laboratoire de Chimie, Ecole Normale Supérieure de Lyon, CNRS, Lyon 69342, France
| | - Nastaran Meftahi
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Aaron Elbourne
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Alessia C G Weiss
- Leibniz-Institut für Polymerforschung e.V., Hohe Straße 6, 01069 Dresden, Germany
| | - Chris F McConville
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia.,Institute for Frontier Materials, Deakin University, Geelong, Victoria 3220, Australia
| | - Agilio Padua
- Laboratoire de Chimie, Ecole Normale Supérieure de Lyon, CNRS, Lyon 69342, France
| | - David A Winkler
- School of Biochemistry and Genetics, La Trobe Institute for Molecular Science, La Trobe University, Bundoora, Victoria 3086, Australia.,Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria 3052, Australia.,School of Pharmacy, University of Nottingham, Nottingham NG7 2QL, U.K
| | | | - Tamar L Greaves
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Tu C Le
- School of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
| | - Quinn A Besford
- Leibniz-Institut für Polymerforschung e.V., Hohe Straße 6, 01069 Dresden, Germany
| | - Andrew J Christofferson
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia.,ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| |
Collapse
|
32
|
Goloviznina K, Gong Z, Padua AAH. The
CL
&Pol polarizable force field for the simulation of ionic liquids and eutectic solvents. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1572] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
| | - Zheng Gong
- Laboratoire de Chimie École Normale Supérieure de Lyon & CNRS Lyon France
| | - Agilio A. H. Padua
- Laboratoire de Chimie École Normale Supérieure de Lyon & CNRS Lyon France
| |
Collapse
|
33
|
Molecular dynamics simulation of condensation phenomenon of nanofluid on different roughness surfaces in the presence of hydrophilic and hydrophobic structures. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116036] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
34
|
Goloviznina K, Gong Z, Costa Gomes MF, Pádua AAH. Extension of the CL&Pol Polarizable Force Field to Electrolytes, Protic Ionic Liquids, and Deep Eutectic Solvents. J Chem Theory Comput 2021; 17:1606-1617. [PMID: 33555860 DOI: 10.1021/acs.jctc.0c01002] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The polarizable CL&Pol force field presented in our previous study, Transferable, Polarizable Force Field for Ionic Liquids (J. Chem. Theory Comput. 2019, 15, 5858, DOI: http://doi.org/10.1021/acs.jctc.9b0068910.1021/acs.jctc.9b00689), is extended to electrolytes, protic ionic liquids (PIL), deep eutectic solvents (DES), and glycols. These systems are problematic in polarizable simulations because they contain either small, highly charged ions or strong hydrogen bonds, which cause trajectory instabilities due to the pull exerted on the induced dipoles. We use a Tang-Toennies (TT) function to dampen, or smear, the interactions between charges and induced dipole at a short range involving small, highly charged atoms (such as hydrogen or lithium), thus preventing the "polarization catastrophe". The new force field gives stable trajectories and is validated through comparison with experimental data on density, viscosity, and ion diffusion coefficients of liquid systems of the above-mentioned classes. The results also shed light on the hydrogen-bonding pattern in ethylammonium nitrate, a PIL, for which the literature contains conflicting views. We describe the implementation of the TT damping function, of the temperature-grouped Nosé-Hoover thermostat for polarizable molecular dynamics (MD) and of the periodic perturbation method for viscosity evaluation from non-equilibrium trajectories in the LAMMPS MD code. The main result of this work is the wider applicability of the CL&Pol polarizable force field to new, important classes of fluids, achieving robust trajectories and a good description of equilibrium and transport properties in challenging systems. The fragment-based approach of CL&Pol will allow ready extension to a wide variety of PILs, DES, and electrolytes.
Collapse
Affiliation(s)
- Kateryna Goloviznina
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| | - Zheng Gong
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| | | | - Agílio A H Pádua
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| |
Collapse
|
35
|
Misra RP, Blankschtein D. Uncovering a Universal Molecular Mechanism of Salt Ion Adsorption at Solid/Water Interfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:722-733. [PMID: 33395299 DOI: 10.1021/acs.langmuir.0c02829] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Solid/water interfaces, in which salt ions come in close proximity to solids, are ubiquitous in nature. Because water is a polar solvent and salt ions are charged, a long-standing puzzle involving solid/water interfaces is how do the electric fields exerted by the salt ions and the interfacial water molecules polarize the charge distribution in the solid and how does this polarization, in turn, influence ion adsorption at any solid/water interface. Here, using state-of-the-art polarizable force fields derived from quantum chemical simulations, we perform all-atomistic molecular dynamics simulations to investigate the adsorption of various ions comprising the well-known Hofmeister series at the graphene/water interface, including comparing with available experimental data. Our findings reveal that, in vacuum, the ionic electric field-induced polarization of graphene results in a significantly large graphene-ion polarization energy, which drives all salt ions to adsorb to graphene. On the contrary, in the presence of water molecules, we show that the ions and the water molecules exert waves of molecular electric fields on graphene which destructively interfere with each other. This remarkable phenomenon is shown to cause a water-mediated screening of more than 85% of the graphene-ion polarization energy. Finally, by investigating superhydrophilic and superhydrophobic model surfaces, we demonstrate that this phenomenon occurs universally at all solid/water interfaces and results in a significant weakening of the ion-solid interactions, such that ion specific effects are governed primarily by a competition between the ion-water and water-water interactions, irrespective of the nature of the solid/water interface.
Collapse
Affiliation(s)
- Rahul Prasanna Misra
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
36
|
Zhu J, Huang J. Methylguanidinium at the Air/Water Interface: A Simulation Study with the Drude Polarizable Force Field. J Phys Chem B 2021; 125:393-405. [PMID: 33373260 DOI: 10.1021/acs.jpcb.0c08556] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Methylguanidinium is an important molecular ion, which also serves as the model compound for the arginine side chain. We studied the structure and dynamics of the methylguanidium ion at the air/water interface by molecular dynamics simulations employing the Drude polarizable force field. We found out that methylguanidinium accumulated at the interface, with a majority adopting tilted conformations. We also demonstrated that methylguanidinium and guanidinium ions had different preference toward the air/water interface. Detailed analysis of induced dipole moments showed how ions adjusted their charge distribution at the interface and revealed how the anisotropy in molecular polarizability impacted the orientation of molecular ions. Our results illustrate the importance of explicitly including the electronic polarization effects in modeling interfacial properties.
Collapse
Affiliation(s)
- Jian Zhu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China.,Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China.,Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
37
|
Hogan A, Space B. Next-Generation Accurate, Transferable, and Polarizable Potentials for Material Simulations. J Chem Theory Comput 2020; 16:7632-7644. [PMID: 33251798 DOI: 10.1021/acs.jctc.0c00837] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
PHAHST (potentials with high accuracy, high speed, and transferability) intermolecular potential energy functions have been developed from first principles for H2, N2, the noble gases, and a metal-organic material, HKUST-1. The potentials are designed from the outset to be transferable to heterogeneous environments including porous materials, interfaces, and material simulations. This is accomplished by theoretically justified choices for all functional forms, parameters, and mixing rules, including explicit polarization in every environment and fitting to high quality electronic structure calculations using methods that are tractable for real systems. The models have been validated in neat systems by comparison to second virial coefficients and bulk pressure-density isotherms. For inhomogeneous applications, our main target, comparisons are presented to previously published experimental studies on the metal-organic material HKUST-1 including adsorption, isosteric heats of adsorption, binding site locations, and binding site energies. A systematic prescription is provided for developing compatible potentials for additional small molecules and materials. The resulting models are recommended for use in complex heterogeneous simulations where existing potentials may be inadequate.
Collapse
Affiliation(s)
- Adam Hogan
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, Florida 33620-5250, United States
| | - Brian Space
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, Florida 33620-5250, United States
| |
Collapse
|
38
|
Sharma S, Emerson MS, Wu F, Wang H, Maginn EJ, Margulis CJ. SEM-Drude Model for the Accurate and Efficient Simulation of MgCl 2-KCl Mixtures in the Condensed Phase. J Phys Chem A 2020; 124:7832-7842. [PMID: 32841019 DOI: 10.1021/acs.jpca.0c06721] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
There is a long history of models that to different extents reproduce structural and dynamical properties of high-temperature molten salts. Whereas rigid ion models can work fairly well for some of the monovalent salts, polarizability is fundamentally important when small divalent or multivalent cations are combined with significantly polarizable anions such as Cl- to form networked liquids that display a first sharp diffraction peak. There are excellent polarizable ion models (PIMs) for these systems, but there has been little success with the less expensive Core-Shell type models, which are often described as unwieldy or difficult to fit. In this article, we present the Sharma-Emerson-Margulis (SEM)-Drude model for MgCl2/KCl mixtures that with the same ingredients used in the latest and most accurate PIM models overcome the aforementioned obstacles at significantly less computational cost; structural and dynamical properties are for all practical purposes very similar to what we obtain from the PIM but typical simulations can be more than 30 times faster. This has allowed us not only to expand our recent studies on the temperature and composition dependence of intermediate range order in MgCl2/KCl mixtures but also to access transport properties that were simply too costly to properly sample in our recently published studies.
Collapse
Affiliation(s)
- Shobha Sharma
- Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| | - Matthew S Emerson
- Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| | - Fei Wu
- Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| | - Haimeng Wang
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Edward J Maginn
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Claudio J Margulis
- Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
39
|
Xu L, Izgorodina EI, Coote ML. Ordered Solvents and Ionic Liquids Can Be Harnessed for Electrostatic Catalysis. J Am Chem Soc 2020; 142:12826-12833. [PMID: 32610899 DOI: 10.1021/jacs.0c05643] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Herein, we employ classical molecular dynamics simulations using the Drude oscillator-based polarizable force field, quantum chemical calculations, and ONIOM multiscale calculations to study (a) how an external field orders the solvent environment in a chemical reaction and then (b) whether in the absence of this same applied field the ordered solvent environment alone can electrostatically catalyze a chemical reaction when compared with the corresponding disordered solvent. Our results show that a 0.2 V/Å external electric field, which is below the threshold for bond breaking of solvent molecules, leads to significant ordering of bulk methanol solvent and the ionic liquid [EMIM][BF4]. Importantly, in the absence of this same field, the ordered solvent lowers the activation energy of the hydrogen-transfer reaction of o-alkylphenyl ketones in excess of 20 kcal/mol when the solvent is methanol and by over 30 kcal/mol for [EMIM][BF4]. Even a 0.1 V/Å external field has effects of ca. 10 and 20 kcal/mol, respectively. This work suggests a possible strategy for scaling electrostatic catalysis by applying a pulsed external field to the reaction medium to maintain solvent ordering while allowing the reaction to proceed largely in the absence of an external field.
Collapse
Affiliation(s)
- Longkun Xu
- ARC Centre of Excellence for Electromaterials Science, Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory 2601, Australia
| | - Ekaterina I Izgorodina
- Monash Computational Chemistry Group, School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, Victoria 3800, Australia
| | - Michelle L Coote
- ARC Centre of Excellence for Electromaterials Science, Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory 2601, Australia
| |
Collapse
|
40
|
Ding Y, Xu Y, Qian C, Chen J, Zhu J, Huang H, Shi Y, Huang J. Predicting partition coefficients of drug-like molecules in the SAMPL6 challenge with Drude polarizable force fields. J Comput Aided Mol Des 2020; 34:421-435. [PMID: 31960252 DOI: 10.1007/s10822-020-00282-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Accepted: 01/08/2020] [Indexed: 12/14/2022]
Abstract
The water-octanol partition coefficient is an important physicochemical property for small molecule drug design. Here, we report our participation in the SAMPL6 logP prediction challenge with free energy perturbation (FEP) calculations in the water phase and in the 1-octanol phase using Drude polarizable force fields. Root mean square error (RMSE) and mean absolute error (MAE) of our prediction are equal to 1.85 and 1.25 logP units. The errors are not evenly distributed. Out of eleven SAMPL6 solutes, FEP/Drude performed very badly on three molecules (deviations all larger than 2 logP units) but good on the remaining eight (deviations all less than 1 logP unit). We find while FEP converges well within one nanosecond in water, simulations in 1-octanol need much longer simulation time and possibly more independent runs for sampling. We also find out that 1-octanol, albeit being a non-polar solvent, still polarizes solute molecules and forms stable hydrogen bonds with them. At the end, we attempt to reweight FEP trajectories with QM/Drude calculations and discuss possible caveats in our simulation setup.
Collapse
Affiliation(s)
- Ye Ding
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - You Xu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Cheng Qian
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jinfeng Chen
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jian Zhu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Houhou Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Yi Shi
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China. .,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.
| |
Collapse
|
41
|
Bernardino K, Goloviznina K, Gomes MC, Pádua AAH, Ribeiro MCC. Ion pair free energy surface as a probe of ionic liquid structure. J Chem Phys 2020; 152:014103. [DOI: 10.1063/1.5128693] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Affiliation(s)
- Kalil Bernardino
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, São Paulo 05508-000, Brazil
| | - Kateryna Goloviznina
- Laboratoire de Chimie, ENS de Lyon, Université de Lyon, CNRS, 69364 Lyon, France
| | | | - Agílio A. H. Pádua
- Laboratoire de Chimie, ENS de Lyon, Université de Lyon, CNRS, 69364 Lyon, France
| | - Mauro C. C. Ribeiro
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, São Paulo 05508-000, Brazil
| |
Collapse
|
42
|
Massaro A, Avila J, Goloviznina K, Rivalta I, Gerbaldi C, Pavone M, Costa Gomes MF, Padua AAH. Sodium diffusion in ionic liquid-based electrolytes for Na-ion batteries: the effect of polarizable force fields. Phys Chem Chem Phys 2020; 22:20114-20122. [DOI: 10.1039/d0cp02760j] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Understanding the transport of sodium ions in ionic liquids is the key to design novel electrolyte materials for sodium-ion batteries.
Collapse
Affiliation(s)
- Arianna Massaro
- Laboratoire de Chimie
- ENS de Lyon
- CNRS
- Université de Lyon
- 69364 Lyon
| | - Jocasta Avila
- Laboratoire de Chimie
- ENS de Lyon
- CNRS
- Université de Lyon
- 69364 Lyon
| | | | - Ivan Rivalta
- Laboratoire de Chimie
- ENS de Lyon
- CNRS
- Université de Lyon
- 69364 Lyon
| | - Claudio Gerbaldi
- GAME Lab
- Department of Applied Science and Technology (DISAT)
- Politecnico di Torino
- 10129 Torino
- Italy
| | - Michele Pavone
- Department of Chemical Sciences
- Università di Napoli “Federico II”
- 80126 Naples
- Italy
| | | | | |
Collapse
|
43
|
Clark R, von Domaros M, McIntosh AJS, Luzar A, Kirchner B, Welton T. Effect of an external electric field on the dynamics and intramolecular structures of ions in an ionic liquid. J Chem Phys 2019; 151:164503. [PMID: 31675893 DOI: 10.1063/1.5129367] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Simulations of the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide in an external electric field have been performed using a Drude particle polarizable force field. The structure of the ions has been analyzed, with close attention paid to the configurations of the ions. The "charge arm" concept is used to explain certain changes of these configurations that would be difficult to rationalize otherwise, e.g., trans → cis isomerization of the bis(trifluoromethylsulfonyl)imide anion and extension of the alkyl chain of the cation. It has also been shown that the ions orient themselves so that their charge arms align with and stretch out along the field, and these effects occur at lower external electric field strengths than cause a change in the inherent diffusion of the ions. The dynamics of the system parallel and perpendicular to the field were analyzed, and it was found that the applied field affected the diffusion normal to the field. This is explained as a secondary effect of a change in the ion cage formed by the surrounding counterions of a given ion in the ionic liquid. The breakdown of the ion cages was rationalized by correlating changes to the inherent diffusion of the ions with other changes to the diffusion and bulk structure of the liquid, as well as considering the average forces on the ions compared to the force the ions would be expected to experience in an electric field. Parallel to the field, a drift was observed at every electric field studied. In electric fields with no changes to the ion cage structure, the relationship between the drift and electric field was found to be linear, becoming nonlinear as the ion cage structure breaks down.
Collapse
Affiliation(s)
- Ryan Clark
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, London W12 0BZ, United Kingdom
| | - Michael von Domaros
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| | - Alastair J S McIntosh
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, London W12 0BZ, United Kingdom
| | - Alenka Luzar
- Department of Chemistry, Virginia Commonwealth University, 1001 West Main St., Richmond, Virginia 23284-2006, USA
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| | - Tom Welton
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, London W12 0BZ, United Kingdom
| |
Collapse
|
44
|
Bernardino K, Lima TA, Ribeiro MCC. Low-Temperature Phase Transitions of the Ionic Liquid 1-Ethyl-3-methylimidazolium Dicyanamide. J Phys Chem B 2019; 123:9418-9427. [DOI: 10.1021/acs.jpcb.9b07654] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kalil Bernardino
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05508-000 São Paulo, São Paulo, Brazil
| | - Thamires A. Lima
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05508-000 São Paulo, São Paulo, Brazil
| | - Mauro C. C. Ribeiro
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05508-000 São Paulo, São Paulo, Brazil
| |
Collapse
|
45
|
Goloviznina K, Canongia Lopes JN, Costa Gomes M, Pádua AAH. Transferable, Polarizable Force Field for Ionic Liquids. J Chem Theory Comput 2019; 15:5858-5871. [DOI: 10.1021/acs.jctc.9b00689] [Citation(s) in RCA: 77] [Impact Index Per Article: 15.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kateryna Goloviznina
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| | - José N. Canongia Lopes
- Centro de Química Estrutural, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
| | | | - Agílio A. H. Pádua
- Laboratoire de Chimie, École Normale Supérieure de Lyon & CNRS, 69364 Lyon, France
| |
Collapse
|
46
|
Červinka C, Klajmon M, Štejfa V. Cohesive Properties of Ionic Liquids Calculated from First Principles. J Chem Theory Comput 2019; 15:5563-5578. [DOI: 10.1021/acs.jctc.9b00625] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Martin Klajmon
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| |
Collapse
|
47
|
Bedrov D, Piquemal JP, Borodin O, MacKerell AD, Roux B, Schröder C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem Rev 2019; 119:7940-7995. [PMID: 31141351 PMCID: PMC6620131 DOI: 10.1021/acs.chemrev.8b00763] [Citation(s) in RCA: 284] [Impact Index Per Article: 56.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2018] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Dmitry Bedrov
- Department
of Materials Science & Engineering, University of Utah, 122 South Central Campus Drive, Room 304, Salt Lake City, Utah 84112, United States
| | - Jean-Philip Piquemal
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, CC137, 4 Place Jussieu, Tour 12-13, 4ème étage, 75252 Paris Cedex 05, France
- Institut
Universitaire de France, 75005, Paris Cedex 05, France
- Department
of Biomedical Engineering, The University
of Texas at Austin, Austin, Texas 78712, United States
| | - Oleg Borodin
- Electrochemistry
Branch, Sensors and Electron Devices Directorate, Army Research Laboratory, 2800 Powder Mill Road, Adelphi, Maryland 20703, United
States
| | - Alexander D. MacKerell
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, United
States
| | - Benoît Roux
- Department
of Biochemistry and Molecular Biology, Gordon Center for Integrative
Science, University of Chicago, 929 57th Street, Chicago, Illinois 60637, United States
| | - Christian Schröder
- Department
of Computational Biological Chemistry, University
of Vienna, Währinger Strasse 17, A-1090 Vienna, Austria
| |
Collapse
|
48
|
Gaudin T, Ma H. A molecular contact theory for simulating polarization: application to dielectric constant prediction. Phys Chem Chem Phys 2019; 21:14846-14857. [PMID: 31232397 DOI: 10.1039/c9cp02358e] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Microscopic polarization in liquids, which is challenging to account for intuitively and quantitatively, can impact the behavior of liquids in numerous ways and thus is ubiquitous in a broad range of domains and applications. To overcome this challenge, in this work, a molecular contact theory was proposed as a proxy to simulate microscopic polarization in liquids. In particular, molecular surfaces from implicit solvation models were used to predict both the dipole moment of individual molecules and mutual orientations arising from contacts between molecules. Then, the calculated dipole moments and orientations were combined in an analytical coupling, which allowed for the prediction of effective (polarized) dipole moments for all distinct species in the liquid. As a proof-of-concept, the model focused on predicting the dielectric constant and was tested on 420 pure liquids, 269 binary organic mixtures (3792 individual compositions) and 46 aqueous mixtures (704 individual compositions). The model proved to be flexible enough to reach an unprecedented satisfactory mean relative error of about 16-22% and a classification accuracy of 84-90% within four meaningful classes of weak, low average, high average and strong dielectric constants. The method also proved to be computationally very efficient, with calculation times ranging from a few seconds to about ten minutes on a personal computer with a single CPU. This success demonstrates that much of the microscopic polarization concept can be satisfactorily described based on a simple molecular contact theory. Moreover, the new model for dielectric constants provides a useful alternative to computationally expensive molecular dynamics simulations for large scale virtual screenings in chemical engineering and material sciences.
Collapse
Affiliation(s)
- Théophile Gaudin
- Key Laboratory of Mesoscopic Chemistry of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, China.
| | - Haibo Ma
- Key Laboratory of Mesoscopic Chemistry of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, China.
| |
Collapse
|
49
|
Seeger ZL, Kobayashi R, Izgorodina EI. Cluster approach to the prediction of thermodynamic and transport properties of ionic liquids. J Chem Phys 2018; 148:193832. [DOI: 10.1063/1.5009791] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Zoe L. Seeger
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, VIC 3800, Australia
| | - Rika Kobayashi
- Australian National University, Leonard Huxley Building 56, Mills Road, Canberra, ACT 2601, Australia
| | | |
Collapse
|
50
|
McDaniel JG. Polarization Effects in Binary [BMIM+][BF4–]/1,2-Dichloroethane, Acetone, Acetonitrile, and Water Electrolytes. J Phys Chem B 2018; 122:4345-4355. [DOI: 10.1021/acs.jpcb.8b01714] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jesse G. McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|