1
|
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
|
2
|
Chen L, Guo Y, Xu Z, Yang X. Multiscale Simulation of the Interaction and Adsorption of Ions on a Hydrophobic Graphene Surface. Chemphyschem 2018; 19:2954-2960. [PMID: 30142233 DOI: 10.1002/cphc.201800428] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2018] [Indexed: 11/07/2022]
Abstract
The adsorption of ions on a graphene surface is very important to control relevant graphene-based processes. In this work, a multiscale simulation was carried out to study the adsorption of Na+ /Cl- ions on graphene by combining quantum mechanics calculations and molecular dynamics (MD) simulations. The interaction energies of the ions with graphene were computed using density functional theory (DFT). It was found that the ions show strong interaction with a graphene cluster and the overwhelming portion of the interaction energy is the ion-π orbital interaction. The large orbital interaction can be ascribed to the two contributions arising from the ion-induced polarization of graphene and the charge transfer between ion and graphene. Their different contribution degrees reveal that the polarization effect plays a main role on the orbital interaction for ion adsorption. Comparatively, for Na/Cl atom adsorption, the charge transfer shows large part to the orbital interaction with weak atom-induced polarization. The obtained interaction energies were applied to develop new interaction potentials between ion and graphene, and then MD simulations were used to study the interfacial adsorption behavior of Na+ /Cl- aqueous solution onto the graphene surface. Due to enhanced ion-π interactions, Na+ /Cl- cooperatively demonstrates a strong ion adsorption layer through direct contact with the hydrophobic graphene surface. Our simulation result presents a new understanding of ion-graphene interactions.
Collapse
Affiliation(s)
- Luohao Chen
- State Key Laboratory of Materials-Oriented Chemical Engineering College of Chemical Engineering, Nanjing Tech University, Nanjing, 210009, China
| | - Yanan Guo
- State Key Laboratory of Materials-Oriented Chemical Engineering College of Chemical Engineering, Nanjing Tech University, Nanjing, 210009, China
| | - Zhijun Xu
- State Key Laboratory of Materials-Oriented Chemical Engineering College of Chemical Engineering, Nanjing Tech University, Nanjing, 210009, China
| | - Xiaoning Yang
- State Key Laboratory of Materials-Oriented Chemical Engineering College of Chemical Engineering, Nanjing Tech University, Nanjing, 210009, China
| |
Collapse
|
3
|
Chaban VV, Prezhdo OV. Ionic Vapor Composition in Critical and Supercritical States of Strongly Interacting Ionic Compounds. J Phys Chem B 2016; 120:4302-9. [DOI: 10.1021/acs.jpcb.6b02405] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Vitaly V. Chaban
- Instituto
de Ciência e Tecnologia, Universidade Federal de São Paulo, 12231-280, São José dos Campos, SP, Brazil
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
4
|
Colherinhas G, Fileti EE, Chaban VV. Can inorganic salts tune electronic properties of graphene quantum dots? Phys Chem Chem Phys 2015; 17:17413-20. [DOI: 10.1039/c5cp02083b] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
In this work, we apply density functional theory to study the effect of neutral ionic clusters adsorbed on the GQD surface. We conclude that both the HOMO and the LUMO of GQDs are very sensitive to the presence of ions and to their distance from the GQD surface. However, the alteration of the band gap itself is modest, as opposed to the case of free ions (recent reports). Our work fosters progress in modulating electronic properties of nanoscale carbonaceous materials.
Collapse
Affiliation(s)
| | - Eudes Eterno Fileti
- Instituto de Ciência e Tecnologia
- Universidade Federal de São Paulo
- São José dos Campos
- Brazil
| | - Vitaly V. Chaban
- Instituto de Ciência e Tecnologia
- Universidade Federal de São Paulo
- São José dos Campos
- Brazil
| |
Collapse
|
5
|
Choi E, McDaniel JG, Schmidt JR, Yethiraj A. First-Principles, Physically Motivated Force Field for the Ionic Liquid [BMIM][BF4]. J Phys Chem Lett 2014; 5:2670-2674. [PMID: 26277961 DOI: 10.1021/jz5010945] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Molecular simulations play an important role in establishing structure-property relations in complex fluids such as room-temperature ionic liquids. Classical force fields are the starting point when large systems or long times are of interest. These force fields must be not only accurate but also transferable. In this work, we report a physically motivated force field for the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]) based on symmetry-adapted perturbation theory. The predictions (from molecular dynamics simulations) of the liquid density, enthalpy of vaporization, diffusion coefficients, viscosity, and conductivity are in excellent agreement with experiment, with no adjustable parameters. The explicit energy decomposition inherent in the force field enables a quantitative analysis of the important physical interactions in these systems. We find that polarization is crucial and there is little evidence of charge transfer. We also argue that the often used procedure of scaling down charges in molecular simulations of ionic liquids is unphysical for [BMIM][BF4]. Because all intermolecular interactions in the force field are parametrized from first-principles, we anticipate good transferability to other ionic liquid systems and physical conditions.
Collapse
Affiliation(s)
- Eunsong Choi
- †Department of Physics, University of Wisconsin, 1150 University Avenue, Madison, Wisconsin 53706, United States
| | - Jesse G McDaniel
- ‡Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - J R Schmidt
- ‡Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- ‡Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|