1
|
Chai Z, Luber S. Grand Canonical Ensemble Approaches in CP2K for Modeling Electrochemistry at Constant Electrode Potentials. J Chem Theory Comput 2024. [PMID: 39240723 DOI: 10.1021/acs.jctc.4c00671] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2024]
Abstract
In electrochemical experiments, the number of electrons of the electrode immersed in the electrolyte is usually variable. Additionally, the numbers of adsorbed substances on the surface of the electrode, the solvent molecules, and counter charge ions in the near-surface region can also vary. Treating electrochemical solid-liquid interfaces with the typical fixed electron number density functional theory (DFT) approach tends to be a challenge. This can be addressed by using grand canonical ensemble approaches. We present the implementation of two grand canonical ensemble approaches in the open-source computational chemistry software CP2K that go beyond the existing canonical ensemble paradigm. The first approach is based on implicit solvent models and explicit atomistic solute (electrode with/without adsorbed species) models, and includes two recent developments: (a) grand canonical self-consistent field (GC-SCF) method (J. Chem. Phys. 2017, 146, 114104) allowing the electron number of the system to fluctuate naturally and accordingly with the experimental electrode potential, (b) planar counter charge (J. Chem. Phys. 2019, 150, 041722, Phys. Rev. B 2003, 68, 245416) salt model completely screening the net charge of the electrode model. In contrast with previous studies, in our implementation, the work function (WF) (absolute electrode potential if the potential drop at the electrolyte-vacuum interface is omitted) is the constrained quantity during an SCF optimization instead of the Fermi energy. The chemical potential of electrons (negative WF) is a natural variable of the grand potential in the GC ensemble of electronic states, and this method can easily achieve stable SCF convergence and obtain an electronic structure that precisely corresponds to a user-specified WF. The second approach referred to as the GC DFT molecular dynamics (DFT-MD) simulation scheme (Phys. Rev. Lett. 2002, 88, 213002, J. Chem. Phys. 2005, 122, 234505, J. Am. Chem. Soc. 2004, 126 (12), 3928-3938) is based on fully explicit modeling the solvent molecules and the ions and is used to calculate the electron chemical potential corresponding to an equilibrium electrochemical half-reaction (M(n+m)+ + ne- ⇌ Mm+) which involves DFT-MD, by allowing the number of electrons to vary during the DFT-MD simulation process. This opens the way for forefront electrochemical calculations in CP2K for a broad range of systems.
Collapse
Affiliation(s)
- Ziwei Chai
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| |
Collapse
|
2
|
Wang Y, Teng C, Begin E, Bussiere M, Bao JL. PW-SMD: A Plane-Wave Implicit Solvation Model Based on Electron Density for Surface Chemistry and Crystalline Systems in Aqueous Solution. J Chem Theory Comput 2024. [PMID: 39024317 DOI: 10.1021/acs.jctc.4c00594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Electron density-based implicit solvation models are a class of techniques for quantifying solvation effects and calculating free energies of solvation without an explicit representation of solvent molecules. Integral to the accuracy of solvation modeling is the proper definition of the solvation shell separating the solute molecule from the solvent environment, allowing for a physical partitioning of the free energies of solvation. Unlike state-of-the-art implicit solvation models for molecular quantum chemistry calculations, e.g., the solvation model based on solute electron density (SMD), solvation models for systems under periodic boundary conditions with plane-wave (PW) basis sets have been limited in their accuracy. Furthermore, a unified implicit solvation model with both homogeneous solution-phase and heterogeneous interfacial structures treated on equal footing is needed. In order to address this challenge, we developed a high-accuracy solvation model for periodic PW calculations that is applicable to molecular, ionic, interfacial, and bulk-phase chemistry. Our model, PW-SMD, is an extension of the SMD molecular solvation model to periodic systems in water. The free energy of solvation is partitioned into the electrostatic and cavity-dispersion-solvent structure (CDS) contributions. The electrostatic contributions of the solvation shell surrounding solute structures are parametrized based on their geometric and physical properties. In addition, the nonelectrostatic contribution to the solvation energy is accounted for by extending the CDS formalism of SMD to incorporate periodic boundary conditions. We validate the accuracy and robustness of our solvation model by comparing predicted solvation free energies against experimental data for molecular and ionic systems, carved-cluster composite energetic models of solvated reaction energies and barriers on surface systems, and deep-learning-accelerated ab initio molecular dynamics (AIMD). Our developed periodic implicit solvation model shows significantly improved accuracy compared to previous work (namely, solvation models in aqueous solution) and can be applied to simulate solvent effects in a wide range of surface and crystalline materials.
Collapse
Affiliation(s)
- Yang Wang
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Chong Teng
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Elijah Begin
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Mason Bussiere
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| |
Collapse
|
3
|
Labat M, Giner E, Jeanmairet G. Coupling molecular density functional theory with converged selected configuration interaction methods to study excited states in aqueous solution. J Chem Phys 2024; 161:014113. [PMID: 38958166 DOI: 10.1063/5.0213426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 06/12/2024] [Indexed: 07/04/2024] Open
Abstract
This paper presents the first implementation of a coupling between advanced wavefunction theories and molecular density functional theory (MDFT). This method enables the modeling of solvent effect into quantum mechanical (QM) calculations by incorporating an electrostatic potential generated by solvent charges into the electronic Hamiltonian. Solvent charges are deduced from the spatially and angularly dependent solvent particle density. Such a density is obtained through the minimization of the functional associated with the molecular mechanics (MM) Hamiltonian describing the interaction between the fluid particles. The introduced QM/MDFT framework belongs to QM/MM family of methods, but its originality lies in the use of MDFT as the MM solver, offering two main advantages. First, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Second, it preserves a molecular-level description lost in macroscopic continuum approaches. The excited state properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. The excitation energies and dipole moments have been compared with experimental data and previous theoretical work. A key finding is that using the Hartree-Fock method to describe the solute allows for predicting the solvent charge around the ground state with sufficient precision for the subsequent SCI calculations of excited states. This significantly reduces the computational cost of the described procedure, paving the way for the study of more complex molecules.
Collapse
Affiliation(s)
- Maxime Labat
- Sorbonne Université, CNRS, Physico-Chimie des électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| | - Emmanuel Giner
- Sorbonne Université, CNRS, Laboratoire de Chimie Théorique, Sorbonne Université, F-75005 Paris, France
| | - Guillaume Jeanmairet
- Sorbonne Université, CNRS, Physico-Chimie des électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
- Réseau sur le Stockage électrochimique de l'énergie (RS2E), FR CNRS 3459, 80039 Amiens Cedex, France
| |
Collapse
|
4
|
Si P, Jayanth A, Andreussi O. Soft-sphere continuum solvation models for nonaqueous solvents. J Comput Chem 2024; 45:719-737. [PMID: 38112395 DOI: 10.1002/jcc.27254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Revised: 10/24/2023] [Accepted: 10/30/2023] [Indexed: 12/21/2023]
Abstract
Solvation effects profoundly influence the characteristics and behavior of chemical systems in liquid solutions. The interaction between solute and solvent molecules intricately impacts solubility, reactivity, stability, and various chemical processes. Continuum solvation models gained prominence in quantum chemistry by implicitly capturing these interactions and enabling efficient investigations of diverse chemical systems in solution. In comparison, continuum solvation models in condensed matter simulation are very recent. Among these, the self-consistent continuum solvation (SCCS) and the soft-sphere continuum solvation models (SSCS) have been among the first to be successfully parameterized and extended to model periodic systems in aqueous solutions and electrolytes. As most continuum approaches, these models depend on a number of parameters that are linked to experimental or theoretical properties of the solvent, or that can be tuned based on reference data. Here, we present a systematic parameterization of the SSCS model for over 100 nonaqueous solvents. We validate the model's efficacy across diverse solvent environments by leveraging experimental solvation-free energies and partition coefficients from comprehensive databases. The average root means square error over all the solvents was calculated as 0.85 kcal/mol which is below the chemical accuracy (1 kcal/mol). Similarly to what has been reported by Hille et al. (J. Chem. Phys. 2019, 150, 041710.) for the SCCS model, a single-parameter model accurately reproduces experimental solvation energies, showcasing the transferability and predictive power of these continuum approaches. Our findings underscore the potential for a unified approach to predict solvation properties, paving the way for enhanced computational studies across various chemical environments.
Collapse
Affiliation(s)
- Pradip Si
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - Ajay Jayanth
- Texas Academy of Math and Science, University of North Texas, Denton, Texas, USA
| | - Oliviero Andreussi
- Department of Chemistry and Biochemistry, Boise State University, Boise, Idaho, USA
| |
Collapse
|
5
|
Putra MH, Bagemihl B, Rau S, Groß A. Prediction of Strong Solvatochromism in a Molecular Photocatalyst. Chemistry 2024; 30:e202302643. [PMID: 37754665 DOI: 10.1002/chem.202302643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2023] [Revised: 09/25/2023] [Accepted: 09/25/2023] [Indexed: 09/28/2023]
Abstract
Based on quantum chemical calculations, we predict strong solvatochromism in a light-driven molecular photocatalyst for hydrogen generation, that is we show that the electronic and optical properties of the photocatalyst strongly depend on the solvent it is dissolved in. Our calculations in particular indicate a solvent-dependent relocation of the highest occupied molecular orbital (HOMO). Ground-state density functional theory and linear response time-dependent density functional theory calculations were applied in order to investigate the influence of implicit solvents on the structural, electronic and optical properties of a molecular photocatalyst. Only at high dielectric constants of the solvent, is the HOMO located at the metal center of the photosensitizer, whereas at low dielectric constants the HOMO is centered at the metal atom of the catalytically active complex. We elucidate the electronic origins of this strong solvatochromic effect and sketch the consequences of these insights for the use of photocatalysts in different environments.
Collapse
Affiliation(s)
| | - Benedikt Bagemihl
- Institute of Inorganic Chemistry I, Materials and Catalysis, Ulm University, 89069, Ulm, Germany
| | - Sven Rau
- Institute of Inorganic Chemistry I, Materials and Catalysis, Ulm University, 89069, Ulm, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069, Ulm, Germany
- Helmholtz Institute Ulm (HIU), Electrochemical Energy Storage, 89069, Ulm, Germany
| |
Collapse
|
6
|
M V, Singh S, Bononi F, Andreussi O, Karmodak N. Thermodynamic and kinetic modeling of electrocatalytic reactions using a first-principles approach. J Chem Phys 2023; 159:111001. [PMID: 37728202 DOI: 10.1063/5.0165835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 08/28/2023] [Indexed: 09/21/2023] Open
Abstract
The computational modeling of electrochemical interfaces and their applications in electrocatalysis has attracted great attention in recent years. While tremendous progress has been made in this area, however, the accurate atomistic descriptions at the electrode/electrolyte interfaces remain a great challenge. The Computational Hydrogen Electrode (CHE) method and continuum modeling of the solvent and electrolyte interactions form the basis for most of these methodological developments. Several posterior corrections have been added to the CHE method to improve its accuracy and widen its applications. The most recently developed grand canonical potential approaches with the embedded diffuse layer models have shown considerable improvement in defining interfacial interactions at electrode/electrolyte interfaces over the state-of-the-art computational models for electrocatalysis. In this Review, we present an overview of these different computational models developed over the years to quantitatively probe the thermodynamics and kinetics of electrochemical reactions in the presence of an electrified catalyst surface under various electrochemical environments. We begin our discussion by giving a brief picture of the different continuum solvation approaches, implemented within the ab initio method to effectively model the solvent and electrolyte interactions. Next, we present the thermodynamic and kinetic modeling approaches to determine the activity and stability of the electrocatalysts. A few applications to these approaches are also discussed. We conclude by giving an outlook on the different machine learning models that have been integrated with the thermodynamic approaches to improve their efficiency and widen their applicability.
Collapse
Affiliation(s)
- Vasanthapandiyan M
- Department of Chemistry, Shiv Nadar Institution of Eminence, Dadri, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| | - Shagun Singh
- Department of Chemistry, Shiv Nadar Institution of Eminence, Dadri, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| | - Fernanda Bononi
- Department of Physics, University of North Texas, Denton, Texas 76203, USA
| | - Oliviero Andreussi
- Department of Chemistry and Biochemistry, Boise State University, Boise, Idaho 83725, USA
| | - Naiwrit Karmodak
- Department of Chemistry, Shiv Nadar Institution of Eminence, Dadri, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| |
Collapse
|
7
|
German E, Gebauer R. The Oxygen Evolution Reaction at MoS 2 Edge Sites: The Role of a Solvent Environment in DFT-Based Molecular Simulations. Molecules 2023; 28:5182. [PMID: 37446844 DOI: 10.3390/molecules28135182] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Revised: 06/23/2023] [Accepted: 06/28/2023] [Indexed: 07/15/2023] Open
Abstract
Density functional theory (DFT) calculations are employed to study the oxygen evolution reaction (OER) on the edges of stripes of monolayer molybdenum disulfide. Experimentally, this material has been shown to evolve oxygen, albeit with low efficiency. Previous DFT studies have traced this low catalytic performance to the unfavourable adsorption energies of some reaction intermediates on the MoS2 edge sites. In this work, we study the effects of the aqueous liquid surrounding the active sites. A computational approach is used, where the solvent is modeled as a continuous medium providing a dielectric embedding of the catalyst and the reaction intermediates. A description at this level of theory can have a profound impact on the studied reactions: the calculated overpotential for the OER is lowered from 1.15 eV to 0.77 eV. It is shown that such variations in the reaction energetics are linked to the polar nature of the adsorbed intermediates, which leads to changes in the calculated electronic charge density when surrounded by water. These results underline the necessity to computationally account for solvation effects, especially in aqueous environments and when highly polar intermediates are present.
Collapse
Affiliation(s)
- Estefania German
- Department of Theoretical, Atomic and Optical Physics, University of Valladolid, 47011 Valladolid, Spain
| | - Ralph Gebauer
- The Abdus Salam International Centre for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
8
|
Zhang KX, Liu ZP. Electrochemical hydrogen evolution on Pt-based catalysts from a theoretical perspective. J Chem Phys 2023; 158:141002. [PMID: 37061480 DOI: 10.1063/5.0142540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
Hydrogen evolution reaction (HER) by splitting water is a key technology toward a clean energy society, where Pt-based catalysts were long known to have the highest activity under acidic electrochemical conditions but suffer from high cost and poor stability. Here, we overview the current status of Pt-catalyzed HER from a theoretical perspective, focusing on the methodology development of electrochemistry simulation, catalytic mechanism, and catalyst stability. Recent developments in theoretical methods for studying electrochemistry are introduced, elaborating on how they describe solid-liquid interface reactions under electrochemical potentials. The HER mechanism, the reaction kinetics, and the reaction sites on Pt are then summarized, which provides an atomic-level picture of Pt catalyst surface dynamics under reaction conditions. Finally, state-of-the-art experimental solutions to improve catalyst stability are also introduced, which illustrates the significance of fundamental understandings in the new catalyst design.
Collapse
Affiliation(s)
- Ke-Xiang Zhang
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Key Laboratory of Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Zhi-Pan Liu
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Key Laboratory of Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
9
|
Hashemi A, Peljo P, Laasonen K. Understanding Electron Transfer Reactions Using Constrained Density Functional Theory: Complications Due to Surface Interactions. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:3398-3407. [PMID: 36865990 PMCID: PMC9969872 DOI: 10.1021/acs.jpcc.2c06537] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 01/25/2023] [Indexed: 06/18/2023]
Abstract
The kinetic rates of electrochemical reactions depend on electrodes and molecules in question. In a flow battery, where the electrolyte molecules are charged and discharged on the electrodes, the efficiency of the electron transfer is of crucial importance for the performance of the device. The purpose of this work is to present a systematic atomic-level computational protocol for studying electron transfer between electrolyte and electrode. The computations are done by using constrained density functional theory (CDFT) to ensure that the electron is either on the electrode or in the electrolyte. The ab initio molecular dynamics (AIMD) is used to simulate the movement of the atoms. We use the Marcus theory to predict electron transfer rates and the combined CDFT-AIMD approach to compute the parameters for the Marcus theory where it is needed. We model the electrode with a single layer of graphene and methylviologen, 4,4'-dimethyldiquat, desalted basic red 5, 2-hydroxy-1,4-naphthaquinone, and 1,1-di(2-ethanol)-4,4-bipyridinium were selected for the electrolyte molecules. All of these molecules undergo consecutive electrochemical reactions with one electron being transferred at each stage. Because of significant electrode-molecule interactions, it is not possible to evaluate outer-sphere ET. This theoretical study contributes toward the development of a realistic-level prediction of electron transfer kinetics suitable for energy storage applications.
Collapse
Affiliation(s)
- Arsalan Hashemi
- Research
Group of Computational Chemistry, Department of Chemistry and Materials
Science, Aalto University, FI-00076 Aalto, Finland
| | - Pekka Peljo
- Research
Group of Battery Materials and Technologies, Department of Mechanical
and Materials Engineering, Faculty of Technology, University of Turku, 20014 Turun Yliopisto, Finland
| | - Kari Laasonen
- Research
Group of Computational Chemistry, Department of Chemistry and Materials
Science, Aalto University, FI-00076 Aalto, Finland
| |
Collapse
|
10
|
Xue W, Li J, Huang H, Zhang W, Mei D. Theoretical Screening of CO 2 Electroreduction over MOF-808-Supported Self-Adaptive Dual-Metal-Site Pairs. Inorg Chem 2023; 62:930-941. [PMID: 36607142 DOI: 10.1021/acs.inorgchem.2c03734] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Electrochemical CO2 reduction to transportation fuels and valuable platform chemicals provides a sustainable avenue for renewable energy storage and realizes an artificially closed carbon loop. However, the rational design of highly active and selective CO2 reduction electrocatalysts remains a challenging task. Herein, a series of metal-organic framework (MOF)-supported flexible, self-adaptive dual-metal-site pairs (DMSPs) including 21 pairwise combinations of six transition metal single sites (MOF-808-EDTA-M1M2, M1/M2 = Fe, Cu, Ni, Pd, Pt, Au) for the CO2 reduction reaction (CO2RR) were theoretically screened using density functional theory calculations. Against the competitive hydrogen evolution reaction, MOF-808-EDTA-FeFe and MOF-808-EDTA-FePt were identified as the promising CO2RR electrocatalysts toward C1 and C2 products. The calculated limiting potential for CO2 electroreduction to C2H6 and C2H5OH over MOF-808-EDTA-FeFe is -0.87 V. Compared with an applied potential of -0.56 eV toward CH4 production over MOF-808-EDTA-FeFe, MOF-808-EDTA-FePt exhibits an even better activity for CO2 reduction to C1 products at a limiting potential of -0.35 V. The present work not only identifies promising candidates for highly selective CO2RR electrocatalysts leading to C1 and C2 products but also provides mechanistic insights into the dynamic nature of DMSPs for stabilizing various reaction intermediates in the CO2RR process.
Collapse
Affiliation(s)
- Wenjuan Xue
- School of Chemical Engineering and Technology, Tiangong University, Tianjin300387, China.,State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin300387, China
| | - Jian Li
- School of Chemical Engineering and Technology, Tiangong University, Tianjin300387, China.,State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin300387, China
| | - Hongliang Huang
- School of Chemical Engineering and Technology, Tiangong University, Tianjin300387, China.,State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin300387, China
| | - Weiwei Zhang
- School of Chemical Engineering and Technology, Tiangong University, Tianjin300387, China.,State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin300387, China
| | - Donghai Mei
- School of Chemical Engineering and Technology, Tiangong University, Tianjin300387, China.,School of Environmental Science and Engineering, Tiangong University, Tianjin300387, China.,State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin300387, China
| |
Collapse
|
11
|
Karibayev M, Kalybekkyzy S, Wang Y, Mentbayeva A. Molecular Modeling in Anion Exchange Membrane Research: A Brief Review of Recent Applications. Molecules 2022; 27:3574. [PMID: 35684512 PMCID: PMC9182285 DOI: 10.3390/molecules27113574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 05/24/2022] [Accepted: 05/30/2022] [Indexed: 12/04/2022] Open
Abstract
Anion Exchange Membrane (AEM) fuel cells have attracted growing interest, due to their encouraging advantages, including high power density and relatively low cost. AEM is a polymer matrix, which conducts hydroxide (OH-) ions, prevents physical contact of electrodes, and has positively charged head groups (mainly quaternary ammonium (QA) groups), covalently bound to the polymer backbone. The chemical instability of the quaternary ammonium (QA)-based head groups, at alkaline pH and elevated temperature, is a significant threshold in AEMFC technology. This review work aims to introduce recent studies on the chemical stability of various QA-based head groups and transportation of OH- ions in AEMFC, via modeling and simulation techniques, at different scales. It starts by introducing the fundamental theories behind AEM-based fuel-cell technology. In the main body of this review, we present selected computational studies that deal with the effects of various parameters on AEMs, via a variety of multi-length and multi-time-scale modeling and simulation methods. Such methods include electronic structure calculations via the quantum Density Functional Theory (DFT), ab initio, classical all-atom Molecular Dynamics (MD) simulations, and coarse-grained MD simulations. The explored processing and structural parameters include temperature, hydration levels, several QA-based head groups, various types of QA-based head groups and backbones, etc. Nowadays, many methods and software packages for molecular and materials modeling are available. Applications of such methods may help to understand the transportation mechanisms of OH- ions, the chemical stability of functional head groups, and many other relevant properties, leading to a performance-based molecular and structure design as well as, ultimately, improved AEM-based fuel cell performances. This contribution aims to introduce those molecular modeling methods and their recent applications to the AEM-based fuel cells research community.
Collapse
Affiliation(s)
- Mirat Karibayev
- Department of Chemical & Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
| | - Sandugash Kalybekkyzy
- Laboratory of Advanced Materials and Systems for Energy Storage, Center for Energy and Advanced Materials Science, National Laboratory Astana, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
| | - Yanwei Wang
- Department of Chemical & Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
- Laboratory of Computational Materials Science for Energy Applications, Center for Energy and Advanced Materials Science, National Laboratory Astana, Nur-Sultan 010000, Kazakhstan
| | - Almagul Mentbayeva
- Department of Chemical & Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
- Laboratory of Advanced Materials and Systems for Energy Storage, Center for Energy and Advanced Materials Science, National Laboratory Astana, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
| |
Collapse
|
12
|
Warburton RE, Soudackov AV, Hammes-Schiffer S. Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer. Chem Rev 2022; 122:10599-10650. [PMID: 35230812 DOI: 10.1021/acs.chemrev.1c00929] [Citation(s) in RCA: 42] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.
Collapse
Affiliation(s)
- Robert E Warburton
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
13
|
Prentice JCA. Efficiently Computing Excitations of Complex Systems: Linear-Scaling Time-Dependent Embedded Mean-Field Theory in Implicit Solvent. J Chem Theory Comput 2022; 18:1542-1554. [PMID: 35133827 PMCID: PMC9082505 DOI: 10.1021/acs.jctc.1c01133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Quantum embedding schemes have the
potential to significantly reduce
the computational cost of first-principles calculations while maintaining
accuracy, particularly for calculations of electronic excitations
in complex systems. In this work, I combine time-dependent embedded
mean field theory (TD-EMFT) with linear-scaling density functional
theory and implicit solvation models, extending previous work within
the ONETEP code. This provides a way to perform multilevel calculations
of electronic excitations on very large systems, where long-range
environmental effects, both quantum and classical in nature, are important.
I demonstrate the power of this method by performing simulations on
a variety of systems, including a molecular dimer, a chromophore in
solution, and a doped molecular crystal. This work paves the way for
high accuracy calculations to be performed on large-scale systems
that were previously beyond the reach of quantum embedding schemes.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, University of Oxford, Parks Road, Oxford OX1 3PH, United Kingdom
| |
Collapse
|
14
|
Yang Y, Peltier CR, Zeng R, Schimmenti R, Li Q, Huang X, Yan Z, Potsi G, Selhorst R, Lu X, Xu W, Tader M, Soudackov AV, Zhang H, Krumov M, Murray E, Xu P, Hitt J, Xu L, Ko HY, Ernst BG, Bundschu C, Luo A, Markovich D, Hu M, He C, Wang H, Fang J, DiStasio RA, Kourkoutis LF, Singer A, Noonan KJT, Xiao L, Zhuang L, Pivovar BS, Zelenay P, Herrero E, Feliu JM, Suntivich J, Giannelis EP, Hammes-Schiffer S, Arias T, Mavrikakis M, Mallouk TE, Brock JD, Muller DA, DiSalvo FJ, Coates GW, Abruña HD. Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies. Chem Rev 2022; 122:6117-6321. [PMID: 35133808 DOI: 10.1021/acs.chemrev.1c00331] [Citation(s) in RCA: 89] [Impact Index Per Article: 44.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.
Collapse
Affiliation(s)
- Yao Yang
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Cheyenne R Peltier
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Rui Zeng
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Roberto Schimmenti
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Qihao Li
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Xin Huang
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | - Zhifei Yan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Georgia Potsi
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Ryan Selhorst
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Xinyao Lu
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Weixuan Xu
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Mariel Tader
- Department of Physics, Cornell University, Ithaca, New York 14853, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Hanguang Zhang
- Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mihail Krumov
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Ellen Murray
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Pengtao Xu
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Jeremy Hitt
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Linxi Xu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Brian G Ernst
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Colin Bundschu
- Department of Physics, Cornell University, Ithaca, New York 14853, United States
| | - Aileen Luo
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Danielle Markovich
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | - Meixue Hu
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Cheng He
- Chemical and Materials Science Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Hongsen Wang
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Jiye Fang
- Department of Chemistry, State University of New York at Binghamton, Binghamton, New York 13902, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Lena F Kourkoutis
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States.,Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, New York 14853, United States
| | - Andrej Singer
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Kevin J T Noonan
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Li Xiao
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Lin Zhuang
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Bryan S Pivovar
- Chemical and Materials Science Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Piotr Zelenay
- Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Enrique Herrero
- Instituto de Electroquímica, Universidad de Alicante, Alicante E-03080, Spain
| | - Juan M Feliu
- Instituto de Electroquímica, Universidad de Alicante, Alicante E-03080, Spain
| | - Jin Suntivich
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States.,Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, New York 14853, United States
| | - Emmanuel P Giannelis
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | | | - Tomás Arias
- Department of Physics, Cornell University, Ithaca, New York 14853, United States
| | - Manos Mavrikakis
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Thomas E Mallouk
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joel D Brock
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | - David A Muller
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States.,Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, New York 14853, United States
| | - Francis J DiSalvo
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Geoffrey W Coates
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Héctor D Abruña
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Center for Alkaline Based Energy Solutions (CABES), Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
15
|
Abstract
Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.
Collapse
Affiliation(s)
- Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany.,Electrochemical Energy Storage, Helmholtz Institute Ulm (HIU), 89069 Ulm, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| |
Collapse
|
16
|
Shao F, Zheng L, Lan J, Zenobi R. Nanoscale Chemical Imaging of Coadsorbed Thiolate Self-Assembled Monolayers on Au(111) by Tip-Enhanced Raman Spectroscopy. Anal Chem 2022; 94:1645-1653. [DOI: 10.1021/acs.analchem.1c03968] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Feng Shao
- Department of Physics and Astronomy, National Graphene Institute, University of Manchester, Manchester M13 9PL, U.K
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
| | - Liqing Zheng
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
| | - Jinggang Lan
- Department of Chemistry, University of Zurich, 8057 Zurich, Switzerland
| | - Renato Zenobi
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
| |
Collapse
|
17
|
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
|
18
|
|
19
|
Nakov S, Sobakinskaya E, Renger T, Kraus J. ARGOS: An adaptive refinement goal-oriented solver for the linearized Poisson-Boltzmann equation. J Comput Chem 2021; 42:1832-1860. [PMID: 34302374 DOI: 10.1002/jcc.26716] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 06/10/2021] [Accepted: 06/15/2021] [Indexed: 11/09/2022]
Abstract
An adaptive finite element solver for the numerical calculation of the electrostatic coupling between molecules in a solvent environment is developed and tested. At the heart of the solver is a goal-oriented a posteriori error estimate for the electrostatic coupling, derived and implemented in the present work, that gives rise to an orders of magnitude improved precision and a shorter computational time as compared to standard finite difference solvers. The accuracy of the new solver ARGOS is evaluated by numerical experiments on a series of problems with analytically known solutions. In addition, the solver is used to calculate electrostatic couplings between two chromophores, linked to polyproline helices of different lengths and between the spike protein of SARS-CoV-2 and the ACE2 receptor. All the calculations are repeated by using the well-known finite difference solvers MEAD and APBS, revealing the advantages of the present finite element solver.
Collapse
Affiliation(s)
- Svetoslav Nakov
- Institute for Theoretical Physics, Johannes Kepler University, Linz, Austria
| | | | - Thomas Renger
- Institute for Theoretical Physics, Johannes Kepler University, Linz, Austria
| | - Johannes Kraus
- Faculty of Mathematics, University of Duisburg-Essen, Essen, Germany
| |
Collapse
|
20
|
Bhandari A, Peng C, Dziedzic J, Anton L, Owen JR, Kramer D, Skylaris CK. Electrochemistry from first-principles in the grand canonical ensemble. J Chem Phys 2021; 155:024114. [PMID: 34266248 DOI: 10.1063/5.0056514] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Progress in electrochemical technologies, such as automotive batteries, supercapacitors, and fuel cells, depends greatly on developing improved charged interfaces between electrodes and electrolytes. The rational development of such interfaces can benefit from the atomistic understanding of the materials involved by first-principles quantum mechanical simulations with Density Functional Theory (DFT). However, such simulations are typically performed on the electrode surface in the absence of its electrolyte environment and at constant charge. We have developed a new hybrid computational method combining DFT and the Poisson-Boltzmann equation (P-BE) capable of simulating experimental electrochemistry under potential control in the presence of a solvent and an electrolyte. The charged electrode is represented quantum-mechanically via linear-scaling DFT, which can model nanoscale systems with thousands of atoms and is neutralized by a counter electrolyte charge via the solution of a modified P-BE. Our approach works with the total free energy of the combined multiscale system in a grand canonical ensemble of electrons subject to a constant electrochemical potential. It is calibrated with respect to the reduction potential of common reference electrodes, such as the standard hydrogen electrode and the Li metal electrode, which is used as a reference electrode in Li-ion batteries. Our new method can be used to predict electrochemical properties under constant potential, and we demonstrate this in exemplar simulations of the differential capacitance of few-layer graphene electrodes and the charging of a graphene electrode coupled to a Li metal electrode at different voltages.
Collapse
Affiliation(s)
- Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Chao Peng
- Engineering Sciences, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Lucian Anton
- Atos UK, 71 High Holborn London WC1V 6EA, United Kingdom
| | - John R Owen
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Denis Kramer
- Engineering Sciences, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
21
|
Oğuz IC, Vassetti D, Labat F. Assessing the performances of different continuum solvation models for the calculation of hydration energies of molecules, polymers and surfaces: a comparison between the SMD, VASPsol and FDPB models. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02799-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
22
|
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
|
23
|
D'Cunha R, Crawford TD. Modeling Complex Solvent Effects on the Optical Rotation of Chiral Molecules: A Combined Molecular Dynamics and Density Functional Theory Study. J Phys Chem A 2021; 125:3095-3108. [PMID: 33829790 DOI: 10.1021/acs.jpca.1c00803] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The challenge of assigning the absolute stereochemical configuration to a chiral compound can be overcome via accurate ab initio predictions of optical rotation, a sensitive molecular property that is further complicated by solvent effects. The solvent's "chiral imprint"-the transfer of the chirality from the solute to the surrounding achiral solvent-is explored here using conformational averaging and time-dependent density-functional theory. These complex solvent effects are taken into account via simple averaging over a molecular dynamics trajectory together with the explicit quantum mechanical consideration of the solvent molecules within the solute's cybotactic region and implicit modeling of the bulk solvent. We consider several axes along which the system's optical rotation varies, including the sampling of the dynamical trajectory, the quality of the one-electron basis set, and the use of continuum solvent models to account for bulk effects.
Collapse
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States.,Molecular Sciences Software Institute, 1880 Pratt Drive, Suite 1100, Blacksburg, Virginia 24060, United States
| |
Collapse
|
24
|
Huang J, Chen S, Eikerling M. Grand-Canonical Model of Electrochemical Double Layers from a Hybrid Density-Potential Functional. J Chem Theory Comput 2021; 17:2417-2430. [PMID: 33787259 DOI: 10.1021/acs.jctc.1c00098] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A hybrid density-potential functional of an electrochemical interface that encompasses major effects in the contacting metal and electrolyte phases is formulated. Variational analysis of this functional yields a grand-canonical model of the electrochemical double layer (EDL). Specifically, metal electrons are described using the Thomas-Fermi-Dirac-Wigner theory of an inhomogeneous electron gas. The electrolyte solution is treated classically at the mean-field level, taking into account electrostatic interactions, ion size effects, and nonlinear solvent polarization. The model uses parametrizable force relations to describe the short-range forces between metal cationic cores, metal electrons, and electrolyte ions and solvent molecules. Therefore, the gap between the metal skeleton and the electrolyte solution, key to properties of the EDL, varies consistently as a function of the electrode potential. Partial charge transfer in the presence of ion specific adsorption is described using an Anderson-Newns type theory. This model is parametrized with density functional theory calculations, compared with experimental data, and then employed to unravel several interfacial properties of fundamental significance in electrochemistry. In particular, a closer approach of the solution phase toward the metal surface, for example, caused by a stronger ion specific adsorption, decreases the potential of zero charge and elevates the double-layer capacitance curve. In addition, the ion specific adsorption can lead to surface depolarization of ions. The present model represents a viable framework to model (reactive) EDLs under the constant potential condition, which can be used to understand multifaceted EDL effects in electrocatalysis.
Collapse
Affiliation(s)
- Jun Huang
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Shengli Chen
- Hubei Key Laboratory of Electrochemical Power Sources, Key Laboratory of Analytical Chemistry for Biology and Medicine (Ministry of Education), Department of Chemistry, Wuhan University, Wuhan 430072, China
| | - Michael Eikerling
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.,Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Engineering, RWTH Aachen University, Aachen 52062, Germany
| |
Collapse
|
25
|
Herbert JM. Dielectric continuum methods for quantum chemistry. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1519] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- John M. Herbert
- Department of Chemistry and Biochemistry The Ohio State University Columbus Ohio USA
| |
Collapse
|
26
|
Mao Y, Loipersberger M, Kron KJ, Derrick JS, Chang CJ, Sharada SM, Head-Gordon M. Consistent inclusion of continuum solvation in energy decomposition analysis: theory and application to molecular CO 2 reduction catalysts. Chem Sci 2020; 12:1398-1414. [PMID: 34163903 PMCID: PMC8179122 DOI: 10.1039/d0sc05327a] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2020] [Accepted: 11/26/2020] [Indexed: 12/13/2022] Open
Abstract
To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to molecular CO2 reduction catalysis. For [FeTPP(CO2-κC)]2- (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO2 binding, -N(CH3)3 + (TMA) and -OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The coulombic interaction between the positively charged TMA group and activated CO2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO2 via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible "pockets" in that electrostatics is unscreened. For the reactant and product complexes associated with the electron transfer from the p-terphenyl radical anion to CO2, we demonstrate that the double terminal substitution of p-terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution.
Collapse
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
| | | | - Kareesa J Kron
- Mork Family Department of Chemical Engineering and Material Science, University of Southern California Los Angeles CA 90089 USA
| | - Jeffrey S Derrick
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
| | - Christopher J Chang
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
- Department of Molecular and Cell Biology, University of California Berkeley Berkeley CA 94720 USA
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Material Science, University of Southern California Los Angeles CA 90089 USA
- Department of Chemistry, University of Southern California Los Angeles CA 90089 USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| |
Collapse
|
27
|
Warburton RE, Hutchison P, Jackson MN, Pegis ML, Surendranath Y, Hammes-Schiffer S. Interfacial Field-Driven Proton-Coupled Electron Transfer at Graphite-Conjugated Organic Acids. J Am Chem Soc 2020; 142:20855-20864. [DOI: 10.1021/jacs.0c10632] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Robert E. Warburton
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Phillips Hutchison
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Megan N. Jackson
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael L. Pegis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Yogesh Surendranath
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
28
|
Bononi FC, Chen Z, Rocca D, Andreussi O, Hullar T, Anastasio C, Donadio D. Bathochromic Shift in the UV–Visible Absorption Spectra of Phenols at Ice Surfaces: Insights from First-Principles Calculations. J Phys Chem A 2020; 124:9288-9298. [DOI: 10.1021/acs.jpca.0c07038] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Fernanda C. Bononi
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| | - Zekun Chen
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| | - Dario Rocca
- Université de Lorraine, CNRS, LPTC, F-54000 Nancy, France
| | - Oliviero Andreussi
- Department of Physics, University of North Texas Denton, Texas 76203, United States
| | - Ted Hullar
- Department of Land, Air and Water Resources, University of California Davis Davis, California 95616-8627, United States
| | - Cort Anastasio
- Department of Land, Air and Water Resources, University of California Davis Davis, California 95616-8627, United States
| | - Davide Donadio
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| |
Collapse
|
29
|
Collinge G, Yuk SF, Nguyen MT, Lee MS, Glezakou VA, Rousseau R. Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations. ACS Catal 2020. [DOI: 10.1021/acscatal.0c01501] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Greg Collinge
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Simuck F. Yuk
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Manh-Thuong Nguyen
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Mal-Soon Lee
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Vassiliki-Alexandra Glezakou
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Roger Rousseau
- Basic & Applied Molecular Foundations, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
30
|
Kühne TD, Iannuzzi M, Del Ben M, Rybkin VV, Seewald P, Stein F, Laino T, Khaliullin RZ, Schütt O, Schiffmann F, Golze D, Wilhelm J, Chulkov S, Bani-Hashemian MH, Weber V, Borštnik U, Taillefumier M, Jakobovits AS, Lazzaro A, Pabst H, Müller T, Schade R, Guidon M, Andermatt S, Holmberg N, Schenter GK, Hehn A, Bussy A, Belleflamme F, Tabacchi G, Glöß A, Lass M, Bethune I, Mundy CJ, Plessl C, Watkins M, VandeVondele J, Krack M, Hutter J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J Chem Phys 2020; 152:194103. [PMID: 33687235 DOI: 10.1063/5.0007045] [Citation(s) in RCA: 986] [Impact Index Per Article: 246.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
Collapse
Affiliation(s)
- Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Mauro Del Ben
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Teodoro Laino
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, CH-801 Sherbrooke St. West, Montreal, Quebec H3A 0B8, Canada
| | - Ole Schütt
- Department of Materials, ETH Zürich, CH-8092 Zürich, Switzerland
| | | | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, Universitätsstraße 31, D-93053 Regensburg, Germany
| | - Sergey Chulkov
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | | | - Valéry Weber
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | | | | | | | | | - Hans Pabst
- Intel Extreme Computing, Software and Systems, Zürich, Switzerland
| | - Tiziano Müller
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Robert Schade
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Manuel Guidon
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Samuel Andermatt
- Integrated Systems Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Nico Holmberg
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Anna Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Gloria Tabacchi
- Department of Science and High Technology, University of Insubria and INSTM, via Valleggio 9, I-22100 Como, Italy
| | - Andreas Glöß
- BASF SE, Carl-Bosch-Straße 38, D-67056 Ludwigshafen am Rhein, Germany
| | - Michael Lass
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Iain Bethune
- Hartree Centre, Science and Technology Facilities Council, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Christian Plessl
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Matt Watkins
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | - Joost VandeVondele
- Swiss National Supercomputing Centre (CSCS), ETH Zürich, Zürich, Switzerland
| | - Matthias Krack
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
31
|
Prentice JCA, Aarons J, Womack JC, Allen AEA, Andrinopoulos L, Anton L, Bell RA, Bhandari A, Bramley GA, Charlton RJ, Clements RJ, Cole DJ, Constantinescu G, Corsetti F, Dubois SMM, Duff KKB, Escartín JM, Greco A, Hill Q, Lee LP, Linscott E, O'Regan DD, Phipps MJS, Ratcliff LE, Serrano ÁR, Tait EW, Teobaldi G, Vitale V, Yeung N, Zuehlsdorff TJ, Dziedzic J, Haynes PD, Hine NDM, Mostofi AA, Payne MC, Skylaris CK. The ONETEP linear-scaling density functional theory program. J Chem Phys 2020; 152:174111. [PMID: 32384832 DOI: 10.1063/5.0004445] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Jolyon Aarons
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Alice E A Allen
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Lampros Andrinopoulos
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Lucian Anton
- UKAEA, Culham Science Centre, Abingdon OX14 3DB, United Kingdom
| | - Robert A Bell
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gabriel A Bramley
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Robert J Charlton
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Rebecca J Clements
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Gabriel Constantinescu
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Fabiano Corsetti
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Simon M-M Dubois
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
| | - Kevin K B Duff
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - José María Escartín
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Andrea Greco
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Quintin Hill
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Louis P Lee
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Edward Linscott
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David D O'Regan
- School of Physics, AMBER, and CRANN Institute, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Maximillian J S Phipps
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Álvaro Ruiz Serrano
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Edward W Tait
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gilberto Teobaldi
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Valerio Vitale
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nelson Yeung
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Peter D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Arash A Mostofi
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
32
|
Schwarz K, Sundararaman R. The electrochemical interface in first-principles calculations. SURFACE SCIENCE REPORTS 2020; 75:10.1016/j.surfrep.2020.100492. [PMID: 34194128 PMCID: PMC8240516 DOI: 10.1016/j.surfrep.2020.100492] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
First-principles predictions play an important role in understanding chemistry at the electrochemical interface. Electronic structure calculations are straightforward for vacuum interfaces, but do not easily account for the interfacial fields and solvation that fundamentally change the nature of electrochemical reactions. Prevalent techniques for first-principles prediction of electrochemical processes range from expensive explicit solvation using ab initio molecular dynamics, through a hierarchy of continuum solvation techniques, to neglecting solvation and interfacial field effects entirely. Currently, no single approach reliably captures all relevant effects of the electrochemical double layer in first-principles calculations. This review systematically lays out the relation between all major approaches to first-principles electrochemistry, including the key approximations and their consequences for accuracy and computational cost. Focusing on ab initio methods for thermodynamic properties of aqueous interfaces, we first outline general considerations for modeling electrochemical interfaces, including solvent and electrolyte dynamics and electrification. We then present the specifics of various explicit and implicit models of the solvent and electrolyte. Finally, we discuss the compromise between computational efficiency and accuracy, and identify key outstanding challenges and future opportunities in the wide range of techniques for first-principles electrochemistry.
Collapse
Affiliation(s)
- Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Dr., Gaithersburg, Maryland 20899, USA
| | - Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th St., Troy, New York 12180, USA
| |
Collapse
|
33
|
Bramley G, Nguyen MT, Glezakou VA, Rousseau R, Skylaris CK. Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study. J Chem Theory Comput 2020; 16:2703-2715. [DOI: 10.1021/acs.jctc.0c00034] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Gabriel Bramley
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K
| | - Manh-Thuong Nguyen
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - Roger Rousseau
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | |
Collapse
|
34
|
Li LF, Li YF, Liu ZP. Oxygen Evolution Activity on NiOOH Catalysts: Four-Coordinated Ni Cation as the Active Site and the Hydroperoxide Mechanism. ACS Catal 2020. [DOI: 10.1021/acscatal.9b04975] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Li-Fen Li
- Collaborative Innovation Center of Chemistry for Energy Material, Key Laboratory of Computational Physical Science (Ministry of Education), Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Ye-Fei Li
- Collaborative Innovation Center of Chemistry for Energy Material, Key Laboratory of Computational Physical Science (Ministry of Education), Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Zhi-Pan Liu
- Collaborative Innovation Center of Chemistry for Energy Material, Key Laboratory of Computational Physical Science (Ministry of Education), Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
35
|
Yang XH, Cuesta A, Cheng J. Computational Ag/AgCl Reference Electrode from Density Functional Theory-Based Molecular Dynamics. J Phys Chem B 2019; 123:10224-10232. [PMID: 31693366 DOI: 10.1021/acs.jpcb.9b06650] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have developed a scheme to compute the standard potential of the Ag/AgCl reference electrode using density functional theory-based molecular dynamics, similar to the computational standard hydrogen electrode (SHE) developed by Cheng, Sulpizi, and Sprik [J. Chem. Phys. 2009, 131, 154504], with which our new computational reference electrode was compared. We have obtained a similar value of the potential of the Ag/AgCl electrode versus SHE to the experiment. The newly developed computational reference electrode will be extended to nonaqueous solvents in the future, where it will be used to predict standard equilibrium potentials to be compared with experimental data.
Collapse
Affiliation(s)
- Xiao-Hui Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering , Xiamen University , Xiamen 361005 , China.,Department of Chemistry , University of Aberdeen , Aberdeen AB24 3UE , U.K
| | - Angel Cuesta
- Department of Chemistry , University of Aberdeen , Aberdeen AB24 3UE , U.K
| | - Jun Cheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering , Xiamen University , Xiamen 361005 , China
| |
Collapse
|
36
|
Abstract
Bilayer graphene consists of two stacked graphene layers bound together by van der Waals interaction. As the molecular analog of bilayer graphene, molecular bilayer graphene (MBLG) can offer useful insights into the structural and functional properties of bilayer graphene. However, synthesis of MBLG, which requires discrete assembly of two graphene fragments, has proved to be challenging. Here, we show the synthesis and characterization of two structurally well-defined MBLGs, both consisting of two π-π stacked nanographene sheets. We find they have excellent stability against variation of concentration, temperature and solvents. The MBLGs show sharp absorption and emission peaks, and further time-resolved spectroscopic studies reveal drastically different lifetimes for the bright and dark Davydov states in these MBLGs.
Collapse
|
37
|
Shan N, Liu B. Elucidating Molecular Interactions in Glycerol Adsorption at the Metal-Water Interface with Density Functional Theory. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:4791-4805. [PMID: 30350699 DOI: 10.1021/acs.langmuir.8b02385] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Glycerol is an extremely versatile platform molecule for chemical and fuel production, as evidenced by successful demonstrations in electrochemical and thermochemical processes, where key catalytic chemistries occur at the solid-liquid interface. Despite the remarkable progress made in enriching the first-principles-based computational tool set to reveal and characterize solvent structures in the past decade, techniques for realistic and efficient molecular-level modeling to study aqueous-phase glycerol chemistry are still far from mature. Many aqueous-phase catalytic systems are deemed too complex for routine modeling because of their highly correlated structures at the heterogeneous solid-liquid interface. This invited feature article merges recent developments in quantum mechanical solvation models and oxygenated hydrocarbon conversion chemistry by revisiting the molecular interactions of adsorbed glycerol and its dehydrogenation intermediates at the water-metal interface. Explicit participation of water through the establishment of water-adsorbate, water-water, and water-metal interactions on Pt(111) was investigated using density functional theory. In periodic models, the adsorption favors networklike structures with adsorbates as nodal points linked by coadsorbed water molecules. We also showed that these adsorption patterns actually preserve the original bond-order-based scaling relationship framework established without the consideration of solvent. This behavior can be exploited to improve computational efficiency for future analysis of catalytic polyol conversions in the aqueous-phase environment.
Collapse
Affiliation(s)
- Nannan Shan
- Department of Chemical Engineering , Kansas State University , Manhattan , Kansas 66506 , United States
| | - Bin Liu
- Department of Chemical Engineering , Kansas State University , Manhattan , Kansas 66506 , United States
| |
Collapse
|
38
|
Affiliation(s)
- Matthew Truscott
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
| | - Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
| |
Collapse
|
39
|
Magnussen OM, Groß A. Toward an Atomic-Scale Understanding of Electrochemical Interface Structure and Dynamics. J Am Chem Soc 2019; 141:4777-4790. [DOI: 10.1021/jacs.8b13188] [Citation(s) in RCA: 145] [Impact Index Per Article: 29.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Olaf M. Magnussen
- Institute of Experimental and Applied Physics, Kiel University, Olshausenstr. 40, 24098 Kiel, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm, Germany
- Helmholtz-Institute Ulm, Helmholtzstr. 11, 89081 Ulm, Germany
| |
Collapse
|
40
|
Andreussi O, Hörmann NG, Nattino F, Fisicaro G, Goedecker S, Marzari N. Solvent-Aware Interfaces in Continuum Solvation. J Chem Theory Comput 2019; 15:1996-2009. [DOI: 10.1021/acs.jctc.8b01174] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Oliviero Andreussi
- Department of Physics, University of North Texas, Denton, Texas 76207, United States
- Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Station 9, CH-1015 Lausanne, Switzerland
| | - Nicolas Georg Hörmann
- Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Station 9, CH-1015 Lausanne, Switzerland
| | - Francesco Nattino
- Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Station 9, CH-1015 Lausanne, Switzerland
| | - Giuseppe Fisicaro
- Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
- Istituto per la Microelettronica e Microsistemi (CNR-IMM), VIII Strada 5, 95121 Catania, Italy
| | - Stefan Goedecker
- Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland
| | - Nicola Marzari
- Theory and Simulations of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Station 9, CH-1015 Lausanne, Switzerland
| |
Collapse
|
41
|
Castañeda Medina A, Schmid R. High Order Compact Multigrid Solver for Implicit Solvation Models. J Chem Theory Comput 2019; 15:1293-1301. [DOI: 10.1021/acs.jctc.8b00774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Arcesio Castañeda Medina
- Computational Materials Chemistry Group, Lehrstuhl für Anorganische Chemie 2, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Rochus Schmid
- Computational Materials Chemistry Group, Lehrstuhl für Anorganische Chemie 2, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
42
|
Labat F, Civalleri B, Dovesi R. Implicit Solvation Using a Generalized Finite-Difference Approach in CRYSTAL: Implementation and Results for Molecules, Polymers, and Surfaces. J Chem Theory Comput 2018; 14:5969-5983. [PMID: 30347161 DOI: 10.1021/acs.jctc.8b00762] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the implementation of an implicit solvation model in the CRYSTAL code. The solvation energy is separated into two components: the electrostatic contribution arising from a self-consistent reaction field treatment obtained within a generalized finite-difference Poisson model, augmented by a nonelectrostatic contribution proportional to the solvent-accessible surface area of the solute. A discontinuous dielectric boundary is used, along with a solvent-excluded surface built from interlocking atom-centered spheres on which apparent surface point charges are mapped. The procedure is general and can be performed at both the Hartree-Fock and density functional theory levels, with pure or hybrid functionals, for systems periodic in 0, 1, and 2 directions, that is, for isolated molecules and extended polymers and surfaces. The Poisson equation resolution and apparent surface charge formalism is first validated on model analytical test cases. The good agreement obtained on solvation free energies is further confirmed by calculations performed on a large test set of 501 neutral molecules, for which a mean unsigned error of 1.3 kcal/mol is obtained when compared to the available experimental data. Importantly, the self-consistent reaction field procedure converges well for all molecules tested. This is further verified for all polymers and surfaces considered. In particular, for periodic systems, results obtained on an infinite glycine chain and on the wettability parameters of SiO2 surfaces are in good agreement with previously published data. The size extensivity of the energetic terms involved in the electrostatic contribution to the solvation energy is also well verified. These encouraging results constitute a first step to take into account complex environments in the CRYSTAL code, potentially allowing for a more accurate modeling of complex processes for both periodic and nonperiodic systems.
Collapse
Affiliation(s)
- Frédéric Labat
- PSL Research University, Chimie Paristech-CNRS , Institut de Recherche de Chimie de Paris , 11 rue P. et M. Curie , 75005 Paris , France
| | - Bartolomeo Civalleri
- Dipartimento di Chimica IFM , Università di Torino and NIS - Nanostructured Interfaces and Surfaces - Centre of Excellence , Via P. Giuria 7 , 10125 Torino , Italy
| | - Roberto Dovesi
- Dipartimento di Chimica IFM , Università di Torino and NIS - Nanostructured Interfaces and Surfaces - Centre of Excellence , Via P. Giuria 7 , 10125 Torino , Italy
| |
Collapse
|
43
|
Howard JC, Crawford TD. Calculating Optical Rotatory Dispersion Spectra in Solution Using a Smooth Dielectric Model. J Phys Chem A 2018; 122:8557-8564. [DOI: 10.1021/acs.jpca.8b07803] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- J. Coleman Howard
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
44
|
Wang H, An W, Liu X, Heath Turner C. Oxygen reduction reaction on Pt(1 1 1), Pt(2 2 1), and Ni/Au1Pt3(2 2 1) surfaces: Probing scaling relationships of reaction energetics and interfacial composition. Chem Eng Sci 2018. [DOI: 10.1016/j.ces.2018.03.054] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
45
|
Cantu DC, Padmaperuma AB, Nguyen MT, Akhade SA, Yoon Y, Wang YG, Lee MS, Glezakou VA, Rousseau R, Lilga MA. A Combined Experimental and Theoretical Study on the Activity and Selectivity of the Electrocatalytic Hydrogenation of Aldehydes. ACS Catal 2018. [DOI: 10.1021/acscatal.8b00858] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- David C. Cantu
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Asanga B. Padmaperuma
- Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Manh-Thuong Nguyen
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Sneha A. Akhade
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Yeohoon Yoon
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Yang-Gang Wang
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Mal-Soon Lee
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Vassiliki-Alexandra Glezakou
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Roger Rousseau
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| | - Michael A. Lilga
- Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland Washington 99352, United States
| |
Collapse
|
46
|
Coons MP, Herbert JM. Quantum chemistry in arbitrary dielectric environments: Theory and implementation of nonequilibrium Poisson boundary conditions and application to compute vertical ionization energies at the air/water interface. J Chem Phys 2018; 148:222834. [DOI: 10.1063/1.5023916] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
47
|
Groß A. Fundamental Challenges for Modeling Electrochemical Energy Storage Systems at the Atomic Scale. Top Curr Chem (Cham) 2018; 376:17. [DOI: 10.1007/s41061-018-0194-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2017] [Accepted: 03/23/2018] [Indexed: 10/17/2022]
|
48
|
Pache D, Schmid R. Molecular Dynamics Investigation of the Dielectric Decrement of Ion Solutions. ChemElectroChem 2018. [DOI: 10.1002/celc.201800158] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Dennis Pache
- Computational Materials Chemistry Research Group, Department of Inorganic Chemistry II; Ruhr University Bochum; Germany
| | - Rochus Schmid
- Computational Materials Chemistry Research Group, Department of Inorganic Chemistry II; Ruhr University Bochum; Germany
| |
Collapse
|
49
|
Womack JC, Anton L, Dziedzic J, Hasnip PJ, Probert MIJ, Skylaris CK. DL_MG: A Parallel Multigrid Poisson and Poisson-Boltzmann Solver for Electronic Structure Calculations in Vacuum and Solution. J Chem Theory Comput 2018; 14:1412-1432. [PMID: 29447447 DOI: 10.1021/acs.jctc.7b01274] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The solution of the Poisson equation is a crucial step in electronic structure calculations, yielding the electrostatic potential-a key component of the quantum mechanical Hamiltonian. In recent decades, theoretical advances and increases in computer performance have made it possible to simulate the electronic structure of extended systems in complex environments. This requires the solution of more complicated variants of the Poisson equation, featuring nonhomogeneous dielectric permittivities, ionic concentrations with nonlinear dependencies, and diverse boundary conditions. The analytic solutions generally used to solve the Poisson equation in vacuum (or with homogeneous permittivity) are not applicable in these circumstances, and numerical methods must be used. In this work, we present DL_MG, a flexible, scalable, and accurate solver library, developed specifically to tackle the challenges of solving the Poisson equation in modern large-scale electronic structure calculations on parallel computers. Our solver is based on the multigrid approach and uses an iterative high-order defect correction method to improve the accuracy of solutions. Using two chemically relevant model systems, we tested the accuracy and computational performance of DL_MG when solving the generalized Poisson and Poisson-Boltzmann equations, demonstrating excellent agreement with analytic solutions and efficient scaling to ∼109 unknowns and 100s of CPU cores. We also applied DL_MG in actual large-scale electronic structure calculations, using the ONETEP linear-scaling electronic structure package to study a 2615 atom protein-ligand complex with routinely available computational resources. In these calculations, the overall execution time with DL_MG was not significantly greater than the time required for calculations using a conventional FFT-based solver.
Collapse
Affiliation(s)
- James C Womack
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , United Kingdom
| | - Lucian Anton
- Cray U.K. Ltd. , Broad Quay House, Prince Street , Bristol BS1 4DJ , United Kingdom
| | - Jacek Dziedzic
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , United Kingdom.,Faculty of Applied Physics and Mathematics , Gdańsk University of Technology , Gdańsk 80-233 , Poland
| | - Phil J Hasnip
- Department of Physics , University of York , Heslington, York YO10 5DD , United Kingdom
| | - Matt I J Probert
- Department of Physics , University of York , Heslington, York YO10 5DD , United Kingdom
| | - Chris-Kriton Skylaris
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , United Kingdom
| |
Collapse
|
50
|
Yoo SH, Todorova M, Neugebauer J. Selective Solvent-Induced Stabilization of Polar Oxide Surfaces in an Electrochemical Environment. PHYSICAL REVIEW LETTERS 2018; 120:066101. [PMID: 29481276 DOI: 10.1103/physrevlett.120.066101] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2017] [Revised: 12/10/2017] [Indexed: 06/08/2023]
Abstract
The impact of an electrochemical environment on the thermodynamic stability of polar oxide surfaces is investigated for the example of ZnO(0001) surfaces immersed in water using density functional theory calculations. We show that solvation effects are highly selective: They have little effect on surfaces showing a metallic character, but largely stabilize semiconducting structures, particularly those that have a high electrostatic penalty in vacuum. The high selectivity is shown to have direct consequences for the surface phase diagram and explains, e.g., why certain surface structures could be observed only in an electrochemical environment.
Collapse
Affiliation(s)
- Su-Hyun Yoo
- Department of Computational Materials Design, Max-Planck-Institut für Eisenforschung GmbH, Max-Planck-Str. 1, D-40237 Düsseldorf, Germany
| | - Mira Todorova
- Department of Computational Materials Design, Max-Planck-Institut für Eisenforschung GmbH, Max-Planck-Str. 1, D-40237 Düsseldorf, Germany
| | - Jörg Neugebauer
- Department of Computational Materials Design, Max-Planck-Institut für Eisenforschung GmbH, Max-Planck-Str. 1, D-40237 Düsseldorf, Germany
| |
Collapse
|