51
|
Seeger ZL, Izgorodina EI. A Systematic Study of DFT Performance for Geometry Optimizations of Ionic Liquid Clusters. J Chem Theory Comput 2020; 16:6735-6753. [PMID: 32865998 DOI: 10.1021/acs.jctc.0c00549] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Clusters of two ion pairs of imidazolium-based ionic liquids were optimized with 43 different levels of theory, including DFT functionals and MP2-based methods combined with varying Dunning's basis sets, and added dispersion corrections. Better preforming DFT functionals were then applied to clusters consisting of four ion pairs. Excellent performance of some DFT functionals for the two ion pair clusters did not always match that of the four ion-paired clusters despite interionic distances remaining constant between the optimized two and four ion-paired clusters of the same ionic liquid. Combinations of DFT functional and basis set such as ωB97X-D/cc-pVDZ, M06-2X/aug-cc-pVDZ, B3LYP-D3/cc-pVTZ, and TPSS-D3/cc-pVTZ gave excellent results for geometry optimization of two ion-paired clusters of imidazolium ionic liquids but gave larger deviations when applied to the four ion-paired clusters of varying ionic liquids. Empirical dispersion corrections were seen to be crucial in correctly capturing correlation effects in the studied ionic liquid clusters, becoming more important in larger clusters. Dunning's double-ζ basis set, cc-pVDZ, is associated with the smallest root mean squared deviations for geometries; however, it also produces the largest deviations in total electronic energies. ωB97X-D and M06-2X produced the best performance with the augmented version of this basis set. The triple-ζ basis set, cc-pVTZ, leads to the best performance of most of the DFT functionals (especially the dispersion-corrected ones) used, whereas its augmented version, aug-cc-pVTZ, was not seen to improve results. The combinations of functional and basis set that gave the best geometry and energetics in both two and four ion-paired clusters were PBE-D3/cc-pVTZ, ωB97X-D/aug-cc-pVDZ, and BLYP-D3/cc-pVTZ. All three combinations are recommended for geometry optimizations of larger clusters of ionic liquids. PBE-D3/cc-pVTZ performed the best with an average deviation of 2.3 kJ mol-1 and a standard deviation of 3.4 kJ mol-1 for total electronic energy when applied to four ion-paired clusters. Geometries optimized with FMO2-SRS-MP2/cc-pVTZ produced total energy within 2.0 kJ mol-1 off the benchmark in two ion-paired clusters, with the cc-pVDZ basis set performing unsurprisingly poorly with the same method. The error increased to 4.8 kJ mol-1 on average in four ion-paired clusters, with the smallest RMSD deviations in geometries when compared to the benchmark ones. This study is the first report that investigated the performance of DFT functionals for two and four ion-paired clusters of a wide range of ionic liquids consisting of commonly used cations such as pyrrolidinium, imidazolium, pyridinium, and ammonium. It also identified the importance of assessing the performance of quantum chemical methods for ionic liquids on a variety of cation-anion combinations.
Collapse
Affiliation(s)
- Zoe L Seeger
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, Victoria 3800, Australia
| | - Ekaterina I Izgorodina
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, Victoria 3800, Australia
| |
Collapse
|
52
|
Philips A, Autschbach J. Quadrupolar NMR Relaxation of Aqueous 127I -, 131Xe, and 133Cs +: A First-Principles Approach from Dynamics to Properties. J Chem Theory Comput 2020; 16:5835-5844. [PMID: 32786904 DOI: 10.1021/acs.jctc.0c00581] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quadrupolar NMR relaxation rates were computed for aqueous 133Cs+, 131Xe, and 127I- via Kohn-Sham (KS) density functional theory-based ab initio molecular dynamics and KS calculations of the electric field gradient (EFG) tensors along the trajectories. The resulting rates are within a factor of 1-3 of the experimental values and can be compared to available results from classical dynamics and EFGs from electrostatic models with corrections via Sternheimer antishielding factors. Relativistic effects are shown to have an enhancing effect on the magnitude of the EFGs. An analysis of the EFGs was carried out in terms of localized molecular orbitals to elucidate contributions from the solvent versus solute polarization and assess the validity of the Sternheimer approximation for these systems.
Collapse
Affiliation(s)
- Adam Philips
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
53
|
Borah S. Spatially resolved hydration shells and dynamics of different sulfur species in water from first-principle molecular dynamics simulations. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113387] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
54
|
Yao Y, Kanai Y. Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional. J Chem Phys 2020; 153:044114. [DOI: 10.1063/5.0012815] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
55
|
Ab Initio Molecular Dynamics Study of Methanol-Water Mixtures under External Electric Fields. Molecules 2020; 25:molecules25153371. [PMID: 32722281 PMCID: PMC7435743 DOI: 10.3390/molecules25153371] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 07/17/2020] [Accepted: 07/20/2020] [Indexed: 11/25/2022] Open
Abstract
Intense electric fields applied on H-bonded systems are able to induce molecular dissociations, proton transfers, and complex chemical reactions. Nevertheless, the effects induced in heterogeneous molecular systems such as methanol-water mixtures are still elusive. Here we report on a series of state-of-the-art ab initio molecular dynamics simulations of liquid methanol-water mixtures at different molar ratios exposed to static electric fields. If, on the one hand, the presence of water increases the proton conductivity of methanol-water mixtures, on the other, it hinders the typical enhancement of the chemical reactivity induced by electric fields. In particular, a sudden increase of the protonic conductivity is recorded when the amount of water exceeds that of methanol in the mixtures, suggesting that important structural changes of the H-bond network occur. By contrast, the field-induced multifaceted chemistry leading to the synthesis of e.g., hydrogen, dimethyl ether, formaldehyde, and methane observed in neat methanol, in 75:25, and equimolar methanol-water mixtures, completely disappears in samples containing an excess of water and in pure water. The presence of water strongly inhibits the chemical reactivity of methanol.
Collapse
|
56
|
Hirakawa T, Bowler DR, Miyazaki T, Morikawa Y, Truflandier LA. Blue moon ensemble simulation of aquation free energy profiles applied to mono and bifunctional platinum anticancer drugs. J Comput Chem 2020; 41:1973-1984. [DOI: 10.1002/jcc.26367] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 05/21/2020] [Accepted: 06/01/2020] [Indexed: 12/31/2022]
Affiliation(s)
- Teruo Hirakawa
- Department of Precision EngineeringGraduate School of Engineering, Osaka University Suita Osaka Japan
- Institut des Sciences Moléculaires (ISM), Université Bordeaux Talence Cedex France
| | - David R. Bowler
- Department of Physics & AstronomyUniversity College London (UCL) London United Kingdom
- London Centre for Nanotechnology, UCL London United Kingdom
- International Centre for Materials Nanoarchitechtonics (WPI‐MANA), National Institute for Materials Science (NIMS) Tsukuba Ibaraki Japan
| | - Tsuyoshi Miyazaki
- International Centre for Materials Nanoarchitechtonics (WPI‐MANA), National Institute for Materials Science (NIMS) Tsukuba Ibaraki Japan
| | - Yoshitada Morikawa
- Department of Precision EngineeringGraduate School of Engineering, Osaka University Suita Osaka Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University Kyoto Japan
- Research Center for Ultra‐Precision Science and TechnologyGraduate School of Engineering, Osaka University Suita Osaka Japan
| | - Lionel A. Truflandier
- Department of Precision EngineeringGraduate School of Engineering, Osaka University Suita Osaka Japan
- Institut des Sciences Moléculaires (ISM), Université Bordeaux Talence Cedex France
| |
Collapse
|
57
|
Borah S. Hydration Properties of H nPO 4n-3 (n = 0-3) From Ab Initio Molecular Dynamics Simulations. J Phys Chem B 2020; 124:5454-5464. [PMID: 32484352 DOI: 10.1021/acs.jpcb.0c01769] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
For a comprehensive and detailed microscopic understanding of the hydration properties of primary aqueous phosphorus species of valence states V (viz., H3PO4, H2PO4-, HPO42-, and PO43-), a series of extensive ab initio molecular dynamics simulations is conducted at ambient temperature. In each of these cases, the spatially resolved, three-dimensional hydration shells are computed, allowing for a direct microscopic visual understanding of the hydration shells around the species. Since these species are excellent agents for the formation of hydrogen bonds (H-bonds) in water, which determine a wide range of their structural, dynamic, and spectroscopic features, a detailed analysis of the qualitative and quantitative aspects of the H-bonds, including their lifetime calculations, is performed. Vibrational density of states (VDOS) is calculated for each of the species in solute phases, resolved for each H-bonding site, and compared against the gas-phase normal modes of H3PO4 for the purpose of understanding the signatures of the peaks in VDOS plots and, in particular, the effects of solvation and H-bonding mechanisms. The results are well in line with available experimental data and other recent computer-aided studies in the literature.
Collapse
Affiliation(s)
- Sangkha Borah
- Okinawa Institute of Science and Technology Graduate University, 1919-1, Okinawa 904-0412, Japan.,Department of Physics, Indian Institute of Technology Guwahati, Guwahati, Assam 781039, India
| |
Collapse
|
58
|
The H2S dimer revisited – Insights from wave-function and density functional theory methods. Ab initio molecular dynamics simulations of liquid H2S. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112821] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
59
|
Heenen HH, Gauthier JA, Kristoffersen HH, Ludwig T, Chan K. Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches. J Chem Phys 2020; 152:144703. [DOI: 10.1063/1.5144912] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Hendrik H. Heenen
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Joseph A. Gauthier
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | | | - Thomas Ludwig
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Karen Chan
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
60
|
Dettori R, Donadio D. Carbon dioxide, bicarbonate and carbonate ions in aqueous solutions under deep Earth conditions. Phys Chem Chem Phys 2020; 22:10717-10725. [PMID: 32103223 DOI: 10.1039/c9cp06904f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
We investigate the effect of pressure, temperature and acidity on the composition of water-rich carbon-bearing fluids under thermodynamic conditions that correspond to the Earth's deep crust and upper mantle. Our first-principles molecular dynamics simulations provide mechanistic insight into the hydration shell of carbon dioxide, bicarbonate and carbonate ions, and into the pathways of the acid/base reactions that convert these carbon species into one another in aqueous solutions. At temperatures of 1000 K and higher, our simulations can sample the chemical equilibrium of these acid/base reactions, thus allowing us to estimate the chemical composition of diluted carbon dioxide and (bi)carbonate ions as a function of acidity and thermodynamic conditions. We find that, especially at the highest temperature, the acidity of the solution is essential to determine the stability domain of CO2vs. HCO3-vs. CO32-.
Collapse
Affiliation(s)
- Riccardo Dettori
- Department of Chemistry, University of California Davis, One Shields Avenue, Davis, California 95616, USA.
| | | |
Collapse
|
61
|
Zhong K, Yu CC, Dodia M, Bonn M, Nagata Y, Ohto T. Vibrational mode frequency correction of liquid water in density functional theory molecular dynamics simulations with van der Waals correction. Phys Chem Chem Phys 2020; 22:12785-12793. [DOI: 10.1039/c9cp06335h] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We develop a frequency correction scheme for the stretch and bending modes of liquid water, which substantially improves the prediction of the vibrational spectra.
Collapse
Affiliation(s)
- Kai Zhong
- Hefei National Laboratory for Physical Sciences at the Microscale
- Collaborative Innovation Center of Chemistry for Energy Materials
- CAS Center for Excellence in Nanoscience
- School of Chemistry and Materials Science
- University of Science and Technology of China
| | - Chun-Chieh Yu
- Max Planck Institute for Polymer Research
- Ackermannweg 10
- 55128 Mainz
- Germany
| | - Mayank Dodia
- Max Planck Institute for Polymer Research
- Ackermannweg 10
- 55128 Mainz
- Germany
| | - Mischa Bonn
- Max Planck Institute for Polymer Research
- Ackermannweg 10
- 55128 Mainz
- Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research
- Ackermannweg 10
- 55128 Mainz
- Germany
| | - Tatsuhiko Ohto
- Graduate School of Engineering Science
- Osaka University
- Osaka 560-8531
- Japan
| |
Collapse
|
62
|
Caratelli C, Hajek J, Meijer EJ, Waroquier M, Van Speybroeck V. Dynamic Interplay between Defective UiO-66 and Protic Solvents in Activated Processes. Chemistry 2019; 25:15315-15325. [PMID: 31461187 PMCID: PMC6916623 DOI: 10.1002/chem.201903178] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2019] [Indexed: 01/24/2023]
Abstract
UiO‐66, composed by Zr‐oxide inorganic bricks [Zr6(μ3‐O)4(μ3‐OH)4] and organic terephthalate linkers, is one of the most studied metal–organic frameworks (MOFs) due to its exceptional thermal, chemical, and mechanical stability. Thanks to its high connectivity, the material can withstand structural deformations during activation processes such as linker exchange, dehydration, and defect formation. These processes do alter the zirconium coordination number in a dynamic way, creating open metal sites for catalysis and thus are able to tune the catalytic properties. In this work, it is shown, by means of first‐principle molecular‐dynamics simulations at operating conditions, how protic solvents may facilitate such changes in the metal coordination. Solvent can induce structural rearrangements in the material that can lead to undercoordinated but also overcoordinated metal sites. This is demonstrated by simulating activation processes along well‐chosen collective variables. Such enhanced MD simulations are able to track the intrinsic dynamics of the framework at realistic conditions.
Collapse
Affiliation(s)
- Chiara Caratelli
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Julianna Hajek
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Evert Jan Meijer
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098, XH, Amsterdam, The Netherlands
| | - Michel Waroquier
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | | |
Collapse
|
63
|
Ko HY, Zhang L, Santra B, Wang H, E W, DiStasio Jr RA, Car R. Isotope effects in liquid water via deep potential molecular dynamics. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1652366] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ, USA
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, PA, USA
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing, People's Republic of China
| | - Weinan E
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Department of Mathematics, Princeton University, Princeton, NJ, USA
| | | | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ, USA
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Department of Physics and Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, NJ, USA
| |
Collapse
|
64
|
Cassone G, Sponer J, Trusso S, Saija F. Ab initio spectroscopy of water under electric fields. Phys Chem Chem Phys 2019; 21:21205-21212. [PMID: 31368466 DOI: 10.1039/c9cp03101d] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Whereas a broad range of literature exists on the spectroscopy of water in disparate conditions, infrared (IR) and Raman spectra of water subjected to electric fields have never extensively been investigated so far. Based on ab initio molecular dynamics simulations, here we present IR and Raman spectra of bulk liquid water under the effect of static electric fields. A contraction of the entire frequency range is recorded upon increasing the field intensity both in the IR and in the Raman spectra. Whilst the OH stretching band is progressively shifted toward lower frequencies - indicating a field-induced strengthening of the H-bond network - all the other bands are up-shifted by the field. Furthermore, an evident modification of the librational mode band appears in all the spectra. Finally, the order-maker action of the field emerges also from the increase of the water orientational tetrahedral order. Upon field exposure, the water structure becomes more "ice like".
Collapse
Affiliation(s)
- Giuseppe Cassone
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61265 Brno, Czech Republic.
| | | | | | | |
Collapse
|
65
|
Apostolidou C. OH radical in water from ab initio molecular dynamics simulation employing hybrid functionals. J Chem Phys 2019. [DOI: 10.1063/1.5107479] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Christina Apostolidou
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
66
|
Daub CD, Halonen L. Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide Salt on the Deprotonation of Formic Acid in Aqueous Solution. J Phys Chem B 2019; 123:6823-6829. [PMID: 31310529 PMCID: PMC6750841 DOI: 10.1021/acs.jpcb.9b04618] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
The deprotonation of formic acid is investigated using metadynamics in tandem with Born-Oppenheimer molecular dynamics simulations. We compare our findings for formic acid in pure water with previous studies before examining formic acid in aqueous solutions of lithium bromide. We carefully consider different definitions for the collective variable(s) used to drive the metadynamics, emphasizing that the variables used must include all of the possible reactive atoms in the system, in this case carboxylate oxygens and water hydrogens. This ensures that all the various possible proton exchange events can be accommodated and the collective variable(s) can distinguish the protonated and deprotonated states, even over rather long ab initio simulation runs (ca. 200-300 ps). Our findings show that the formic acid deprotonation barrier and the free energy of the deprotonated state are higher in concentrated lithium bromide, in agreement with the available experimental data for acids in salt solution. We show that the presence of Br- in proximity to the formic acid hydroxyl group effectively inhibits deprotonation. Our study extends previous work on acid deprotonation in pure water and at air-water interfaces to more complex multicomponent systems of importance in atmospheric and marine chemistry.
Collapse
Affiliation(s)
- Christopher D Daub
- Department of Chemistry , University of Helsinki , P.O. Box 55, Helsinki FIN-00014 , Finland
| | - Lauri Halonen
- Department of Chemistry , University of Helsinki , P.O. Box 55, Helsinki FIN-00014 , Finland
| |
Collapse
|
67
|
Choudhary A, Chandra A. An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute. J Chem Phys 2019; 151:044508. [PMID: 31370512 DOI: 10.1063/1.5094570] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.
Collapse
Affiliation(s)
- Ashu Choudhary
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
68
|
Dodia M, Ohto T, Imoto S, Nagata Y. Structure and Dynamics of Water at the Water-Air Interface Using First-Principles Molecular Dynamics Simulations. II. NonLocal vs Empirical van der Waals Corrections. J Chem Theory Comput 2019; 15:3836-3843. [PMID: 31074989 PMCID: PMC6750744 DOI: 10.1021/acs.jctc.9b00253] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
van der Waals (vdW) correction schemes
have been recognized to
be essential for an accurate description of liquid water in first-principles
molecular dynamics simulation. The description of the structure and
dynamics of water is governed by the type of the vdW corrections.
So far, two vdW correction schemes have been often used: empirical
vdW corrections and nonlocal vdW corrections. In this paper, we assess
the influence of the empirical vs nonlocal vdW correction schemes
on the structure and dynamics of water at the water–air interface.
Since the structure of water at the water–air interface is
established by a delicate balance of hydrogen bond formation and breaking,
the simulation at the water–air interface provides a unique
platform to testify as to the heterogeneous interaction of water.
We used the metrics [Ohto et al. J. Chem. Theory Comput., 2019, 15, 595−60230468702] which
are directly connected with the sum-frequency generation spectroscopic
measurement. We find that the overall performance of nonlocal vdW
methods is either similar or worse compared to the empirical vdW methods.
We also investigated the performance of the optB88-DRSLL functional,
which showed slightly less accuracy than the revPBE-D3 method. We
conclude that the revPBE-D3 method shows the best performance for
describing the interfacial water.
Collapse
Affiliation(s)
- Mayank Dodia
- Max Planck Institute for Polymer Research, Ackermannweg 10 , 55128 Mainz , Germany
| | - Tatsuhiko Ohto
- Graduate School of Engineering Science , Osaka University , 1-3 Machikaneyama , Toyonaka, Osaka 560-8531 , Japan
| | - Sho Imoto
- Max Planck Institute for Polymer Research, Ackermannweg 10 , 55128 Mainz , Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10 , 55128 Mainz , Germany
| |
Collapse
|
69
|
Biswas S, Mallik BS. Heterogeneous Occupancy and Vibrational Dynamics of Spatially Patterned Water Molecules. J Phys Chem B 2019; 123:4278-4290. [PMID: 31018092 DOI: 10.1021/acs.jpcb.9b00271] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We performed first-principles molecular dynamics simulations of relatively dilute aqueous solutions of sulfate and thiosulfate dianions to analyze the structure, dynamics, and vibrational spectral properties of water molecules around the solute, especially the spatially patterned solvent molecules in the first solvation layer and the extended layers. This study also involves the investigation of dynamics of dangling OH groups in these layers and their role in patterning the water molecules around the dianions. Structural evaluation of the systems is carried out by radial distribution functions, number integrals, and spatial distribution functions. The lifetime of dangling OH groups inside the solvation shell is compared more to that of the bulk. By constructing the O-H groups in three ensembles (S1, S2, and S3) around the anion, we show that the frequency distribution of OH modes in the S1 ensemble show red-shifting for both sulfate and thiosulfate. The O-H groups in the S2 ensemble of the sulfate-water system show red-shifting by 10 cm-1, while in the case of thiosulfate-water, these O-H groups show blue-shifting by 8 cm-1. The water molecules in S1 and S2 subensembles have slower dynamics compared to those in the bulk (S3). The dynamics of various kinds of hydrogen bonds were characterized by hydrogen bond population correlation functions. The spectral diffusion of solvation shell O-H modes was performed through a frequency-time correlation function. We find a significant amount of orientational retardation of water molecules in the S1 layer and moderate retardation in the S2 layer as compared to that in the bulk, S3 layer. All these findings, the red shift of the OH stretching frequency in S1 and S2 layers, slowing down of the orientational dynamics of OH vectors in S1 and S2 layers, and less diffusivity of water in S1 and S2 layers, show the long-range kosmotropic effect of multivalent sulfate and thiosulfate oxyanions. Due to the long-range effect, heterogeneous occupancy of water molecules is observed, and the water molecules are found to arrange in a patterned manner in the vicinity of anions with varied local density.
Collapse
Affiliation(s)
- Sohag Biswas
- Department of Chemistry , Indian Institute of Technology Hyderabad , Kandi, Sangareddy 502285 , Telangana , India
| | - Bhabani S Mallik
- Department of Chemistry , Indian Institute of Technology Hyderabad , Kandi, Sangareddy 502285 , Telangana , India
| |
Collapse
|
70
|
DelloStritto MJ, Piontek SM, Klein ML, Borguet E. Effect of Functional and Electron Correlation on the Structure and Spectroscopy of the Al 2O 3(001)-H 2O Interface. J Phys Chem Lett 2019; 10:2031-2036. [PMID: 30977654 DOI: 10.1021/acs.jpclett.9b00016] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Oxide-water interfaces are ubiquitous, with many applications in industry and the environment, yet there is a great deal of controversy over their properties and microscopic structure. This controversy stems, in part, from the unique H-bond networks formed at different surface terminations and mineral compositions. Density functional theory simulations of these interfaces require an accurate description of both the oxide mineral and water in diverse H-bond environments. Thus, herein we simulate the Al2O3(001)-H2O interface using the PBE, PBE-TS, RPBE, SCAN, and HSE06-TS functionals to determine how calculated interfacial properties depend on the choice of functional. We find that the structure of the first few layers of water at the surface is determined by electron correlation in a way that cannot be approximated using semiemipirical van der Waals corrections. Of the functionals investigated, we find that SCAN yields the most accurate interfacial structure, dynamics, and sum frequency generation spectrum. Furthermore, SCAN leads to a reduction in the order of the 2D H-bond network of water at the alumina surface predicted by GGA functionals, leading to a significant decrease in the anisotropy of the diffusion coefficient at the surface. We emphasize the importance of using a functional which accurately describes electron correlation for more complex oxides, such as transition-metal oxides, where electron correlation may play an even greater role in determining the structure and dynamics of the oxide-water interface.
Collapse
|
71
|
Zhou L, Xu J, Xu L, Wu X. Importance of van der Waals effects on the hydration of metal ions from the Hofmeister series. J Chem Phys 2019; 150:124505. [PMID: 30927898 DOI: 10.1063/1.5086939] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The van der Waals (vdW) interaction plays a crucial role in the description of liquid water. Based on ab initio molecular dynamics simulations, including the non-local and fully self-consistent density-dependent implementation of the Tkatchenko-Scheffler dispersion correction, we systematically studied the aqueous solutions of metal ions (K+, Na+, and Ca2+) from the Hofmeister series. Similar to liquid water, the vdW interactions strengthen the attractions among water molecules in the long-range, leading to the hydrogen bond networks softened in all the ion solutions. However, the degree that the hydration structure is revised by the vdW interactions is distinct for different ions, depending on the strength of short-range interactions between the hydrated ion and surrounding water molecules. Such revisions by the vdW interactions are important for the understanding of biological functionalities of ion channels.
Collapse
Affiliation(s)
- Liying Zhou
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Jianhang Xu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Limei Xu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
72
|
Biswas S, Mallik BS. Vibration Spectral Dynamics of Weakly Coordinating Water Molecules near an Anion: FPMD Simulations of an Aqueous Solution of Tetrafluoroborate. J Phys Chem B 2019; 123:2135-2146. [PMID: 30759344 DOI: 10.1021/acs.jpcb.9b00069] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The extent to which the ions affect the nearby water molecules will decide the structure-making or breaking nature of those ions in aqueous solutions. The effects of a weakly coordinating anion on the structure, dynamics, and vibrational properties of water molecules are not so significant as compared to an anion capable of making strong ion-water hydrogen bonds. The present work deals with the first-principles molecular dynamics study of an aqueous solution of such a weakly coordinating anion, tetrafluoroborate (BF4-), using dispersion-corrected DFT-based first-principles molecular dynamics (FPMD) simulations. Various structural, dynamical, and spectral properties, such as radial distribution functions (RDFs), rotational dynamics, vibrational density of states (VDOS), hydrogen bond as well as dangling OH autocorrelation functions, and residence dynamics, were calculated to investigate the effects of the anion on nearby water molecules. The process of spectral diffusion was assessed through a time series wavelet transformation of trajectories obtained from FPMD simulations. The first ion-water solvation shell extends up to 5.5 Å, containing around 20 water molecules. The lifetime of the ion-water hydrogen bond is found to be 1.19 ps, whereas the water-water hydrogen bond lifetime is found to be 1.13 ps. Inside the solvation shell, the persistence time of dangling OH chromophores and the average frequency of OH modes inside the solvation shell are found to be more compared to bulk. Three time scales are found for solvation shell OH modes from the frequency-frequency correlation function. A very short time scale is found for the intact ion-water interaction; the short time scale is for the ion-water hydrogen bond, and the long time scale is for escape dynamics of water molecules from the ion solvation shell. From the mean squared displacement, it is found that solvation water molecules diffuse slower than the bulk. However, solvation shell water molecules show faster relaxation from the analysis of rotational anisotropy. Within the longer time scale of spectral diffusion, this process (which is related to various dynamics of the molecules) is not yet complete, as compared to fast anisotropic decay. This fact is similar to the experimental finding of spectral diffusion and anisotropy time scales in the aqueous solution of borohydride anion. The calculated results are also compared with available experimental data wherever possible.
Collapse
Affiliation(s)
- Sohag Biswas
- Department of Chemistry , Indian Institute of Technology Hyderabad , Kandi, Sangareddy , 502 285 Telangana , India
| | - Bhabani S Mallik
- Department of Chemistry , Indian Institute of Technology Hyderabad , Kandi, Sangareddy , 502 285 Telangana , India
| |
Collapse
|
73
|
Dziedzic J, Head-Gordon T, Head-Gordon M, Skylaris CK. Mutually polarizable QM/MM model with in situ optimized localized basis functions. J Chem Phys 2019; 150:074103. [PMID: 30795653 DOI: 10.1063/1.5080384] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
We extend our recently developed quantum-mechanical/molecular mechanics (QM/MM) approach [Dziedzic et al., J. Chem. Phys. 145, 124106 (2016)] to enable in situ optimization of the localized orbitals. The quantum subsystem is described with onetep linear-scaling density functional theory and the classical subsystem - with the AMOEBA polarizable force field. The two subsystems interact via multipolar electrostatics and are fully mutually polarizable. A total energy minimization scheme is employed for the Hamiltonian of the coupled QM/MM system. We demonstrate that, compared to simpler models using fixed basis sets, the additional flexibility offered by in situ optimized basis functions improves the accuracy of the QM/MM interface, but also poses new challenges, making the QM subsystem more prone to overpolarization and unphysical charge transfer due to increased charge penetration. We show how these issues can be efficiently solved by replacing the classical repulsive van der Waals term for QM/MM interactions with an interaction of the electronic density with a fixed, repulsive MM potential that mimics Pauli repulsion, together with a modest increase in the damping of QM/MM polarization. We validate our method, with particular attention paid to the hydrogen bond, in tests on water-ion pairs, the water dimer, first solvation shells of neutral and charged species, and solute-solvent interaction energies. As a proof of principle, we determine suitable repulsive potential parameters for water, K+, and Cl-. The mechanisms we employed to counteract the unphysical overpolarization of the QM subsystem are demonstrated to be adequate, and our approach is robust. We find that the inclusion of explicit polarization in the MM part of QM/MM improves agreement with fully QM calculations. Our model permits the use of minimal size QM regions and, remarkably, yields good energetics across the well-balanced QM/MM interface.
Collapse
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Teresa Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
74
|
Subasinghege Don V, David R, Du P, Milet A, Kumar R. Interfacial Water at Graphene Oxide Surface: Ordered or Disordered? J Phys Chem B 2019; 123:1636-1649. [DOI: 10.1021/acs.jpcb.8b10987] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Visal Subasinghege Don
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Rolf David
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
- Université Grenoble Alpes, CNRS, DCM, 38000 Grenoble, France
| | - Pu Du
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Anne Milet
- Université Grenoble Alpes, CNRS, DCM, 38000 Grenoble, France
| | - Revati Kumar
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
75
|
Arnold C, Inhester L, Carbajo S, Welsch R, Santra R. Simulated XUV photoelectron spectra of THz-pumped liquid water. J Chem Phys 2019; 150:044505. [PMID: 30709301 DOI: 10.1063/1.5054272] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Highly intense, sub-picosecond terahertz (THz) pulses can be used to induce ultrafast temperature jumps (T-jumps) in liquid water. A supercritical state of gas-like water with liquid density is established, and the accompanying structural changes are expected to give rise to time-dependent chemical shifts. We investigate the possibility of using extreme ultraviolet photoelectron spectroscopy as a probe for ultrafast dynamics induced by sub-picosecond THz pulses of varying intensities and frequencies. To this end, we use ab initio methods to calculate photoionization cross sections and photoelectron energies of (H2O)20 clusters embedded in an aqueous environment represented by point charges. The cluster geometries are sampled from ab initio molecular dynamics simulations modeling the THz-water interactions. We find that the peaks in the valence photoelectron spectrum are shifted by up to 0.4 eV after the pump pulse and that they are broadened with respect to unheated water. The shifts can be connected to structural changes caused by the heating, but due to saturation effects they are not sensitive enough to serve as a thermometer for T-jumped water.
Collapse
Affiliation(s)
- Caroline Arnold
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Ludger Inhester
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Sergio Carbajo
- SLAC National Accelerator Laboratory, Stanford University, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Robin Santra
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
76
|
Daub CD, Hänninen V, Halonen L. Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide on the Structure of the Aqueous Solution-Air Interface. J Phys Chem B 2019; 123:729-737. [PMID: 30605330 PMCID: PMC6727360 DOI: 10.1021/acs.jpcb.8b10552] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We
present the results of ab initio molecular dynamics simulations
of the solution–air interface of aqueous lithium bromide (LiBr).
We find that, in agreement with the experimental data and previous
simulation results with empirical polarizable force field models,
Br– anions prefer to accumulate just below the first
molecular water layer near the interface, whereas Li+ cations
remain deeply buried several molecular layers from the interface,
even at very high concentration. The separation of ions has a profound
effect on the average orientation of water molecules in the vicinity
of the interface. We also find that the hydration number of Li+ cations in the center of the slab Nc,Li+–H2O ≈ 4.7 ±
0.3, regardless of the salt concentration. This estimate is consistent
with the recent experimental neutron scattering data, confirming that
results from nonpolarizable empirical models, which consistently predict
tetrahedral coordination of Li+ to four solvent molecules,
are incorrect. Consequently, disruption of the hydrogen bond network
caused by Li+ may be overestimated in nonpolarizable empirical
models. Overall, our results suggest that empirical models, in particular
nonpolarizable models, may not capture all of the properties of the
solution–air interface necessary to fully understand the interfacial
chemistry.
Collapse
Affiliation(s)
- Christopher D Daub
- Department of Chemistry , University of Helsinki , P.O. Box 55, Helsinki 00014 , Finland
| | - Vesa Hänninen
- Department of Chemistry , University of Helsinki , P.O. Box 55, Helsinki 00014 , Finland
| | - Lauri Halonen
- Department of Chemistry , University of Helsinki , P.O. Box 55, Helsinki 00014 , Finland
| |
Collapse
|
77
|
Philips A, Autschbach J. Proton NMR relaxation from molecular dynamics: intramolecular and intermolecular contributions in water and acetonitrile. Phys Chem Chem Phys 2019; 21:26621-26629. [DOI: 10.1039/c9cp04976b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Accurate 1H NMR relaxation rates for protons in pure water and acetonitrile are computed via ab initio and force field molecular dynamics. Dipole–dipole and spin-rotation mechanisms are considered.
Collapse
Affiliation(s)
- Adam Philips
- Department of Chemistry
- University at Buffalo
- State University of New York
- Buffalo
- USA
| | - Jochen Autschbach
- Department of Chemistry
- University at Buffalo
- State University of New York
- Buffalo
- USA
| |
Collapse
|
78
|
Milne AW, Jorge M. Polarization Corrections and the Hydration Free Energy of Water. J Chem Theory Comput 2018; 15:1065-1078. [DOI: 10.1021/acs.jctc.8b01115] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Andrew W. Milne
- Department of Chemical and Process Engineering, University of Strathclyde, Glasgow G1 1XQ, United Kingdom
| | - Miguel Jorge
- Department of Chemical and Process Engineering, University of Strathclyde, Glasgow G1 1XQ, United Kingdom
| |
Collapse
|
79
|
Philips A, Marchenko A, Ducati LC, Autschbach J. Quadrupolar 14N NMR Relaxation from Force-Field and Ab Initio Molecular Dynamics in Different Solvents. J Chem Theory Comput 2018; 15:509-519. [PMID: 30462503 DOI: 10.1021/acs.jctc.8b00807] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quadrupolar NMR spin relaxation rates and corresponding line widths were computed for the quadrupolar nucleus 14N for neat acetonitrile as well as for 1-methyl-1,3-imidazole and 1-methyl-1,3,4-triazole in different solvents. Molecular dynamics (MD) was performed with forces from the Kohn-Sham (KS) theory (ab initio MD) and force-field molecular mechanics (classical MD), followed by KS electric field gradient (EFG) calculations. For acetonitrile the agreement of the 14N line width with experiment is very good. Relative line widths for the azole nitrogens are improved over simpler approximations used previously in conjunction with single-point calculations at the multiconfigurational self-consistent field level. Overall, the NMR line widths are computed within a factor of 2 of the experimental values, giving access to reasonable estimates both of the dynamic EFG variance in the solvated systems as well as the associated correlation times that determine the relaxation rates.
Collapse
Affiliation(s)
- Adam Philips
- Department of Chemistry , University at Buffalo, State University of New York , Buffalo , New York 14260-3000 , United States
| | - Alex Marchenko
- Department of Chemistry , University at Buffalo, State University of New York , Buffalo , New York 14260-3000 , United States
| | - Lucas C Ducati
- Department of Fundamental Chemistry Institute of Chemistry , University of São Paulo , Av. Prof. Lineu Prestes 748 , São Paulo , SP 05508-000 , Brazil
| | - Jochen Autschbach
- Department of Chemistry , University at Buffalo, State University of New York , Buffalo , New York 14260-3000 , United States
| |
Collapse
|
80
|
Ohto T, Dodia M, Imoto S, Nagata Y. Structure and Dynamics of Water at the Water–Air Interface Using First-Principles Molecular Dynamics Simulations within Generalized Gradient Approximation. J Chem Theory Comput 2018; 15:595-602. [DOI: 10.1021/acs.jctc.8b00567] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tatsuhiko Ohto
- Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Mayank Dodia
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Sho Imoto
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
81
|
Xie WJ, Cha S, Ohto T, Mizukami W, Mao Y, Wagner M, Bonn M, Hunger J, Nagata Y. Large Hydrogen-Bond Mismatch between TMAO and Urea Promotes Their Hydrophobic Association. Chem 2018. [DOI: 10.1016/j.chempr.2018.08.020] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
82
|
König G, Pickard FC, Huang J, Thiel W, MacKerell AD, Brooks BR, York DM. A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules 2018; 23:E2695. [PMID: 30347691 PMCID: PMC6222909 DOI: 10.3390/molecules23102695] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 12/01/2022] Open
Abstract
Maintaining a proper balance between specific intermolecular interactions and non-specific solvent interactions is of critical importance in molecular simulations, especially when predicting binding affinities or reaction rates in the condensed phase. The most rigorous metric for characterizing solvent affinity are solvation free energies, which correspond to a transfer from the gas phase into solution. Due to the drastic change of the electrostatic environment during this process, it is also a stringent test of polarization response in the model. Here, we employ both the CHARMM fixed charge and polarizable force fields to predict hydration free energies of twelve simple solutes. The resulting classical ensembles are then reweighted to obtain QM/MM hydration free energies using a variety of QM methods, including MP2, Hartree⁻Fock, density functional methods (BLYP, B3LYP, M06-2X) and semi-empirical methods (OM2 and AM1 ). Our simulations test the compatibility of quantum-mechanical methods with molecular-mechanical water models and solute Lennard⁻Jones parameters. In all cases, the resulting QM/MM hydration free energies were inferior to purely classical results, with the QM/MM Drude force field predictions being only marginally better than the QM/MM fixed charge results. In addition, the QM/MM results for different quantum methods are highly divergent, with almost inverted trends for polarizable and fixed charge water models. While this does not necessarily imply deficiencies in the QM models themselves, it underscores the need to develop consistent and balanced QM/MM interactions. Both the QM and the MM component of a QM/MM simulation have to match, in order to avoid artifacts due to biased solute⁻solvent interactions. Finally, we discuss strategies to improve the convergence and efficiency of multi-scale free energy simulations by automatically adapting the molecular-mechanics force field to the target quantum method.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
- School of Life Sciences, Westlake University, 18 Shilongshan Street, Hangzhou 310024, China.
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander D MacKerell
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
| |
Collapse
|
83
|
Ruiz Pestana L, Marsalek O, Markland TE, Head-Gordon T. The Quest for Accurate Liquid Water Properties from First Principles. J Phys Chem Lett 2018; 9:5009-5016. [PMID: 30118601 DOI: 10.1021/acs.jpclett.8b02400] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of density functional theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, display accuracy comparable to a computationally expensive dispersion-corrected hybrid functional, revPBE0-D3. We show that the meta-GGA DFT functional reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared spectra with comparable accuracy of a much more expensive hybrid functional. This work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei and that high-temperature simulations are a poor proxy for properly treating NQEs.
Collapse
Affiliation(s)
| | - Ondrej Marsalek
- Charles University , Faculty of Mathematics and Physics , Ke Karlovu 3 , 121 16 Prague 2 , Czech Republic
| | - Thomas E Markland
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | | |
Collapse
|
84
|
Tsimpanogiannis IN, Moultos OA, Franco LFM, Spera MBDM, Erdős M, Economou IG. Self-diffusion coefficient of bulk and confined water: a critical review of classical molecular simulation studies. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1511903] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Luís F. M. Franco
- School of Chemical Engineering, University of Campinas, Campinas, Brazil
| | | | - Máté Erdős
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
85
|
Lesnicki D, Sulpizi M. A Microscopic Interpretation of Pump-Probe Vibrational Spectroscopy Using Ab Initio Molecular Dynamics. J Phys Chem B 2018; 122:6604-6609. [PMID: 29799755 DOI: 10.1021/acs.jpcb.8b04159] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
What happens when extra vibrational energy is added to water? Using nonequilibrium molecular dynamics simulations, also including the full electronic structure, and novel descriptors, based on projected vibrational density of states, we are able to follow the flow of excess vibrational energy from the excited stretching and bending modes. We find that the energy relaxation, mostly mediated by a stretching-stretching coupling in the first solvation shell, is highly heterogeneous and strongly depends on the local environment, where a strong hydrogen bond network can transport energy with a time scale of 200 fs, whereas a weaker network can slow down the transport by a factor 2-3.
Collapse
Affiliation(s)
- Dominika Lesnicki
- Institute of Physics , Johannes Gutenberg University Mainz , Staudingerweg 7 , 55099 Mainz , Germany
| | - Marialore Sulpizi
- Institute of Physics , Johannes Gutenberg University Mainz , Staudingerweg 7 , 55099 Mainz , Germany
| |
Collapse
|
86
|
Singraber A, Morawietz T, Behler J, Dellago C. Density anomaly of water at negative pressures from first principles. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:254005. [PMID: 29762140 DOI: 10.1088/1361-648x/aac4f4] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Using molecular dynamics simulations based on ab initio trained high-dimensional neural network potentials, we study the equation of state of liquid water at negative pressures. From density isobars computed for various pressures down to p = -230 MPa we determine the line of density maxima for two potentials based on the BLYP and the RPBE functionals, respectively. In both cases, dispersion corrections are included to account for non-local long-range correlations that give rise to van der Waals forces. We have followed the density maximum down to negative pressures close to the spinodal instability. For both functionals, the temperature of maximum density increases with decreasing pressure under moderate stretching, but changes slope at [Formula: see text] MPa and [Formula: see text] MPa for BLYP and RPBE, respectively. Our calculations confirm qualitatively the retracing shape of the line of density maxima found for empirical water models, indicating that the spinodal line maintains a positive slope even at strongly negative pressures.
Collapse
Affiliation(s)
- Andreas Singraber
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | | | | | | |
Collapse
|
87
|
Agosta L, Brandt EG, Lyubartsev AP. Diffusion and reaction pathways of water near fully hydrated TiO 2 surfaces from ab initio molecular dynamics. J Chem Phys 2018; 147:024704. [PMID: 28711052 DOI: 10.1063/1.4991381] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Ab initio molecular dynamics simulations are reported for water-embedded TiO2 surfaces to determine the diffusive and reactive behavior at full hydration. A three-domain model is developed for six surfaces [rutile (110), (100), and (001), and anatase (101), (100), and (001)] which describes waters as "hard" (irreversibly bound to the surface), "soft" (with reduced mobility but orientation freedom near the surface), or "bulk." The model explains previous experimental data and provides a detailed picture of water diffusion near TiO2 surfaces. Water reactivity is analyzed with a graph-theoretic approach that reveals a number of reaction pathways on TiO2 which occur at full hydration, in addition to direct water splitting. Hydronium (H3O+) is identified to be a key intermediate state, which facilitates water dissociation by proton hopping between intact and dissociated waters near the surfaces. These discoveries significantly improve the understanding of nanoscale water dynamics and reactivity at TiO2 interfaces under ambient conditions.
Collapse
Affiliation(s)
- Lorenzo Agosta
- Department of Materials and Environmental Chemistry, Stockholm University, 10691 Stockholm, Sweden
| | - Erik G Brandt
- Department of Materials and Environmental Chemistry, Stockholm University, 10691 Stockholm, Sweden
| | - Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, 10691 Stockholm, Sweden
| |
Collapse
|
88
|
Jiang T, Simko S, Bulo RE. Accurate Quantum Mechanics/Molecular Mechanics Simulation of Aqueous Solutions with Tailored Molecular Mechanics Models. J Chem Theory Comput 2018; 14:3943-3954. [DOI: 10.1021/acs.jctc.7b01218] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tao Jiang
- Laboratoire de chimie, École Normale Supérieure de lyon, 46 allée d’italie, 69007 Lyon, France
| | - Stanislav Simko
- Inorganic Chemistry and Catalysis Group, Debye Institute for Nanomaterials Science, Utrecht University, Universiteitsweg 99, 3584 CG, Utrecht, The Netherlands
| | - Rosa E. Bulo
- Inorganic Chemistry and Catalysis Group, Debye Institute for Nanomaterials Science, Utrecht University, Universiteitsweg 99, 3584 CG, Utrecht, The Netherlands
| |
Collapse
|
89
|
Liberatore E, Meli R, Rothlisberger U. A Versatile Multiple Time Step Scheme for Efficient ab Initio Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:2834-2842. [PMID: 29624388 DOI: 10.1021/acs.jctc.7b01189] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present here our implementation of a time-reversible, multiple time step (MTS) method for full QM and hybrid QM/MM Born-Oppenheimer molecular dynamics simulations. The method relies on a fully flexible combination of electronic structure methods, from density functional theory to wave function-based quantum chemistry methods, to evaluate the nuclear forces in the reference and in the correction steps. The possibility of combining different electronic structure methods is based on the observation that exchange and correlation terms only contribute to low frequency modes of nuclear forces. We show how a pair of low/high level electronic structure methods that individually would lead to very different system properties can be efficiently combined in the reference and correction steps of this MTS scheme. The current MTS implementation makes it possible to perform highly accurate ab initio molecular dynamics simulations at reduced computational cost. Stable and accurate trajectories were obtained with time steps of several femtoseconds, similar to and even exceeding the ones usually adopted in classical molecular dynamics, in particular when using a generalized Langevin stochastic thermostat. Compared to the standard Velocity Verlet integration, the present MTS scheme allows for a 5- to 6-fold overall speedup, at an unaltered level of accuracy.
Collapse
Affiliation(s)
- Elisa Liberatore
- École Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne , Switzerland
| | - Rocco Meli
- École Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne , Switzerland
| | - Ursula Rothlisberger
- École Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne , Switzerland
| |
Collapse
|
90
|
Karmakar A, Chandra A. Water under Supercritical Conditions: Hydrogen Bonds, Polarity, and Vibrational Frequency Fluctuations from Ab Initio Simulations with a Dispersion Corrected Density Functional. ACS OMEGA 2018; 3:3453-3462. [PMID: 31458597 PMCID: PMC6641502 DOI: 10.1021/acsomega.7b02036] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Accepted: 03/09/2018] [Indexed: 05/27/2023]
Abstract
We have studied the effects of dispersion interactions on the dynamics of vibrational frequency fluctuations, hydrogen bonds, and free OD modes in supercritical heavy water at three different densities by means of ab initio molecular dynamics simulations. The vibrational spectral diffusion, as described by the frequency fluctuations, is studied through calculations of frequency time correlation of stretch modes of deuterated water, and its relations to the dynamics of hydrogen bonds and free OD modes are described. In addition, some of the other dynamical, structural, and electronic properties such as diffusion, rotational relaxation, radial distribution functions, hydrogen bond and coordination numbers, and dipole moments are also investigated from the perspectives of their variation with inclusion of dispersion interactions at varying density of the solvent. Although some changes in the structural properties are found on inclusion of dispersion corrections, no significant difference in the fluctuation dynamics of OD stretching frequencies and also in other dynamical quantities of supercritical water are found because of dispersion effects. The dynamics of water molecules under supercritical conditions is very fast compared to the corresponding dynamics under ambient conditions. The large thermal effects at such a high temperature seem to take over any relatively minor changes that might be introduced by weak dispersion interaction.
Collapse
|
91
|
Zhang X, Kumar R, Kuroda DG. Acetate ion and its interesting solvation shell structure and dynamics. J Chem Phys 2018. [DOI: 10.1063/1.5019363] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiaoliu Zhang
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Revati Kumar
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Daniel G. Kuroda
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| |
Collapse
|
92
|
Biswas S, Chakraborty D, Mallik BS. Interstitial Voids and Resultant Density of Liquid Water: A First-Principles Molecular Dynamics Study. ACS OMEGA 2018; 3:2010-2017. [PMID: 31458509 PMCID: PMC6641304 DOI: 10.1021/acsomega.7b01996] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Accepted: 02/07/2018] [Indexed: 05/31/2023]
Abstract
Many anomalous properties of water can be explained on the basis of the coexistence of more than one density states: high-density water (HDW) and low-density water (LDW). We investigated these two phases of water molecules through first-principles molecular dynamics simulations using density functional theory (DFT) in conjunction with various van der Waals-corrected exchange and correlation functionals. Different density regions were found to exist due to the difference in short-range and long-range forces present in DFT potentials. These density regions were identified and analyzed on the basis of the distribution of molecules and voids present. We defined a local structure index to distinguish and find the probability of occurrence of these different states. HDW and LDW arise due to the presence of "interstitial water" molecules in between the first and second coordination shells. The population of interstitial water molecules is found to affect the overall dynamics of the system as they change the hydrogen bond pattern.
Collapse
Affiliation(s)
- Sohag Biswas
- Department
of Chemistry, Indian Institute of Technology
Hyderabad, Kandi, 502285 Sangareddy, Telangana, India
| | - Debashree Chakraborty
- Department
of Chemistry, National Institute of Technology
Karnataka, Surathkal, Mangalore 575025, Karnataka, India
| | - Bhabani S. Mallik
- Department
of Chemistry, Indian Institute of Technology
Hyderabad, Kandi, 502285 Sangareddy, Telangana, India
| |
Collapse
|
93
|
Ayala R, Martínez JM, R Pappalardo R, Refson K, Sánchez Marcos E. Effect of Basicity on the Hydrolysis of the Bi(III) Aqua Ion in Solution: An Ab Initio Molecular Dynamics Study. J Phys Chem A 2018; 122:1905-1915. [PMID: 29376356 DOI: 10.1021/acs.jpca.7b12402] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hydrolysis of the Bi(III) aqua ion under a range of solution conditions has been studied by means of ab initio molecular dynamics simulations. While the Bi(III) aqua ion is stable in pure water, there is an increasing degree of hydrolysis with the number of hydroxide anions in the medium. This is accompanied by a monotonic decrease of the total coordination number to an asymptotic value of ∼6, reached under extreme basicity conditions. Comparison of the simulated Bi(III) hydrolyzed species with the experimental species distribution at different degrees of basicity suggests that, at the PBE/DFT level of theory here employed, liquid water shows an overly acidic character. Predictions of theoretical EXAFS and XANES spectra were generated from the AIMD trajectories for different Bi hydrolyzed species, [Bi(HO)m(H2O)n]3-m+, m = 0-3 and n = 7-2. Comparison with available experimental spectra is presented. Spectral features joined to the degree of hydrolysis and hydration are analyzed.
Collapse
Affiliation(s)
- Regla Ayala
- Departamento Química Inorgánica, Universidad de Sevilla , 41012 Seville, Spain
| | | | | | - Keith Refson
- Royal Holloway University of London , Egham Hill, Egham, Surrey TW20 0EX, United Kingdom.,ISIS Facility, Science and Technology Facility Council , Harwell Science Campus, Didcot, OX11 0QX, United Kingdom
| | | |
Collapse
|
94
|
Liu J, He X, Zhang JZH. Structure of liquid water - a dynamical mixture of tetrahedral and 'ring-and-chain' like structures. Phys Chem Chem Phys 2018; 19:11931-11936. [PMID: 28440370 DOI: 10.1039/c7cp00667e] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
The nature of the dynamical hydrogen-bond network of liquid water under ambient conditions has challenged both experimental and theoretical researchers for decades and remains a topic of intense debate. In this work, we addressed the structural issue of the hydrogen-bond network of liquid water based on an accurate ab initio molecular dynamics simulation. The present work showed clearly that liquid water is neither accurately described by a static picture of mostly tetrahedral water molecules nor dominated by "ring-and-chain" like structures. Instead, the structure of water is a dynamical mixture of tetrahedral and 'ring-and-chain' like structures with a slight bias toward the former. On average, each water molecule forms about three hydrogen bonds with the surrounding water molecules. The present accurate ab initio molecular dynamics simulation of liquid water was made possible by using a fragment-based second-order Møller-Plesset perturbation theory (MP2) with a large basis set to treat a large body of water molecules. This level of ab initio theory is sufficiently accurate for describing water interactions, and the simulated structural and dynamical properties of liquid water, including radial distribution functions, diffusion coefficient, dipole moment, etc., are uniformly in excellent agreement with experimental observations.
Collapse
Affiliation(s)
- Jinfeng Liu
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.
| | | | | |
Collapse
|
95
|
Yao Y, Kanai Y. Free Energy Profile of NaCl in Water: First-Principles Molecular Dynamics with SCAN and ωB97X-V Exchange–Correlation Functionals. J Chem Theory Comput 2018; 14:884-893. [DOI: 10.1021/acs.jctc.7b00846] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
96
|
Wang H, Xu H, Jia W, Ren S. Functionalized carbon black nanoparticles used for separation of emulsified oil from oily wastewater. J DISPER SCI TECHNOL 2018. [DOI: 10.1080/01932691.2017.1298040] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Huanjiang Wang
- State Key Laboratory of Solid Lubrication, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou, P. R. China
- University of Chinese Academy of Sciences, Beijing, P. R. China
| | - Haiyan Xu
- State Key Laboratory of Solid Lubrication, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou, P. R. China
- University of Chinese Academy of Sciences, Beijing, P. R. China
| | - Weihong Jia
- State Key Laboratory of Solid Lubrication, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou, P. R. China
- University of Chinese Academy of Sciences, Beijing, P. R. China
| | - Sili Ren
- State Key Laboratory of Solid Lubrication, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou, P. R. China
| |
Collapse
|
97
|
Dolgonos GA, Loboda OA, Boese AD. Development of Embedded and Performance of Density Functional Methods for Molecular Crystals. J Phys Chem A 2018; 122:708-713. [PMID: 29265819 DOI: 10.1021/acs.jpca.7b12467] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report an alternative quantum mechanical:quantum mechanical (QM:QM) method to the currently used periodic density functional calculations including dispersion and investigate its performance with respect to main structural and energetic properties of the X23 set of molecular crystals. By setting the goal of reproducing reference periodic BLYP+D3 values and by embedding BLYP+D3 into DFTB, we obtain results similar to those of periodic BLYP+D3-typically within 1-2% in lattice energies and ∼0.4% in cell volumes. The accuracy of this QM:QM method in comparison to DFTB+D and DFT+D for the X23 set of molecular crystals is discussed.
Collapse
Affiliation(s)
- Grygoriy A Dolgonos
- Institute of Chemistry, University of Graz , Heinrichstrasse 28/IV, 8010 Graz, Austria
| | - Oleksandr A Loboda
- Institute of Chemistry, University of Graz , Heinrichstrasse 28/IV, 8010 Graz, Austria
| | - A Daniel Boese
- Institute of Chemistry, University of Graz , Heinrichstrasse 28/IV, 8010 Graz, Austria
| |
Collapse
|
98
|
Lawson Daku LM. Spin-state dependence of the structural and vibrational properties of solvated iron(ii) polypyridyl complexes from AIMD simulations: aqueous [Fe(bpy)3]Cl2, a case study. Phys Chem Chem Phys 2018; 20:6236-6253. [DOI: 10.1039/c7cp07862e] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
LS and HS IR spectra of aqueous [Fe(bpy)3]2+ and corresponding HS–LS difference IR spectrum as obtained from state-of-the-art ab initio molecular dynamics simulations applied to the determination of the structural and vibrational properties of the solvated complex.
Collapse
|
99
|
Karmakar A, Chandra A. Dynamics of vibrational spectral diffusion in water: Effects of dispersion interactions, temperature, density, system size and fictitious orbital mass. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2017.11.027] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
100
|
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
|