1
|
Bhattacharjee R, Miró P. Aqueous Speciation of Tetravalent Actinides in the Presence of Chloride and Nitrate Ligands. Inorg Chem 2022; 61:14718-14725. [PMID: 36050286 DOI: 10.1021/acs.inorgchem.2c02064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Speciation of hexachloride tetravalent uranium, neptunium, and plutonium species in aqueous media has been investigated using density functional theory in the presence of inner sphere ligands such as chloride, nitrate, and solvent molecules. All possible structures with the formula [AnIV(Cl)x(H2O)y(NO3)z]4-x-z (An = U, Np, and Pu; x = 0-6, y = 0-8, and z = 0-6) were considered to explore the speciation chemical space of each actinide. The nature of the mixed-ligand complexes present in solution is controlled by the concentration of free ligands in solution. A low chloride concentration is suitable to drive the speciation away from the highly thermodynamically stable hexachloride species. Furthermore, the formation of dimeric species can proceed through both olation and oxolation mechanisms. Oxolation is preferred for monomers that contain fewer water ligands, while olation becomes favorable for complexes with more water ligands.
Collapse
Affiliation(s)
- Rameswar Bhattacharjee
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Pere Miró
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| |
Collapse
|
2
|
Asthagiri DN, Paulaitis ME, Pratt LR. Thermodynamics of Hydration from the Perspective of the Molecular Quasichemical Theory of Solutions. J Phys Chem B 2021; 125:8294-8304. [PMID: 34313434 DOI: 10.1021/acs.jpcb.1c04182] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
The quasichemical organization of the potential distribution theorem, molecular quasichemical theory (QCT), enables practical calculations and also provides a conceptual framework for molecular hydration phenomena. QCT can be viewed from multiple perspectives: (a) as a way to regularize an ill-conditioned statistical thermodynamic problem; (b) as an introduction of and emphasis on the neighborship characteristics of a solute of interest; or (c) as a way to include accurate electronic structure descriptions of near-neighbor interactions in defensible statistical thermodynamics by clearly defining neighborship clusters. The theory has been applied to solutes of a wide range of chemical complexity, ranging from ions that interact with water with both long-ranged and chemically intricate short-ranged interactions, to solutes that interact with water solely through traditional van der Waals interations, and including water itself. The solutes range in variety from monatomic ions to chemically heterogeneous macromolecules. A notable feature of QCT is that, in applying the theory to this range of solutes, the theory itself provides guidance on the necessary approximations and simplifications that can facilitate the calculations. In this Perspective, we develop these ideas and document them with examples that reveal the insights that can be extracted using the QCT formulation.
Collapse
Affiliation(s)
- Dilipkumar N Asthagiri
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, United States
| | - Michael E Paulaitis
- Center for Nanomedicine, Johns Hopkins School of Medicine, Baltimore, Maryland 21231, United States
| | - Lawrence R Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
3
|
Absolute ion hydration free energy scale and the surface potential of water via quantum simulation. Proc Natl Acad Sci U S A 2020; 117:30151-30158. [PMID: 33203676 DOI: 10.1073/pnas.2017214117] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
With a goal of determining an absolute free energy scale for ion hydration, quasi-chemical theory and ab initio quantum mechanical simulations are employed to obtain an accurate value for the bulk hydration free energy of the Na+ ion. The free energy is partitioned into three parts: 1) the inner-shell or chemical contribution that includes direct interactions of the ion with nearby waters, 2) the packing free energy that is the work to produce a cavity of size λ in water, and 3) the long-range contribution that involves all interactions outside the inner shell. The interfacial potential contribution to the free energy resides in the long-range term. By averaging cation and anion data for that contribution, cumulant terms of all odd orders in the electrostatic potential are removed. The computed total is then the bulk hydration free energy. Comparison with the experimentally derived real hydration free energy produces an effective surface potential of water in the range -0.4 to -0.5 V. The result is consistent with a variety of experiments concerning acid-base chemistry, ion distributions near hydrophobic interfaces, and electric fields near the surface of water droplets.
Collapse
|
4
|
Muralidharan A, Pratt LR, Chaudhari MI, Rempe SB. Quasi-Chemical Theory with Cluster Sampling from Ab Initio Molecular Dynamics: Fluoride (F -) Anion Hydration. J Phys Chem A 2018; 122:9806-9812. [PMID: 30475612 DOI: 10.1021/acs.jpca.8b08474] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Accurate predictions of the hydration free energy for anions typically has been more challenging than that for cations. Hydrogen bond donation to the anion in hydrated clusters such as F(H2O) n - can lead to delicate structures. Consequently, the energy landscape contains many local minima, even for small clusters, and these minima present a challenge for computational optimization. Utilization of cluster experimental results for the free energies of gas-phase clusters shows that even though anharmonic effects are interesting they need not be of troublesome magnitudes for careful applications of quasi-chemical theory to ion hydration. Energy-optimized cluster structures for anions can leave the central ion highly exposed, and application of implicit solvation models to these structures can incur more serious errors than those for metal cations. Utilizing cluster structures sampled from ab initio molecular dynamics simulations substantially fixes those issues.
Collapse
Affiliation(s)
- A Muralidharan
- Department of Chemical and Biomolecular Engineering , Tulane University , New Orleans , Louisiana 70118 , United States
| | - L R Pratt
- Department of Chemical and Biomolecular Engineering , Tulane University , New Orleans , Louisiana 70118 , United States
| | - M I Chaudhari
- Center for Biological and Engineering Sciences , Sandia National Laboratories , Albuquerque , New Mexico 87185 , United States
| | - S B Rempe
- Center for Biological and Engineering Sciences , Sandia National Laboratories , Albuquerque , New Mexico 87185 , United States
| |
Collapse
|
5
|
König G, Pickard FC, Huang J, Thiel W, MacKerell AD, Brooks BR, York DM. A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules 2018; 23:E2695. [PMID: 30347691 PMCID: PMC6222909 DOI: 10.3390/molecules23102695] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 12/01/2022] Open
Abstract
Maintaining a proper balance between specific intermolecular interactions and non-specific solvent interactions is of critical importance in molecular simulations, especially when predicting binding affinities or reaction rates in the condensed phase. The most rigorous metric for characterizing solvent affinity are solvation free energies, which correspond to a transfer from the gas phase into solution. Due to the drastic change of the electrostatic environment during this process, it is also a stringent test of polarization response in the model. Here, we employ both the CHARMM fixed charge and polarizable force fields to predict hydration free energies of twelve simple solutes. The resulting classical ensembles are then reweighted to obtain QM/MM hydration free energies using a variety of QM methods, including MP2, Hartree⁻Fock, density functional methods (BLYP, B3LYP, M06-2X) and semi-empirical methods (OM2 and AM1 ). Our simulations test the compatibility of quantum-mechanical methods with molecular-mechanical water models and solute Lennard⁻Jones parameters. In all cases, the resulting QM/MM hydration free energies were inferior to purely classical results, with the QM/MM Drude force field predictions being only marginally better than the QM/MM fixed charge results. In addition, the QM/MM results for different quantum methods are highly divergent, with almost inverted trends for polarizable and fixed charge water models. While this does not necessarily imply deficiencies in the QM models themselves, it underscores the need to develop consistent and balanced QM/MM interactions. Both the QM and the MM component of a QM/MM simulation have to match, in order to avoid artifacts due to biased solute⁻solvent interactions. Finally, we discuss strategies to improve the convergence and efficiency of multi-scale free energy simulations by automatically adapting the molecular-mechanics force field to the target quantum method.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
- School of Life Sciences, Westlake University, 18 Shilongshan Street, Hangzhou 310024, China.
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander D MacKerell
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
| |
Collapse
|
6
|
Bansal A, Chapman WG, Asthagiri D. Quasichemical theory and the description of associating fluids relative to a reference: Multiple bonding of a single site solute. J Chem Phys 2017; 147:124505. [DOI: 10.1063/1.4997663] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Artee Bansal
- Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005-1827, USA
| | - Walter G. Chapman
- Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005-1827, USA
| | - D. Asthagiri
- Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005-1827, USA
| |
Collapse
|
7
|
Duignan TT, Baer MD, Schenter GK, Mundy CJ. Real single ion solvation free energies with quantum mechanical simulation. Chem Sci 2017; 8:6131-6140. [PMID: 28989643 PMCID: PMC5625628 DOI: 10.1039/c7sc02138k] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Accepted: 05/26/2017] [Indexed: 01/11/2023] Open
Abstract
Single ion solvation free energies are one of the most important properties of electrolyte solutions and yet there is ongoing debate about what these values are. Only the values for neutral ion pairs are known. Here, we use DFT interaction potentials with molecular dynamics simulation (DFT-MD) combined with a modified version of the quasi-chemical theory (QCT) to calculate these energies for the lithium and fluoride ions. A method to correct for the error in the DFT functional is developed and very good agreement with the experimental value for the lithium fluoride pair is obtained. Moreover, this method partitions the energies into physically intuitive terms such as surface potential, cavity and charging energies which are amenable to descriptions with reduced models. Our research suggests that lithium's solvation free energy is dominated by the free energetics of a charged hard sphere, whereas fluoride exhibits significant quantum mechanical behavior that cannot be simply described with a reduced model.
Collapse
Affiliation(s)
- Timothy T Duignan
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| | - Marcel D Baer
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| | - Gregory K Schenter
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| | - Christopher J Mundy
- Physical Science Division , Pacific Northwest National Laboratory , P.O. Box 999 , Richland , Washington 99352 , USA . ; Tel: +1 509 3756940
| |
Collapse
|
8
|
Asthagiri D, Valiya Parambathu A, Ballal D, Chapman WG. Electrostatic and induction effects in the solubility of water in alkanes. J Chem Phys 2017; 147:074506. [DOI: 10.1063/1.4997916] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- D. Asthagiri
- Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005-1827, USA
| | | | - Deepti Ballal
- Materials Science and Engineering, Ames Laboratory, Ames, Iowa 50011, USA
| | - Walter G. Chapman
- Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005-1827, USA
| |
Collapse
|
9
|
Brandt EG, Agosta L, Lyubartsev AP. Reactive wetting properties of TiO2 nanoparticles predicted by ab initio molecular dynamics simulations. NANOSCALE 2016; 8:13385-13398. [PMID: 27341183 DOI: 10.1039/c6nr02791a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Small-sized wet TiO2 nanoparticles have been investigated by ab initio molecular dynamics simulations. Chemical and physical adsorption of water on the TiO2-water interface was studied as a function of water content, ranging from dry nanoparticles to wet nanoparticles with monolayer coverage of water. The surface reactivity was shown to be a concave function of water content and driven by surface defects. The local coordination number at the defect was identified as the key factor to decide whether water adsorption proceeds through dissociation or physisorption on the surface. A consistent picture of TiO2 nanoparticle wetting at the microscopic level emerges, which corroborates existing experimental data and gives further insight into the molecular mechanisms behind nanoparticle wetting. These calculations will facilitate the engineering of metal oxide nanoparticles with a controlled catalytic water activity.
Collapse
Affiliation(s)
- Erik G Brandt
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden.
| | | | | |
Collapse
|
10
|
Shevkunov SV. Hydration of Cl– ion in a planar nanopore with hydrophilic walls. 1. Molecular structure. COLLOID JOURNAL 2016. [DOI: 10.1134/s1061933x15060186] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
11
|
Brorsen KR, Willow SY, Xantheas SS, Gordon MS. The Melting Temperature of Liquid Water with the Effective Fragment Potential. J Phys Chem Lett 2015; 6:3555-3559. [PMID: 26722723 DOI: 10.1021/acs.jpclett.5b01702] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The direct simulation of the solid-liquid water interface with the effective fragment potential (EFP) via the constant enthalpy and pressure (NPH) ensemble was used to estimate the melting temperature (T(m)) of ice-I(h). Initial configurations and velocities, taken from equilibrated constant pressure and temperature (NPT) simulations at P = 1 atm and T = 305 K, 325 K and 399 K, respectively, yielded corresponding T(m) values of 378 ± 16 K, 382 ± 14 K and 384 ± 15 K. These estimates are consistently higher than experiment, albeit to the same degree as previously reported estimates using density functional theory (DFT)-based Born-Oppenheimer simulations with the Becke-Lee-Yang-Parr functional plus dispersion corrections (BLYP-D).
Collapse
Affiliation(s)
- Kurt R Brorsen
- Ames Laboratory, US-DOE and Department of Chemistry, Iowa State University , Ames, Iowa 50011, United States
| | - Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Sotiris S Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory , 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| | - Mark S Gordon
- Ames Laboratory, US-DOE and Department of Chemistry, Iowa State University , Ames, Iowa 50011, United States
| |
Collapse
|
12
|
Iftimie R, Tremblay MH, Thomas V, Hétu S, de Lasalle F, Rivard U. Moderately strong phenols dissociate by forming an ion-pair kinetic intermediate. J Phys Chem A 2013; 117:13976-87. [PMID: 24299203 DOI: 10.1021/jp410858d] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We show computational evidence that ground-state moderately strong hydroxyarenes (Ar-OH, pKa ∼ 0) dissociate by forming an ion-pair intermediate that lives for 3-5 ps. The concentration of this intermediate is approximately 2 times smaller than that of the un-ionized acid at pH ∼ 0.6 and is characterized by average C-O bond lengths (1.30 Å) that are intermediate between those of un-ionized (1.29 Å) and fully dissociated (1.34 Å) species. During the lifetime of the ion-pair intermediate the excess proton fluctuates between the oxygen atom of the phenolic moiety and those of water molecules in the first and second solvation shells on a subpicosecond time scale (∼100-300 fs).
Collapse
Affiliation(s)
- Radu Iftimie
- Département de Chimie, Université de Montréal , CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | | | | | | | | | | |
Collapse
|
13
|
Hong QJ, van de Walle A. Solid-liquid coexistence in small systems: A statistical method to calculate melting temperatures. J Chem Phys 2013; 139:094114. [DOI: 10.1063/1.4819792] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
14
|
You X, Chaudhari MI, Pratt LR, Pesika N, Aritakula KM, Rick SW. Interfaces of propylene carbonate. J Chem Phys 2013; 138:114708. [DOI: 10.1063/1.4794792] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Xinli You
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, USA.
| | | | | | | | | | | |
Collapse
|
15
|
Avendaño C, Lafitte T, Adjiman CS, Galindo A, Müller EA, Jackson G. SAFT-γ force field for the simulation of molecular fluids: 2. Coarse-grained models of greenhouse gases, refrigerants, and long alkanes. J Phys Chem B 2013; 117:2717-33. [PMID: 23311931 DOI: 10.1021/jp306442b] [Citation(s) in RCA: 110] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In the first paper of this series [C. Avendaño, T. Lafitte, A. Galindo, C. S. Adjiman, G. Jackson, and E. A. Müller, J. Phys. Chem. B2011, 115, 11154] we introduced the SAFT-γ force field for molecular simulation of fluids. In our approach, a molecular-based equation of state (EoS) is used to obtain coarse-grained (CG) intermolecular potentials that can then be employed in molecular simulation over a wide range of thermodynamic conditions of the fluid. The macroscopic experimental data for the vapor-liquid equilibria (saturated liquid density and vapor pressure) of a given system are represented with the SAFT-VR Mie EoS and used to estimate effective intermolecular parameters that provide a good description of the thermodynamic properties by exploring a wide parameter space for models based on the Mie (generalized Lennard-Jones) potential. This methodology was first used to develop a simple single-segment CG Mie model of carbon dioxide (CO2) which allows for a reliable representation of the fluid-phase equilibria (for which the model was parametrized), as well as an accurate prediction of other properties such as the enthalpy of vaporization, interfacial tension, supercritical density, and second-derivative thermodynamic properties (thermal expansivity, isothermal compressibility, heat capacity, Joule-Thomson coefficient, and speed of sound). In our current paper, the methodology is further applied and extended to develop effective SAFT-γ CG Mie force fields for some important greenhouse gases including carbon tetrafluoride (CF4) and sulfur hexafluoride (SF6), modeled as simple spherical molecules, and for long linear alkanes including n-decane (n-C10H22) and n-eicosane (n-C20H42), modeled as homonuclear chains of spherical Mie segments. We also apply the SAFT-γ methodology to obtain a CG homonuclear two-segment Mie intermolecular potential for the more challenging polar and asymmetric compound 2,3,3,3-tetrafluoro-1-propene (HFO-1234yf), a novel replacement refrigerant with promising properties. The description of the fluid-phase behavior and the prediction of the other thermophysical properties obtained by molecular simulation using our SAFT-γ CG Mie force fields are found to be of comparable quality (and sometimes superior) to that obtained using the more sophisticated all-atom (AA) and united-atom (UA) models commonly employed in the field. We should emphasize that though the focus of our current work is on simple homonuclear models, the SAFT-γ methodology is based on a group contribution methodology which is naturally suited to the development of more sophisticated heteronuclear models.
Collapse
Affiliation(s)
- Carlos Avendaño
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | | | | | | | | | | |
Collapse
|
16
|
Pascal TA, Schärf D, Jung Y, Kühne TD. On the absolute thermodynamics of water from computer simulations: A comparison of first-principles molecular dynamics, reactive and empirical force fields. J Chem Phys 2012; 137:244507. [DOI: 10.1063/1.4771974] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
|
17
|
Weber V, Asthagiri D. Regularizing Binding Energy Distributions and the Hydration Free Energy of Protein Cytochrome C from All-Atom Simulations. J Chem Theory Comput 2012; 8:3409-15. [DOI: 10.1021/ct300505b] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
| | - D. Asthagiri
- Department of Chemical and Biomolecular
Engineering, Johns Hopkins University, Baltimore, Maryland 21218,
United States
| |
Collapse
|
18
|
Leung K, Criscenti LJ. Predicting the acidity constant of a goethite hydroxyl group from first principles. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2012; 24:124105. [PMID: 22395040 DOI: 10.1088/0953-8984/24/12/124105] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Accurate predictions of the acid-base behavior of hydroxyl groups at mineral surfaces are critical for understanding the trapping of toxic and radioactive ions in soil samples. In this work, we apply ab initio molecular dynamics (AIMD) simulations and potential-of-mean-force techniques to calculate the pK(a) of a doubly protonated oxygen atom bonded to a single Fe atom (Fe(I)OH(2)) on the goethite (101) surface. Using formic acid as a reference system, pK(a) = 7.0 is predicted, suggesting that isolated, positively charged groups of this type are marginally stable at neutral pH. Similarities and differences between AIMD and the more empirical multi-site complexation methodology are highlighted, particularly with respect to the treatment of hydrogen bonding with water and proton sharing among surface hydroxyl groups. We also highlight the importance of an electronic structure method that can accurately predict transition metal ion properties for goethite pK(a) calculations.
Collapse
Affiliation(s)
- Kevin Leung
- Sandia National Laboratories, MS 1415 and 0754, Albuquerque, NM 87185, USA.
| | | |
Collapse
|
19
|
|
20
|
Weber V, Merchant S, Asthagiri D. Communication: Regularizing binding energy distributions and thermodynamics of hydration: Theory and application to water modeled with classical and ab initio simulations. J Chem Phys 2011; 135:181101. [DOI: 10.1063/1.3660205] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Valéry Weber
- IBM Research Division, Zurich Research Laboratory, 8803 Ruschlikon, Switzerland
| | - Safir Merchant
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - D. Asthagiri
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
21
|
Zhang C, Spanu L, Galli G. Entropy of Liquid Water from Ab Initio Molecular Dynamics. J Phys Chem B 2011; 115:14190-5. [DOI: 10.1021/jp204981y] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Cui Zhang
- Department of Chemistry and ‡Department of Physics, University of California, Davis, California 95616, United States
| | - Leonardo Spanu
- Department of Chemistry and ‡Department of Physics, University of California, Davis, California 95616, United States
| | - Giulia Galli
- Department of Chemistry and ‡Department of Physics, University of California, Davis, California 95616, United States
| |
Collapse
|
22
|
Baer MD, Mundy CJ, McGrath MJ, Kuo IFW, Siepmann JI, Tobias DJ. Re-examining the properties of the aqueous vapor–liquid interface using dispersion corrected density functional theory. J Chem Phys 2011; 135:124712. [DOI: 10.1063/1.3633239] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
McGrath MJ, Kuo IFW, Ghogomu JN, Mundy CJ, Siepmann JI. Vapor–Liquid Coexistence Curves for Methanol and Methane Using Dispersion-Corrected Density Functional Theory. J Phys Chem B 2011; 115:11688-92. [DOI: 10.1021/jp205072v] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matthew J. McGrath
- Department of Physics, University of Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - I.-F. Will Kuo
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, California 94550, United States
| | - Julius N. Ghogomu
- Department of Chemistry, University of Dschang, B.P. 67, Dschang, Cameroon
- Departments of Chemistry and Chemical Engineering and Materials Science and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Christopher J. Mundy
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Washington 99352, United States
| | - J. Ilja Siepmann
- Departments of Chemistry and Chemical Engineering and Materials Science and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
24
|
Abstract
Monovalent ion hydration entropies are analyzed via energetic partitioning of the potential distribution theorem free energy. Extensive molecular dynamics simulations and free energy calculations are performed over a range of temperatures to determine the electrostatic and van der Waals components of the entropy. The far-field electrostatic contribution is negative and small in magnitude, and it does not vary significantly as a function of ion size, consistent with dielectric models. The local electrostatic contribution, however, varies widely as a function of ion size; the sign yields a direct indication of the kosmotropic (strongly hydrated) or chaotropic (weakly hydrated) nature of the ion hydration. The results provide a thermodynamic signature for specific ion effects in hydration and are consistent with experiments that suggest minimal perturbations of water structure outside the first hydration shell. The hydration entropies are also examined in relation to the corresponding entropies for the isoelectronic rare gas pairs; an inverse correlation is observed, as expected from thermodynamic hydration data.
Collapse
Affiliation(s)
- Thomas L Beck
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221-0172, United States.
| |
Collapse
|
25
|
Merchant S, Shah JK, Asthagiri D. Water coordination structures and the excess free energy of the liquid. J Chem Phys 2011; 134:124514. [DOI: 10.1063/1.3572058] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
|
26
|
Maurer P, Thomas V, Iftimie R. A computational study of ultrafast acid dissociation and acid–base neutralization reactions. II. The relationship between the coordination state of solvent molecules and concerted versus sequential acid dissociation. J Chem Phys 2011; 134:094505. [DOI: 10.1063/1.3554654] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
|
27
|
McGrath MJ, Kuo IFW, Siepmann JI. Liquid structures of water, methanol, and hydrogen fluoride at ambient conditions from first principles molecular dynamics simulations with a dispersion corrected density functional. Phys Chem Chem Phys 2011; 13:19943-50. [DOI: 10.1039/c1cp21890e] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
Watkins M, Berkowitz ML, Shluger AL. Role of water in atomic resolution AFM in solutions. Phys Chem Chem Phys 2011; 13:12584-94. [DOI: 10.1039/c1cp21021a] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|