1
|
Depew DD, Vaghjiani GL, Parmar SM, Wang JJ. Liquid Structure and Hydrogen Bonding in Aqueous Hydroxylammonium Nitrate. J Phys Chem B 2024; 128:824-840. [PMID: 38194505 DOI: 10.1021/acs.jpcb.3c05623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2024]
Abstract
Hydroxylammonium nitrate (HAN) has emerged as a promising component in ionic liquid-based spacecraft propellants. However, the physicochemical and structural properties of aqueous HAN have been largely overlooked. The purpose of this study is to investigate the hydrogen bonding in aqueous HAN and understand its implications on these properties and the proton transfer mechanism as a function of concentration. Classical polarizable molecular dynamics simulations have been employed with the APPLE&P force field to analyze the geometry of individual hydrogen bonds and the overall hydrogen-bonding network in various concentrations of aqueous HAN. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) indicate the structural arrangement of the species and their hydrogen bonds. Projections of water density and the orientation of its electric dipole moment near the ions provide insight into the hydrogen-bonding network. The incorporation of water into the hydrogen-bonding network at high ion concentrations occurs via interstitial accommodation around the ions immediately outside the first solvation shell. While ion pairs are observed at all concentrations considered, the frequency of Ha···On hydrogen bonds increases substantially with the ion concentration. The findings contribute to a better fundamental understanding of HAN and the precursors of reactivity, crucial to the development of "green" spacecraft propellants.
Collapse
Affiliation(s)
- Daniel D Depew
- Department of Astronautical Engineering, University of Southern California, Los Angeles, California 90089, United States
| | - Ghanshyam L Vaghjiani
- Aerospace Systems Directorate, Air Force Research Laboratory, AFRL/RQRS, Edwards Air Force Base, California 93524, United States
| | - Shehan M Parmar
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Joseph J Wang
- Department of Astronautical Engineering, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
2
|
Chen YL, Nie HN, Dong H, Gao Q, Peng X, Wang N, Chen X, Liu QZ, Li JK, Xu XB, Xue YL. Revealing the mechanism underlying the effects of γ-aminobutyric acid-dioscorin interactions on dioscorin structure and emulsifying properties by molecular dynamic simulations. Food Res Int 2023; 171:112982. [PMID: 37330840 DOI: 10.1016/j.foodres.2023.112982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Revised: 04/04/2023] [Accepted: 05/15/2023] [Indexed: 06/19/2023]
Abstract
Many studies have shown that γ-aminobutyric acid (GABA) exhibits various beneficial biological activities, including gut-modulating, neuro-stimulating, and cardio-protecting activities. Naturally, GABA exists in small amounts in yam, which is primarily synthesized by the decarboxylation of L-glutamic acid in the presence of glutamate decarboxylase. Dioscorin, the major tuber storage protein of yam, has been shown to have good solubility and emulsifying activity. However, how GABA interacts with dioscorin and affects their properties has yet to be clarified. In this research, the physicochemical and emulsifying properties of GABA-fortified dioscorin, which was dried by spray drying and freeze drying, were studied. As results, the freeze-dried (FD) dioscorin produced more stable emulsions, while the spray-dried (SD) dioscorin adsorbed more rapidly to oil/water (O/W) interface. The fluorescence spectroscopy, ultraviolet spectroscopy and circular dichroism spectroscopy showed that GABA changed the structure of dioscorin, by exposing its hydrophobic groups. The addition of GABA significantly promoted the adsorption of dioscorin to the O/W interface and prevented droplets coalescence. The results of molecular dynamics simulation (MD) showed that GABA destroyed the H-bond network between dioscorin and water, increased surface hydrophobicity and finally improved the emulsifying properties of dioscorin.
Collapse
Affiliation(s)
- Yun-Long Chen
- College of Light Industry, Liaoning University, Shenyang 110036, China
| | - Hao-Nan Nie
- College of Light Industry, Liaoning University, Shenyang 110036, China
| | - Hui Dong
- College of Light Industry, Liaoning University, Shenyang 110036, China
| | - Qi Gao
- College of Light Industry, Liaoning University, Shenyang 110036, China; Department of Regional Economic Development, Party School of Liaoning Provincial Party Committee, Shenyang 110161, China
| | - Xue Peng
- College of Light Industry, Liaoning University, Shenyang 110036, China
| | - Ning Wang
- College of Light Industry, Liaoning University, Shenyang 110036, China
| | - Xueling Chen
- Institute for Farm Products Processing and Nuclear-agricultural Technology, Hubei Academy of Agricultural Sciences, Wuhan 430064, China
| | | | - Jiang-Kuo Li
- Tianjin Academy of Agricultural Sciences, National Engineering and Technology Research Center for Preservation of Agricultural Products (Tianjin), Tianjin 300384, China
| | - Xiang-Bin Xu
- College of Food Science and Engineering, Hainan University, Haikou 570228, China
| | - You-Lin Xue
- College of Light Industry, Liaoning University, Shenyang 110036, China.
| |
Collapse
|
3
|
de Oliveira PMC, de Souza JIR, da Silva JAB, Longo RL. Temperature Dependence of Hydrogen Bond Networks of Liquid Water: Thermodynamic Properties and Structural Heterogeneity from Topological Descriptors. J Phys Chem B 2023; 127:2250-2257. [PMID: 36877152 DOI: 10.1021/acs.jpcb.2c08873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
Topological analyses of hydrogen bond networks were performed based on the complex network and island statistics of liquid water at different temperatures. The influence of temperature on the liquid water structures and the topological properties of the hydrogen bond networks was investigated by Metropolis Monte Carlo simulations with the TIP4P/2005 potential model. The bilinear behavior of the second peak in the radial distribution function with the temperature was properly reproduced by these simulations. The average connectivity also displayed a bilinear behavior consistent with being a local descriptor. The semiglobal average path length (or geodesic distance) descriptor showed an unprecedented trimodal distribution, whose areas were dependent on the temperature. Considering equilibrium between these three sets of networks, standard enthalpy and entropy of equilibrium were determined for the first time, providing new insights into the structural heterogeneities of liquid water with interesting perspectives for modeling these quantitative properties of hydrogen bond networks.
Collapse
Affiliation(s)
- Paulo M C de Oliveira
- Núcleo Interdisciplinar de Ciências Exatas e da Natureza, Campus do Agreste. Universidade Federal de Pernambuco, 55.014-900 Caruaru, Pernambuco, Brazil
| | - Jéssica I R de Souza
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50740-560 Recife, Pernambuco, Brazil
| | - Juliana A B da Silva
- Núcleo Interdisciplinar de Ciências Exatas e da Natureza, Campus do Agreste. Universidade Federal de Pernambuco, 55.014-900 Caruaru, Pernambuco, Brazil.,Programa de Pós-Graduação em Química, Universidade Federal Rural de Pernambuco, 52.171-900 Recife, Pernambuco, Brazil
| | - Ricardo L Longo
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50740-560 Recife, Pernambuco, Brazil.,Departamento de Química Fundamental, Universidade Federal de Pernambuco, 50740-540 Recife, Pernambuco, Brazil
| |
Collapse
|
4
|
Biswas A, Mallik BS. Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations. J Phys Chem A 2022; 126:5134-5147. [PMID: 35900106 DOI: 10.1021/acs.jpca.2c03387] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.
Collapse
Affiliation(s)
- Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502284, Telangana, India
| | - Bhabani S Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502284, Telangana, India
| |
Collapse
|
5
|
Biswas A, Mallik BS. 2D IR spectra of the intrinsic vibrational probes of ionic liquid from dispersion corrected DFT-MD simulations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118390] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
6
|
Biswas A, Mallik BS. Dynamics of Ionic Liquid through Intrinsic Vibrational Probes Using the Dispersion-Corrected DFT Functionals. J Phys Chem B 2021; 125:6994-7008. [PMID: 34142827 DOI: 10.1021/acs.jpcb.1c04960] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
First principles molecular dynamics simulations have been utilized to study the spectral properties of the protic ionic liquid, methylammonium formate (MAF). All simulations were performed using density functional theory (DFT) and various van der Waals-corrected exchange-correlation functionals. We calculated the vibrational stretch frequency distributions, determined the time-frequency correlations of the intrinsic vibrational probes, the N-H and C-O modes in MAF, and the frequency-structure correlations. We also estimated the average hydrogen-bond lifetimes and orientation dynamics to capture the ultrafast spectral response. The spectroscopic signature of the N-H stretching vibrations using the Becke-Lee-Yang-Parr (BLYP) and Perdew-Burke-Ernzerhof (PBE) functionals displays a spectral shift in the lower frequency side, suggesting stronger hydrogen-bonding interactions represented by the gradient approximation functionals than the van der Waals (vdW)-corrected simulations. The carboxylate frequency profiles with the dispersion-corrected representations are almost similar without a significant difference in the normalized distributions. Besides, the COO stretching frequencies at the peak maxima positions of the PBE functionals exhibit a lesser deviation from the experimental data. Spectral diffusion dynamics of the intrinsic vibrational probes on the cationic and anionic sites of the ionic liquid proceed through a short time relaxation of the intact hydrogen bonds followed by an intermediate time constant and a longer time decay indicating the switchover of hydrogen bonds. Dispersion-corrected atom-centered one-electron potential (DCACP) correction added to the BLYP system slows down the picosecond time scales of frequency correlation and the time constants of rotational motion, lengthening the overall system dynamics. The observed trends in the time-dependent decays of frequency fluctuations and the orientation autocorrelation functions correlate with the structural interactions in liquid MAF and hydrogen-bond dynamics. In this study, we examine the predictions made by different density functional treatments comparing the results of the uncorrected BLYP and PBE representations with the semiempirical vdW methods of Grimme and matching our calculated data with the experimental observations.
Collapse
Affiliation(s)
- Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502285, Telangana, India
| | - Bhabani S Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502285, Telangana, India
| |
Collapse
|
7
|
Biswas A, Dasari S, Mallik BS. Cohesiveness and Nondiffusive Rotational Jump Dynamics of Protic Ionic Liquid from Dispersion-Corrected FPMD Simulations. J Phys Chem B 2020; 124:10752-10765. [DOI: 10.1021/acs.jpcb.0c05866] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502285, Telangana, India
| | - Sathish Dasari
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502285, Telangana, India
| | - Bhabani S. Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502285, Telangana, India
| |
Collapse
|
8
|
Lynch C, Rao S, Sansom MSP. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem Rev 2020; 120:10298-10335. [PMID: 32841020 PMCID: PMC7517714 DOI: 10.1021/acs.chemrev.9b00830] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Indexed: 12/18/2022]
Abstract
This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
Collapse
Affiliation(s)
- Charlotte
I. Lynch
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Shanlin Rao
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| |
Collapse
|
9
|
Dohn AO, Jónsson EÖ, Jónsson H. Polarizable Embedding with a Transferable H 2O Potential Function II: Application to (H 2O) n Clusters and Liquid Water. J Chem Theory Comput 2019; 15:6578-6587. [PMID: 31692344 DOI: 10.1021/acs.jctc.9b00778] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The incorporation of polarization in multiscale quantum-mechanics/molecular-mechanics (QM/MM) simulations is important for a variety of applications, for example, charge-transfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H2O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H2O)n clusters shows that the interquartile range of the distributions of the QM/MM binding energies is never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme. Optimized hexamer geometries as well as the liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.
Collapse
Affiliation(s)
- Asmus Ougaard Dohn
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| | - Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| |
Collapse
|
10
|
Cheng B, Engel EA, Behler J, Dellago C, Ceriotti M. Ab initio thermodynamics of liquid and solid water. Proc Natl Acad Sci U S A 2019; 116:1110-1115. [PMID: 30610171 PMCID: PMC6347673 DOI: 10.1073/pnas.1815117116] [Citation(s) in RCA: 152] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.
Collapse
Affiliation(s)
- Bingqing Cheng
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland;
| | - Edgar A Engel
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Universität Göttingen, 37073 Göttingen, Germany
| | | | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
11
|
Dohn AO, Jónsson EÖ, Levi G, Mortensen JJ, Lopez-Acevedo O, Thygesen KS, Jacobsen KW, Ulstrup J, Henriksen NE, Møller KB, Jónsson H. Grid-Based Projector Augmented Wave (GPAW) Implementation of Quantum Mechanics/Molecular Mechanics (QM/MM) Electrostatic Embedding and Application to a Solvated Diplatinum Complex. J Chem Theory Comput 2017; 13:6010-6022. [PMID: 29083921 DOI: 10.1021/acs.jctc.7b00621] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
A multiscale density functional theory-quantum mechanics/molecular mechanics (DFT-QM/MM) scheme is presented, based on an efficient electrostatic coupling between the electronic density obtained from a grid-based projector augmented wave (GPAW) implementation of density functional theory and a classical potential energy function. The scheme is implemented in a general fashion and can be used with various choices for the descriptions of the QM or MM regions. Tests on H2O clusters, ranging from dimer to decamer show that no systematic energy errors are introduced by the coupling that exceeds the differences in the QM and MM descriptions. Over 1 ns of liquid water, Born-Oppenheimer QM/MM molecular dynamics (MD) are sampled combining 10 parallel simulations, showing consistent liquid water structure over the QM/MM border. The method is applied in extensive parallel MD simulations of an aqueous solution of the diplatinum [Pt2(P2O5H2)4]4- complex (PtPOP), spanning a total time period of roughly half a nanosecond. An average Pt-Pt distance deviating only 0.01 Å from experimental results, and a ground-state Pt-Pt oscillation frequency deviating by <2% from experimental results were obtained. The simulations highlight a remarkable harmonicity of the Pt-Pt oscillation, while also showing clear signs of Pt-H hydrogen bonding and directional coordination of water molecules along the Pt-Pt axis of the complex.
Collapse
Affiliation(s)
- A O Dohn
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - E Ö Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - G Levi
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J J Mortensen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - O Lopez-Acevedo
- Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| | - K S Thygesen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K W Jacobsen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J Ulstrup
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - N E Henriksen
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K B Møller
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - H Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland.,Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| |
Collapse
|
12
|
Abstract
Water challenges our fundamental understanding of emergent materials properties from a molecular perspective. It exhibits a uniquely rich phenomenology including dramatic variations in behavior over the wide temperature range of the liquid into water's crystalline phases and amorphous states. We show that many-body responses arising from water's electronic structure are essential mechanisms harnessed by the molecule to encode for the distinguishing features of its condensed states. We treat the complete set of these many-body responses nonperturbatively within a coarse-grained electronic structure derived exclusively from single-molecule properties. Such a "strong coupling" approach generates interaction terms of all symmetries to all orders, thereby enabling unique transferability to diverse local environments such as those encountered along the coexistence curve. The symmetries of local motifs that can potentially emerge are not known a priori. Consequently, electronic responses unfiltered by artificial truncation are then required to embody the terms that tip the balance to the correct set of structures. Therefore, our fully responsive molecular model produces, a simple, accurate, and intuitive picture of water's complexity and its molecular origin, predicting water's signature physical properties from ice, through liquid-vapor coexistence, to the critical point.
Collapse
|
13
|
Molecular dynamics simulation study of solvation effects of water and trifluoroethanol on gamma-aminobutyric acid (GABA). J Mol Liq 2014. [DOI: 10.1016/j.molliq.2014.04.022] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
14
|
|
15
|
Thermal stability of water up to super-critical states: Application of the singular value decomposition and grund functions. J Mol Liq 2013. [DOI: 10.1016/j.molliq.2013.07.013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Jones A, Cipcigan F, Sokhan VP, Crain J, Martyna GJ. Electronically coarse-grained model for water. PHYSICAL REVIEW LETTERS 2013; 110:227801. [PMID: 23767748 DOI: 10.1103/physrevlett.110.227801] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2012] [Indexed: 06/02/2023]
Abstract
We introduce an electronically coarse-grained description of water representing all long range, many-body electronic responses via an embedded quantum oscillator. Leading-order response coefficients and gas phase electrostatic moments are exactly reproduced. Molecular dynamics, using electronic path integral sampling, shows that this framework is sufficient for a realistic liquid to emerge naturally with transferability extending further to nonambient state points and to the free water surface. The model allows the strength of many-body dispersion and polarization to be adjusted independently and these are found to have significant effects on the condensed phase.
Collapse
Affiliation(s)
- A Jones
- School of Physics and Astronomy, University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, United Kingdom
| | | | | | | | | |
Collapse
|
17
|
Khaliullin RZ, Kühne TD. Microscopic properties of liquid water from combined ab initio molecular dynamics and energy decomposition studies. Phys Chem Chem Phys 2013; 15:15746-66. [DOI: 10.1039/c3cp51039e] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
18
|
Ma Z, Zhang Y, Tuckerman ME. Ab initio molecular dynamics study of water at constant pressure using converged basis sets and empirical dispersion corrections. J Chem Phys 2012; 137:044506. [PMID: 22852630 DOI: 10.1063/1.4736712] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It is generally believed that studies of liquid water using the generalized gradient approximation to density functional theory require dispersion corrections in order to obtain reasonably accurate structural and dynamical properties. Here, we report on an ab initio molecular dynamics study of water in the isothermal-isobaric ensemble using a converged discrete variable representation basis set and an empirical dispersion correction due to Grimme [J. Comp. Chem. 27, 1787 (2006)]. At 300 K and an applied pressure of 1 bar, the density obtained without dispersion corrections is approximately 0.92 g/cm(3) while that obtained with dispersion corrections is 1.07 g/cm(3), indicating that the empirical dispersion correction overestimates the density by almost as much as it is underestimated without the correction for this converged basis. Radial distribution functions exhibit a loss of structure in the second solvation shell. Comparison of our results with other studies using the same empirical correction suggests the cause of the discrepancy: the Grimme dispersion correction is parameterized for use with a particular basis set; this parameterization is sensitive to this choice and, therefore, is not transferable to other basis sets.
Collapse
Affiliation(s)
- Zhonghua Ma
- Department of Chemistry, New York University, New York, New York 10003, USA.
| | | | | |
Collapse
|
19
|
Lin IC, Seitsonen AP, Tavernelli I, Rothlisberger U. Structure and Dynamics of Liquid Water from ab Initio Molecular Dynamics—Comparison of BLYP, PBE, and revPBE Density Functionals with and without van der Waals Corrections. J Chem Theory Comput 2012; 8:3902-10. [DOI: 10.1021/ct3001848] [Citation(s) in RCA: 216] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- I-Chun Lin
- Department of Chemistry,
New
York University, 100 Washington Square East, New York, New York 10003,
United States
| | - Ari P. Seitsonen
- Physikalisch-Chemisches
Institut,
University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Ivano Tavernelli
- Laboratory of Computational Chemistry
and Biochemistry, Ecole Polytechnique Fédérale de Lausanne,
CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry
and Biochemistry, Ecole Polytechnique Fédérale de Lausanne,
CH-1015 Lausanne, Switzerland
| |
Collapse
|
20
|
Lin IC, Seitsonen AP, Tavernelli I, Rothlisberger U. Structure and Dynamics of Liquid Water from ab Initio Molecular Dynamics—Comparison of BLYP, PBE, and revPBE Density Functionals with and without van der Waals Corrections. J Chem Theory Comput 2012. [DOI: 10.10.1021/ct3001848] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Affiliation(s)
- I-Chun Lin
- Department of Chemistry,
New
York University, 100 Washington Square East, New York, New York 10003,
United States
| | - Ari P. Seitsonen
- Physikalisch-Chemisches
Institut,
University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Ivano Tavernelli
- Laboratory of Computational Chemistry
and Biochemistry, Ecole Polytechnique Fédérale de Lausanne,
CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry
and Biochemistry, Ecole Polytechnique Fédérale de Lausanne,
CH-1015 Lausanne, Switzerland
| |
Collapse
|
21
|
|
22
|
Temelso B, Archer KA, Shields GC. Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections. J Phys Chem A 2011; 115:12034-46. [DOI: 10.1021/jp2069489] [Citation(s) in RCA: 256] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Berhane Temelso
- Dean’s Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Kaye A. Archer
- Dean’s Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - George C. Shields
- Dean’s Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
23
|
Paesani F. Hydrogen bond dynamics in heavy water studied with quantum dynamical simulations. Phys Chem Chem Phys 2011; 13:19865-75. [PMID: 21892511 DOI: 10.1039/c1cp21863h] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The structure and dynamics of the hydrogen-bond network in heavy water (D(2)O) is studied as a function of the temperature using quantum dynamical simulations. Our approach combines an ab initio-based representation of the water interactions with an explicit quantum treatment of the molecular motion. A direct connection between the calculated linear and nonlinear vibrational spectra and the underlying molecular dynamics is made, which provides new insights into the rearrangement of the hydrogen-bond network in heavy water. A comparison with previous calculations on liquid H(2)O suggests that tunneling does not effectively contribute to the dynamics of the water hydrogen-bond network above the melting point. However, the effects of nuclear quantization are not negligible at all temperatures and become increasingly important near the melting point, which is in agreement with recent experimental analysis of the structural properties of liquid water as well as of the proton momentum distribution in supercooled water.
Collapse
Affiliation(s)
- Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, La Jolla, San Diego, CA 92093, USA.
| |
Collapse
|
24
|
Medina J, Prosmiti R, Villarreal P, Delgado-Barrio G, Winter G, González B, Alemán J, Collado C. Molecular dynamics simulations of rigid and flexible water models: Temperature dependence of viscosity. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.07.001] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
25
|
Shields RM, Temelso B, Archer KA, Morrell TE, Shields GC. Accurate predictions of water cluster formation, (H₂O)(n=2-10). J Phys Chem A 2011; 114:11725-37. [PMID: 20882961 DOI: 10.1021/jp104865w] [Citation(s) in RCA: 161] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2-10 water molecules. This approach reveals new low energy conformers for (H(2)O)(n=7,9,10). The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.
Collapse
Affiliation(s)
- Robert M Shields
- Dean's Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, USA
| | | | | | | | | |
Collapse
|
26
|
Shimazaki T, Kubo M. Analytical energy gradient of Gaussian and Fourier transform (GFT) method for periodic condensed systems. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2010.12.089] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
27
|
da Silva JAB, Moreira FGB, dos Santos VML, Longo RL. Hydrogen bond networks in water and methanol with varying interaction strengths. Phys Chem Chem Phys 2011; 13:593-603. [DOI: 10.1039/c0cp01204a] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
Gebauer JS, Treuel L. Influence of individual ionic components on the agglomeration kinetics of silver nanoparticles. J Colloid Interface Sci 2010; 354:546-54. [PMID: 21146829 DOI: 10.1016/j.jcis.2010.11.016] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2010] [Revised: 10/25/2010] [Accepted: 11/09/2010] [Indexed: 10/18/2022]
Abstract
The precise characteristic of the agglomeration behavior of colloidal suspensions is of paramount interest to many current studies in nanoscience. This work seeks to elucidate the influence that differently charged salts have on the agglomeration state of a Lee-Meisel-type silver colloid. Moreover, we investigate the influence of the chemical nature of individual ions on their potential to induce agglomeration. Raman spectroscopy and surface-enhanced Raman spectroscopy are used to give insights into mechanistic aspects of the agglomeration process and to assess the differences in the influence of different salts on the agglomeration behavior. Finally, we demonstrate the potential of the measurement procedure used in this work to determine the elementary charge on colloidal NPs.
Collapse
Affiliation(s)
- J S Gebauer
- Institute for Physical Chemistry and Centre for Water and Environmental Research, University of Duisburg-Essen, Universitaetsstrasse 5-7, 45117 Essen, Germany
| | | |
Collapse
|
29
|
|
30
|
Probing the Equilibrium Size and Hydrogen Bonding Structure in Aqueous Aerosol Droplets. ACTA ACUST UNITED AC 2010. [DOI: 10.1524/zpch.2010.6147] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Abstract
Comparative measurements on two or more optically trapped aerosol droplets can allow the hygroscopic properties of aerosol to be investigated in exquisite detail. In this paper we consider two methods for assessing the dependence on relative humidity of the water activity and solute concentration in the condensed phase. We demonstrate in this paper that using a control droplet as a microprobe of relative humidity in the vicinity of a second droplet of interest can allow the RH to be determined with an accuracy of between 0.01 to 0.1 %. Not only is such a probe highly responsive, but the accuracy of the RH determination improves with increase in GF or RH, counter to the performance of macroscopic RH probes. We also demonstrate that spontaneous Raman scattering recorded from excitation of the OH stretching vibration in supersaturated solution droplets of sodium chloride shows the considerable perturbation induced by the chloride ion on the hydrogen bonding network. Chloride is recognised as a hydrogen bonding structure breaker and these observations are shown to be consistent with measurements on more dilute solutions.
Collapse
|
31
|
Jalili S, Akhavan M. Molecular dynamics simulation study of association in trifluoroethanol/water mixtures. J Comput Chem 2010; 31:286-94. [PMID: 19475633 DOI: 10.1002/jcc.21316] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Mixtures of Trifluoroethanol (TFE) and water with different proportions are studied using molecular dynamics simulations. The radial and spatial distribution functions, as well as the size distribution of TFE clusters are obtained from the trajectories. The variation of radial and spatial distribution functions with composition show that the addition of TFE enhances the water structure, but the hydrogen bonds between TFE molecules are broken as TFE is diluted with water. The TFE-rich solutions have stronger TFE-water hydrogen bonds. The clustering of TFE molecules in low concentration region is attributed to the hydrophobic interactions between CF(3) groups. The distribution of cluster sizes in solution supports these conclusions.
Collapse
Affiliation(s)
- S Jalili
- Department of Chemistry and Nanosciences Research Center, K. N. Toosi University of Technology, Tehran, Iran.
| | | |
Collapse
|
32
|
Ranaghan KE, Mulholland AJ. Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods. INT REV PHYS CHEM 2010. [DOI: 10.1080/01442350903495417] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
33
|
Furmanchuk A, Isayev O, Shishkin OV, Gorb L, Leszczynski J. Hydration of nucleic acid bases: a Car–Parrinello molecular dynamics approach. Phys Chem Chem Phys 2010; 12:3363-75. [DOI: 10.1039/b923930h] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
Berkelbach TC, Lee HS, Tuckerman ME. Concerted hydrogen-bond dynamics in the transport mechanism of the hydrated proton: a first-principles molecular dynamics study. PHYSICAL REVIEW LETTERS 2009; 103:238302. [PMID: 20366181 DOI: 10.1103/physrevlett.103.238302] [Citation(s) in RCA: 152] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/17/2009] [Indexed: 05/07/2023]
Abstract
First-principles molecular dynamics calculations performed in a fully converged basis set are used to reveal new details about the mechanism of the anomalous proton-transport process in water, a fundamental question dating back over 200 years. By separating actual structural diffusion from simple rattling events, wherein a proton shuttles forth and back in a hydrogen bond, it is found that the former are driven by a concerted mechanism in which hydronium begins to accept a hydrogen bond from a donor water molecule while the proton-receiving water molecule simultaneously loses one of its acceptor hydrogen bonds. The kinetics of the process are found to be in good agreement with recent experiments.
Collapse
|
35
|
Bakker HJ, Skinner JL. Vibrational Spectroscopy as a Probe of Structure and Dynamics in Liquid Water. Chem Rev 2009; 110:1498-517. [DOI: 10.1021/cr9001879] [Citation(s) in RCA: 586] [Impact Index Per Article: 39.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- H. J. Bakker
- FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands, and Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - J. L. Skinner
- FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands, and Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| |
Collapse
|
36
|
Paesani F, Xantheas SS, Voth GA. Infrared Spectroscopy and Hydrogen-Bond Dynamics of Liquid Water from Centroid Molecular Dynamics with an Ab Initio-Based Force Field. J Phys Chem B 2009; 113:13118-30. [DOI: 10.1021/jp907648y] [Citation(s) in RCA: 114] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Francesco Paesani
- Center for Biophysical Modeling and Simulation and Department of Chemistry, University of Utah, 315 South 1400 East Room 2020, Salt Lake City, Utah 84112-0850, and Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352
| | - Sotiris S. Xantheas
- Center for Biophysical Modeling and Simulation and Department of Chemistry, University of Utah, 315 South 1400 East Room 2020, Salt Lake City, Utah 84112-0850, and Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352
| | - Gregory A. Voth
- Center for Biophysical Modeling and Simulation and Department of Chemistry, University of Utah, 315 South 1400 East Room 2020, Salt Lake City, Utah 84112-0850, and Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352
| |
Collapse
|
37
|
Tulip PR, Bates SP. Peptide aggregation and solvent electrostriction in a simple zwitterionic dipeptide via molecular dynamics simulations. J Chem Phys 2009; 131:015103. [DOI: 10.1063/1.3160682] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- P. R. Tulip
- Scottish Universities Physics Alliance (SUPA), School of Physics and Astronomy, The University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, United Kingdom
| | - S. P. Bates
- Scottish Universities Physics Alliance (SUPA), School of Physics and Astronomy, The University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, United Kingdom
| |
Collapse
|
38
|
Liu LM, Krack M, Michaelides A. Interfacial water: A first principles molecular dynamics study of a nanoscale water film on salt. J Chem Phys 2009; 130:234702. [DOI: 10.1063/1.3152845] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
39
|
Paesani F, Voth GA. The properties of water: insights from quantum simulations. J Phys Chem B 2009; 113:5702-19. [PMID: 19385690 DOI: 10.1021/jp810590c] [Citation(s) in RCA: 189] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The properties of water play a central role in many phenomena of relevance to different areas of science, including physics, chemistry, biology, geology, and climate research. Although well studied for decades, the behavior of water under different conditions and in different environments still remains mysterious and often surprising. In this article, various efforts aimed at providing a comprehensive representation of the water properties at a molecular level through computer modeling and simulation will be described. In particular, the unique role played by the hydrogen-bond network will be examined, first in liquid water, then in the solvation of model biological compounds, and finally in ice, especially highlighting the important effects related to the quantization of the nuclear motion.
Collapse
Affiliation(s)
- Francesco Paesani
- Center for Biophysical Modeling and Simulation and Department of Chemistry, University of Utah, 315 South 1400 East Room 2020, Salt Lake City, Utah 84112-0850, USA
| | | |
Collapse
|
40
|
Schmidt DA, Scipioni R, Boero M. Water Solvation Properties: An Experimental and Theoretical Investigation of Salt Solutions at Finite Dilution. J Phys Chem A 2009; 113:7725-9. [DOI: 10.1021/jp9016932] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Diedrich A. Schmidt
- International Center for Young Scientists, National Institute for Materials Science, Tsukuba, 305-0044, Japan, Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany, Institut de Physique et Chimie des Matériaux de Strasbourg, UMR 7504 CNRS-UDS, 23 rue du Loess, BP 43, F-67034 Strasbourg, France, and CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Roberto Scipioni
- International Center for Young Scientists, National Institute for Materials Science, Tsukuba, 305-0044, Japan, Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany, Institut de Physique et Chimie des Matériaux de Strasbourg, UMR 7504 CNRS-UDS, 23 rue du Loess, BP 43, F-67034 Strasbourg, France, and CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Mauro Boero
- International Center for Young Scientists, National Institute for Materials Science, Tsukuba, 305-0044, Japan, Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany, Institut de Physique et Chimie des Matériaux de Strasbourg, UMR 7504 CNRS-UDS, 23 rue du Loess, BP 43, F-67034 Strasbourg, France, and CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
41
|
Lin IC, Seitsonen AP, Coutinho-Neto MD, Tavernelli I, Rothlisberger U. Importance of van der Waals interactions in liquid water. J Phys Chem B 2009; 113:1127-31. [PMID: 19123911 DOI: 10.1021/jp806376e] [Citation(s) in RCA: 154] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present ab initio molecular dynamics studies on liquid water using density functional theory in conjunction with either dispersion-corrected atom-centered potentials or empirical van der Waals corrections. Our results show that improving the description of van der Waals interactions in DFT-GGA leads to a softening of liquid water's structure with higher mobility. The results obtained with dispersion-corrected atom-centered potentials are especially encouraging. In particular, the radial distribution functions are in better agreement with experiment, and the self-diffusion coefficient increases by more than three-fold compared with the one predicted by the BLYP functional. This work demonstrates that van der Waals interactions are essential in fine-tuning both structural and dynamical properties of liquid water.
Collapse
Affiliation(s)
- I-Chun Lin
- Ecole Polytechnique Federale de Lausanne, CH-1015 Lausanne, Switzerland
| | | | | | | | | |
Collapse
|
42
|
Markovitch O, Agmon N. Reversible geminate recombination of hydrogen-bonded water molecule pair. J Chem Phys 2009; 129:084505. [PMID: 19044833 DOI: 10.1063/1.2968608] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The (history independent) autocorrelation function for a hydrogen-bonded water molecule pair, calculated from classical molecular dynamics trajectories of liquid water, exhibits a t(-3/2) asymptotic tail. Its whole time dependence agrees quantitatively with the solution for reversible diffusion-influenced geminate recombination derived by Agmon and Weiss [J. Chem. Phys. 91, 6937 (1989)]. Agreement with diffusion theory is independent of the precise definition of the bound state. Given the water self-diffusion constant, this theory enables us to determine the dissociation and bimolecular recombination rate parameters for a water dimer. (The theory is indispensable for obtaining the bimolecular rate coefficient.) Interestingly, the activation energies obtained from the temperature dependence of these rate coefficients are similar, rather than differing by the hydrogen-bond (HB) strength. This suggests that recombination requires displacing another water molecule, which meanwhile occupied the binding site. Because these activation energies are about twice the HB strength, cleavage of two HBs may be required to allow pair separation. The autocorrelation function without the HB angular restriction yields a recombination rate coefficient that is larger than that for rebinding to all four tetrahedral water sites (with angular restrictions), suggesting the additional participation of interstitial sites. Following dissociation, the probability of the pair to be unbound but within the reaction sphere rises more slowly than expected, possibly because binding to the interstitial sites delays pair separation. An extended diffusion model, which includes an additional binding site, can account for this behavior.
Collapse
Affiliation(s)
- Omer Markovitch
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | | |
Collapse
|
43
|
Kühne TD, Krack M, Parrinello M. Static and Dynamical Properties of Liquid Water from First Principles by a Novel Car−Parrinello-like Approach. J Chem Theory Comput 2009; 5:235-41. [DOI: 10.1021/ct800417q] [Citation(s) in RCA: 175] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Thomas D. Kühne
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Matthias Krack
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| |
Collapse
|
44
|
Tu Y, Fang H. Anomalies of liquid water at low temperature due to two types of hydrogen bonds. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:016707. [PMID: 19257166 DOI: 10.1103/physreve.79.016707] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2008] [Revised: 01/08/2009] [Indexed: 05/27/2023]
Abstract
It is believed that water anomalies are related to the hydrogen bonds between water molecules, but the relationship still remains a challenge. Here, by introducing two types of hydrogen bond--strong and weak, we successfully reproduce many liquid water anomalies, particularly the thermodynamic anomalies at supercooled temperature. We find that the exchange of strong and weak H bonds, which enhance the competition between the open and collapsed structures of liquid water, might be the key to those anomalies of supercooled liquid water. Our study is helpful in understanding the hydrogen bond structure in liquid water.
Collapse
Affiliation(s)
- Yusong Tu
- Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P. O. Box 800-204, Shanghai 201800, China
| | | |
Collapse
|
45
|
Khaliullin R, Bell A, Head-Gordon M. Electron Donation in the Water-Water Hydrogen Bond. Chemistry 2008; 15:851-5. [DOI: 10.1002/chem.200802107] [Citation(s) in RCA: 122] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
46
|
Douady J, Calvo F, Spiegelman F. Structure, stability, and infrared spectroscopy of (H2O)nNH4+ clusters: A theoretical study at zero and finite temperature. J Chem Phys 2008; 129:154305. [DOI: 10.1063/1.2987304] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
Nair NN, Schreiner E, Marx D. Peptide Synthesis in Aqueous Environments: The Role of Extreme Conditions on Amino Acid Activation. J Am Chem Soc 2008; 130:14148-60. [DOI: 10.1021/ja802370c] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Nisanth N. Nair
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Eduard Schreiner
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
48
|
Abstract
A path-integral Car-Parrinello molecular dynamics simulation of liquid water and ice is performed. It is found that the inclusion of nuclear quantum effects systematically improves the agreement of first-principles simulations of liquid water with experiment. In addition, the proton momentum distribution is computed utilizing a recently developed open path-integral molecular dynamics methodology. It is shown that these results are in good agreement with experimental data.
Collapse
Affiliation(s)
- Joseph A Morrone
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
49
|
High resolution X-ray emission spectroscopy of liquid water: The observation of two structural motifs. Chem Phys Lett 2008. [DOI: 10.1016/j.cplett.2008.04.077] [Citation(s) in RCA: 298] [Impact Index Per Article: 18.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
50
|
Kumar R, Skinner JL. Water Simulation Model with Explicit Three-Molecule Interactions. J Phys Chem B 2008; 112:8311-8. [DOI: 10.1021/jp8009468] [Citation(s) in RCA: 94] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- R. Kumar
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - J. L. Skinner
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| |
Collapse
|