1
|
Wang X, Wang Y, Kuang Y, Le JB. Understanding the Effects of Electrode Material, Single Crystal Facet, and Electrolyte Ion on the Helmholtz Capacitance of Metal/Aqueous Solution Interfaces. J Phys Chem Lett 2023; 14:7833-7839. [PMID: 37624858 DOI: 10.1021/acs.jpclett.3c02108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/27/2023]
Abstract
The comprehensive interpretation of the measured differential Helmholtz capacitance curve is vital for advancing our understanding of the interfacial structure. While several possible physical effects contributing to the Helmholtz capacitance have been proposed theoretically, combining those factors to explain the experimentally observed potential-dependent capacitance profile remains a significant challenge. In this study, we employ ab initio molecular dynamics simulations to model various metal/solution interfaces. Our investigation primarily emphasizes the substantial effect of water chemisorption on the potential-dependent behavior of the Helmholtz capacitance. Additionally, we identify other critical factors that profoundly impact the Helmholtz capacitance: (1) Ions with low hydration energy hinder the availability of surface sites for water adsorption, resulting in a diminished enhancement of capacitance from water chemisorption. (2) Using large-sized ions leads to an expansion of the Helmholtz layer, causing a decrease in the Helmholtz capacitance. (3) Metal surfaces with higher affinity for water attract water adsorption at lower potentials, resulting in a lower peak potential for the differential Helmholtz capacitance curve.
Collapse
Affiliation(s)
- Xue Wang
- Key Laboratory of Advanced Fuel Cells and Electrolyzers Technology of Zhejiang Province, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo 315201, China
| | - Ying Wang
- Institute of Materials, Ningbo University of Technology, Ningbo 315211, China
| | - Yongbo Kuang
- Key Laboratory of Advanced Fuel Cells and Electrolyzers Technology of Zhejiang Province, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo 315201, China
| | - Jia-Bo Le
- Key Laboratory of Advanced Fuel Cells and Electrolyzers Technology of Zhejiang Province, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo 315201, China
| |
Collapse
|
2
|
Shin SJ, Gittins JW, Golomb MJ, Forse AC, Walsh A. Microscopic Origin of Electrochemical Capacitance in Metal-Organic Frameworks. J Am Chem Soc 2023; 145:14529-14538. [PMID: 37341453 PMCID: PMC10326873 DOI: 10.1021/jacs.3c04625] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2023] [Indexed: 06/22/2023]
Abstract
Electroconductive metal-organic frameworks (MOFs) have emerged as high-performance electrode materials for supercapacitors, but the fundamental understanding of the underlying chemical processes is limited. Here, the electrochemical interface of Cu3(HHTP)2 (HHTP = 2,3,6,7,10,11-hexahydroxytriphenylene) with an organic electrolyte is investigated using a multiscale quantum-mechanics/molecular-mechanics (QM/MM) procedure and experimental electrochemical measurements. Our simulations reproduce the observed capacitance values and reveals the polarization phenomena of the nanoporous framework. We find that excess charges mainly form on the organic ligand, and cation-dominated charging mechanisms give rise to greater capacitance. The spatially confined electric double-layer structure is further manipulated by changing the ligand from HHTP to HITP (HITP = 2,3,6,7,10,11-hexaiminotriphenylene). This minimal change to the electrode framework not only increases the capacitance but also increases the self-diffusion coefficients of in-pore electrolytes. The performance of MOF-based supercapacitors can be systematically controlled by modifying the ligating group.
Collapse
Affiliation(s)
- Seung-Jae Shin
- Department
of Materials Science and Engineering, Yonsei
University, Seoul 03722, Korea
| | - Jamie W. Gittins
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, U.K.
| | - Matthias J. Golomb
- Thomas
Young Centre and Department of Materials, Imperial College London, London SW7 2AZ, U.K.
| | - Alexander C. Forse
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, U.K.
| | - Aron Walsh
- Thomas
Young Centre and Department of Materials, Imperial College London, London SW7 2AZ, U.K.
- Department
of Physics, Ewha Womans University, Seoul 03760, Korea
| |
Collapse
|
3
|
Oh H, Shin SJ, Choi E, Yamagishi H, Ohta T, Yabuuchi N, Jung HG, Kim H, Byon HR. Anion-Induced Interfacial Liquid Layers on LiCoO 2 in Salt-in-Water Lithium-Ion Batteries. JACS AU 2023; 3:1392-1402. [PMID: 37234123 PMCID: PMC10206596 DOI: 10.1021/jacsau.3c00061] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/08/2023] [Revised: 04/17/2023] [Accepted: 04/18/2023] [Indexed: 05/27/2023]
Abstract
The incompatibility of lithium intercalation electrodes with water has impeded the development of aqueous Li-ion batteries. The key challenge is protons which are generated by water dissociation and deform the electrode structures through intercalation. Distinct from previous approaches utilizing large amounts of electrolyte salts or artificial solid-protective films, we developed liquid-phase protective layers on LiCoO2 (LCO) using a moderate concentration of 0.5∼3 mol kg-1 lithium sulfate. Sulfate ion strengthened the hydrogen-bond network and easily formed ion pairs with Li+, showing strong kosmotropic and hard base characteristics. Our quantum mechanics/molecular mechanics (QM/MM) simulations revealed that sulfate ion paired with Li+ helped stabilize the LCO surface and reduced the density of free water in the interface region below the point of zero charge (PZC) potential. In addition, in situ electrochemical surface-enhanced infrared absorption spectroscopy (SEIRAS) proved the appearance of inner-sphere sulfate complexes above the PZC potential, serving as the protective layers of LCO. The role of anions in stabilizing LCO was correlated with kosmotropic strength (sulfate > nitrate > perchlorate > bistriflimide (TFSI-)) and explained better galvanostatic cyclability in LCO cells.
Collapse
Affiliation(s)
- Hyunjeong Oh
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
- Advanced
Battery Center, KAIST Institute for NanoCentury, Daejeon 34141, Republic of Korea
| | - Seung-Jae Shin
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Eunjin Choi
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
- Advanced
Battery Center, KAIST Institute for NanoCentury, Daejeon 34141, Republic of Korea
| | - Hirona Yamagishi
- Sychrotron
Radiation (SR) Center, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Toshiaki Ohta
- Sychrotron
Radiation (SR) Center, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Naoaki Yabuuchi
- Department
of Chemistry and Life Science, Yokohama
National University, 79-5 Tokiwadai, Hodogaya-ku, Yokohama, Kanagawa 240-8501, Japan
| | - Hun-Gi Jung
- Energy
Storage Research Center, Clean Energy Research Division, Korea Institute of Science and Technology (KIST), Seoul 02792, Republic of Korea
| | - Hyungjun Kim
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Hye Ryung Byon
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
- Advanced
Battery Center, KAIST Institute for NanoCentury, Daejeon 34141, Republic of Korea
| |
Collapse
|
4
|
Abidi N, Steinmann SN. An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces. ACS APPLIED MATERIALS & INTERFACES 2023; 15:25009-25017. [PMID: 37163568 DOI: 10.1021/acsami.3c01430] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Atomistic modeling of electrified interfaces remains a major issue for detailed insights in electrocatalysis, corrosion, electrodeposition, batteries, and related devices such as pseudocapacitors. In these domains, the use of grand-canonical density functional theory (GC-DFT) in combination with implicit solvation models has become popular. GC-DFT can be conveniently applied not only to metallic surfaces but also to semiconducting oxides and sulfides and is, furthermore, sufficiently robust to achieve a consistent description of reaction pathways. However, the accuracy of implicit solvation models for solvation effects at interfaces is in general unknown. One promising way to overcome the limitations of implicit solvents is going toward hybrid quantum mechanical (QM)/molecular mechanics (MM) models. For capturing the electrochemical potential dependence, the key quantity is the capacitance, i.e., the relation between the surface charge and the electrochemical potential. In order to retrieve the electrochemical potential from a QM/MM hybrid scheme, an electrostatic embedding is required. Furthermore, the charge of the surface and of the solvent regions has to be strictly opposite in order to consistently simulate charge-neutral unit cells in MM and in QM. To achieve such a QM/MM scheme, we present the implementation of electrostatic embedding in the VASP code. This scheme is broadly applicable to any neutral or charged solid/liquid interface. Here, we demonstrate its use in the context of GC-DFT for the hydrogen evolution reaction (HER) over a noble-metal-free electrocatalyst, MoS2. We investigate the effect of electrostatic embedding compared to the implicit solvent model for three contrasting active sites on MoS2: (i) the sulfur vacancy defect, which is rather apolar; (ii) a Mo antisite defect, where the active site is a surface bound highly polar OH group; and (iii) a reconstructed edge site, which is generally believed to be responsible for most of the catalytic activity. According to our results, the electrostatic embedding leads to almost indistinguishable results compared to the implicit solvent for the apolar system but has a significant effect on polar sites. This demonstrates the reliability of the hybrid QM/MM, electrostatically embedded solvation model for electrified interfaces.
Collapse
Affiliation(s)
- Nawras Abidi
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| | - Stephan N Steinmann
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| |
Collapse
|
5
|
Bramley GA, Beynon OT, Stishenko PV, Logsdail AJ. The application of QM/MM simulations in heterogeneous catalysis. Phys Chem Chem Phys 2023; 25:6562-6585. [PMID: 36810655 DOI: 10.1039/d2cp04537k] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
The QM/MM simulation method is provenly efficient for the simulation of biological systems, where an interplay of extensive environment and delicate local interactions drives a process of interest through a funnel on a complex energy landscape. Recent advances in quantum chemistry and force-field methods present opportunities for the adoption of QM/MM to simulate heterogeneous catalytic processes, and their related systems, where similar intricacies exist on the energy landscape. Herein, the fundamental theoretical considerations for performing QM/MM simulations, and the practical considerations for setting up QM/MM simulations of catalytic systems, are introduced; then, areas of heterogeneous catalysis are explored where QM/MM methods have been most fruitfully applied. The discussion includes simulations performed for adsorption processes in solvent at metallic interfaces, reaction mechanisms within zeolitic systems, nanoparticles, and defect chemistry within ionic solids. We conclude with a perspective on the current state of the field and areas where future opportunities for development and application exist.
Collapse
Affiliation(s)
- Gabriel Adrian Bramley
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | - Owain Tomos Beynon
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | | | - Andrew James Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| |
Collapse
|
6
|
Takahashi K, Nakano H, Sato H. Accelerated constant-voltage quantum mechanical/molecular mechanical method for molecular systems at electrochemical interfaces. J Chem Phys 2022; 157:234107. [PMID: 36550044 DOI: 10.1063/5.0128358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022] Open
Abstract
The structure and electronic properties of a molecule at an electrochemical interface are changed by interactions with the electrode surface and the electrolyte solution, which can be significantly modulated by an applied voltage. We present an efficient self-consistent quantum mechanics/molecular mechanics (QM/MM) approach to study a physisorbed molecule at a metal electrode-electrolyte interface under the constant-voltage condition. The approach employs a classical polarizable double electrode model, which enables us to study the QM/MM system in the constant-voltage ensemble. A mean-field embedding approximation is further introduced in order to overcome the difficulties associated with statistical sampling of the electrolyte configurations. The results of applying the method to a test system indicate that the adsorbed molecule is no less or slightly more polarized at the interface than in the bulk electrolyte solution. The geometry of the horizontally adsorbed molecule is modulated by their electrostatic interactions with the polarizable electrode surfaces and also the interactions with cations attracted toward the interface when the adsorbate is reduced. We also demonstrate that the approach can be used to quantitatively evaluate the reorganization energy of a one electron reduction reaction of a molecule in an electrochemical cell.
Collapse
Affiliation(s)
- Ken Takahashi
- Department of Molecular Engineering, Kyoto University, Kyoto Daigaku Katsura, Kyoto 615-8246, Japan
| | - Hiroshi Nakano
- CD-FMat, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba Central 2, Tsukuba, Ibaraki 305-8568, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Kyoto University, Kyoto Daigaku Katsura, Kyoto 615-8246, Japan
| |
Collapse
|
7
|
Kim MK, Shin S, Lee JM, Park YB, Kim YM, Kim HJ, Choi JW. Cationic Additive with a Rigid Solvation Shell for High‐Performance Zinc Ion Batteries. Angew Chem Int Ed Engl 2022; 61:e202211589. [DOI: 10.1002/anie.202211589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2022] [Indexed: 11/06/2022]
Affiliation(s)
- Minkwan K. Kim
- School of Chemical and Biological Engineering and Institute of Chemical Process Seoul National University 1 Gwanak-ro, Gwanak-gu Seoul 08826 Republic of Korea
| | - Seung‐Jae Shin
- Department of Chemistry Korea Advanced Institute of Science and Technology Daejeon 34141 Republic of Korea
| | - Jimin M. Lee
- School of Chemical and Biological Engineering and Institute of Chemical Process Seoul National University 1 Gwanak-ro, Gwanak-gu Seoul 08826 Republic of Korea
| | - Youngbin B. Park
- School of Chemical and Biological Engineering and Institute of Chemical Process Seoul National University 1 Gwanak-ro, Gwanak-gu Seoul 08826 Republic of Korea
| | - Yangmoon M. Kim
- School of Chemical and Biological Engineering and Institute of Chemical Process Seoul National University 1 Gwanak-ro, Gwanak-gu Seoul 08826 Republic of Korea
| | - Hyungjun J. Kim
- Department of Chemistry Korea Advanced Institute of Science and Technology Daejeon 34141 Republic of Korea
| | - Jang Wook Choi
- School of Chemical and Biological Engineering and Institute of Chemical Process Seoul National University 1 Gwanak-ro, Gwanak-gu Seoul 08826 Republic of Korea
| |
Collapse
|
8
|
Tran B, Milner ST, Janik MJ. Kinetics of Acid-Catalyzed Dehydration of Alcohols in Mixed Solvent Modeled by Multiscale DFT/MD. ACS Catal 2022. [DOI: 10.1021/acscatal.2c03978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bolton Tran
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania16802, United States
| | - Scott T. Milner
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania16802, United States
| | - Michael J. Janik
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania16802, United States
| |
Collapse
|
9
|
Shin SJ, Choi H, Ringe S, Won DH, Oh HS, Kim DH, Lee T, Nam DH, Kim H, Choi CH. A unifying mechanism for cation effect modulating C1 and C2 productions from CO 2 electroreduction. Nat Commun 2022; 13:5482. [PMID: 36123326 PMCID: PMC9485141 DOI: 10.1038/s41467-022-33199-8] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Accepted: 09/07/2022] [Indexed: 02/05/2023] Open
Abstract
Electrocatalysis, whose reaction venue locates at the catalyst-electrolyte interface, is controlled by the electron transfer across the electric double layer, envisaging a mechanistic link between the electron transfer rate and the electric double layer structure. A fine example is in the CO2 reduction reaction, of which rate shows a strong dependence on the alkali metal cation (M+) identity, but there is yet to be a unified molecular picture for that. Using quantum-mechanics-based atom-scale simulation, we herein scrutinize the M+-coupling capability to possible intermediates, and establish H+- and M+-associated ET mechanisms for CH4 and CO/C2H4 formations, respectively. These theoretical scenarios are successfully underpinned by Nernstian shifts of polarization curves with the H+ or M+ concentrations and the first-order kinetics of CO/C2H4 formation on the electrode surface charge density. Our finding further rationalizes the merit of using Nafion-coated electrode for enhanced C2 production in terms of enhanced surface charge density.
Collapse
Affiliation(s)
- Seung-Jae Shin
- grid.37172.300000 0001 2292 0500Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, 34141 Republic of Korea
| | - Hansol Choi
- grid.61221.360000 0001 1033 9831School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju, 61005 Republic of Korea
| | - Stefan Ringe
- grid.222754.40000 0001 0840 2678Department of Chemistry, Korea University, Seoul, 02841 Republic of Korea
| | - Da Hye Won
- grid.35541.360000000121053345Clean Energy Research Center, Korea Institute of Science and Technology, Seoul, 02792 Republic of Korea
| | - Hyung-Suk Oh
- grid.35541.360000000121053345Clean Energy Research Center, Korea Institute of Science and Technology, Seoul, 02792 Republic of Korea
| | - Dong Hyun Kim
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 Republic of Korea
| | - Taemin Lee
- grid.417736.00000 0004 0438 6721Department of Energy Science and Engineering, Daegu Gyeongbuk Institute of Science and Technology, Daegu, 42988 Republic of Korea
| | - Dae-Hyun Nam
- grid.417736.00000 0004 0438 6721Department of Energy Science and Engineering, Daegu Gyeongbuk Institute of Science and Technology, Daegu, 42988 Republic of Korea
| | - Hyungjun Kim
- grid.37172.300000 0001 2292 0500Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, 34141 Republic of Korea
| | - Chang Hyuck Choi
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 Republic of Korea
| |
Collapse
|
10
|
Paik D, Lee H, Kim H, Choi JM. Thermodynamics of π–π Interactions of Benzene and Phenol in Water. Int J Mol Sci 2022; 23:ijms23179811. [PMID: 36077201 PMCID: PMC9456204 DOI: 10.3390/ijms23179811] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 08/20/2022] [Accepted: 08/22/2022] [Indexed: 11/16/2022] Open
Abstract
The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.
Collapse
Affiliation(s)
- Dooam Paik
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
| | - Hankyul Lee
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
| | - Hyungjun Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
- Correspondence: (H.K.); (J.-M.C.)
| | - Jeong-Mo Choi
- Department of Chemistry, Pusan National University, Busan 46241, Korea
- Correspondence: (H.K.); (J.-M.C.)
| |
Collapse
|
11
|
Jang T, Paik D, Shin S, Kim H. Density functional theory in classical explicit solvents: Mean‐field
QM
/
MM
method for simulating solid–liquid interfaces. B KOREAN CHEM SOC 2022. [DOI: 10.1002/bkcs.12485] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Taehwan Jang
- Department of Chemistry Korea Advanced Institute of Science and Technology (KAIST) Daejeon Republic of Korea
| | - Dooam Paik
- Department of Chemistry Korea Advanced Institute of Science and Technology (KAIST) Daejeon Republic of Korea
| | - Seung‐Jae Shin
- Department of Chemistry Korea Advanced Institute of Science and Technology (KAIST) Daejeon Republic of Korea
| | - Hyungjun Kim
- Department of Chemistry Korea Advanced Institute of Science and Technology (KAIST) Daejeon Republic of Korea
| |
Collapse
|
12
|
Shin SJ, Kim DH, Bae G, Ringe S, Choi H, Lim HK, Choi CH, Kim H. On the importance of the electric double layer structure in aqueous electrocatalysis. Nat Commun 2022; 13:174. [PMID: 35013347 PMCID: PMC8748683 DOI: 10.1038/s41467-021-27909-x] [Citation(s) in RCA: 55] [Impact Index Per Article: 27.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2021] [Accepted: 12/07/2021] [Indexed: 11/09/2022] Open
Abstract
To design electrochemical interfaces for efficient electric-chemical energy interconversion, it is critical to reveal the electric double layer (EDL) structure and relate it with electrochemical activity; nonetheless, this has been a long-standing challenge. Of particular, no molecular-level theories have fully explained the characteristic two peaks arising in the potential-dependence of the EDL capacitance, which is sensitively dependent on the EDL structure. We herein demonstrate that our first-principles-based molecular simulation reproduces the experimental capacitance peaks. The origin of two peaks emerging at anodic and cathodic potentials is unveiled to be an electrosorption of ions and a structural phase transition, respectively. We further find a cation complexation gradually modifies the EDL structure and the field strength, which linearly scales the carbon dioxide reduction activity. This study deciphers the complex structural response of the EDL and highlights its catalytic importance, which bridges the mechanistic gap between the EDL structure and electrocatalysis.
Collapse
Affiliation(s)
- Seung-Jae Shin
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, 34141, Republic of Korea
| | - Dong Hyun Kim
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju, 61005, Republic of Korea
| | - Geunsu Bae
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju, 61005, Republic of Korea
| | - Stefan Ringe
- Department of Energy Science and Engineering, Daegu Gyeongbuk Institute of Science and Technology, Daegu, 42988, Republic of Korea.,Energy Science and Engineering Research Center, Daegu Gyeongbuk Institute of Science and Technology (DGIST), Daegu, 42988, Republic of Korea
| | - Hansol Choi
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju, 61005, Republic of Korea
| | - Hyung-Kyu Lim
- Division of Chemical Engineering and Bioengineering, Kangwon National University, Chuncheon, Gangwon-do, 24341, Republic of Korea
| | - Chang Hyuck Choi
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju, 61005, Republic of Korea.
| | - Hyungjun Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, 34141, Republic of Korea.
| |
Collapse
|
13
|
Ringe S, Hörmann NG, Oberhofer H, Reuter K. Implicit Solvation Methods for Catalysis at Electrified Interfaces. Chem Rev 2021; 122:10777-10820. [PMID: 34928131 PMCID: PMC9227731 DOI: 10.1021/acs.chemrev.1c00675] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Implicit solvation
is an effective, highly coarse-grained approach
in atomic-scale simulations to account for a surrounding liquid electrolyte
on the level of a continuous polarizable medium. Originating in molecular
chemistry with finite solutes, implicit solvation techniques are now
increasingly used in the context of first-principles modeling of electrochemistry
and electrocatalysis at extended (often metallic) electrodes. The
prevalent ansatz to model the latter electrodes and the reactive surface
chemistry at them through slabs in periodic boundary condition supercells
brings its specific challenges. Foremost this concerns the difficulty
of describing the entire double layer forming at the electrified solid–liquid
interface (SLI) within supercell sizes tractable by commonly employed
density functional theory (DFT). We review liquid solvation methodology
from this specific application angle, highlighting in particular its
use in the widespread ab initio thermodynamics approach
to surface catalysis. Notably, implicit solvation can be employed
to mimic a polarization of the electrode’s electronic density
under the applied potential and the concomitant capacitive charging
of the entire double layer beyond the limitations of the employed
DFT supercell. Most critical for continuing advances of this effective
methodology for the SLI context is the lack of pertinent (experimental
or high-level theoretical) reference data needed for parametrization.
Collapse
Affiliation(s)
- Stefan Ringe
- Department of Energy Science and Engineering, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea.,Energy Science & Engineering Research Center, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea
| | - Nicolas G Hörmann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany.,Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany.,Chair for Theoretical Physics VII and Bavarian Center for Battery Technology (BayBatt), University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
14
|
Kim M, Gould T, Izgorodina EI, Rocca D, Lebègue S. Establishing the accuracy of density functional approaches for the description of noncovalent interactions in ionic liquids. Phys Chem Chem Phys 2021; 23:25558-25564. [PMID: 34782901 DOI: 10.1039/d1cp03888e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We test a number of dispersion corrected versatile Generalized Gradient Approximation (GGA) and meta-GGA functionals for their ability to predict the interactions of ionic liquids, and show that most can achieve energies within 1 kcal mol-1 of benchmarks. This compares favorably with an accurate dispersion corrected hybrid, ωB97X-V. Our tests also reveal that PBE (Perdew-Burke-Ernzerhof GGA) calculations using the plane-wave projector augmented wave method and Gaussian Type Orbitals (GTOs) differ by less than 0.6 kJ mol-1 for ionic liquids, despite ions being difficult to evaluate in periodic cells - thus revealing that GTO benchmarks may be used also for plane-wave codes. Finally, the relatively high success of explicit van der Waals density functionals, compared to elemental and ionic dispersion models, suggests that improvements are required for low-cost dispersion correction models of ions.
Collapse
Affiliation(s)
- Minho Kim
- Université de Lorraine and CNRS, LPCT, UMR 7019, Vandoeuvre-lès-Nancy 54506, France.
| | - Tim Gould
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan, Queensland 4111, Australia
| | | | - Dario Rocca
- Université de Lorraine and CNRS, LPCT, UMR 7019, Vandoeuvre-lès-Nancy 54506, France.
| | - Sébastien Lebègue
- Université de Lorraine and CNRS, LPCT, UMR 7019, Vandoeuvre-lès-Nancy 54506, France.
| |
Collapse
|
15
|
Le JB, Yang XH, Zhuang YB, Jia M, Cheng J. Recent Progress toward Ab Initio Modeling of Electrocatalysis. J Phys Chem Lett 2021; 12:8924-8931. [PMID: 34499508 DOI: 10.1021/acs.jpclett.1c02086] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Electrode potential is the key factor for controlling electrocatalytic reactions at electrochemical interfaces, and moreover, it is also known that the pH and solutes (e.g., cations) of the solution have prominent effects on electrocatalysis. Understanding these effects requires microscopic information on the electrochemical interfaces, in which theoretical simulations can play an important role. This Perspective summarizes the recent progress in method development for modeling electrochemical interfaces, including different methods for describing the electrolytes at the interfaces and different schemes for charging up the electrode surfaces. In the final section, we provide an outlook for future development in modeling methods and their applications to electrocatalysis.
Collapse
Affiliation(s)
- Jia-Bo Le
- Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo 315201, China
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Xiao-Hui Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yong-Bin Zhuang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Mei Jia
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Jun Cheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
16
|
Yang X, Zhuang Y, Zhu J, Le J, Cheng J. Recent progress on multiscale modeling of electrochemistry. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1559] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Xiao‐Hui Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yong‐Bin Zhuang
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Jia‐Xin Zhu
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Jia‐Bo Le
- Ningbo Institute of Materials Technology and Engineering Chinese Academy of Sciences Ningbo China
| | - Jun Cheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM, College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| |
Collapse
|
17
|
Li S, Pokuri BSS, Ryno SM, Nkansah A, De'Vine C, Ganapathysubramanian B, Risko C. Determination of the Free Energies of Mixing of Organic Solutions through a Combined Molecular Dynamics and Bayesian Statistics Approach. J Chem Inf Model 2020; 60:1424-1431. [PMID: 31935097 DOI: 10.1021/acs.jcim.9b01113] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
As new generations of thin-film semiconductors are moving toward solution-based processing, the development of printing formulations will require information pertaining to the free energies of mixing of complex mixtures. From the standpoint of in silico material design, this move necessitates the development of methods that can accurately and quickly evaluate these formulations in order to maximize processing speed and reproducibility. Here, we make use of molecular dynamics (MD) simulations, in combination with the two-phase thermodynamic (2PT) model, to explore the free energy of mixing surfaces for a series of halogenated solvents and high-boiling point solvent additives used in the development of thin-film organic semiconductors. Although the combined methods generally show good agreement with available experimental data, the computational cost to traverse the free-energy landscape is considerable. Hence, we demonstrate how a Bayesian optimization scheme, coupled with the MD and 2PT approaches, can drastically reduce the number of simulations required, in turn shrinking both the computational cost and time.
Collapse
Affiliation(s)
- Shi Li
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | | | - Sean M Ryno
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | - Asare Nkansah
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | - Camron De'Vine
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | | | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
18
|
Kim J, Kim H, Lee WJ, Ruqia B, Baik H, Oh HS, Paek SM, Lim HK, Choi CH, Choi SI. Theoretical and Experimental Understanding of Hydrogen Evolution Reaction Kinetics in Alkaline Electrolytes with Pt-Based Core–Shell Nanocrystals. J Am Chem Soc 2019; 141:18256-18263. [DOI: 10.1021/jacs.9b09229] [Citation(s) in RCA: 71] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jeonghyeon Kim
- Department of Chemistry and Green-Nano Materials Research Center, Kyungpook National University, Daegu 41566, Korea
| | - Haesol Kim
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju 61005, Korea
| | - Won-Jae Lee
- Department of Chemistry and Green-Nano Materials Research Center, Kyungpook National University, Daegu 41566, Korea
| | - Bibi Ruqia
- Department of Chemistry and Green-Nano Materials Research Center, Kyungpook National University, Daegu 41566, Korea
| | - Hionsuck Baik
- Korea Basic Science Institute (KBSI), Seoul 02841, Korea
| | - Hyung-Suk Oh
- Clean Energy Research Center, Korea Institute of Science and Technology (KIST), Seoul 02792, Korea
| | - Seung-Min Paek
- Department of Chemistry and Green-Nano Materials Research Center, Kyungpook National University, Daegu 41566, Korea
| | - Hyung-Kyu Lim
- Division of Chemical Engineering and Bioengineering, Kangwon National University, Chuncheon, Gangwon-do 24341, Korea
| | - Chang Hyuck Choi
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology, Gwangju 61005, Korea
| | - Sang-Il Choi
- Department of Chemistry and Green-Nano Materials Research Center, Kyungpook National University, Daegu 41566, Korea
| |
Collapse
|
19
|
Structure, Dynamics, and Wettability of Water at Metal Interfaces. Sci Rep 2019; 9:14805. [PMID: 31616006 PMCID: PMC6794260 DOI: 10.1038/s41598-019-51323-5] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Accepted: 09/27/2019] [Indexed: 11/24/2022] Open
Abstract
The water/metal interface often governs important chemophysical processes in various technologies. Therefore, from scientific and engineering perspectives, the detailed molecular-level elucidation of the water/metal interface is of high priority, but the related research is limited. In experiments, the surface-science techniques, which can provide full structural details of the surface, are not easy to directly apply to the interfacial systems under ambient conditions, and the well-defined facets cannot be entirely free from contamination at the contact with water. To answer long-standing debates regarding the wettability, structure, and dynamics of water at metal interfaces, we here develop reliable first-principles-based multiscale simulations. Using the state-of-the-art simulations, we find that the clean metal surfaces are actually superhydrophilic and yield zero contact angles. Furthermore, we disclose an inadequacy of widespread ice-like bilayer model of the water adlayers on metal surfaces from both averaged structural and dynamic points of view. Our findings on the nature of water on metal surfaces provide new molecular level perspectives on the tuning and design of water/metal interfaces that are at the heart of many energy applications.
Collapse
|
20
|
Lee H, Dehez F, Chipot C, Lim HK, Kim H. Enthalpy-Entropy Interplay in π-Stacking Interaction of Benzene Dimer in Water. J Chem Theory Comput 2019; 15:1538-1545. [PMID: 30721623 DOI: 10.1021/acs.jctc.8b00880] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Aromatic groups can engage in an interesting class of noncovalent interactions termed π-π interactions, which play a pivotal role in stabilizing a variety of molecular architectures, including nucleic acids, proteins, and supramolecular assemblies. When the aromatic compounds interact with each other in an aqueous environment, their association is facilitated by the hydrophobic effect-the trend of nonpolar solutes to aggregate in a polar solution. To develop an in-depth understanding of hydrophobic association, we investigate in the present work π-π interactions in water, employing as a paradigm the benzene dimer. Using DFT-CES, a mean-field QM/MM method recently developed by our group, we describe the benzene solute at a quantum-mechanical level. Full consideration of detailed solute-electron density enables an optimal description of the solute-solvent interactions, leading to an accurate prediction of hydration free energies. In π-stacking of benzene, we find an entropic stabilization associated with the shrinkage of the solvent-excluded volume, which agrees with the theory of hydrophobic effect at subnanoscales. However, at the equilibrium binding distance of the benzene dimer, we find that the entropic stabilization nearly cancels out due to the enthalpic cost required for dewetting the internal space. Such an enthalpy-entropy compensation leads the association free energy to be predominantly dictated by the solute-solute interaction enthalpy. The present work offers new insight into the mechanistic role of water and the primary thermodynamic driving force of hydrophobic association.
Collapse
Affiliation(s)
- Hankyul Lee
- Department of Chemistry , Korea Advanced Institute of Science and Technology (KAIST) , Yuseong-gu, Daejeon 34141 , Korea
| | - François Dehez
- LPCT , UMR 7019 Université de Lorraine CNRS , Vandœuvre-lès-Nancy F-54500 , France.,Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign, Vandœuvre-lès-Nancy F-54506 , France
| | - Christophe Chipot
- LPCT , UMR 7019 Université de Lorraine CNRS , Vandœuvre-lès-Nancy F-54500 , France.,Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign, Vandœuvre-lès-Nancy F-54506 , France.,Department of Physics , University of Illinois at Urbana-Champaign , 1110 West Green Street , Urbana , Illinois 61801 , United States
| | - Hyung-Kyu Lim
- Division of Chemical Engineering and Bioengineering , Kangwon National University , Chuncheon , Gangwon-do 24341 , Korea
| | - Hyungjun Kim
- Department of Chemistry , Korea Advanced Institute of Science and Technology (KAIST) , Yuseong-gu, Daejeon 34141 , Korea
| |
Collapse
|
21
|
Zhang P, Shen L, Yang W. Solvation Free Energy Calculations with Quantum Mechanics/Molecular Mechanics and Machine Learning Models. J Phys Chem B 2019; 123:901-908. [PMID: 30557020 PMCID: PMC6448400 DOI: 10.1021/acs.jpcb.8b11905] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
For exploration of chemical and biological systems, the combined quantum mechanics and molecular mechanics (QM/MM) and machine learning (ML) models have been developed recently to achieve high accuracy and efficiency for molecular dynamics (MD) simulations. Despite its success on reaction free energy calculations, how to identify new configurations on insufficiently sampled regions during MD and how to update the current ML models with the growing database on the fly are both very important but still challenging. In this article, we apply the QM/MM ML method to solvation free energy calculations and address these two challenges. We employ three approaches to detect new data points and introduce the gradient boosting algorithm to reoptimize efficiently the ML model during ML-based MD sampling. The solvation free energy calculations on several typical organic molecules demonstrate that our developed method provides a systematic, robust, and efficient way to explore new chemistry using ML-based QM/MM MD simulations.
Collapse
Affiliation(s)
- Pan Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Lin Shen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry and Department of Physics, Duke University, Durham, NC 27708, United States
- Key laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment, South China Normal University, Guangzhou 510006, P.R.China
| |
Collapse
|
22
|
Lee H, Lim HK, Kim H. Hydration Thermodynamics of Non-Polar Aromatic Hydrocarbons: Comparison of Implicit and Explicit Solvation Models. Molecules 2018; 23:molecules23112927. [PMID: 30423973 PMCID: PMC6278349 DOI: 10.3390/molecules23112927] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Revised: 11/05/2018] [Accepted: 11/07/2018] [Indexed: 11/29/2022] Open
Abstract
The precise description of solute-water interactions is essential to understand the chemo-physical nature in hydration processes. Such a hydration thermodynamics for various solutes has been explored by means of explicit or implicit solvation methods. Using the Poisson-Boltzmann solvation model, the implicit models are well designed to reasonably predict the hydration free energies of polar solutes. The implicit model, however, is known to have shortcomings in estimating those for non-polar aromatic compounds. To investigate a cause of error, we employed a novel systematic framework of quantum-mechanical/molecular-mechanical (QM/MM) coupling protocol in explicit solvation manner, termed DFT-CES, based on the grid-based mean-field treatment. With the aid of DFT-CES, we delved into multiple energy parts, thereby comparing DFT-CES and PB models component-by-component. By applying the modified PB model to estimate the hydration free energies of non-polar solutes, we find a possibility to improve the predictability of PB models. We expect that this study could shed light on providing an accurate route to study the hydration thermodynamics for various solute compounds.
Collapse
Affiliation(s)
- Hankyul Lee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Yuseong-gu, Daejeon 34141, Korea.
| | - Hyung-Kyu Lim
- Division of Chemical Engineering and Bioengineering, Kangwon National University, Chuncheon, Gangwon-do 24341, Korea.
| | - Hyungjun Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Yuseong-gu, Daejeon 34141, Korea.
| |
Collapse
|
23
|
Steinmann SN, Ferreira De Morais R, Götz AW, Fleurat-Lessard P, Iannuzzi M, Sautet P, Michel C. Force Field for Water over Pt(111): Development, Assessment, and Comparison. J Chem Theory Comput 2018; 14:3238-3251. [PMID: 29660272 DOI: 10.1021/acs.jctc.7b01177] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Metal/water interfaces are key in many natural and industrial processes, such as corrosion, atmospheric, or environmental chemistry. Even today, the only practical approach to simulate large interfaces between a metal and water is to perform force-field simulations. In this work, we propose a novel force field, GAL17, to describe the interaction of water and a Pt(111) surface. GAL17 builds on three terms: (i) a standard Lennard-Jones potential for the bonding interaction between the surface and water, (ii) a Gaussian term to improve the surface corrugation, and (iii) two terms describing the angular dependence of the interaction energy. The 12 parameters of this force field are fitted against a set of 210 adsorption geometries of water on Pt(111). The performance of GAL17 is compared to several other approaches that have not been validated against extensive first-principles computations yet. Their respective accuracy is evaluated on an extended set of 802 adsorption geometries of H2O on Pt(111), 52 geometries derived from icelike layers, and an MD simulation of an interface between a c(4 × 6) Pt(111) surface and a water layer of 14 Å thickness. The newly developed GAL17 force field provides a significant improvement over previously existing force fields for Pt(111)/H2O interactions. Its well-balanced performance suggests that it is an ideal candidate to generate relevant geometries for the metal/water interface, paving the way to a representative sampling of the equilibrium distribution at the interface and to predict solvation free energies at the solid/liquid interface.
Collapse
Affiliation(s)
- Stephan N Steinmann
- Univ Lyon, Ecole Normale Supérieure de Lyon , CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182 , 46 allée d'Italie , F-69364 Lyon , France
| | - Rodrigo Ferreira De Morais
- Univ Lyon, Ecole Normale Supérieure de Lyon , CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182 , 46 allée d'Italie , F-69364 Lyon , France
| | - Andreas W Götz
- San Diego Supercomputer Center , University of California San Diego , La Jolla , California 92093 , United States
| | - Paul Fleurat-Lessard
- Institut de Chimie Moléculaire de l'Université de Bourgogne (ICMUB, UMR 6302, CNRS) , Université de Bourgogne Franche-Comté , 9 Avenue Alain Savary , 21078 Dijon , France
| | - Marcella Iannuzzi
- Institut für Chemie , University of Zurich , Winterthurerstrasse 190 , CH-8057 Zurich , Switzerland
| | | | - Carine Michel
- Univ Lyon, Ecole Normale Supérieure de Lyon , CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182 , 46 allée d'Italie , F-69364 Lyon , France
| |
Collapse
|
24
|
Gim S, Lim HK, Kim H. Multiscale Simulation Method for Quantitative Prediction of Surface Wettability at the Atomistic Level. J Phys Chem Lett 2018; 9:1750-1758. [PMID: 29558139 DOI: 10.1021/acs.jpclett.8b00466] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The solid-liquid interface is of great interest because of its highly heterogeneous character and its ubiquity in various applications. The most fundamental physical variable determining the strength of the solid-liquid interface is the solid-liquid interfacial tension, which is usually measured according to the contact angle. However, an accurate experimental measurement and a reliable theoretical prediction of the contact angle remain lacking because of many practical issues. Here, we propose a first-principles-based simulation approach to quantitatively predict the contact angle of an ideally clean surface using our recently developed multiscale simulation method of density functional theory in classical explicit solvents (DFT-CES). Using this approach, we simulate the surface wettability of a graphene and graphite surface, resulting in a reliable contact angle value that is comparable to the experimental data. From our simulation results, we find that the surface wettability is dominantly affected by the strength of the solid-liquid van der Waal's interaction. However, we further elucidate that there exists a secondary contribution from the change of water-water interaction, which is manifested by the change of liquid structure and dynamics of interfacial water layer. We expect that our proposed method can be used to quantitatively predict and understand the intriguing wetting phenomena at an atomistic level and can eventually be utilized to design a surface with a controlled hydrophobic(philic)ity.
Collapse
Affiliation(s)
- Suji Gim
- Department of Chemistry and Graduate School of EEWS , Korea Advanced Institute of Science and Technology (KAIST) , Yuseong-gu, Daejeon 34141 , Korea
| | - Hyung-Kyu Lim
- Department of Chemical Engineering , Kangwon National University , Chuncheon , Gangwon-do 24341 , Korea
| | - Hyungjun Kim
- Department of Chemistry and Graduate School of EEWS , Korea Advanced Institute of Science and Technology (KAIST) , Yuseong-gu, Daejeon 34141 , Korea
| |
Collapse
|
25
|
Lim HK, Kwon Y, Kim HS, Jeon J, Kim YH, Lim JA, Kim BS, Choi J, Kim H. Insight into the Microenvironments of the Metal–Ionic Liquid Interface during Electrochemical CO2 Reduction. ACS Catal 2018. [DOI: 10.1021/acscatal.7b03777] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hyung-Kyu Lim
- Department
of Chemical Engineering, Kangwon National University, Chuncheon, Gangwon-do 24341, Republic of Korea
| | - Youngkook Kwon
- Carbon
Resource Institute, Korea Research Institute of Chemical Technology, 141 Gajeong-ro, Yuseong-gu, Daejeon 34114, Republic of Korea
- Advanced Materials and Chemical Engineering, University of Science & Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Republic of Korea
| | - Han Seul Kim
- Graduate
School of EEWS, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Jiwon Jeon
- Graduate
School of EEWS, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Yong-Hoon Kim
- Graduate
School of EEWS, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Jung-Ae Lim
- Carbon
Resource Institute, Korea Research Institute of Chemical Technology, 141 Gajeong-ro, Yuseong-gu, Daejeon 34114, Republic of Korea
| | - Beom-Sik Kim
- Carbon
Resource Institute, Korea Research Institute of Chemical Technology, 141 Gajeong-ro, Yuseong-gu, Daejeon 34114, Republic of Korea
- Advanced Materials and Chemical Engineering, University of Science & Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Republic of Korea
| | - Jina Choi
- Carbon
Resource Institute, Korea Research Institute of Chemical Technology, 141 Gajeong-ro, Yuseong-gu, Daejeon 34114, Republic of Korea
| | - Hyungjun Kim
- Graduate
School of EEWS, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
- Department
of Chemistry, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| |
Collapse
|
26
|
Dohn AO, Jónsson EÖ, Levi G, Mortensen JJ, Lopez-Acevedo O, Thygesen KS, Jacobsen KW, Ulstrup J, Henriksen NE, Møller KB, Jónsson H. Grid-Based Projector Augmented Wave (GPAW) Implementation of Quantum Mechanics/Molecular Mechanics (QM/MM) Electrostatic Embedding and Application to a Solvated Diplatinum Complex. J Chem Theory Comput 2017; 13:6010-6022. [PMID: 29083921 DOI: 10.1021/acs.jctc.7b00621] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
A multiscale density functional theory-quantum mechanics/molecular mechanics (DFT-QM/MM) scheme is presented, based on an efficient electrostatic coupling between the electronic density obtained from a grid-based projector augmented wave (GPAW) implementation of density functional theory and a classical potential energy function. The scheme is implemented in a general fashion and can be used with various choices for the descriptions of the QM or MM regions. Tests on H2O clusters, ranging from dimer to decamer show that no systematic energy errors are introduced by the coupling that exceeds the differences in the QM and MM descriptions. Over 1 ns of liquid water, Born-Oppenheimer QM/MM molecular dynamics (MD) are sampled combining 10 parallel simulations, showing consistent liquid water structure over the QM/MM border. The method is applied in extensive parallel MD simulations of an aqueous solution of the diplatinum [Pt2(P2O5H2)4]4- complex (PtPOP), spanning a total time period of roughly half a nanosecond. An average Pt-Pt distance deviating only 0.01 Å from experimental results, and a ground-state Pt-Pt oscillation frequency deviating by <2% from experimental results were obtained. The simulations highlight a remarkable harmonicity of the Pt-Pt oscillation, while also showing clear signs of Pt-H hydrogen bonding and directional coordination of water molecules along the Pt-Pt axis of the complex.
Collapse
Affiliation(s)
- A O Dohn
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - E Ö Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - G Levi
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J J Mortensen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - O Lopez-Acevedo
- Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| | - K S Thygesen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K W Jacobsen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J Ulstrup
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - N E Henriksen
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K B Møller
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - H Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland.,Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| |
Collapse
|
27
|
Steinmann SN, Sautet P, Michel C. Solvation free energies for periodic surfaces: comparison of implicit and explicit solvation models. Phys Chem Chem Phys 2016; 18:31850-31861. [DOI: 10.1039/c6cp04094b] [Citation(s) in RCA: 61] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A strategy based on molecular mechanics free energy of perturbation, seeded by quantum mechanics, is presented to take solvation energies into account in the context of periodic, solid–liquid interfaces.
Collapse
Affiliation(s)
- Stephan N. Steinmann
- Univ Lyon
- Ecole Normale Supérieure de Lyon
- CNRS Université Lyon 1
- Laboratoire de Chimie UMR 5182
- Lyon
| | - Philippe Sautet
- Univ Lyon
- Ecole Normale Supérieure de Lyon
- CNRS Université Lyon 1
- Laboratoire de Chimie UMR 5182
- Lyon
| | - Carine Michel
- Univ Lyon
- Ecole Normale Supérieure de Lyon
- CNRS Université Lyon 1
- Laboratoire de Chimie UMR 5182
- Lyon
| |
Collapse
|