1
|
Yan S, Wang B, Lin H. Reshaping the QM Region On-the-Fly: Adaptive-Shape QM/MM Dynamic Simulations of a Hydrated Proton in Bulk Water. J Chem Theory Comput 2024; 20:3462-3472. [PMID: 38671391 DOI: 10.1021/acs.jctc.4c00164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2024]
Abstract
Adaptive quantum mechanics/molecular mechanics (QM/MM) reclassifies on-the-fly a molecule or molecular fragment as QM or MM during dynamics simulations without abrupt changes in the energy or forces. Notably, the permuted adaptive-partitioning (PAP) algorithms have been applied to simulate a hydrated proton, with a mobile QM zone anchored at a pseudoatom called a proton indicator. The position of the proton indicator approximates the location of the delocalized excess proton, yielding a smooth trajectory of the proton diffusing via the Grotthuss mechanism in aqueous solutions. The mobile QM zone, which has been taken to be a sphere with a preset radius, follows the proton wherever it goes. Although the simulations are successful, the use of a spherical QM zone has one disadvantage: A large preset radius must be utilized to minimize the chance of missing water molecules that are important to proton translocation. A large radius leads to a large QM zone, which is computationally expensive. In this work, we report a new way to set up the QM zone, where one includes only the water molecules important to proton transfer. The importance of a given water molecule is quantified by its "weight" that depends on its relation to the reaction path of proton transfer. The weight varies smoothly, ensuring that a water molecule gradually appears in or disappears from the QM zone without abrupt changes, as required by the PAP method. Consequently, the shape of the QM zone evolves on-the-fly, keeping the QM zone as small as possible and as large as necessary. Test simulations demonstrate that the new algorithm significantly improves the computation efficiency while maintaining the proper descriptions of proton transfer in bulk water.
Collapse
Affiliation(s)
- Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Hai Lin
- Department of Chemistry, CB 194, University of Colorado Denver, Denver, P.O. Box 173364, Colorado 80217, United States
| |
Collapse
|
2
|
Stoppelman JP, McDaniel JG. Proton Transport in [BMIM+][BF4–]/Water Mixtures Near the Percolation Threshold. J Phys Chem B 2020; 124:5957-5970. [DOI: 10.1021/acs.jpcb.0c02487] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- John P. Stoppelman
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia30332-0400, United States
| | - Jesse G. McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia30332-0400, United States
| |
Collapse
|
3
|
Schran C, Behler J, Marx D. Automated Fitting of Neural Network Potentials at Coupled Cluster Accuracy: Protonated Water Clusters as Testing Ground. J Chem Theory Comput 2019; 16:88-99. [DOI: 10.1021/acs.jctc.9b00805] [Citation(s) in RCA: 56] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr−Universität Bochum, 44780 Bochum, Germany
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr−Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
4
|
Hofer TS, Kilchert FM, Tanjung BA. An effective partial charge model for bulk and surface properties of cubic ZrO 2, Y 2O 3 and yttrium-stabilised zirconia. Phys Chem Chem Phys 2019; 21:25635-25648. [PMID: 31720638 DOI: 10.1039/c9cp04307a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
In this work a newly parametrised Coulomb plus Buckingham potential formulation for cubic ZrO2, Y2O3 and yttrium-stabilised zirconia (YSZ) is presented. The density and pair distributions obtained for neat ZrO2 and Y2O3 under ambient conditions are in excellent agreement with experimental data, while the vibrational power spectra are highly similar compared to those obtained via ab initio molecular dynamics simulations at the PBEsol level. In addition, it is shown that the use of effective partial charges has several advantages compared to interaction potentials employing the oxidation states in the evaluation of the coulombic interactions: (i) the diffusion coefficient and the associated activation energy of oxygen ions evaluated for YSZn (n = 4 to 12) display the best agreement with experimental data; (ii) no unphysical reorganisation of the interface and the bulk are observed in simulations of the (110) and (111) surfaces of cubic ZrO2 and Y2O3, while due to the strong coulombic contributions in the case of the tested full-charge models a pronounced restructuring of the interface and the bulk is observed in the ZrO2 case, and (iii) the use of effective partial charges ensures compatibility with existing solvent models and force-fields for the treatment of molecular compounds.
Collapse
Affiliation(s)
- Thomas S Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria.
| | - Franziska M Kilchert
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria.
| | - Bagas A Tanjung
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria.
| |
Collapse
|
5
|
Bresnahan CG, David R, Milet A, Kumar R. Ion Pairing in HCl-Water Clusters: From Electronic Structure Investigations to Multiconfigurational Force-Field Development. J Phys Chem A 2019; 123:9371-9381. [PMID: 31589444 DOI: 10.1021/acs.jpca.9b07775] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the bulk, condensed-phase HCl exists as a dissociated Cl- ion and a proton that is delocalized over solvating water molecules. However, in the gas phase, HCl is covalent, and even on the introduction of hydrating water molecules, the HCl covalent state dominates small clusters and is relevant at larger clusters including 21 water molecules. Electronic structure calculations (at the MP2 level) and ab initio metadynamics simulations (at the DFT level) have been carried out on HCl-(H2O)n clusters with n = 2-22 to investigate distinct solvation environments in clusters from covalent HCl structure, to contact ion pairs and solvent-separated ion pairs. The data were further used to train and validate a multiconfigurational force-field for HCl-water clusters that incorporates covalent HCl states into the MS-EVB3.2 formalism. Additionally, the many-body interaction of the Cl- ion with water and the excess proton was modeled by the introduction of two geometric three-body terms that incorporates the dominant many-body interaction in an efficient noniterative manner.
Collapse
Affiliation(s)
- Caitlin G Bresnahan
- Department of Chemistry , 232 Choppin Hall , Louisiana State University , Baton Rouge , Louisiana 70803 , United States
| | - Rolf David
- Department of Chemistry , 232 Choppin Hall , Louisiana State University , Baton Rouge , Louisiana 70803 , United States.,Univ. Grenoble Alpes, CNRS, DCM , 38000 Grenoble , France
| | - Anne Milet
- Univ. Grenoble Alpes, CNRS, DCM , 38000 Grenoble , France
| | - Revati Kumar
- Department of Chemistry , 232 Choppin Hall , Louisiana State University , Baton Rouge , Louisiana 70803 , United States
| |
Collapse
|
6
|
Berthias F, Feketeová L, Abdoul-Carime H, Calvo F, Farizon B, Farizon M, Märk TD. Maxwell-Boltzmann versus non-ergodic events in the velocity distribution of water molecules evaporated from protonated water nanodroplets. J Chem Phys 2018; 149:084308. [PMID: 30193492 DOI: 10.1063/1.5037281] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Measurement of velocity distributions of evaporated water monomers from small mass- and energy-selected protonated water clusters allows probing the extent of thermalization after excitation of these ultimately small nanodroplets. Electronic excitation of a molecule in the cluster is here induced by a single collision with an argon atom in the keV energy range. The measured velocity distributions of the departing neutral molecules exhibit bimodal shapes with a lower-velocity part consistent with a complete redistribution of the deposited energy in the entire cluster and a higher-velocity contribution corresponding to evaporation before complete energy redistribution. Statistical molecular dynamics calculations reproduce the bimodal shape of the velocity distributions by assuming an initial spreading of the excitation energy among all modes, thereby reproducing the lower velocity contribution of the distribution. By contrast, assuming the deposited energy to be initially localized among the modes of a single molecule leads to calculated distributions with two components whose shape is in accordance with the experimental results. The characteristics and the relative abundance of these two contributions in the velocity distributions obtained are presented and discussed as a function of the number of molecules (n = 2-10) in the ionized nanodroplet H+(H2O) n .
Collapse
Affiliation(s)
- F Berthias
- Institut de Physique Nucléaire de Lyon; CNRS/IN2P3, UMR5822, Université de Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne, France
| | - L Feketeová
- Institut de Physique Nucléaire de Lyon; CNRS/IN2P3, UMR5822, Université de Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne, France
| | - H Abdoul-Carime
- Institut de Physique Nucléaire de Lyon; CNRS/IN2P3, UMR5822, Université de Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne, France
| | - F Calvo
- LIPhy, Université Grenoble 1 and CNRS, UMR 5588, 140, Ave. de la Physique, 38402 Saint-Martin-d'Hères, France
| | - B Farizon
- Institut de Physique Nucléaire de Lyon; CNRS/IN2P3, UMR5822, Université de Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne, France
| | - M Farizon
- Institut de Physique Nucléaire de Lyon; CNRS/IN2P3, UMR5822, Université de Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne, France
| | - T D Märk
- Institut für Ionenphysik und Angewandte Physik, Leopold Franzens Universität, Technikerstrasse 25, A-6020 Innsbruck, Austria
| |
Collapse
|
7
|
Hofer TS, Wiedemair MJ. Towards a dissociative SPC-like water model II. The impact of Lennard-Jones and Buckingham non-coulombic forces. Phys Chem Chem Phys 2018; 20:28523-28534. [DOI: 10.1039/c8cp04957b] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dissociative water potential by Garofalini and coworkers has been re-formulated in the framework of the widely employed Lennard-Jones and Buckingham potentials, enhancing the transferability of the model to third party simulation programs.
Collapse
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- Center for Chemistry and Biomedicine
- University of Innsbruck
| | - Martin J. Wiedemair
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- Center for Chemistry and Biomedicine
- University of Innsbruck
| |
Collapse
|
8
|
Wick CD. A comparison of sodium and hydrogen halides at the air-water interface. J Chem Phys 2017; 147:161703. [PMID: 29096503 PMCID: PMC5648577 DOI: 10.1063/1.4984114] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2017] [Accepted: 05/12/2017] [Indexed: 11/14/2022] Open
Abstract
New molecular models, parameterized to ab initio calculations, were developed to describe HBr and HI at the air-water interface. These were used to compare how the air-water interface influenced dissociation of NaX and HX, with X being Cl, Br, or I, and also their propensity for the interface. The polarizable multistate empirical valence bond method, which explicitly describes proton sharing, was used to model HX. Results showed that the air-water interface suppressed HX dissociation from a contact ion pair to a solvent separated to a greater degree than NaX dissociation. Furthermore, HX had a greater propensity for the interface than NaX, which was a consequence of the hydronium ion having a greatest interfacial activity of all species studied. As a consequence of this, the average configuration of dissociated HX, while in both contact ion and solvent separated ion pairs near the air-water interface, is with the dissociated hydrogen oriented more towards the air than the X atom.
Collapse
Affiliation(s)
- Collin D Wick
- Louisiana Tech University, Ruston, Louisiana 71270, USA
| |
Collapse
|
9
|
Affiliation(s)
- Jesse G. McDaniel
- Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
10
|
Duster AW, Lin H. Restrained Proton Indicator in Combined Quantum-Mechanics/Molecular-Mechanics Dynamics Simulations of Proton Transfer through a Carbon Nanotube. J Phys Chem B 2017; 121:8585-8592. [PMID: 28820594 DOI: 10.1021/acs.jpcb.7b06657] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Recently, a collective variable "proton indicator" was purposed for tracking an excess proton solvated in bulk water in molecular dynamics simulations. In this work, we demonstrate the feasibility of utilizing the position of this proton indicator as a reaction coordinate to model an excess proton migrating through a hydrophobic carbon nanotube in combined quantum-mechanics/molecular-mechanics simulations. Our results indicate that applying a harmonic restraint to the proton indicator in the bulk solvent near the nanotube pore entrance leads to the recruitment of water molecules into the pore. This is consistent with an earlier study that employed a multistate empirical valence bond potential and a different representation (center of excess charge) of the proton. We attribute this water recruitment to the delocalized nature of the solvated proton, which prefers to be in high-dielectric bulk solvent. While water recruitment into the pore is considered an artifact in the present simulations (because of the artificially imposed restraint on the proton), if the proton were naturally restrained, it could assist in building water wires prior to proton transfer through the pore. The potential of mean force for a proton translocation through the water-filled pore was computed by umbrella sampling, where the bias potentials were applied to the proton indicator. The free energy curve and barrier heights agree reasonably with those in the literature. The results suggest that the proton indicator can be used as a reaction coordinate in simulations of proton transport in confined environments.
Collapse
Affiliation(s)
- Adam W Duster
- Chemistry Department, CB 194, University of Colorado Denver , Denver, Colorado 80217, United States
| | - Hai Lin
- Chemistry Department, CB 194, University of Colorado Denver , Denver, Colorado 80217, United States
| |
Collapse
|
11
|
Duster AW, Wang C, Garza CM, Miller DE, Lin H. Adaptive quantum/molecular mechanics: what have we learned, where are we, and where do we go from here? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1310] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Adam W. Duster
- Chemistry Department University of Colorado Denver Denver CO USA
| | - Chun‐Hung Wang
- Chemistry Department University of Colorado Denver Denver CO USA
| | | | | | - Hai Lin
- Chemistry Department University of Colorado Denver Denver CO USA
| |
Collapse
|
12
|
Hori Y, Ida T, Mizuno M. Potential energy construction in the diabatic picture for quantum mechanical rate constants of intermolecular proton transfer. Phys Chem Chem Phys 2017. [DOI: 10.1039/c7cp03024j] [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/21/2022]
Abstract
We propose a simple method for potential construction in the diabatic picture and the estimation of thermal rate constants for intermolecular proton transfer reactions using quantum dynamics simulations carried out on the constructed potentials.
Collapse
Affiliation(s)
- Yuta Hori
- Chemistry Course
- Division of Material Chemistry
- Graduate School of Natural Science and Technology
- Kanazawa University
- Kanazawa 920-1192
| | - Tomonori Ida
- Chemistry Course
- Division of Material Chemistry
- Graduate School of Natural Science and Technology
- Kanazawa University
- Kanazawa 920-1192
| | - Motohiro Mizuno
- Chemistry Course
- Division of Material Chemistry
- Graduate School of Natural Science and Technology
- Kanazawa University
- Kanazawa 920-1192
| |
Collapse
|
13
|
Wiedemair MJ, Hofer TS. Towards a dissociative SPC-like water model – probing the impact of intramolecular Coulombic contributions. Phys Chem Chem Phys 2017; 19:31910-31920. [DOI: 10.1039/c7cp06191a] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A modification of the dissociative Garofalini water model towards an SPC-like Coulombic formulation proved to enhance accuracy and transferability of this successful force field approach.
Collapse
Affiliation(s)
- Martin J. Wiedemair
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- 6020 Innsbruck
| | - Thomas S. Hofer
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- 6020 Innsbruck
| |
Collapse
|
14
|
Abdoul‐Carime H, Berthias F, Feketeová L, Marciante M, Calvo F, Forquet V, Chermette H, Farizon B, Farizon M, Märk TD. Velocity of a Molecule Evaporated from a Water Nanodroplet: Maxwell–Boltzmann Statistics versus Non‐Ergodic Events. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/ange.201505890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Hassan Abdoul‐Carime
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Francis Berthias
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Linda Feketeová
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Mathieu Marciante
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institute Lumière Matière, CNRS, 69622 Villeurbanne (France)
| | - Florent Calvo
- Université Grenoble 1, CNRS, LIPhy UMR 5588, 38041 Grenoble (France)
| | - Valérian Forquet
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut des Sciences Analytiques, CNRS UMR 5280, 43 boulevard 11 novembre, 69622 Villeurbanne Cedex (France)
| | - Henry Chermette
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut des Sciences Analytiques, CNRS UMR 5280, 43 boulevard 11 novembre, 69622 Villeurbanne Cedex (France)
| | - Bernadette Farizon
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Michel Farizon
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Tilmann D. Märk
- Institut für Ionenphysik und Angewandte Physik, Leopold Franzens Universität, 6020 Innsbruck (Austria)
| |
Collapse
|
15
|
Abdoul-Carime H, Berthias F, Feketeová L, Marciante M, Calvo F, Forquet V, Chermette H, Farizon B, Farizon M, Märk TD. Velocity of a Molecule Evaporated from a Water Nanodroplet: Maxwell-Boltzmann Statistics versus Non-Ergodic Events. Angew Chem Int Ed Engl 2015; 54:14685-9. [PMID: 26473406 PMCID: PMC4678507 DOI: 10.1002/anie.201505890] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2015] [Revised: 09/07/2015] [Indexed: 11/22/2022]
Abstract
The velocity of a molecule evaporated from a mass-selected protonated water nanodroplet is measured by velocity map imaging in combination with a recently developed mass spectrometry technique. The measured velocity distributions allow probing statistical energy redistribution in ultimately small water nanodroplets after ultrafast electronic excitation. As the droplet size increases, the velocity distribution rapidly approaches the behavior expected for macroscopic droplets. However, a distinct high-velocity contribution provides evidence of molecular evaporation before complete energy redistribution, corresponding to non-ergodic events.
Collapse
Affiliation(s)
- Hassan Abdoul-Carime
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Francis Berthias
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Linda Feketeová
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Mathieu Marciante
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institute Lumière Matière, CNRS, 69622 Villeurbanne (France)
| | - Florent Calvo
- Université Grenoble 1, CNRS, LIPhy UMR 5588, 38041 Grenoble (France)
| | - Valérian Forquet
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut des Sciences Analytiques, CNRS UMR 5280, 43 boulevard 11 novembre, 69622 Villeurbanne Cedex (France)
| | - Henry Chermette
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut des Sciences Analytiques, CNRS UMR 5280, 43 boulevard 11 novembre, 69622 Villeurbanne Cedex (France)
| | - Bernadette Farizon
- Université de Lyon, 69003 Lyon (France)
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France)
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France)
| | - Michel Farizon
- Université de Lyon, 69003 Lyon (France).
- Université Claude Bernard Lyon1, 69622 Villeurbanne (France).
- Institut de Physique Nucléaire de Lyon, CNRS/IN2P3 UMR 5822, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France).
| | - Tilmann D Märk
- Institut für Ionenphysik und Angewandte Physik, Leopold Franzens Universität, 6020 Innsbruck (Austria)
| |
Collapse
|
16
|
Biswas R, Tse YLS, Tokmakoff A, Voth GA. Role of Presolvation and Anharmonicity in Aqueous Phase Hydrated Proton Solvation and Transport. J Phys Chem B 2015; 120:1793-804. [PMID: 26575795 DOI: 10.1021/acs.jpcb.5b09466] [Citation(s) in RCA: 58] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Results from condensed phase ab initio molecular dynamics (AIMD) simulations suggest a proton transfer reaction is facilitated by "presolvation" in which the hydronium is transiently solvated by four water molecules, similar to the typical solvation structure of water, by accepting a weak hydrogen bond from the fourth water molecule. A new version 3.2 multistate empirical valence bond (MS-EVB 3.2) model for the hydrated excess proton incorporating this presolvation behavior is therefore developed. The classical MS-EVB simulations show similar structural properties as those of the previous model but with significantly improved diffusive behavior. The inclusion of nuclear quantum effects in the MS-EVB also provides an even better description of the proton diffusion rate. To quantify the influence of anharmonicity, a second model (aMS-EVB 3.2) is developed using the anharmonic aSPC/Fw water model, which provides similar structural properties but improved spectroscopic responses at high frequencies.
Collapse
Affiliation(s)
- Rajib Biswas
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| | - Ying-Lung Steve Tse
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| |
Collapse
|
17
|
Wiedemair MJ, Hitzenberger M, Hofer TS. Tuning the reactivity of a dissociative force field: proton transfer properties of aqueous H3O+ and their dependence on the three-body interaction. Phys Chem Chem Phys 2015; 17:10934-43. [DOI: 10.1039/c4cp05607h] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Selective adjustment of the three-body interaction of a dissociative water potential results in a significant improvement in the description of proton transport properties.
Collapse
Affiliation(s)
- Martin J. Wiedemair
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- Innrain 80-82
| | - Manuel Hitzenberger
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- Innrain 80-82
| | - Thomas S. Hofer
- Theoretical Chemistry Division
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- Innrain 80-82
| |
Collapse
|
18
|
Hofer TS. Probing Proton Transfer Reactions in Molecular Dynamics—A Crucial Prerequisite for QM/MM Simulations Using Dissociative Models. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2015. [DOI: 10.1007/978-3-319-21626-3_4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
|
19
|
Kumar R, Knight C, Wick CD, Chen B. Bringing Reactivity to the Aggregation-Volume-Bias Monte Carlo Based Simulation Framework: Water Nucleation Induced by a Reactive Proton. J Phys Chem B 2014; 119:9068-75. [DOI: 10.1021/jp508749x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Revati Kumar
- Department
of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
- Center for Computation and Technology, Louisiana State University, Baton
Rouge, Louisiana 70803, United States
| | - Chris Knight
- Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Collin D. Wick
- College
of Engineering and Science, Louisiana Tech University, Ruston, Louisiana 71272, United States
| | - Bin Chen
- Department
of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
20
|
Canaval LR, Lutz OMD, Weiss AKH, Huck CW, Hofer TS. A Dissociative Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulation and Infrared Experiments Reveal Characteristics of the Strongly Hydrolytic Arsenic(III). Inorg Chem 2014; 53:11861-70. [DOI: 10.1021/ic4031156] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Lorenz R. Canaval
- Theoretical Chemistry Division, Institute of General,
Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain
80-82, A-6020 Innsbruck, Austria
| | - Oliver M. D. Lutz
- Institute for Analytical Chemistry and
Radiochemistry, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Alexander K. H. Weiss
- Theoretical Chemistry Division, Institute of General,
Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain
80-82, A-6020 Innsbruck, Austria
| | - Christian W. Huck
- Institute for Analytical Chemistry and
Radiochemistry, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General,
Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain
80-82, A-6020 Innsbruck, Austria
| |
Collapse
|
21
|
Wiebe J, Spohr E. Double layer effects in a model of proton discharge on charged electrodes. BEILSTEIN JOURNAL OF NANOTECHNOLOGY 2014; 5:973-982. [PMID: 25161833 PMCID: PMC4143088 DOI: 10.3762/bjnano.5.111] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2014] [Accepted: 06/10/2014] [Indexed: 06/03/2023]
Abstract
We report first results on double layer effects on proton discharge reactions from aqueous solutions to charged platinum electrodes. We have extended a recently developed combined proton transfer/proton discharge model on the basis of empirical valence bond theory to include specifically adsorbed sodium cations and chloride anions. For each of four studied systems 800-1000 trajectories of a discharging proton were integrated by molecular dynamics simulations until discharge occurred. The results show significant influences of ion presence on the average behavior of protons prior to the discharge event. Rationalization of the observed behavior cannot be based solely on the electrochemical potential (or surface charge) but needs to resort to the molecular details of the double layer structure.
Collapse
Affiliation(s)
- Johannes Wiebe
- Faculty of Chemistry, Duisburg-Essen University, D-45117 Essen, Germany
| | - Eckhard Spohr
- Faculty of Chemistry and Center for Computational Sciences and Simulation (CCSS), Universität Duisburg-Essen, D-45117 Essen, Germany
| |
Collapse
|
22
|
Hofer TS. Perspectives for hybrid ab initio/molecular mechanical simulations of solutions: from complex chemistry to proton-transfer reactions and interfaces. PURE APPL CHEM 2014. [DOI: 10.1515/pac-2014-5019] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
As a consequence of the ongoing development of enhanced computational resources, theoretical chemistry has become an increasingly valuable field for the investigation of a variety of chemical systems. Simulations employing a hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) technique have been shown to be a particularly promising approach, whenever ultrafast (i.e., picosecond) dynamical properties are to be studied, which are in many cases difficult to access via experimental techniques. Details of the quantum mechanical charge field (QMCF) ansatz, an advanced QM/MM protocol, are discussed and simulation results for various systems ranging from simple ionic hydrates to solvated organic molecules and coordination complexes in solution are presented. A particularly challenging application is the description of proton-transfer reactions in chemical simulations, which is a prerequisite to study acidified and basic systems. The methodical requirements for a combination of the QMCF methodology with a dissociative potential model for the description of the solvent are discussed. Furthermore, the possible extension of QM/MM approaches to solid/liquid interfaces is outlined.
Collapse
|
23
|
Marciante M, Calvo F. Modelling infrared action spectra of protonated water clusters. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.844804] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- M. Marciante
- ILM, CNRS UMR 5306, Université Lyon 1, Bâtiment Kastler, 43 Bd du 11 Novembre 1918, F69622, Villeurbanne Cedex, France
| | - F. Calvo
- ILM, CNRS UMR 5306, Université Lyon 1, Bâtiment Kastler, 43 Bd du 11 Novembre 1918, F69622, Villeurbanne Cedex, France
| |
Collapse
|
24
|
Wick CD. HCl Accommodation, Dissociation, and Propensity for the Surface of Water. J Phys Chem A 2013; 117:12459-67. [DOI: 10.1021/jp4084212] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Collin D. Wick
- Louisiana Tech University, P.O. Box 10348, Ruston, Louisiana 71270, United States
| |
Collapse
|
25
|
Dreyer J, Zhang C, Ippoliti E, Carloni P. Role of the Membrane Dipole Potential for Proton Transport in Gramicidin A Embedded in a DMPC Bilayer. J Chem Theory Comput 2013; 9:3826-31. [PMID: 26584128 DOI: 10.1021/ct400374n] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The membrane potential at the water/phospholipid interfaces may play a key role for proton conduction of gramicidin A (gA). Here we address this issue by Density Functional Theory-based molecular dynamics and metadynamics simulations. The calculations, performed on gA embedded in a solvated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) model membrane environment (about 2,000 atoms), indicate that (i) the membrane dipole potential rises at the channel mouth by ∼0.4 V. A similar value has been measured for gA embedded in a DMPC monolayer; (ii) the calculated free energy barrier is located at the channel entrance, consistent with experiments comparing gA proton conduction in different bilayers. The electronic structures of the proton ligands (water molecules and peptide units) are similar to those in the bulk solvent. Based on these results, we suggest an important role of the membrane dipole potential for the free energy barrier of proton permeation of gA. This may provide a rationale for the large increase in the rate of proton conduction under application of a transmembrane voltage, as observed experimentally. Our calculations might suggest also a role for proton desolvation for the permeation process. This role has already emerged from EVB calculations on gA embedded in a model membrane.
Collapse
Affiliation(s)
- Jens Dreyer
- Computational Biophysics, German Research School for Simulation Sciences, Joint venture of RWTH Aachen University and Forschungszentrum Jülich , Germany, D-52425 Jülich, Germany.,IAS-5, Computational Biomedicine, Institute for Advanced Simulation , Forschungszentrum Jülich, D-52425 Jülich, Germany
| | - Chao Zhang
- Computational Biophysics, German Research School for Simulation Sciences, Joint venture of RWTH Aachen University and Forschungszentrum Jülich , Germany, D-52425 Jülich, Germany.,IAS-5, Computational Biomedicine, Institute for Advanced Simulation , Forschungszentrum Jülich, D-52425 Jülich, Germany
| | - Emiliano Ippoliti
- Computational Biophysics, German Research School for Simulation Sciences, Joint venture of RWTH Aachen University and Forschungszentrum Jülich , Germany, D-52425 Jülich, Germany.,IAS-5, Computational Biomedicine, Institute for Advanced Simulation , Forschungszentrum Jülich, D-52425 Jülich, Germany
| | - Paolo Carloni
- Computational Biophysics, German Research School for Simulation Sciences, Joint venture of RWTH Aachen University and Forschungszentrum Jülich , Germany, D-52425 Jülich, Germany.,IAS-5, Computational Biomedicine, Institute for Advanced Simulation , Forschungszentrum Jülich, D-52425 Jülich, Germany
| |
Collapse
|
26
|
Probing the temperature dependence of proton transfer to charged platinum electrodes by reactive molecular dynamics trajectory studies. Electrochim Acta 2013. [DOI: 10.1016/j.electacta.2013.01.146] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
27
|
Asthana A, Wheeler DR. A polarizable reactive force field for water to enable molecular dynamics simulations of proton transport. J Chem Phys 2013; 138:174502. [DOI: 10.1063/1.4798457] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
|
28
|
Lockwood GK, Garofalini SH. Lifetimes of excess protons in water using a dissociative water potential. J Phys Chem B 2013; 117:4089-97. [PMID: 23565831 DOI: 10.1021/jp310300x] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Molecular dynamics simulations using a dissociative water potential were applied to study transport of excess protons in water and determine the applicability of this potential to describe such behavior. While originally developed for gas-phase molecules and bulk liquid water, the potential is transferrable to nanoconfinement and interface scenarios. Applied here, it shows proton behavior consistent with ab initio calculations and empirical models specifically designed to describe proton transport. Both Eigen and Zundel complexes are observed in the simulations showing the Eigen-Zundel-Eigen-type mechanism. In addition to reproducing the short-time rattling of the excess proton between the two oxygens of Zundel complexes, a picosecond-scale lifetime was also found. These longer-lived H3O(+) ions are caused by the rapid conversion of the local solvation structure around the transferring proton from a Zundel-like form to an Eigen-like form following the transfer, effectively severing the path along which the proton can rattle. The migration of H(+) over long times (>100 ps) deviates from the conventional short-time multiexponentially decaying lifetime autocorrelation model and follows the t(-3/2) power-law behavior. The potential function employed here matches many of the features of proton transport observed in ab initio molecular dynamics simulations as well as the highly developed empirical valence bond models, yet is computationally very efficient, enabling longer time and larger systems to be studied.
Collapse
Affiliation(s)
- Glenn K Lockwood
- Interfacial Molecular Science Laboratory, Department of Materials Science and Engineering, Rutgers University, 607 Taylor Rd., Piscataway, New Jersey 08855, USA
| | | |
Collapse
|
29
|
Abstract
Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data directly from condensed phase ab initio simulations to develop reactive molecular dynamics models that do not require predefined empirical functions. Instead, the interactions used in the reactive model are expressed as linear combinations of interpolating functions that are optimized by using a linear least-squares algorithm. One notable benefit of the procedure outlined here is the capability to minimize the number of parameters requiring nonlinear optimization. The method presented can be generally applied to multiscale problems and is demonstrated by generating reactive models for the hydrated excess proton and hydroxide ion based directly on condensed phase ab initio molecular dynamics simulations. The resulting models faithfully reproduce the water-ion structural properties and diffusion constants from the ab initio simulations. Additionally, the free energy profiles for proton transfer, which is sensitive to the structural diffusion of both ions in water, are reproduced. The high fidelity of these models to ab initio simulations will permit accurate modeling of general chemical reactions in condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus facilitating a proper statistical sampling of the coupling to slow, large-scale motions of the system.
Collapse
Affiliation(s)
- Chris Knight
- Computing, Environment, and Life Sciences, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | | | | |
Collapse
|
30
|
Kale S, Herzfeld J. Natural polarizability and flexibility via explicit valency: the case of water. J Chem Phys 2012; 136:084109. [PMID: 22380034 DOI: 10.1063/1.3688228] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
As the dominant physiological solvent, water drives the folding of biological macromolecules, influences conformational changes, determines the ionization states of surface groups, actively participates in catalytic events, and provides "wires" for long-range proton transfer. Elucidation of all these roles calls for atomistic simulations. However, currently available methods do not lend themselves to efficient simulation of proton transfer events, or even polarizability and flexibility. Here, we report that an explicit account of valency can provide a unified description for the polarizability, flexibility, and dissociability of water in one intuitive and efficient setting. We call this approach LEWIS, after the chemical theory that inspires the use of valence electron pairs. In this paper, we provide details of the method, the choice of the training set, and predictions for the neat ambient liquid, with emphasis on structure, dynamics, and polarization. LEWIS water provides a good description of bulk properties, and dipolar and quadrupolar responses.
Collapse
Affiliation(s)
- Seyit Kale
- Graduate Program in Biophysics and Structural Biology, Brandeis University, Waltham, Massachusetts 02454-9110, USA
| | | |
Collapse
|
31
|
Bellucci MA, Coker DF. Molecular dynamics of excited state intramolecular proton transfer: 3-hydroxyflavone in solution. J Chem Phys 2012; 136:194505. [DOI: 10.1063/1.4707736] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Hofer TS, Hitzenberger M, Randolf BR. Combining a Dissociative Water Model with a Hybrid QM/MM Approach—A Simulation Strategy for the Study of Proton Transfer Reactions in Solution. J Chem Theory Comput 2012; 8:3586-95. [DOI: 10.1021/ct300062k] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General,
Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Manuel Hitzenberger
- Theoretical Chemistry Division, Institute of General,
Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Bernhard R. Randolf
- Theoretical Chemistry Division, Institute of General,
Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| |
Collapse
|
33
|
Park K, Lin W, Paesani F. A Refined MS-EVB Model for Proton Transport in Aqueous Environments. J Phys Chem B 2011; 116:343-52. [DOI: 10.1021/jp208946p] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Kyoyeon Park
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Wei Lin
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| |
Collapse
|
34
|
Wilhelm F, Schmickler W, Nazmutdinov R, Spohr E. Modeling proton transfer to charged silver electrodes. Electrochim Acta 2011. [DOI: 10.1016/j.electacta.2011.04.036] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
35
|
Somtua T, Tongraar A. Correlation effects on the structure and dynamics of the H3O+ hydrate: B3LYP/MM and MP2/MM MD simulations. Phys Chem Chem Phys 2011; 13:16190-6. [PMID: 21847469 DOI: 10.1039/c1cp20823c] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations, namely B3LYP/MM and MP2/MM, have been performed to investigate the possible influence of electron correlation on the structure and dynamics of the H(3)O(+) hydrate. In comparison to the previously published HF/MM results, both B3LYP/MM and MP2/MM simulations clearly reveal stronger H(3)O(+)-water hydrogen bond interactions, which are reflected in a slightly greater compactness of the H(3)O(+) hydrate. However, the B3LYP/MM simulation, although providing structural details very close to the MP2/MM data, shows an artificially slow dynamic nature of some first shell water molecules as a consequence of the formation of a long-lived H(3)O(+)···H(2)O hydrogen bonding structure.
Collapse
Affiliation(s)
- Thanawat Somtua
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
| | | |
Collapse
|
36
|
Donnini S, Tegeler F, Groenhof G, Grubmüller H. Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics. J Chem Theory Comput 2011; 7:1962-1978. [PMID: 21687785 PMCID: PMC3114466 DOI: 10.1021/ct200061r] [Citation(s) in RCA: 142] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2011] [Indexed: 12/22/2022]
Abstract
pH is an important parameter in condensed-phase systems, because it determines the protonation state of titratable groups and thus influences the structure, dynamics, and function of molecules in solution. In most force field simulation protocols, however, the protonation state of a system (rather than its pH) is kept fixed and cannot adapt to changes of the local environment. Here, we present a method, implemented within the MD package GROMACS, for constant pH molecular dynamics simulations in explicit solvent that is based on the λ-dynamics approach. In the latter, the dynamics of the titration coordinate λ, which interpolates between the protonated and deprotonated states, is driven by generalized forces between the protonated and deprotonated states. The hydration free energy, as a function of pH, is included to facilitate constant pH simulations. The protonation states of titratable groups are allowed to change dynamically during a simulation, thus reproducing average protonation probabilities at a certain pH. The accuracy of the method is tested against titration curves of single amino acids and a dipeptide in explicit solvent.
Collapse
Affiliation(s)
- Serena Donnini
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | | | - Gerrit Groenhof
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| |
Collapse
|
37
|
Agostini F, Vuilleumier R, Ciccotti G. Infrared spectroscopy of small protonated water clusters at room temperature: An effective modes analysis. J Chem Phys 2011; 134:084302. [DOI: 10.1063/1.3533229] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Agostini F, Vuilleumier R, Ciccotti G. Infrared spectroscopy and effective modes analysis of the protonated water dimer H+(H2O)2 at room temperature under H/D substitution. J Chem Phys 2011; 134:084303. [DOI: 10.1063/1.3521273] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
39
|
Abstract
We have performed the multistate empirical valence bond (MS-EVB) molecular dynamics simulations of a dilute hydrofluoric acid solution at ambient temperature to study the hydration structure associated with its weak acidity. The developed MS-EVB model showed reasonable agreement with experimental and previous ab initio molecular dynamics and reference interaction site model self-consistent field simulations for the free energy and structural properties. The local tetrahedral and translational order parameters around the fluorine atom significantly increase in the transition and product states of the HF dissociation reaction. This indicates that the angular and translational rearrangements of the hydrogen-bond topology are necessary especially around the fluorine atom. At the transition state of the proton transfer, the tetrahedral order parameters are very large, whereas the translational order parameters are not. This suggests that for the proton transfer to occur the large angular rearrangements of the hydrogen-bond topology are more necessary than the translational ones.
Collapse
Affiliation(s)
- Tatsuya Joutsuka
- Department of Chemistry, Kyoto University, Sakyo-ku, Kyoto, Japan.
| | | |
Collapse
|
40
|
Bucher D, Gray-Weale A, Kuyucak S. Ab Initio Study of Water Polarization in the Hydration Shell of Aqueous Hydroxide: Comparison between Polarizable and Nonpolarizable Water Models. J Chem Theory Comput 2010; 6:2888-95. [DOI: 10.1021/ct1003719] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Denis Bucher
- School of Physics, University of Sydney, NSW 2006, Australia
| | | | - Serdar Kuyucak
- School of Physics, University of Sydney, NSW 2006, Australia
| |
Collapse
|
41
|
Wick CD, Dang LX. The behavior of NaOH at the air-water interface: A computational study. J Chem Phys 2010; 133:024705. [DOI: 10.1063/1.3455332] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
42
|
Wilhelm F, Schmickler W, Spohr E. Proton transfer to charged platinum electrodes. A molecular dynamics trajectory study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2010; 22:175001. [PMID: 21393659 DOI: 10.1088/0953-8984/22/17/175001] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
A recently developed empirical valence bond (EVB) model for proton transfer on Pt(111) electrodes (Wilhelm et al 2008 J. Phys. Chem. C 112 10814) has been applied in molecular dynamics (MD) simulations of a water film in contact with a charged Pt surface. A total of seven negative surface charge densities σ between -7.5 and -18.9 µC cm(-2) were investigated. For each value of σ, between 30 and 84 initial conditions of a solvated proton within a water slab were sampled, and the trajectories were integrated until discharge of a proton occurred on the charged surfaces. We have calculated the mean rates for discharge and for adsorption of solvated protons within the adsorbed water layer in contact with the metal electrode as a function of surface charge density. For the less negative values of σ we observe a Tafel-like exponential increase of discharge rate with decreasing σ. At the more negative values this exponential increase levels off and the discharge process is apparently transport limited. Mechanistically, the Tafel regime corresponds to a stepwise proton transfer: first, a proton is transferred from the bulk into the contact water layer, which is followed by transfer of a proton to the charged surface and concomitant discharge. At the more negative surface charge densities the proton transfer into the contact water layer and the transfer of another proton to the surface and its discharge occur almost simultaneously.
Collapse
Affiliation(s)
- Florian Wilhelm
- Institute of Theoretical Chemistry, Ulm University, D-89069 Ulm, Germany
| | | | | |
Collapse
|
43
|
Choi TH, Jordan KD. Application of the SCC-DFTB Method to H+(H2O)6, H+(H2O)21, and H+(H2O)22. J Phys Chem B 2010; 114:6932-6. [DOI: 10.1021/jp912289e] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tae Hoon Choi
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| | - Kenneth D. Jordan
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| |
Collapse
|
44
|
Ufimtsev IS, Kalinichev AG, Martinez TJ, Kirkpatrick RJ. A multistate empirical valence bond model for solvation and transport simulations of OH- in aqueous solutions. Phys Chem Chem Phys 2009; 11:9420-30. [PMID: 19830325 DOI: 10.1039/b907859b] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
We describe a new multistate empirical valence bond (MS-EVB) model of OH(-) in aqueous solutions. This model is based on the recently proposed "charged ring" parameterization for the intermolecular interaction of hydroxyl ion with water [Ufimtsev, et al., Chem. Phys. Lett., 2007, 442, 128] and is suitable for classical molecular simulations of OH(-) solvation and transport. The model reproduces the hydration structure of OH(-)(aq) in good agreement with experimental data and the results of ab initio molecular dynamics simulations. It also accurately captures the major structural, energetic, and dynamic aspects of the proton transfer processes involving OH(-) (aq). The model predicts an approximately two-fold increase of the OH(-) mobility due to proton exchange reactions.
Collapse
Affiliation(s)
- Ivan S Ufimtsev
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | | | | | | |
Collapse
|
45
|
Kumar R, Christie RA, Jordan KD. A Modified MSEVB Force Field for Protonated Water Clusters. J Phys Chem B 2008; 113:4111-8. [DOI: 10.1021/jp8066475] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- R. Kumar
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| | - R. A. Christie
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| | - K. D. Jordan
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| |
Collapse
|
46
|
Mitchell-Koch KR, Thompson WH. Infrared spectra of a model phenol-amine proton transfer complex in nanoconfined CH3Cl. J Phys Chem B 2008; 112:7448-59. [PMID: 18517239 DOI: 10.1021/jp076714e] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The vibrational spectra of a model phenol-amine proton transfer complex dissolved in CH3Cl solvent confined in a 12 A radius spherical hydrophobic cavity were calculated using mixed quantum-classical molecular dynamics simulations. The reaction free energy of the proton transfer complex was varied in order to explore the contributions to the vibrational absorption band from product and reactant species. The vibrational spectra of the model proton transfer complex resulted in motionally narrowed spectral linewidths with two distinct peaks for products and reactants in cases where the system undergoes chemical exchange. It was found that the n=1 and n=2 vibrational excited states combine to form diabatic states such that the spectra have contributions from both n=0 --> n=1 and n=0 --> n=2 transitions. A strong relationship between the instantaneous vibrational frequency and a collective solvent coordinate was found that assists in understanding the origin of the spectral features.
Collapse
|
47
|
Kurtz I, Kraut J, Ornekian V, Nguyen MK. Acid-base analysis: a critique of the Stewart and bicarbonate-centered approaches. Am J Physiol Renal Physiol 2008; 294:F1009-31. [PMID: 18184741 DOI: 10.1152/ajprenal.00475.2007] [Citation(s) in RCA: 103] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
When approaching the analysis of disorders of acid-base balance, physical chemists, physiologists, and clinicians, tend to focus on different aspects of the relevant phenomenology. The physical chemist focuses on a quantitative understanding of proton hydration and aqueous proton transfer reactions that alter the acidity of a given solution. The physiologist focuses on molecular, cellular, and whole organ transport processes that modulate the acidity of a given body fluid compartment. The clinician emphasizes the diagnosis, clinical causes, and most appropriate treatment of acid-base disturbances. Historically, two different conceptual frameworks have evolved among clinicians and physiologists for interpreting acid-base phenomena. The traditional or bicarbonate-centered framework relies quantitatively on the Henderson-Hasselbalch equation, whereas the Stewart or strong ion approach utilizes either the original Stewart equation or its simplified version derived by Constable. In this review, the concepts underlying the bicarbonate-centered and Stewart formulations are analyzed in detail, emphasizing the differences in how each approach characterizes acid-base phenomenology at the molecular level, tissue level, and in the clinical realm. A quantitative comparison of the equations that are currently used in the literature to calculate H+concentration ([H+]) is included to clear up some of the misconceptions that currently exist in this area. Our analysis demonstrates that while the principle of electroneutrality plays a central role in the strong ion formulation, electroneutrality mechanistically does not dictate a specific [H+], and the strong ion and bicarbonate-centered approaches are quantitatively identical even in the presence of nonbicarbonate buffers. Finally, our analysis indicates that the bicarbonate-centered approach utilizing the Henderson-Hasselbalch equation is a mechanistic formulation that reflects the underlying acid-base phenomenology.
Collapse
|
48
|
Dahlke EE, Orthmeyer MA, Truhlar DG. Assessment of Multicoefficient Correlation Methods, Second-Order Møller−Plesset Perturbation Theory, and Density Functional Theory for H3O+(H2O)n (n = 1−5) and OH-(H2O)n (n = 1−4). J Phys Chem B 2008; 112:2372-81. [DOI: 10.1021/jp075823q] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Erin E. Dahlke
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Michelle A. Orthmeyer
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
49
|
Ishiyama T, Morita A. Molecular Dynamics Analysis of Interfacial Structures and Sum Frequency Generation Spectra of Aqueous Hydrogen Halide Solutions. J Phys Chem A 2007; 111:9277-85. [PMID: 17705456 DOI: 10.1021/jp072997z] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular dynamics simulation for gas/liquid interfaces of aqueous hydrochloric (HCl) and hydroiodic (HI) acid solutions is performed to calculate and analyze their sum frequency generation (SFG) spectra. The present MD simulation supports the strong preference of hydronium ions at the topmost surface layer and a consequent formation of ionic double layers by the hydronium and halide ions near the interface. Accordingly, the orientational order of surface water in the double layers is reversed in the acid solutions from that in the salt (NaCl or NaI). The calculated SFG spectra of the O-H stretching region reproduce the experimental spectra of ssp and sps polarizations well. In the ssp spectra, the strong enhancement in the hydrogen-bonding region for the acid solutions is elucidated by two mechanisms, ordered orientation of water in the double layer and symmetric OH stretching of the surface hydronium ions. In the sps spectra, reversed orientation of surface water is evidenced in the spectral line shapes, which are quite different from those of the salt solutions.
Collapse
Affiliation(s)
- Tatsuya Ishiyama
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | | |
Collapse
|
50
|
Abstract
A molecular dynamics (MD) simulation based on a combined ab initio quantum mechanics/molecular mechanics (QM/MM) method has been performed to investigate the solvation structure and dynamics of H3O+ in water. The QM region is a sphere around the central H3O+ ion, and contains about 6-8 water molecules. It is treated at the Hartree-Fock (HF) level, while the rest of the system is described by means of classical pair potentials. The Eigen complex (H9O4+) is found to be the most prevalent species in the aqueous solution, partly due to the selection scheme of the center of the QM region. The QM/MM results show that the Eigen complex frequently converts back and forth into the Zundel (H5O2+) structure. Besides the three nearest-neighbor water molecules directly hydrogen-bonded to H3O+, other neighbor waters, such as a fourth water molecule which interacts preferentially with the oxygen atom of the hydronium ion, are found occasionally near the ion. Analyses of the water exchange processes and the mean residence times of water molecules in the ion's hydration shell indicate that such next-nearest neighbor water molecules participate in the rearrangement of the hydrogen bond network during fluctuative formation of the Zundel ion and, thus, contribute to the Grotthuss transport of the proton.
Collapse
Affiliation(s)
- Pathumwadee Intharathep
- School of Chemistry, Institute of Science, Suranaree University of Technology, 30000 Nakhon Ratchasima, Thailand
| | | | | |
Collapse
|