1
|
Tang F, Shi K, Wu X. Exploring the impact of ions on oxygen K-edge X-ray absorption spectroscopy in NaCl solution using the GW-Bethe-Salpeter-equation approach. J Chem Phys 2023; 159:174501. [PMID: 37909453 DOI: 10.1063/5.0167999] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2023] [Accepted: 10/13/2023] [Indexed: 11/03/2023] Open
Abstract
X-ray absorption spectroscopy (XAS) is a powerful experimental tool to probe the local structure in materials with the core hole excitations. Here, the oxygen K-edge XAS spectra of the NaCl solution and pure water are computed by using a recently developed GW-Bethe-Salpeter equation approach, based on configurations modeled by path-integral molecular dynamics with the deep-learning technique. The neural network is trained on ab initio data obtained with strongly constrained and appropriately normed density functional theory. The observed changes in the XAS features of the NaCl solution, compared to those of pure water, are in good agreement between experimental and theoretical results. We provided detailed explanations for these spectral changes that occur when NaCl is solvated in pure water. Specifically, the presence of solvating ion pairs leads to localization of electron-hole excitons. Our theoretical XAS results support the theory that the effects of the solvating ions on the H-bond network are mainly confined within the first hydration shell of ions, however beyond the shell the arrangement of water molecules remains to be comparable to that observed in pure water.
Collapse
Affiliation(s)
- Fujie Tang
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Kefeng Shi
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
2
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
3
|
Bajaj A, Duan C, Nandy A, Taylor MG, Kulik HJ. Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry. J Chem Phys 2022; 156:184112. [DOI: 10.1063/5.0089460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d1 MOPs to be suitable for nearly eliminating all energetic delocalization error and static correlation error. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization error and static correlation error within a non-empirical framework.
Collapse
Affiliation(s)
- Akash Bajaj
- Massachusetts Institute of Technology, United States of America
| | - Chenru Duan
- Massachusetts Institute of Technology, United States of America
| | - Aditya Nandy
- Massachusetts Institute of Technology, United States of America
| | | | - Heather J. Kulik
- Dept of Chemical Engineering, Massachusetts Institute of Technology, United States of America
| |
Collapse
|
4
|
Bajaj A, Kulik HJ. Molecular DFT+U: A Transferable, Low-Cost Approach to Eliminate Delocalization Error. J Phys Chem Lett 2021; 12:3633-3640. [PMID: 33826346 DOI: 10.1021/acs.jpclett.1c00796] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
While density functional theory (DFT) is widely applied for its combination of cost and accuracy, corrections (e.g., DFT+U) that improve it are often needed to tackle correlated transition-metal chemistry. In principle, the functional form of DFT+U, consisting of a set of localized atomic orbitals (AOs) and a quadratic energy penalty for deviation from integer occupations of those AOs, enables the recovery of the exact conditions of piecewise linearity and the derivative discontinuity. Nevertheless, for practical transition-metal complexes, where both atomic states and ligand orbitals participate in bonding, standard DFT+U can fail to eliminate delocalization error (DE). Here, we show that by introducing an alternative valence-state (i.e., molecular orbital or MO) basis to the DFT+U approach, we recover exact conditions in cases for which standard DFT+U corrections have no error-reducing effect. This MO-based DFT+U also eliminates DE where standard AO-based DFT+U is already successful. We demonstrate the transferability of our approach on representative transition-metal complexes with a range of ligand field strengths, electron configurations (i.e., from Sc to Zn), and spin states.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
5
|
Marin TW, Janik I, Bartels DM, Chipman DM. Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum. Phys Chem Chem Phys 2021; 23:9109-9120. [PMID: 33885094 DOI: 10.1039/d1cp00930c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl-, Na+, and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl-.
Collapse
Affiliation(s)
- Timothy W Marin
- Department of Physical Sciences, Benedictine University, 5700 College Rd, Lisle, IL 60532, USA
| | | | | | | |
Collapse
|
6
|
Hirakawa T, Bowler DR, Miyazaki T, Morikawa Y, Truflandier LA. Blue moon ensemble simulation of aquation free energy profiles applied to mono and bifunctional platinum anticancer drugs. J Comput Chem 2020; 41:1973-1984. [DOI: 10.1002/jcc.26367] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 05/21/2020] [Accepted: 06/01/2020] [Indexed: 12/31/2022]
Affiliation(s)
- Teruo Hirakawa
- Department of Precision EngineeringGraduate School of Engineering, Osaka University Suita Osaka Japan
- Institut des Sciences Moléculaires (ISM), Université Bordeaux Talence Cedex France
| | - David R. Bowler
- Department of Physics & AstronomyUniversity College London (UCL) London United Kingdom
- London Centre for Nanotechnology, UCL London United Kingdom
- International Centre for Materials Nanoarchitechtonics (WPI‐MANA), National Institute for Materials Science (NIMS) Tsukuba Ibaraki Japan
| | - Tsuyoshi Miyazaki
- International Centre for Materials Nanoarchitechtonics (WPI‐MANA), National Institute for Materials Science (NIMS) Tsukuba Ibaraki Japan
| | - Yoshitada Morikawa
- Department of Precision EngineeringGraduate School of Engineering, Osaka University Suita Osaka Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University Kyoto Japan
- Research Center for Ultra‐Precision Science and TechnologyGraduate School of Engineering, Osaka University Suita Osaka Japan
| | - Lionel A. Truflandier
- Department of Precision EngineeringGraduate School of Engineering, Osaka University Suita Osaka Japan
- Institut des Sciences Moléculaires (ISM), Université Bordeaux Talence Cedex France
| |
Collapse
|
7
|
Ko HY, Jia J, Santra B, Wu X, Car R, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance. J Chem Theory Comput 2020; 16:3757-3785. [PMID: 32045232 DOI: 10.1021/acs.jctc.9b01167] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Princeton University, Princeton, New Jersey 08544, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
8
|
Rozsa V, Pham TA, Galli G. Molecular polarizabilities as fingerprints of perturbations to water by ions and confinement. J Chem Phys 2020; 152:124501. [DOI: 10.1063/1.5143317] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Viktor Rozsa
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - Tuan Anh Pham
- Quantum Simulations Group and Laboratory for Energy Applications for the Future (LEAF), Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Giulia Galli
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, USA
- Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
9
|
DelloStritto M, Xu J, Wu X, Klein ML. Aqueous solvation of the chloride ion revisited with density functional theory: impact of correlation and exchange approximations. Phys Chem Chem Phys 2020; 22:10666-10675. [DOI: 10.1039/c9cp06821j] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Aqueous chloride is simulated using PBE-D3, PBE0-D3, and SCAN to investigate the impact of exchange and correlation approximations; we find the exact exchange fraction strongly impacts the energetics and polarizability of solvated chloride.
Collapse
Affiliation(s)
- Mark DelloStritto
- Institute for Computational Molecular Science
- Temple University SERC
- Philadelphia
- USA
| | - Jianhang Xu
- Department of Physics
- Temple University SERC
- Philadelphia
- USA
| | - Xifan Wu
- Department of Physics
- Temple University SERC
- Philadelphia
- USA
| | - Michael L. Klein
- Institute for Computational Molecular Science
- Temple University SERC
- Philadelphia
- USA
| |
Collapse
|
10
|
Awoonor-Williams E, Isley WC, Dale SG, Johnson ER, Yu H, Becke AD, Roux B, Rowley CN. Quantum Chemical Methods for Modeling Covalent Modification of Biological Thiols. J Comput Chem 2019; 41:427-438. [PMID: 31512279 DOI: 10.1002/jcc.26064] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Revised: 07/24/2019] [Accepted: 08/16/2019] [Indexed: 02/06/2023]
Abstract
Targeted covalent inhibitor drugs require computational methods that go beyond simple molecular-mechanical force fields in order to model the chemical reactions that occur when they bind to their targets. Here, several semiempirical and density-functional theory (DFT) methods are assessed for their ability to describe the potential energy surface and reaction energies of the covalent modification of a thiol by an electrophile. Functionals such as PBE and B3LYP fail to predict a stable enolate intermediate. This is largely due to delocalization error, which spuriously stabilizes the prereaction complex, in which excess electron density is transferred from the thiolate to the electrophile. Functionals with a high-exact exchange component, range-separated DFT functionals, and variationally optimized exact exchange (i.e., the LC-B05minV functional) correct this issue to various degrees. The large gradient behavior of the exchange enhancement factor is also found to significantly affect the results, leading to the improved performance of PBE0. While ωB97X-D and M06-2X were reasonably accurate, no method provided quantitative accuracy for all three electrophiles, making this a very strenuous test of functional performance. Additionally, one drawback of M06-2X was that molecular dynamics (MD) simulations using this functional were only stable if a fine integration grid was used. The low-cost semiempirical methods, PM3, AM1, and PM7, provide a qualitatively correct description of the reaction mechanism, although the energetics is not quantitatively reliable. As a proof of concept, the potential of mean force for the addition of methylthiolate to methylvinyl ketone was calculated using quantum mechanical/molecular mechanical MD in an explicit polarizable aqueous solvent. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ernest Awoonor-Williams
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador, Canada
| | - William C Isley
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois
| | - Stephen G Dale
- Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia
| | - Erin R Johnson
- Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Haibo Yu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong, New South Wales, Australia
| | - Axel D Becke
- Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador, Canada
| |
Collapse
|
11
|
Muchová E, Slavíček P. Beyond Koopmans' theorem: electron binding energies in disordered materials. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:043001. [PMID: 30524069 DOI: 10.1088/1361-648x/aaf130] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The topical review focuses on calculating ionization energies (IE), or electronic polarons in quasi-particle terminology, in large disordered systems, e.g. for a solute dissolved in a molecular solvent. The simplest estimate of the ionization energy is provided by one-electron energies in the Hartree-Fock theory, but the calculated quantities are not accurate. Density functional theory as many-body theory provides a principal opportunity for calculating one-electron energies including correlation and relaxation effects, i.e. the true energies of electronic polarons. We argue that such a principal possibility materializes within the concept of optimally tuned range-separated hybrid functionals (OT-RSH). We describe various schemes for optimal tuning. Importantly, the OT-RSH scheme is investigated for systems capped with dielectric continuum, providing a consistent picture on the QM/dielectric boundary. Finally, some limitations and open issues of the OT-RSH approach are addressed.
Collapse
Affiliation(s)
- Eva Muchová
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Czech Republic
| | | |
Collapse
|
12
|
Bouchafra Y, Shee A, Réal F, Vallet V, Severo Pereira Gomes A. Predictive Simulations of Ionization Energies of Solvated Halide Ions with Relativistic Embedded Equation of Motion Coupled Cluster Theory. PHYSICAL REVIEW LETTERS 2018; 121:266001. [PMID: 30636145 DOI: 10.1103/physrevlett.121.266001] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2018] [Indexed: 06/09/2023]
Abstract
A subsystem approach for obtaining electron binding energies in the valence region is presented and applied to the case of halide ions (X^{-},X=F-At) in water. This approach is based on electronic structure calculations combining the relativistic equation-of-motion coupled cluster method for electron detachment and density functional theory via the frozen density embedding approach, using structures from classical molecular dynamics with polarizable force fields for discrete systems (in our study, droplets containing the anion and 50 water molecules). Our results indicate that one can accurately capture both the large solvent effect observed for the halides and the splitting of their ionization signals due to the increasingly large spin-orbit coupling of the p_{3/2}-p_{1/2} manifold across the series, at an affordable computational cost. Furthermore, owing to the quantum mechanical treatment of both solute and solvent electron binding energies of semiquantitative quality are also obtained for (bulk) water as by-products of the calculations for the halogens (in droplets).
Collapse
Affiliation(s)
- Yassine Bouchafra
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Avijit Shee
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Florent Réal
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Valérie Vallet
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| |
Collapse
|
13
|
Awoonor-Williams E, Rowley CN. The hydration structure of methylthiolate from QM/MM molecular dynamics. J Chem Phys 2018; 149:045103. [PMID: 30068187 DOI: 10.1063/1.5038010] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Thiols are widely present in biological systems, most notably as the side chain of cysteine amino acids in proteins. Thiols can be deprotonated to form a thiolate which affords a diverse range of enzymatic activity and modes for chemical modification of proteins. Parameters for modeling thiolates using molecular mechanical force fields have not yet been validated, in part due to the lack of structural data on thiolate solvation. Here, the CHARMM36 and Amber models for thiolates in aqueous solutions are assessed using free energy perturbation and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. The hydration structure of methylthiolate was calculated from 1 ns of QM/MM MD (PBE0-D3/def2-TZVP//TIP3P), which shows that the water-S- distances are approximately 2 Å with a coordination number near 6. The CHARMM thiolate parameters predict a thiolate S radius close to the QM/MM value and predict a hydration Gibbs energy of -329.2 kJ/mol, close to the experimental value of -318 kJ/mol. The cysteine thiolate model in the Amber force field underestimates the thiolate radius by 0.2 Å and overestimates the thiolate hydration energy by 119 kJ/mol because it uses the same Lennard-Jones parameters for thiolates as for thiols. A recent Drude polarizable model for methylthiolate with optimized thiolate parameters also performs well. SAPT2+ [Symmetry Adapted Perturbation Theory (SAPT)] analysis indicates that exchange repulsion is larger for the methylthiolate, consistent with it having a more diffuse electron density distribution in comparison with the parent thiol. These data demonstrate that it is important to define distinct non-bonded parameters for the protonated/deprotonated states of amino acid side chains in molecular mechanical force fields.
Collapse
Affiliation(s)
- Ernest Awoonor-Williams
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador A1C 5S7, Canada
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador A1C 5S7, Canada
| |
Collapse
|
14
|
Gaiduk AP, Gustafson J, Gygi F, Galli G. First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional. J Phys Chem Lett 2018; 9:3068-3073. [PMID: 29768015 DOI: 10.1021/acs.jpclett.8b01017] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.
Collapse
Affiliation(s)
- Alex P Gaiduk
- Institute for Molecular Engineering , The University of Chicago , Chicago , Illinois 60637 , United States
- Materials Science Division , Argonne National Laboratory , Argonne , Illinois 60439 , United States
| | - Jeffrey Gustafson
- Department of Chemistry , The University of Chicago , Chicago , Illinois 60637 , United States
| | - François Gygi
- Department of Computer Science , University of California , Davis , California 95616 , United States
| | - Giulia Galli
- Institute for Molecular Engineering , The University of Chicago , Chicago , Illinois 60637 , United States
- Materials Science Division , Argonne National Laboratory , Argonne , Illinois 60439 , United States
- Department of Chemistry , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
15
|
Philips A, Marchenko A, Truflandier LA, Autschbach J. Quadrupolar NMR Relaxation from ab Initio Molecular Dynamics: Improved Sampling and Cluster Models versus Periodic Calculations. J Chem Theory Comput 2017; 13:4397-4409. [PMID: 28719202 DOI: 10.1021/acs.jctc.7b00584] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quadrupolar NMR relaxation rates are computed for 17O and 2H nuclei of liquid water, and of 23Na+, and 35Cl- in aqueous solution via Kohn-Sham (KS) density functional theory ab initio molecular dynamics (aiMD) and subsequent KS electric field gradient (EFG) calculations along the trajectories. The calculated relaxation rates are within about a factor of 2 of experimental results and improved over previous aiMD simulations. The relaxation rates are assessed with regard to the lengths of the simulations as well as configurational sampling. The latter is found to be the more limiting factor in obtaining good statistical sampling and is improved by averaging over many equivalent nuclei of a system or over several independent trajectories. Further, full periodic plane-wave basis calculations of the EFGs are compared with molecular-cluster atomic-orbital basis calculations. The two methods deliver comparable results with nonhybrid functionals. With the molecular-cluster approach, a larger variety of electronic structure methods is available. For chloride, the EFG computations benefit from using a hybrid KS functional.
Collapse
Affiliation(s)
- Adam Philips
- Department of Chemistry University at Buffalo, State University of New York , Buffalo, New York 14260-3000, United States
| | - Alex Marchenko
- Department of Chemistry University at Buffalo, State University of New York , Buffalo, New York 14260-3000, United States
| | - Lionel A Truflandier
- CNRS UMR 5255, Institut des Sciences Moléculaires Université Bordeaux , 351 cours de la Libération, 33405 Talence cedex, France
| | - Jochen Autschbach
- Department of Chemistry University at Buffalo, State University of New York , Buffalo, New York 14260-3000, United States
| |
Collapse
|
16
|
Brawand NP, Govoni M, Vörös M, Galli G. Performance and Self-Consistency of the Generalized Dielectric Dependent Hybrid Functional. J Chem Theory Comput 2017; 13:3318-3325. [DOI: 10.1021/acs.jctc.7b00368] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Nicholas P. Brawand
- Institute
for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Govoni
- Institute
for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Argonne National
Laboratory, Lemont, Illinois 60439, United States
| | - Márton Vörös
- Institute
for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Argonne National
Laboratory, Lemont, Illinois 60439, United States
| | - Giulia Galli
- Institute
for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Argonne National
Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
17
|
Mao Y, Shao Y, Dziedzic J, Skylaris CK, Head-Gordon T, Head-Gordon M. Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis. J Chem Theory Comput 2017; 13:1963-1979. [PMID: 28430427 DOI: 10.1021/acs.jctc.7b00089] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH3, NH4+, CN-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.
Collapse
Affiliation(s)
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma , Norman, Oklahoma 73019, United States
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K.,Faculty of Applied Physics and Mathematics, Gdańsk University of Technology , Gdańsk 80-233, Poland
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K
| | | | - Martin Head-Gordon
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
18
|
Pham TA, Ping Y, Galli G. Modelling heterogeneous interfaces for solar water splitting. NATURE MATERIALS 2017; 16:401-408. [PMID: 28068314 DOI: 10.1038/nmat4803] [Citation(s) in RCA: 113] [Impact Index Per Article: 16.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Accepted: 10/18/2016] [Indexed: 05/17/2023]
Abstract
The generation of hydrogen from water and sunlight offers a promising approach for producing scalable and sustainable carbon-free energy. The key of a successful solar-to-fuel technology is the design of efficient, long-lasting and low-cost photoelectrochemical cells, which are responsible for absorbing sunlight and driving water splitting reactions. To this end, a detailed understanding and control of heterogeneous interfaces between photoabsorbers, electrolytes and catalysts present in photoelectrochemical cells is essential. Here we review recent progress and open challenges in predicting physicochemical properties of heterogeneous interfaces for solar water splitting applications using first-principles-based approaches, and highlights the key role of these calculations in interpreting increasingly complex experiments.
Collapse
Affiliation(s)
- Tuan Anh Pham
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Yuan Ping
- Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, California 95064, USA
| | - Giulia Galli
- The Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
19
|
Pham TA, Ogitsu T, Lau EY, Schwegler E. Structure and dynamics of aqueous solutions from PBE-based first-principles molecular dynamics simulations. J Chem Phys 2016; 145:154501. [DOI: 10.1063/1.4964865] [Citation(s) in RCA: 82] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Tuan Anh Pham
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Tadashi Ogitsu
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Edmond Y. Lau
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Eric Schwegler
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| |
Collapse
|
20
|
Antalek M, Pace E, Hedman B, Hodgson KO, Chillemi G, Benfatto M, Sarangi R, Frank P. Solvation structure of the halides from x-ray absorption spectroscopy. J Chem Phys 2016; 145:044318. [PMID: 27475372 PMCID: PMC4967075 DOI: 10.1063/1.4959589] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2016] [Accepted: 07/11/2016] [Indexed: 11/14/2022] Open
Abstract
Three-dimensional models for the aqueous solvation structures of chloride, bromide, and iodide are reported. K-edge extended X-ray absorption fine structure (EXAFS) and Minuit X-ray absorption near edge (MXAN) analyses found well-defined single shell solvation spheres for bromide and iodide. However, dissolved chloride proved structurally distinct, with two solvation shells needed to explain its strikingly different X-ray absorption near edge structure (XANES) spectrum. Final solvation models were as follows: iodide, 8 water molecules at 3.60 ± 0.13 Å and bromide, 8 water molecules at 3.40 ± 0.14 Å, while chloride solvation included 7 water molecules at 3.15 ± 0.10 Å, and a second shell of 7 water molecules at 4.14 ± 0.30 Å. Each of the three derived solvation shells is approximately uniformly disposed about the halides, with no global asymmetry. Time-dependent density functional theory calculations simulating the chloride XANES spectra following from alternative solvation spheres revealed surprising sensitivity of the electronic state to 6-, 7-, or 8-coordination, implying a strongly bounded phase space for the correct structure during an MXAN fit. MXAN analysis further showed that the asymmetric solvation predicted from molecular dynamics simulations using halide polarization can play no significant part in bulk solvation. Classical molecular dynamics used to explore chloride solvation found a 7-water solvation shell at 3.12 (-0.04/+0.3) Å, supporting the experimental result. These experiments provide the first fully three-dimensional structures presenting to atomic resolution the aqueous solvation spheres of the larger halide ions.
Collapse
Affiliation(s)
- Matthew Antalek
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| | - Elisabetta Pace
- Laboratori Nazionali di Frascati-INFN, P.O. Box 13, 00044 Frascati, Italy
| | - Britt Hedman
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| | - Keith O Hodgson
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Giovanni Chillemi
- CINECA, SCAI-SuperComputing Applications and Innovation Department, Via dei Tizii 6, 00185 Roma, Italy
| | - Maurizio Benfatto
- Laboratori Nazionali di Frascati-INFN, P.O. Box 13, 00044 Frascati, Italy
| | - Ritimukta Sarangi
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| | - Patrick Frank
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| |
Collapse
|
21
|
Gaiduk AP, Govoni M, Seidel R, Skone JH, Winter B, Galli G. Photoelectron Spectra of Aqueous Solutions from First Principles. J Am Chem Soc 2016; 138:6912-5. [PMID: 27105336 DOI: 10.1021/jacs.6b00225] [Citation(s) in RCA: 56] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a combined computational and experimental study of the photoelectron spectrum of a simple aqueous solution of NaCl. Measurements were conducted on microjets, and first-principles calculations were performed using hybrid functionals and many-body perturbation theory at the G0W0 level, starting with wave functions computed in ab initio molecular dynamics simulations. We show excellent agreement between theory and experiments for the positions of both the solute and solvent excitation energies on an absolute energy scale and for peak intensities. The best comparison was obtained using wave functions obtained with dielectric-dependent self-consistent and range-separated hybrid functionals. Our computational protocol opens the way to accurate, predictive calculations of the electronic properties of electrolytes, of interest to a variety of energy problems.
Collapse
Affiliation(s)
- Alex P Gaiduk
- Institute for Molecular Engineering, The University of Chicago , Chicago, Illinois 60637, United States
| | - Marco Govoni
- Institute for Molecular Engineering, The University of Chicago , Chicago, Illinois 60637, United States.,Materials Science Division, Argonne National Laboratory , Argonne, Illinois 60439, United States
| | - Robert Seidel
- Methods for Material Development, Helmholtz-Zentrum Berlin für Materialien und Energie , D-12489 Berlin, Germany
| | - Jonathan H Skone
- Institute for Molecular Engineering, The University of Chicago , Chicago, Illinois 60637, United States
| | - Bernd Winter
- Methods for Material Development, Helmholtz-Zentrum Berlin für Materialien und Energie , D-12489 Berlin, Germany
| | - Giulia Galli
- Institute for Molecular Engineering, The University of Chicago , Chicago, Illinois 60637, United States.,Materials Science Division, Argonne National Laboratory , Argonne, Illinois 60439, United States
| |
Collapse
|
22
|
Blase X, Boulanger P, Bruneval F, Fernandez-Serra M, Duchemin I. GW and Bethe-Salpeter study of small water clusters. J Chem Phys 2016; 144:034109. [DOI: 10.1063/1.4940139] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | | | - Fabien Bruneval
- CEA, DEN, Service de Recherches de Métallurgie Physique, F-91191 Gif-sur-Yvette, France
| | - Marivi Fernandez-Serra
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA
- Institute for Advanced Computational Sciences, Stony Brook University, Stony Brook, New York 11794-3800, USA
| | - Ivan Duchemin
- INAC, SP2M/L_Sim, CEA/UJF Cedex 09, 38054 Grenoble, France
| |
Collapse
|
23
|
Gaiduk AP, Gygi F, Galli G. Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals. J Phys Chem Lett 2015; 6:2902-2908. [PMID: 26267178 DOI: 10.1021/acs.jpclett.5b00901] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water. Approximate inclusion of dispersion interactions on computed molecular-dynamics trajectories led to a substantial improvement of the PBE0 results for the density of liquid water, which, however, resulted to be slightly lower than that of ice.
Collapse
Affiliation(s)
- Alex P Gaiduk
- †Institute for Molecular Engineering, The University of Chicago, 5747 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - François Gygi
- ‡Department of Computer Science, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Giulia Galli
- †Institute for Molecular Engineering, The University of Chicago, 5747 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
24
|
Opalka D, Pham TA, Sprik M, Galli G. Electronic Energy Levels and Band Alignment for Aqueous Phenol and Phenolate from First Principles. J Phys Chem B 2015; 119:9651-60. [DOI: 10.1021/acs.jpcb.5b04189] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daniel Opalka
- Max Planck Institute
for Solid State Research, Heisenbergstraße
1, 70569 Stuttgart, Germany
| | - Tuan Anh Pham
- Lawrence Livermore
National Laboratory, Livermore, California 94551, United States
| | - Michiel Sprik
- Department
of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Giulia Galli
- The
Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
25
|
Bankura A, Santra B, DiStasio RA, Swartz CW, Klein ML, Wu X. A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1059959] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Arindam Bankura
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Charles W. Swartz
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Xifan Wu
- Department of Physics, Temple University , Philadelphia, PA, USA
| |
Collapse
|
26
|
Tse YLS, Knight C, Voth GA. An analysis of hydrated proton diffusion in ab initio molecular dynamics. J Chem Phys 2015; 142:014104. [DOI: 10.1063/1.4905077] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ying-Lung Steve Tse
- Department of Chemistry, James Franck Institute, and Computation Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Chris Knight
- Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, and Computation Institute, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
27
|
Pham TA, Lee D, Schwegler E, Galli G. Interfacial Effects on the Band Edges of Functionalized Si Surfaces in Liquid Water. J Am Chem Soc 2014; 136:17071-7. [DOI: 10.1021/ja5079865] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tuan Anh Pham
- Department
of Chemistry, University of California, Davis, California 95616, United States
- Lawrence
Livermore
National Laboratory, Livermore, California 94551, United States
| | - Donghwa Lee
- Lawrence
Livermore
National Laboratory, Livermore, California 94551, United States
| | - Eric Schwegler
- Lawrence
Livermore
National Laboratory, Livermore, California 94551, United States
| | - Giulia Galli
- The
Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
28
|
Liu X, Cheng J, Sprik M. Aqueous Transition-Metal Cations as Impurities in a Wide Gap Oxide: The Cu2+/Cu+ and Ag2+/Ag+ Redox Couples Revisited. J Phys Chem B 2014; 119:1152-63. [DOI: 10.1021/jp506691h] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Xiandong Liu
- Department
of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
- State Key Laboratory
for Mineral Deposits Research, School of Earth Sciences and Engineering, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Jun Cheng
- Department
of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
- Department
of Chemistry, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
| | - Michiel Sprik
- Department
of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
29
|
Wan Q, Spanu L, Gygi F, Galli G. Electronic Structure of Aqueous Sulfuric Acid from First-Principles Simulations with Hybrid Functionals. J Phys Chem Lett 2014; 5:2562-2567. [PMID: 26277943 DOI: 10.1021/jz501168p] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Affiliation(s)
- Quan Wan
- †Department of Chemistry, University of California Davis, One Shields Avenue, Davis, California 95616, United States
- ‡The Institute for Molecular Engineering, The University of Chicago, 5801 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Leonardo Spanu
- †Department of Chemistry, University of California Davis, One Shields Avenue, Davis, California 95616, United States
| | - Francois Gygi
- §Department of Computer Science, University of California Davis, One Shields Avenue, Davis, California 95616, United States
| | - Giulia Galli
- ‡The Institute for Molecular Engineering, The University of Chicago, 5801 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
30
|
Opalka D, Pham TA, Sprik M, Galli G. The ionization potential of aqueous hydroxide computed using many-body perturbation theory. J Chem Phys 2014; 141:034501. [DOI: 10.1063/1.4887259] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniel Opalka
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Tuan Anh Pham
- Department of Chemistry, University of California, Davis, California 95616, USA
- Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Michiel Sprik
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Giulia Galli
- The Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
31
|
Baer MD, Fulton JL, Balasubramanian M, Schenter GK, Mundy CJ. Persistent Ion Pairing in Aqueous Hydrochloric Acid. J Phys Chem B 2014; 118:7211-20. [DOI: 10.1021/jp501091h] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Marcel D. Baer
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - John L. Fulton
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - Gregory K. Schenter
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Christopher J. Mundy
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
32
|
Gaiduk AP, Zhang C, Gygi F, Galli G. Structural and electronic properties of aqueous NaCl solutions from ab initio molecular dynamics simulations with hybrid density functionals. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.04.037] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
33
|
Kuechler ER, York DM. Quantum mechanical study of solvent effects in a prototype SN2 reaction in solution: Cl- attack on CH3Cl. J Chem Phys 2014; 140:054109. [PMID: 24511924 PMCID: PMC3977776 DOI: 10.1063/1.4863344] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2013] [Accepted: 01/13/2014] [Indexed: 11/15/2022] Open
Abstract
The nucleophilic attack of a chloride ion on methyl chloride is an important prototype SN2 reaction in organic chemistry that is known to be sensitive to the effects of the surrounding solvent. Herein, we develop a highly accurate Specific Reaction Parameter (SRP) model based on the Austin Model 1 Hamiltonian for chlorine to study the effects of solvation into an aqueous environment on the reaction mechanism. To accomplish this task, we apply high-level quantum mechanical calculations to study the reaction in the gas phase and combined quantum mechanical/molecular mechanical simulations with TIP3P and TIP4P-ew water models and the resulting free energy profiles are compared with those determined from simulations using other fast semi-empirical quantum models. Both gas phase and solution results with the SRP model agree very well with experiment and provide insight into the specific role of solvent on the reaction coordinate. Overall, the newly parameterized SRP Hamiltonian is able to reproduce both the gas phase and solution phase barriers, suggesting it is an accurate and robust model for simulations in the aqueous phase at greatly reduced computational cost relative to comparably accurate ab initio and density functional models.
Collapse
Affiliation(s)
- Erich R Kuechler
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, USA
| | - Darrin M York
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, USA
| |
Collapse
|
34
|
Holmberg N, Chen JC, Foster AS, Laasonen K. Dissolution of NaCl nanocrystals: an ab initio molecular dynamics study. Phys Chem Chem Phys 2014; 16:17437-46. [DOI: 10.1039/c4cp00635f] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
NaCl nanocrystal dissolution was investigated in atomistic detail revealing a difference in the solvation of two different ionic species.
Collapse
Affiliation(s)
- Nico Holmberg
- Department of Chemistry
- Aalto University
- FI-00076 Aalto, Finland
| | - Jian-Cheng Chen
- Department of Applied Physics
- Aalto University
- FI-00076 Aalto, Finland
- COMP Centre of Excellence in Computational Nanoscience
- Aalto University
| | - Adam S. Foster
- Department of Applied Physics
- Aalto University
- FI-00076 Aalto, Finland
- COMP Centre of Excellence in Computational Nanoscience
- Aalto University
| | - Kari Laasonen
- Department of Chemistry
- Aalto University
- FI-00076 Aalto, Finland
- COMP Centre of Excellence in Computational Nanoscience
- Aalto University
| |
Collapse
|
35
|
Lucking M, Sun YY, West D, Zhang S. Absolute redox potential of liquid water: a first-principles theory. Chem Sci 2014. [DOI: 10.1039/c3sc52287c] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The absolute redox potential of liquid water has been calculated by a first-principles method with good agreement with experimental data.
Collapse
Affiliation(s)
- Michael Lucking
- Department of Physics, Applied Physics & Astronomy
- Rensselaer Polytechnic Institute
- Troy, USA
| | - Yi-Yang Sun
- Department of Physics, Applied Physics & Astronomy
- Rensselaer Polytechnic Institute
- Troy, USA
| | - Damien West
- Department of Physics, Applied Physics & Astronomy
- Rensselaer Polytechnic Institute
- Troy, USA
| | - Shengbai Zhang
- Department of Physics, Applied Physics & Astronomy
- Rensselaer Polytechnic Institute
- Troy, USA
| |
Collapse
|
36
|
Swartz CW, Wu X. Ab initio studies of ionization potentials of hydrated hydroxide and hydronium. PHYSICAL REVIEW LETTERS 2013; 111:087801. [PMID: 24010478 DOI: 10.1103/physrevlett.111.087801] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2013] [Indexed: 06/02/2023]
Abstract
The ionization potential distributions of hydrated hydroxide and hydronium are computed with the many-body approach for electron excitations with configurations generated by ab initio molecular dynamics. The experimental features are well reproduced and found to be closely related to the molecular excitations. In the stable configurations, the ionization potential is mainly perturbed by solvent water molecules within the first solvation shell. On the other hand, electron excitation is delocalized on both proton receiving and donating complex during proton transfer, which shifts the excitation energies and broadens the spectra for both hydrated ions.
Collapse
Affiliation(s)
- Charles W Swartz
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | | |
Collapse
|