1
|
Shi T, Wang Z, Aldossary A, Liu Y, Li XS, Head-Gordon M. Local Second Order Mo̷ller-Plesset Theory with a Single Threshold Using Orthogonal Virtual Orbitals: A Distributed Memory Implementation. J Chem Theory Comput 2024. [PMID: 39221855 DOI: 10.1021/acs.jctc.4c01016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
In order to alleviate the computational burden associated with superlinear compute scalings with molecular size in electron correlation methods, researchers have developed local correlation methods that wisely treat relatively small contributions as zeros but still yield accurate energy approximation. Such local correlation techniques can also be combined with parallel computing resources to obtain further efficiency and scalability. This work focuses on the distributed memory parallel implementation of a local correlation method for second order Mo̷ller-Plesset (MP2) theory. This method also only has a single threshold to control the dropping of terms and accuracy of different computing kernels in the algorithm. The process partitioning strategy and distributed parallel implementation with the message passing interface (MPI) are discussed. In particular, the algorithm relies on a fixed sparsity pattern matrix multiplication and a corresponding distributed conjugate gradient solver, which exhibits almost linear scaling in both strong and weak scaling analyses. Numerical experiments on a range of molecules, including linear chains and molecules with 2 and 3-dimensional characters, are reported. For example, with only 32 MPI ranks, this MP2 implementation can calculate the correlation energy of vancomycin in def2-TZVP basis within 0.003% accuracy (10-6.5 threshold) in half an hour, where the same problem is unfeasible to solve with sequential or pure shared memory implementations.
Collapse
Affiliation(s)
- Tianyi Shi
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Zhenling Wang
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | | | - Yang Liu
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Xiaoye S Li
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Mackie CJ, Lu W, Liang J, Kostko O, Bandyopadhyay B, Gupta I, Ahmed M, Head-Gordon M. Magic Numbers and Stabilities of Photoionized Water Clusters: Computational and Experimental Characterization of the Nanosolvated Hydronium Ion. J Phys Chem A 2023. [PMID: 37441795 DOI: 10.1021/acs.jpca.3c02230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
The stability and distributions of small water clusters generated in a supersonic beam expansion are interrogated by tunable vacuum ultraviolet (VUV) radiation generated at a synchrotron. Time-of-flight mass spectrometry reveals enhanced population of various protonated water clusters (H+(H2O)n) based upon ionization energy and photoionization distance from source, suggesting there are "magic" numbers below the traditional n = 21 that predominates in the literature. These intensity distributions suggest that VUV threshold photoionization (11.0-11.5 eV) of neutral water clusters close to the nozzle exit leads to a different nonequilibrium state compared to a skimmed molecular beam. This results in the appearance of a new magic number at 14. Metadynamics conformer searches coupled with modern density functional calculations are used to identify the global minimum energy structures of protonated water clusters between n = 2 and 21, as well as the manifold of low-lying metastable minima. New lowest energy structures are reported for the cases of n = 5, 6, 11, 12, 16, and 18, and special stability is identified by several measures. These theoretical results are in agreement with the experiments performed in this work in that n = 14 is shown to exhibit additional stability, based on the computed second-order stabilization energy relative to most cluster sizes, though not to the extent of the well-known n = 21 cluster. Other cluster sizes that show some additional energetic stability are n = 7, 9, 12, 17, and 19. To gain insight into the balance between ion-water and water-water interactions as a function of the cluster size, an analysis of the effective two-body interactions (which sum exactly to the total interaction energy) was performed. This analysis reveals a crossover as a function of cluster size between a water-hydronium-dominated regime for small clusters and a water-water-dominated regime for larger clusters around n = 17.
Collapse
Affiliation(s)
- Cameron J Mackie
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Wenchao Lu
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jiashu Liang
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Oleg Kostko
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Biswajit Bandyopadhyay
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Ishan Gupta
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Musahid Ahmed
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
3
|
Herman KM, Xantheas SS. An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25. Phys Chem Chem Phys 2023; 25:7120-7143. [PMID: 36853239 DOI: 10.1039/d2cp03241d] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2-11, 16-17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (De's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2-25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83-119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (De's) throughout the entire cluster range (n = 2-25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ∼4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6-3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA.
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA. .,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, WA, 99352, USA.
| |
Collapse
|
4
|
Tokmachev AM. Statistical Equivalence of Quantum Chemical Methods for Energy Distribution in Water Clusters. Chemphyschem 2023; 24:e202300019. [PMID: 36787108 DOI: 10.1002/cphc.202300019] [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: 01/09/2023] [Revised: 02/13/2023] [Accepted: 02/14/2023] [Indexed: 02/15/2023]
Abstract
Methods of quantum chemistry are instrumental in understanding molecular structures and properties. However, the results demonstrate significant variability, which is difficult to predict and rationalize. The fundamental question is whether some molecular systems exhibit properties invariant with respect to the computational method. The idea explored here is that collective properties of statistical ensembles should be more robust than characteristics of individual molecules and their arbitrary sets. This effect is demonstrated for the complete set of hydrogen-bond topologies of the dodecahedral water cluster (H2 O)20 . Non-Gaussian energy distributions produced by various methods have the same functional form despite strong differences in mean values and standard deviations. The conclusion is tested on methods of different complexity and origin employing a number of criteria. A linear mapping between the energies produced by different methods is discussed. The significance of the results is in establishing a collective equivalence property of quantum chemical methods.
Collapse
Affiliation(s)
- Andrey M Tokmachev
- National Research Center "Kurchatov Institute", Kurchatov Sq. 1, 123182, Moscow, Russia
| |
Collapse
|
5
|
Abella D, Franzese G, Hernández-Rojas J. Many-Body Contributions in Water Nanoclusters. ACS NANO 2023; 17:1959-1964. [PMID: 36695562 PMCID: PMC10781035 DOI: 10.1021/acsnano.2c06077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Accepted: 01/20/2023] [Indexed: 06/17/2023]
Abstract
Many-body interactions in water are known to be important but difficult to treat in atomistic models and often are included only as a correction. Polarizable models treat them explicitly, with long-range many-body potentials, within their classical approximation. However, their calculation is computationally expensive. Here, we evaluate how relevant the contributions to the many-body interaction associated with different coordination shells are. We calculate the global energy minimum, and the corresponding configuration, for nanoclusters of up to 20 water molecules. We find that including the first coordination shell, i.e., the five-body term of the central molecule, is enough to approximate within 5% the global energy minimum and its structure. We show that this result is valid for three different polarizable models, the Dang-Chang, the MB-pol, and the Kozack-Jordan potentials. This result suggests a strategy to develop many-body potentials for water that are reliable and, at the same time, computationally efficient.
Collapse
Affiliation(s)
- David Abella
- Instituto
de Física Interdisciplinar y Sistemas Complejos IFISC (CSIC-UIB), Campus UIB, 07122 Palma de Mallorca, Spain
- Secció
de Física Estadística i Interdisciplinària, Departament
de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
| | - Giancarlo Franzese
- Secció
de Física Estadística i Interdisciplinària, Departament
de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
- Institut
de Nanociència i Nanotecnologia, Universitat de Barcelona, 08028 Barcelona, Spain
| | - Javier Hernández-Rojas
- Departamento
de Física e IUdEA, Universidad de
La Laguna, 38205 La Laguna, Tenerife, Spain
| |
Collapse
|
6
|
Shape-shifters among water clusters. Struct Chem 2022. [DOI: 10.1007/s11224-022-02107-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
7
|
Mauger N, Plé T, Lagardère L, Huppert S, Piquemal JP. Improving Condensed-Phase Water Dynamics with Explicit Nuclear Quantum Effects: The Polarizable Q-AMOEBA Force Field. J Phys Chem B 2022; 126:8813-8826. [PMID: 36270033 DOI: 10.1021/acs.jpcb.2c04454] [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/2023]
Abstract
We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14-7 potential to model van der Waals interactions. It has been obtained by using a ForceBalance fitting strategy including high-level quantum chemistry reference energies and selected condensed-phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas-phase and condensed-phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H-O-H angle compared to its gas-phase equilibrium value, isotope effects, and so on. Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins, and nucleic acids opening the possibility for the large-scale study of the importance of NQEs in biophysics.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| |
Collapse
|
8
|
Heindel JP, Kirov MV, Xantheas SS. Hydrogen Bond Arrangements in (H2O)20, 24, 28 Clathrate Hydrate Cages: Optimization and Many-Body Analysis. J Chem Phys 2022; 157:094301. [DOI: 10.1063/5.0095335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We provide a detailed study of the hydrogen bonding arrangements, relative stability, residual entropy and an analysis of the many-body effects in the (H2O)20 (D-cage), (H2O)24 (T-cage) and (H2O)28 (H-cage) hollow cages making up structures I (sI) and II (sII) of clathrate hydrate lattices. Based on the enumeration of the possible hydrogen bonding networks for a fixed oxygen atom scaffold, the residual entropy () of these three gas phase cages was estimated at 0.75482, 0.7544 and 0.75417, where is the number of molecules and is Boltzmann's constant. A previously identified descriptor of enhanced stability based on the relative arrangement and connectivity of nearest-neighbor fragments on the polyhedral water cluster (Strong-Weak-Effective-Bond (SWEB) model), also applies to the larger hollow cages. The three cages contain a maximum of 7, 9 and 11 such preferable arrangements of trans nearest dimer pairs with one "free" OH bond on the donor molecule ( t1d dimers). The Many-Body Expansion (MBE) up to the 4-body suggests that the many-body terms vary nearly linearly with the cluster binding energy. Using a hierarchical approach of screening the relative stability of networks starting from optimizations with the TIP4P, TTM2.1-F and MB-pol classical potentials, subsequently refining at more accurate levels of electronic structure theory (DFT and MP2) and finally correcting for zero-point energy, we were able to identify a group of 4 low-lying isomers of the (H2O)24 T-cage, two of which are antisymmetric and the other two form a pair of antipode configurations.
Collapse
Affiliation(s)
| | - Mikhail V. Kirov
- Tyumen Scientific Center, Earth Cryosphere Institute Tyumen Scientific Centre of the Siberian Branch of the Russian Academy of Sciences, Russia
| | - Sotiris S. Xantheas
- Pacific Northwest National Laboratory, United States of America
- University of Washington
| |
Collapse
|
9
|
Naseem-Khan S, Lagardère L, Narth C, Cisneros GA, Ren P, Gresh N, Piquemal JP. Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:3607-3621. [PMID: 35575306 PMCID: PMC10851344 DOI: 10.1021/acs.jctc.2c00029] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We present the extension of the Sum of Interactions Between Fragments Ab initio Computed (SIBFA) many-body polarizable force field to condensed-phase molecular dynamics (MD) simulations. The quantum-inspired SIBFA procedure is grounded on simplified integrals obtained from localized molecular orbital theory and achieves full separability of its intermolecular potential. It embodies long-range multipolar electrostatics (up to quadrupole) coupled to a short-range penetration correction (up to charge-quadrupole), exchange repulsion, many-body polarization, many-body charge transfer/delocalization, exchange dispersion, and dispersion (up to C10). This enables the reproduction of all energy contributions of ab initio symmetry-adapted perturbation theory (SAPT(DFT)) gas-phase reference computations. The SIBFA approach has been integrated within the Tinker-HP massively parallel MD package. To do so, all SIBFA energy gradients have been derived and the approach has been extended to enable periodic boundary conditions simulations using smooth particle mesh Ewald. This novel implementation also notably includes a computationally tractable simplification of the many-body charge transfer/delocalization contribution. As a proof of concept, we perform a first computational experiment defining a water model fitted on a limited set of SAPT(DFT) data. SIBFA is shown to enable a satisfactory reproduction of both gas-phase energetic contributions and condensed-phase properties highlighting the importance of its physically motivated functional form.
Collapse
Affiliation(s)
- Sehr Naseem-Khan
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Louis Lagardère
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- IP2CT, FR 2622, CNRS, Sorbonne Université, 75005 Paris, France
| | | | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Nohad Gresh
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
| | - Jean-Philip Piquemal
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
10
|
Jong GC, Pak GI, Pak NJ, Jong RJ, Jon Y. Classification of small water clusters (SWCs) by (K, J, nO0)-X notation and study on conductive-like screening model (COSMO) densities of SWCs. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113613] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
11
|
Shee J, Loipersberger M, Rettig A, Lee J, Head-Gordon M. Regularized Second-Order Møller-Plesset Theory: A More Accurate Alternative to Conventional MP2 for Noncovalent Interactions and Transition Metal Thermochemistry for the Same Computational Cost. J Phys Chem Lett 2021; 12:12084-12097. [PMID: 34910484 PMCID: PMC10037552 DOI: 10.1021/acs.jpclett.1c03468] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Second-order Møller-Plesset theory (MP2) notoriously breaks down for π-driven dispersion interactions and dative bonds in transition metal complexes. Herein, we investigate three physically justified forms of single-parameter, energy-gap dependent regularization which can yield high and transferable accuracy for a variety of noncovalent interactions (including S22, S66, and L7 test sets) and (mostly closed shell) transition metal thermochemistry. Regularization serves to damp overestimated pairwise additive contributions, renormalizing first-order amplitudes such that the effects of higher-order correlations are incorporated. The optimal parameter values for the noncovalent and transition metal sets are 1.1, 0.7, and 0.4 for κ, σ, and σ2 regularizers, respectively. However, such regularization slightly degrades the accuracy of conventional MP2 for some small-molecule test sets, most of which have relatively large average frontier energy gaps. Our results suggest that appropriately regularized MP2 models may improve double hybrid density functionals, at no additional cost over conventional MP2.
Collapse
Affiliation(s)
- James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthias Loipersberger
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
12
|
Gudkovskikh SV, Kirov MV. Thermal stability of water polyhedra with square faces. J Mol Model 2021; 27:366. [PMID: 34850292 DOI: 10.1007/s00894-021-04996-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 11/24/2021] [Indexed: 11/30/2022]
Abstract
The ability to form numerous crystalline modifications of ice and gas hydrate frameworks is a characteristic feature of water. In fact, this structural variety is much wider due to the proton disorder. Configurations with different arrangements of hydrogen atoms (protons) in hydrogen bonds are not equivalent in their properties. Polyhedral water clusters are convenient objects for studying the effect of proton disorder on the properties of ice-like systems. It was previously established that the stability of water polyhedra is determined by the competition of two factors. The geometric factor gives preference to tetrahedrally coordinated structures with a large number of pentagonal faces. The topological factor takes into account the number of energetically most favorable types of H-bonds. This number increases with the number of square faces. It was found that tetrahedrally coordinated structures are not the most stable. However, these calculations were performed without taking thermal effects into account (Kirov M. V., J Phys Chem A, 124:4463 - 4470, 2020). The purpose of the present article is to study the structural stability of various water polyhedra at different temperatures. In the course of modeling, using the Amoeba force field, the advantage of configurations with a large number of square faces is demonstrated. The structure and energetics of surface defects are studied. Several very stable structures of unusual shape were found, including polyhedra which contain 4-coordinated molecules and polyhedra whose O-H groups are directed to the cluster center. The comparative analysis of cluster stability includes the temperature intervals of melting-like transitions.
Collapse
Affiliation(s)
- Sergey V Gudkovskikh
- RAS, Earth Cryosphere Inst Tyumen SC SB, Malygina Str 86, Tyumen, 625026, Russia
| | - Mikhail V Kirov
- RAS, Earth Cryosphere Inst Tyumen SC SB, Malygina Str 86, Tyumen, 625026, Russia.
| |
Collapse
|
13
|
|
14
|
Kananenka AA, Skinner JL. Unusually strong hydrogen bond cooperativity in particular (H 2O) 20 clusters. Phys Chem Chem Phys 2020; 22:18124-18131. [PMID: 32761035 DOI: 10.1039/d0cp02343d] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Drawing upon an intuitive charge-transfer-based picture of hydrogen bonding, we demonstrate that cooperativity effects acting in concert can lead to unusually strong hydrogen bonds in neutral water clusters. The structure, vibrational, and NMR properties of a (H2O)20 pentagonal dodecahedron cluster containing such a strong hydrogen bond were studied using second-order perturbation theory and density functional theory. The hydrogen bond length was found to be shorter than 2.50 Å. A large redshift of over 2000 cm-1 with respect to the isolated water molecule was predicted for the OH stretching frequency of the donor water molecule. A large downfield shift to 13.5 ppm of the isotropic part of the 1H magnetic shielding tensor together with an unusually large shielding anisotropy of 49.9 ppm was obtained. The hydrogen bond energy was calculated using symmetry-adapted perturbation theory and was found to be more than three times stronger than a typical hydrogen bond in liquid water.
Collapse
Affiliation(s)
- Alexei A Kananenka
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA.
| | - J L Skinner
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
15
|
Bilbrey JA, Heindel JP, Schram M, Bandyopadhyay P, Xantheas SS, Choudhury S. A look inside the black box: Using graph-theoretical descriptors to interpret a Continuous-Filter Convolutional Neural Network (CF-CNN) trained on the global and local minimum energy structures of neutral water clusters. J Chem Phys 2020; 153:024302. [DOI: 10.1063/5.0009933] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jenna A. Bilbrey
- Computing and Analytics Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Joseph P. Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Malachi Schram
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Sotiris S. Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Sutanay Choudhury
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
16
|
Menzel S, Millan S, Höfert SP, Nuhnen A, Gökpinar S, Schmitz A, Janiak C. Increase of network hydrophilicity from sql to lvt supramolecular isomers of Cu-MOFs with the bifunctional 4-(3,5-dimethyl-1H-pyrazol-4-yl)benzoate linker. Dalton Trans 2020; 49:12854-12864. [DOI: 10.1039/d0dt02642e] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A slight difference in the H-bonding of the linker pyrazole-NH group changes the framework hydrophilicity drastically.
Collapse
Affiliation(s)
- Saskia Menzel
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| | - Simon Millan
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| | - Simon-Patrick Höfert
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| | - Alexander Nuhnen
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| | - Serkan Gökpinar
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| | - Alexa Schmitz
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| | - Christoph Janiak
- Institut für Anorganische Chemie und Strukturchemie
- Heinrich-Heine Universität Düsseldorf
- 40204 Düsseldorf
- Germany
| |
Collapse
|
17
|
Rakshit A, Bandyopadhyay P, Heindel JP, Xantheas SS. Atlas of putative minima and low-lying energy networks of water clusters n = 3-25. J Chem Phys 2019; 151:214307. [PMID: 31822087 DOI: 10.1063/1.5128378] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a database consisting of the putative minima and ∼3.2 × 106 local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as low-lying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation (MP2)/aug-cc-pVTZ level of theory. In total, we have identified 3 138 303 networks corresponding to local minima of the clusters n = 3-25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1-3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but not limited to neural networks and machine and deep learning) to understand the properties of water, nature's most important substance.
Collapse
Affiliation(s)
- Avijit Rakshit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
18
|
What is the hydrophobic interaction contribution to the stabilization of micro-hydrated complexes of trimethylamine oxide (TMAO)? A joint DFT-D, QTAIM, and MESP study. J Mol Model 2019; 25:363. [DOI: 10.1007/s00894-019-4217-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Accepted: 09/30/2019] [Indexed: 11/27/2022]
|
19
|
Das AK, Urban L, Leven I, Loipersberger M, Aldossary A, Head-Gordon M, Head-Gordon T. Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis. J Chem Theory Comput 2019; 15:5001-5013. [DOI: 10.1021/acs.jctc.9b00478] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
20
|
Liu C, Piquemal JP, Ren P. AMOEBA+ Classical Potential for Modeling Molecular Interactions. J Chem Theory Comput 2019; 15:4122-4139. [PMID: 31136175 DOI: 10.1021/acs.jctc.9b00261] [Citation(s) in RCA: 76] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Classical potentials based on isotropic and additive atomic charges have been widely used to model molecules in computers for the past few decades. The crude approximations in the underlying physics are hindering both their accuracy and transferability across chemical and physical environments. Here we present a new classical potential, AMOEBA+, to capture essential intermolecular forces, including permanent electrostatics, repulsion, dispersion, many-body polarization, short-range charge penetration, and charge transfer, by extending the polarizable multipole-based AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) model. For a set of common organic molecules, we show that AMOEBA+ with general parameters can reproduce both quantum mechanical interactions and energy decompositions according to Symmetry-Adapted Perturbation Theory (SAPT). Additionally, a new water model based on the AMOEBA+ framework captures various liquid-phase properties in molecular dynamics simulations while remaining consistent with SAPT energy decompositions, utilizing both ab initio data and experimental liquid properties. Our results demonstrate that it is possible to improve the physical basis of classical force fields to advance their accuracy and general applicability.
Collapse
Affiliation(s)
- Chengwen Liu
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States.,Laboratoire de Chimie Théorique , Sorbonne Université, UMR7616 CNRS , Paris 75252 , France.,Institut Universitaire de France , Paris Cedex 05, 75005 , France
| | - Pengyu Ren
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
21
|
Yang L, Ji H, Liu X, Lu W. Ring-Stacking Water Clusters: Morphology and Stabilities. ChemistryOpen 2019; 8:210-218. [PMID: 30815330 PMCID: PMC6376211 DOI: 10.1002/open.201800284] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 01/16/2019] [Indexed: 12/12/2022] Open
Abstract
The structures and interaction energies of water clusters with ring stacking motifs are studied by using ab initio calculations. The structures of the water clusters are constructed by stacking either single rings or multi-rings of tetramer, pentamer, and hexamer. We found that, in the single-ring-stacking motif, the most stable isomers exhibit an alternative clockwise-anticlockwise stacking pattern. We also show that four-layer single-ring-stacking isomers are not energetically favorable in comparison with those of two-layer multi-ring-stacking isomers. The relative stability of the isomers is also analyzed in terms of H-bond strength and elastic distortions of the water molecules.
Collapse
Affiliation(s)
- Liu Yang
- Center for Quantum Sciences and School of PhysicsNortheast Normal UniversityChangchun130117China
| | - Hanyang Ji
- Center for Quantum Sciences and School of PhysicsNortheast Normal UniversityChangchun130117China
| | - Xiaojie Liu
- Center for Quantum Sciences and School of PhysicsNortheast Normal UniversityChangchun130117China
- Center for Advanced Optoelectronic Functional Materials Research and Key Laboratory of UV Light-Emitting Materials and Technology of Ministry of EducationsNortheast Normal UniversityChangchun130024China
| | - Wen‐Cai Lu
- State Key Laboratory of Theoretical and Computational Chemistry Institute of Theoretical ChemistryJilin University ChangchunJilin130021PR China
| |
Collapse
|
22
|
Morgante P, Peverati R. ACCDB: A collection of chemistry databases for broad computational purposes. J Comput Chem 2018; 40:839-848. [DOI: 10.1002/jcc.25761] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 10/09/2018] [Accepted: 11/11/2018] [Indexed: 01/04/2023]
Affiliation(s)
- Pierpaolo Morgante
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| | - Roberto Peverati
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| |
Collapse
|
23
|
Vibrational, energetic-dynamical and dissociation properties of water clusters in static electric fields: Non-equilibrium molecular-dynamics insights. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.08.061] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
24
|
Mardirossian N, Head-Gordon M. Survival of the most transferable at the top of Jacob's ladder: Defining and testing the ωB97M(2) double hybrid density functional. J Chem Phys 2018; 148:241736. [PMID: 29960332 PMCID: PMC5991970 DOI: 10.1063/1.5025226] [Citation(s) in RCA: 117] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 05/09/2018] [Indexed: 01/19/2023] Open
Abstract
A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resulting functional, ωB97M(2), is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 data points clearly surpasses that of all of the tested density functionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
25
|
Demerdash O, Wang L, Head‐Gordon T. Advanced models for water simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1355] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Omar Demerdash
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California, Davis Davis CA USA
| | - Teresa Head‐Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
- Department of Bioengineering University of California Berkeley CA USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley CA USA
| |
Collapse
|
26
|
Hamm P, Fanourgakis GS, Xantheas SS. A surprisingly simple correlation between the classical and quantum structural networks in liquid water. J Chem Phys 2017; 147:064506. [DOI: 10.1063/1.4993166] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Peter Hamm
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| | - George S. Fanourgakis
- Environmental Chemical Processes Laboratory, Department of Chemistry, University of Crete, P.O. Box 2208, 70013 Heraklion, Greece
| | - Sotiris S. Xantheas
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
27
|
Witte J, Neaton JB, Head-Gordon M. Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C. J Chem Phys 2017. [DOI: 10.1063/1.4986962] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
|
28
|
Mardirossian N, Head-Gordon M. Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1333644] [Citation(s) in RCA: 709] [Impact Index Per Article: 101.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
29
|
Manna D, Kesharwani MK, Sylvetsky N, Martin JML. Conventional and Explicitly Correlated ab Initio Benchmark Study on Water Clusters: Revision of the BEGDB and WATER27 Data Sets. J Chem Theory Comput 2017; 13:3136-3152. [DOI: 10.1021/acs.jctc.6b01046] [Citation(s) in RCA: 69] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Debashree Manna
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Manoj K. Kesharwani
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| |
Collapse
|
30
|
Schmitz G, Hättig C. Accuracy of Explicitly Correlated Local PNO-CCSD(T). J Chem Theory Comput 2017; 13:2623-2633. [DOI: 10.1021/acs.jctc.7b00180] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gunnar Schmitz
- Department
of Chemistry, Aarhus Universitet, Aarhus 8000, Denmark
| | - Christof Hättig
- Arbeitsgruppe
Quantenchemie, Ruhr-Universität, Bochum 44780, Germany
| |
Collapse
|
31
|
Ruiz Pestana L, Mardirossian N, Head-Gordon M, Head-Gordon T. Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals. Chem Sci 2017; 8:3554-3565. [PMID: 30155200 PMCID: PMC6092720 DOI: 10.1039/c6sc04711d] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2016] [Accepted: 02/24/2017] [Indexed: 01/23/2023] Open
Abstract
We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.
Collapse
Affiliation(s)
- Luis Ruiz Pestana
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , USA .
| | - Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry , Department of Chemistry , University of California , Berkeley , USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry , Department of Chemistry , University of California , Berkeley , USA
| | - Teresa Head-Gordon
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , USA .
- Kenneth S. Pitzer Center for Theoretical Chemistry , Department of Chemistry , University of California , Berkeley , USA
- Departments of Chemistry , Bioengineering , Chemical and Biomolecular Engineering , University of California , Berkeley , USA
| |
Collapse
|
32
|
Marques JMC, Pereira FB, Llanio-Trujillo JL, Abreu PE, Albertí M, Aguilar A, Pirani F, Bartolomei M. A global optimization perspective on molecular clusters. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:rsta.2016.0198. [PMID: 28320902 PMCID: PMC5360898 DOI: 10.1098/rsta.2016.0198] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 09/12/2016] [Indexed: 06/06/2023]
Abstract
Although there is a long history behind the idea of chemical structure, this is a key concept that continues to challenge chemists. Chemical structure is fundamental to understanding most of the properties of matter and its knowledge for complex systems requires the use of state-of-the-art techniques, either experimental or theoretical. From the theoretical view point, one needs to establish the interaction potential among the atoms or molecules of the system, which contains all the information regarding the energy landscape, and employ optimization algorithms to discover the relevant stationary points. In particular, global optimization methods are of major importance to search for the low-energy structures of molecular aggregates. We review the application of global optimization techniques to several molecular clusters; some new results are also reported. Emphasis is given to evolutionary algorithms and their application in the study of the microsolvation of alkali-metal and Ca2+ ions with various types of solvents.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.
Collapse
Affiliation(s)
- J M C Marques
- CQC, Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - F B Pereira
- Instituto Superior de Engenharia de Coimbra, Quinta da Nora, 3030-199 Coimbra, Portugal
- Centro de Informática e Sistemas da Universidade de Coimbra (CISUC), 3030-290 Coimbra, Portugal
| | - J L Llanio-Trujillo
- CQC, Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - P E Abreu
- CQC, Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - M Albertí
- IQTCUB, Departament de Ciència de Materials i Química Física, Universitat de Barcelona, 08028 Barcelona, Spain
| | - A Aguilar
- IQTCUB, Departament de Ciència de Materials i Química Física, Universitat de Barcelona, 08028 Barcelona, Spain
| | - F Pirani
- Dipartimento di Chimica, Biologia e Biotecnologie, Univertità di Perugia, 06123 Perugia, Italy
| | - M Bartolomei
- Consejo Superior de Investigaciones Científicas (IFF-CSIC), Instituto de Física Fundamental, Serrano 123, Madrid 28006, Spain
| |
Collapse
|
33
|
Witte J, Mardirossian N, Neaton JB, Head-Gordon M. Assessing DFT-D3 Damping Functions Across Widely Used Density Functionals: Can We Do Better? J Chem Theory Comput 2017; 13:2043-2052. [DOI: 10.1021/acs.jctc.7b00176] [Citation(s) in RCA: 56] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | | | - Jeffrey B. Neaton
- Kavli Energy
Nanosciences
Institute at Berkeley, Berkeley, California 94720, United States
| | | |
Collapse
|
34
|
A theoretical investigation of water–solute interactions: from facial parallel to guest–host structures. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2074-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
35
|
Structure and energetic characteristics of methane hydrates. From single cage to triple cage: A DFT-D study. J Mol Struct 2017. [DOI: 10.1016/j.molstruc.2016.10.093] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
36
|
Mardirossian N, Head-Gordon M. How Accurate Are the Minnesota Density Functionals for Noncovalent Interactions, Isomerization Energies, Thermochemistry, and Barrier Heights Involving Molecules Composed of Main-Group Elements? J Chem Theory Comput 2016; 12:4303-25. [PMID: 27537680 DOI: 10.1021/acs.jctc.6b00637] [Citation(s) in RCA: 286] [Impact Index Per Article: 35.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The 14 Minnesota density functionals published between the years 2005 and early 2016 are benchmarked on a comprehensive database of 4986 data points (84 data sets) involving molecules composed of main-group elements. The database includes noncovalent interactions, isomerization energies, thermochemistry, and barrier heights, as well as equilibrium bond lengths and equilibrium binding energies of noncovalent dimers. Additionally, the sensitivity of the Minnesota density functionals to the choice of basis set and integration grid is explored for both noncovalent interactions and thermochemistry. Overall, the main strength of the hybrid Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., M06-2X), barrier heights (e.g., M08-HX, M08-SO, MN15), and systems heavily characterized by self-interaction error (e.g., M06-2X, M08-HX, M08-SO, MN15), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-2X is recommended from the 10 hybrid Minnesota functionals). Similarly, the main strength of the local Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., MN15-L), barrier heights (e.g., MN12-L), and systems heavily characterized by self-interaction error (e.g., MN12-L and MN15-L), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-L is clearly the best from the four local Minnesota functionals). As an overall guide, M06-2X and MN15 are perhaps the most broadly useful hybrid Minnesota functionals, while M06-L and MN15-L are perhaps the most broadly useful local Minnesota functionals, although each has different strengths and weaknesses.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California , Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory , Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California , Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory , Berkeley, California 94720, United States
| |
Collapse
|
37
|
Mezei PD, Ruzsinszky A, Csonka GI. Application of a Dual-Hybrid Direct Random Phase Approximation to Water Clusters. J Chem Theory Comput 2016; 12:4222-32. [DOI: 10.1021/acs.jctc.6b00323] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Pál D. Mezei
- Department
of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, H-1521 Budapest, Hungary
| | - Adrienn Ruzsinszky
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Gábor I. Csonka
- Department
of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, H-1521 Budapest, Hungary
| |
Collapse
|
38
|
Singh G, Nandi A, Gadre SR. Breaking the bottleneck: Use of molecular tailoring approach for the estimation of binding energies at MP2/CBS limit for large water clusters. J Chem Phys 2016; 144:104102. [PMID: 26979676 DOI: 10.1063/1.4943115] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
A pragmatic method based on the molecular tailoring approach (MTA) for estimating the complete basis set (CBS) limit at Møller-Plesset second order perturbation (MP2) theory accurately for large molecular clusters with limited computational resources is developed. It is applied to water clusters, (H2O)n (n = 7, 8, 10, 16, 17, and 25) optimized employing aug-cc-pVDZ (aVDZ) basis-set. Binding energies (BEs) of these clusters are estimated at the MP2/aug-cc-pVNZ (aVNZ) [N = T, Q, and 5 (whenever possible)] levels of theory employing grafted MTA (GMTA) methodology and are found to lie within 0.2 kcal/mol of the corresponding full calculation MP2 BE, wherever available. The results are extrapolated to CBS limit using a three point formula. The GMTA-MP2 calculations are feasible on off-the-shelf hardware and show around 50%-65% saving of computational time. The methodology has a potential for application to molecular clusters containing ∼100 atoms.
Collapse
Affiliation(s)
- Gurmeet Singh
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Apurba Nandi
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
39
|
Mardirossian N, Head-Gordon M. ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation. J Chem Phys 2016; 144:214110. [DOI: 10.1063/1.4952647] [Citation(s) in RCA: 370] [Impact Index Per Article: 46.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
40
|
Dral P, Wu X, Spörkel L, Koslowski A, Thiel W. Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks for Ground-State Properties. J Chem Theory Comput 2016; 12:1097-120. [PMID: 26771261 PMCID: PMC4785506 DOI: 10.1021/acs.jctc.5b01047] [Citation(s) in RCA: 62] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2015] [Indexed: 11/30/2022]
Abstract
The semiempirical orthogonalization-corrected OMx methods (OM1, OM2, and OM3) go beyond the standard MNDO model by including additional interactions in the electronic structure calculation. When augmented with empirical dispersion corrections, the resulting OMx-Dn approaches offer a fast and robust treatment of noncovalent interactions. Here we evaluate the performance of the OMx and OMx-Dn methods for a variety of ground-state properties using a large and diverse collection of benchmark sets from the literature, with a total of 13035 original and derived reference data. Extensive comparisons are made with the results from established semiempirical methods (MNDO, AM1, PM3, PM6, and PM7) that also use the NDDO (neglect of diatomic differential overlap) integral approximation. Statistical evaluations show that the OMx and OMx-Dn methods outperform the other methods for most of the benchmark sets.
Collapse
Affiliation(s)
- Pavlo
O. Dral
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Xin Wu
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Lasse Spörkel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Axel Koslowski
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| |
Collapse
|
41
|
|
42
|
Affiliation(s)
- Imre Bakó
- Institute
of Organic Chemistry,
Research Centre for Natural Sciences, Hungarian Academy of Sciences, P.O. Box 286, Budapest, 1051 Hungary
| | - István Mayer
- Institute
of Organic Chemistry,
Research Centre for Natural Sciences, Hungarian Academy of Sciences, P.O. Box 286, Budapest, 1051 Hungary
| |
Collapse
|
43
|
Pezzotti G, Puppulin L, La Rosa A, Boffelli M, Zhu W, McEntire BJ, Hosogi S, Nakahari T, Marunaka Y. Effect of pH and monovalent cations on the Raman spectrum of water: Basics revisited and application to measure concentration gradients at water/solid interface in Si3N4 biomaterial. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.10.010] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
44
|
Pradzynski CC, Dierking CW, Zurheide F, Forck RM, Buck U, Zeuch T, Xantheas SS. Infrared detection of (H2O)20 isomers of exceptional stability: a drop-like and a face-sharing pentagonal prism cluster. Phys Chem Chem Phys 2015; 16:26691-6. [PMID: 25231162 DOI: 10.1039/c4cp03642e] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Water clusters with internally solvated water molecules are widespread models that mimic the local environment of the condensed phase. The appearance of stable (H2O)n cluster isomers having a fully coordinated interior molecule has been theoretically predicted to occur around the n = 20 size range. However, our current knowledge about the size regime in which those structures become energetically more stable has remained hypothetical from simulations in lieu of the absence of precisely size-resolved experimental measurements. Here we report size and isomer selective infrared (IR) spectra of (H2O)20 clusters tagged with a sodium atom by employing IR excitation modulated photoionization spectroscopy. The observed absorption patterns in the OH stretching region are consistent with the theoretically predicted spectra of two structurally distinct isomers of exceptional stability: a drop-like cluster with a fully coordinated (interior) water molecule and an edge-sharing pentagonal prism cluster in which all atoms are on the surface. The drop-like structure is the first experimentally detected water cluster exhibiting the local connectivity found in liquid water.
Collapse
Affiliation(s)
- Christoph C Pradzynski
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany.
| | | | | | | | | | | | | |
Collapse
|
45
|
Sahu N, Gadre SR, Rakshit A, Bandyopadhyay P, Miliordos E, Xantheas SS. Low energy isomers of (H2O)25 from a hierarchical method based on Monte Carlo temperature basin paving and molecular tailoring approaches benchmarked by MP2 calculations. J Chem Phys 2015; 141:164304. [PMID: 25362296 DOI: 10.1063/1.4897535] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report new global minimum candidate structures for the (H2O)25 cluster that are lower in energy than the ones reported previously and correspond to hydrogen bonded networks with 42 hydrogen bonds and an interior, fully coordinated water molecule. These were obtained as a result of a hierarchical approach based on initial Monte Carlo Temperature Basin Paving sampling of the cluster's Potential Energy Surface with the Effective Fragment Potential, subsequent geometry optimization using the Molecular Tailoring Approach with the fragments treated at the second order Møller-Plesset (MP2) perturbation (MTA-MP2) and final refinement of the entire cluster at the MP2 level of theory. The MTA-MP2 optimized cluster geometries, constructed from the fragments, were found to be within <0.5 kcal/mol from the minimum geometries obtained from the MP2 optimization of the entire (H2O)25 cluster. In addition, the grafting of the MTA-MP2 energies yields electronic energies that are within <0.3 kcal/mol from the MP2 energies of the entire cluster while preserving their energy rank order. Finally, the MTA-MP2 approach was found to reproduce the MP2 harmonic vibrational frequencies, constructed from the fragments, quite accurately when compared to the MP2 ones of the entire cluster in both the HOH bending and the OH stretching regions of the spectra.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Avijit Rakshit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Evangelos Miliordos
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| | - Sotiris S Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| |
Collapse
|
46
|
Sahu N, Khire SS, Gadre SR. Structures, energetics and vibrational spectra of (H2O)32clusters: a journey from model potentials to correlated theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1062150] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
47
|
Qi R, Wang LP, Wang Q, Pande VS, Ren P. United polarizable multipole water model for molecular mechanics simulation. J Chem Phys 2015; 143:014504. [PMID: 26156485 PMCID: PMC4499046 DOI: 10.1063/1.4923338] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2015] [Accepted: 06/21/2015] [Indexed: 11/14/2022] Open
Abstract
We report the development of a united AMOEBA (uAMOEBA) polarizable water model, which is computationally 3-5 times more efficient than the three-site AMOEBA03 model in molecular dynamics simulations while providing comparable accuracy for gas-phase and liquid properties. In this coarse-grained polarizable water model, both electrostatic (permanent and induced) and van der Waals representations have been reduced to a single site located at the oxygen atom. The permanent charge distribution is described via the molecular dipole and quadrupole moments and the many-body polarization via an isotropic molecular polarizability, all located at the oxygen center. Similarly, a single van der Waals interaction site is used for each water molecule. Hydrogen atoms are retained only for the purpose of defining local frames for the molecular multipole moments and intramolecular vibrational modes. The parameters have been derived based on a combination of ab initio quantum mechanical and experimental data set containing gas-phase cluster structures and energies, and liquid thermodynamic properties. For validation, additional properties including dimer interaction energy, liquid structures, self-diffusion coefficient, and shear viscosity have been evaluated. The results demonstrate good transferability from the gas to the liquid phase over a wide range of temperatures, and from nonpolar to polar environments, due to the presence of molecular polarizability. The water coordination, hydrogen-bonding structure, and dynamic properties given by uAMOEBA are similar to those derived from the all-atom AMOEBA03 model and experiments. Thus, the current model is an accurate and efficient alternative for modeling water.
Collapse
Affiliation(s)
- Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Lee-Ping Wang
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Qiantao Wang
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
48
|
Gimelshein N, Gimelshein S, Pradzynski CC, Zeuch T, Buck U. The temperature and size distribution of large water clusters from a non-equilibrium model. J Chem Phys 2015; 142:244305. [DOI: 10.1063/1.4922312] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - S. Gimelshein
- University of Southern California, Los Angeles, California 90089, USA
| | - C. C. Pradzynski
- Institut für Physikalische Chemie, Universität Göttingen, Tammanstr. 6, D-37077 Göttingen, Germany
| | - T. Zeuch
- Institut für Physikalische Chemie, Universität Göttingen, Tammanstr. 6, D-37077 Göttingen, Germany
| | - U. Buck
- Max-Planck-Institut für Dynamik und Selbstorganisation, Am Faßberg 17, D-37077 Göttingen, Germany
| |
Collapse
|
49
|
Miliordos E, Xantheas SS. An accurate and efficient computational protocol for obtaining the complete basis set limits of the binding energies of water clusters at the MP2 and CCSD(T) levels of theory: Application to (H2O)m, m = 2-6, 8, 11, 16, and 17. J Chem Phys 2015; 142:234303. [DOI: 10.1063/1.4922262] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Evangelos Miliordos
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| | - Sotiris S. Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| |
Collapse
|
50
|
Vítek A, Arismendi-Arrieta DJ, Rodríguez-Cantano R, Prosmiti R, Villarreal P, Kalus R, Delgado-Barrio G. Computational investigations of the thermodynamic properties of size-selected water and Ar-water clusters: high-pressure transitions. Phys Chem Chem Phys 2015; 17:8792-801. [PMID: 25745673 DOI: 10.1039/c4cp04862h] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Classical parallel-tempering Monte Carlo simulations in the isothermal-isobaric ensemble were carried out for the (H2O)20 and Ar(H2O)20 clusters, over a wide range of temperatures (30-1000 K) and pressures (3 kPa-10 GPa) in order to study their thermodynamic properties and structural changes. The TIP4P/ice water model is employed for the water-water interactions, while both semiempirical and ab initio-based potentials are used to model the interaction between the rare-gas atoms and the water molecules. Temperature-pressure phase diagrams for these cluster systems were constructed by employing a two-dimensional multiple-histogram method. Structural changes were detected by analyzing the heat capacity landscape and the Pearson correlation coefficient profile for the interaction energy and volume. Those at high pressure correspond to solid-to-solid transitions and are found to be related to clathrate-like cages around the Ar atom. It is also shown that the formation and thermodynamic stability of such structures are determined by the intermolecular interaction between the rare-gas atoms and the host water molecules.
Collapse
Affiliation(s)
- Aleš Vítek
- IT4Innovations National Supercomputing Center, VSB - Technical University of Ostrava, 70833 Ostrava, Czech Republic
| | | | | | | | | | | | | |
Collapse
|