1
|
Rodriguez L, Natalizio M, Sode O. Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes. J Phys Chem A 2024; 128:4199-4205. [PMID: 38770817 DOI: 10.1021/acs.jpca.4c00639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.
Collapse
Affiliation(s)
- Larry Rodriguez
- Department of Chemistry and Biochemistry, California State University, Los Angeles, 5151 State University Drive, Los Angeles, California 90032, United States
| | - Michael Natalizio
- Department of Chemistry and Biochemistry, California State University, Los Angeles, 5151 State University Drive, Los Angeles, California 90032, United States
| | - Olaseni Sode
- Department of Chemistry and Biochemistry, California State University, Los Angeles, 5151 State University Drive, Los Angeles, California 90032, United States
| |
Collapse
|
2
|
Riera M, Knight C, Bull-Vulpe EF, Zhu X, Agnew H, Smith DGA, Simmonett AC, Paesani F. MBX: A many-body energy and force calculator for data-driven many-body simulations. J Chem Phys 2023; 159:054802. [PMID: 37526156 PMCID: PMC10550339 DOI: 10.1063/5.0156036] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 07/11/2023] [Indexed: 08/02/2023] Open
Abstract
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Argonne National Laboratory, Computational Science Division, Lemont, Illinois 60439, USA
| | - Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Xuanyu Zhu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | | | - Andrew C. Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
3
|
Gartner T, Lauzin C, McKellar ARW, Moazzen-Ahmadi N. Infrared Spectra of the Water-CO 2 Complex in the 4.3-3.6 μm Region and Determination of the Ground State Tunneling Splitting for HDO-CO 2. J Phys Chem A 2023; 127:3668-3674. [PMID: 37067811 DOI: 10.1021/acs.jpca.3c01664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2023]
Abstract
Spectra of water─CO2 dimers are studied using a tunable mid-infrared source to probe a pulsed slit jet supersonic expansion. H2O-CO2 and D2O-CO2 are observed in the CO2 ν3 fundamental region (≈2350 cm-1), D2O-CO2 is also observed in the D2O ν3 fundamental region (≈2790 cm-1), and HDO-CO2 is observed in the HDO O-D stretch fundamental region (≈2720 cm-1), all for the first time in these regions. Analysis of the spectra yields excited state rotational parameters and vibrational shifts. They also yield the first experimental values of the ground state internal rotation tunneling splittings for D2O-CO2 (0.003 cm-1) and HDO-CO2 (0.0234 cm-1). The latter value is a direct determination made possible by the reduced symmetry of HDO-CO2. These results provide stringent and easily interpreted tests for theoretical water-CO2 potential energy surface calculations.
Collapse
Affiliation(s)
- Travis Gartner
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary T2N 1N4, Alberta, Canada
| | - Clément Lauzin
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve B-1348, Belgium
| | - A R W McKellar
- National Research Council of Canada, Ottawa K1A 0R6, Ontario, Canada
| | - Nasser Moazzen-Ahmadi
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary T2N 1N4, Alberta, Canada
| |
Collapse
|
4
|
Houston PL, Qu C, Yu Q, Conte R, Nandi A, Li JK, Bowman JM. PESPIP: Software to fit complex molecular and many-body potential energy surfaces with permutationally invariant polynomials. J Chem Phys 2023; 158:044109. [PMID: 36725524 DOI: 10.1063/5.0134442] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We wish to describe a potential energy surface by using a basis of permutationally invariant polynomials whose coefficients will be determined by numerical regression so as to smoothly fit a dataset of electronic energies as well as, perhaps, gradients. The polynomials will be powers of transformed internuclear distances, usually either Morse variables, exp(-ri,j/λ), where λ is a constant range hyperparameter, or reciprocals of the distances, 1/ri,j. The question we address is how to create the most efficient basis, including (a) which polynomials to keep or discard, (b) how many polynomials will be needed, (c) how to make sure the polynomials correctly reproduce the zero interaction at a large distance, (d) how to ensure special symmetries, and (e) how to calculate gradients efficiently. This article discusses how these questions can be answered by using a set of programs to choose and manipulate the polynomials as well as to write efficient Fortran programs for the calculation of energies and gradients. A user-friendly interface for access to monomial symmetrization approach results is also described. The software for these programs is now publicly available.
Collapse
Affiliation(s)
- Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA and Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Jeffrey K Li
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
5
|
Bull-Vulpe EF, Riera M, Bore SL, Paesani F. Data-Driven Many-Body Potential Energy Functions for Generic Molecules: Linear Alkanes as a Proof-of-Concept Application. J Chem Theory Comput 2022. [PMID: 36113028 DOI: 10.1021/acs.jctc.2c00645] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a generalization of the many-body energy (MB-nrg) theoretical/computational framework that enables the development of data-driven potential energy functions (PEFs) for generic covalently bonded molecules, with arbitrary quantum mechanical accuracy. The "nearsightedness of electronic matter" is exploited to define monomers as "natural building blocks" on the basis of their distinct chemical identity. The energy of generic molecules is then expressed as a sum of individual many-body energies of incrementally larger subsystems. The MB-nrg PEFs represent the low-order n-body energies, with n = 1-4, using permutationally invariant polynomials derived from electronic structure data carried out at an arbitrary quantum mechanical level of theory, while all higher-order n-body terms (n > 4) are represented by a classical many-body polarization term. As a proof-of-concept application of the general MB-nrg framework, we present MB-nrg PEFs for linear alkanes. The MB-nrg PEFs are shown to accurately reproduce reference energies, harmonic frequencies, and potential energy scans of alkanes, independently of their length. Since, by construction, the MB-nrg framework introduced here can be applied to generic covalently bonded molecules, we envision future computer simulations of complex molecular systems using data-driven MB-nrg PEFs, with arbitrary quantum mechanical accuracy.
Collapse
Affiliation(s)
- Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Sigbjørn L. Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
6
|
Yue S, Riera M, Ghosh R, Panagiotopoulos AZ, Paesani F. Transferability of data-driven, many-body models for CO2 simulations in the vapor and liquid phases. J Chem Phys 2022; 156:104503. [DOI: 10.1063/5.0080061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shuwen Yue
- Princeton University, United States of America
| | - Marc Riera
- Chemistry and Biochemistry, University of California San Diego Department of Chemistry and Biochemistry, United States of America
| | - Raja Ghosh
- University of California San Diego, United States of America
| | | | | |
Collapse
|
7
|
Felker PM, Bacic Z. Intermolecular rovibrational states of the H2O-CO2 and D2O-CO2 van der Waals complexes. J Chem Phys 2022; 156:064301. [DOI: 10.1063/5.0083754] [Citation(s) in RCA: 2] [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 M. Felker
- Department of Chemistry and Biochemistry, UCLA, United States of America
| | - Zlatko Bacic
- Department of Chemistry, New York University Department of Chemistry, United States of America
| |
Collapse
|
8
|
Lambros E, Dasgupta S, Palos E, Swee S, Hu J, Paesani F. General Many-Body Framework for Data-Driven Potentials with Arbitrary Quantum Mechanical Accuracy: Water as a Case Study. J Chem Theory Comput 2021; 17:5635-5650. [PMID: 34370954 DOI: 10.1021/acs.jctc.1c00541] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a general framework for the development of data-driven many-body (MB) potential energy functions (MB-QM PEFs) that represent the interactions between small molecules at an arbitrary quantum-mechanical (QM) level of theory. As a demonstration, a family of MB-QM PEFs for water is rigorously derived from density functionals belonging to different rungs across Jacob's ladder of approximations within density functional theory (MB-DFT) and from Møller-Plesset perturbation theory (MB-MP2). Through a systematic analysis of individual MB contributions to the interaction energies of water clusters, we demonstrate that all MB-QM PEFs preserve the same accuracy as the corresponding ab initio calculations, with the exception of those derived from density functionals within the generalized gradient approximation (GGA). The differences between the DFT and MB-DFT results are traced back to density-driven errors that prevent GGA functionals from accurately representing the underlying molecular interactions for different cluster sizes and hydrogen-bonding arrangements. We show that this shortcoming may be overcome, within the MB formalism, by using density-corrected functionals (DC-DFT) that provide a more consistent representation of each individual MB contribution. This is demonstrated through the development of a MB-DFT PEF derived from DC-PBE-D3 data, which more accurately reproduce the corresponding ab initio results.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Steven Swee
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Jie Hu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
9
|
Heindel JP, Herman KM, Aprà E, Xantheas SS. Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH 4, CO 2, and H 2S in (H 2O) 20 Cages. J Phys Chem Lett 2021; 12:7574-7582. [PMID: 34347487 DOI: 10.1021/acs.jpclett.1c01884] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely, the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II, and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with an increasing basis set to the complete basis set (CBS) limit, necessitating the use of large basis sets up to the aug-cc-pV5Z and subsequent correction for basis set superposition error (BSSE). For the host hollow (H2O)20 cages, we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Møller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20, and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions, and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350 208 cores) of the "Cori" supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 min (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path toward parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as density functional theory (DFT) and MP2 including its spin-biased variants.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Edoardo Aprà
- William R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| |
Collapse
|
10
|
Cabrera-Ramírez A, Arismendi-Arrieta DJ, Valdés Á, Prosmiti R. Structural Stability of the CO 2 @sI Hydrate: a Bottom-Up Quantum Chemistry Approach on the Guest-Cage and Inter-Cage Interactions. Chemphyschem 2020; 21:2618-2628. [PMID: 33001534 DOI: 10.1002/cphc.202000753] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Revised: 09/30/2020] [Indexed: 12/27/2022]
Abstract
Through reliable first-principles computations, we have demonstrated the impact of CO2 molecules enclathration on the stability of sI clathrate hydrates. Given the delicate balance between the interaction energy components (van der Waals, hydrogen bonds) present on such systems, we follow a systematic bottom-up approach starting from the individual 512 and 512 62 sI cages, up to all existing combinations of two-adjacent sI crystal cages to evaluate how such clathrate-like models perform on the evaluation of the guest-host and first-neighbors inter-cage effects, respectively. Interaction and binding energies of the CO2 occupation of the sI cages were computed using DF-MP2 and different DFT/DFT-D electronic structure methodologies. The performance of selected DFT functionals, together with various semi-classical dispersion corrections schemes, were validated by comparison with reference ab initio DF-MP2 data, as well as experimental data from x-ray and neutron diffraction studies available. Our investigation confirms that the inclusion of the CO2 in the cage/s is an energetically favorable process, with the CO2 molecule preferring to occupy the large 512 62 sI cages compared to the 512 ones. Further, the present results conclude on the rigidity of the water cages arrangements, showing the importance of the inter-cage couplings in the cluster models under study. In particular, the guest-cage interaction is the key factor for the preferential orientation of the captured CO2 molecules in the sI cages, while the inter-cage interactions seems to cause minor distortions with the CO2 guest neighbors interactions do not extending beyond the large 512 62 sI cages. Such findings on these clathrate-like model systems are in accord with experimental observations, drawing a direct relevance to the structural stability of CO2 @sI clathrates.
Collapse
Affiliation(s)
| | - Daniel J Arismendi-Arrieta
- Donostia International Physics Center (DIPC), Paseo, Manuel de Lardizabal 4, 20018, Donostia-San Sebastián, Spain
| | - Álvaro Valdés
- Escuela de Física, Universidad Nacional de Colombia, Sede, A. A., 3840, Medellín, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), Serrano 123, Madrid, Spain
| |
Collapse
|
11
|
Ollitrault PJ, Baiardi A, Reiher M, Tavernelli I. Hardware efficient quantum algorithms for vibrational structure calculations. Chem Sci 2020; 11:6842-6855. [PMID: 32874524 PMCID: PMC7448527 DOI: 10.1039/d0sc01908a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Accepted: 06/02/2020] [Indexed: 11/21/2022] Open
Abstract
We introduce a framework for the calculation of ground and excited state energies of bosonic systems suitable for near-term quantum devices and apply it to molecular vibrational anharmonic Hamiltonians. Our method supports generic reference modal bases and Hamiltonian representations, including the ones that are routinely used in classical vibrational structure calculations. We test different parametrizations of the vibrational wavefunction, which can be encoded in quantum hardware, based either on heuristic circuits or on the bosonic Unitary Coupled Cluster Ansatz. In particular, we define a novel compact heuristic circuit and demonstrate that it provides a good compromise in terms of circuit depth, optimization costs, and accuracy. We evaluate the requirements, number of qubits and circuit depth, for the calculation of vibrational energies on quantum hardware and compare them with state-of-the-art classical vibrational structure algorithms for molecules with up to seven atoms.
Collapse
Affiliation(s)
- Pauline J Ollitrault
- IBM Quantum , IBM Research - Zurich , Säumerstrasse 4 , 8803 Rüschlikon , Switzerland .
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Alberto Baiardi
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Markus Reiher
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Ivano Tavernelli
- IBM Quantum , IBM Research - Zurich , Säumerstrasse 4 , 8803 Rüschlikon , Switzerland .
| |
Collapse
|
12
|
Yanes-Rodríguez R, Arismendi-Arrieta DJ, Prosmiti R. He Inclusion in Ice-like and Clathrate-like Frameworks: A Benchmark Quantum Chemistry Study of Guest-Host Interactions. J Chem Inf Model 2020; 60:3043-3056. [PMID: 32469514 DOI: 10.1021/acs.jcim.0c00349] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Energetics and structural properties of selected type and size He@hydrate frameworks, e.g., from regular structured ice channels to clathrate-like cages, are presented from first-principles quantum chemistry methods. The scarcity of information on He@hydrates makes such complexes challenging targets, while their computational study entails an interesting and arduous task. Some of them have been synthesized in the laboratory, which motivates further investigations on their stability. Hence, the main focus is to examine the performance and accuracy of different wave function-based electronic structure methods, such as MP2, CCSD(T), their explicitly correlated (F12) and domain-based local pair-natural orbital (DLPNO) analogs, as well as modern and conventional density functional theory (DFT) approaches, and analytical model potentials available. Different structures are considered, starting from the "simplest system" formed by a noble gas atom (such as He) and one water molecule, followed by the study of the "fundamental units" present in all ice-like and clathrate-like frameworks (such as pentamers and hexamers) and finally the description of interactions in the "building blocks" of three-dimensional (3D) ice channels (e.g., horizontal and perpendicular ice II and Ih) and clathrate-like cages, such as the 512 present in the most common sI, sII, and sH clathrate-hydrate structures. The idea is to provide well-converged DLPNO-CCSD(T) and DFMP2/CBS reference datasets that in turn are used to validate how DFT functionals (in total, 29 approaches from generalized-gradient approximation (GGA), meta-GGA, to hybrid and range-separated functionals, including dispersion correction treatments, were checked) and analytical semiempirical/ab initio-based potentials perform compared with high-level alternatives. Within all tested approaches, those best-performing were identified and classified. Most of the DFT/DFT-D functionals, as well as available analytical pairwise model potentials, face difficulties in describing both hydrogen-bonded water frameworks and dispersion bound He-water interactions. Including dispersion corrections yields an overall well-balanced performance for LCωPBE-D3BJ and PBE0-D4 functionals. Such benchmark datasets can benefit research into the development of new cheminformatics models, as can serve to guide and cross-check methodologies, lending increased predicted power to future molecular simulations for investigating the role of structures and phase transitions from nanoscale clusters to macroscopic crystalline structures.
Collapse
Affiliation(s)
| | - Daniel J Arismendi-Arrieta
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.,Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, Gipuzkoa, 20018 Donostia-San Sebastián, Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
13
|
Riera M, Yeh EP, Paesani F. Data-Driven Many-Body Models for Molecular Fluids: CO2/H2O Mixtures as a Case Study. J Chem Theory Comput 2020; 16:2246-2257. [DOI: 10.1021/acs.jctc.9b01175] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Eric P. Yeh
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
14
|
Lauzin C, Imbreckx AC, Foldes T, Vanfleteren T, Moazzen-Ahmadi N, Herman M. High-resolution spectroscopic study of the H2O–CO2 van der Waals complex in the 2OH overtone range. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1706776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- C. Lauzin
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain, Louvain-la-Neuve, Belgium
| | - A. C. Imbreckx
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain, Louvain-la-Neuve, Belgium
| | - T. Foldes
- Service de Chimie quantique et Photophysique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - T. Vanfleteren
- Service de Chimie quantique et Photophysique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - N. Moazzen-Ahmadi
- Department of Physics and Astronomy, University of Calgary, Calgary, Canada
| | - M. Herman
- Service de Chimie quantique et Photophysique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| |
Collapse
|
15
|
Xia GJ, Liu J, Liu ZF. Structural inhomogeneity as a factor promoting the homogenous catalysis of CO 2 hydrogenation by (PMe 3) 4RuH 2. Phys Chem Chem Phys 2019; 21:19252-19268. [PMID: 31441925 DOI: 10.1039/c9cp03288f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
During homogenous catalysis by organometallic complexes, the dissociation of a ligand to produce an unsaturated site on the metal center is often invoked as the first step of activation, especially when photo-excitation is involved. In this theoretical study, we demonstrated that under mild conditions, a thermodynamically unstable yet dynamically favorable active intermediate could be produced by the inhomogeneity of the solvent distribution around the catalyst rather than by ligand dissociation. This occurred at the end of the first catalytic cycle when the product was eliminated. The empty site was immediately filled by one of the additive molecules aggregated around the reaction center even when the intermediate complex was unstable, producing a transient and more active catalyst. This process accounted for the accelerated reaction rate observed in the landmark CO2 hydrogenation catalyzed by (PMe3)4RuH2 in supercritical CO2 when H2O, MeOH, or HNMe2 was added. This also suggests a new way to exploit the structural inhomogeneity around an organometallic complex for the design of superior catalysts.
Collapse
Affiliation(s)
- Guang-Jie Xia
- Department of Chemistry and Centre for Scientific Modeling and Computation Chinese University of Hong Kong, Shatin, Hong Kong, China.
| | - Jianwen Liu
- College of Materials Science and Engineering, Shenzhen University, Shenzhen, China
| | - Zhi-Feng Liu
- Department of Chemistry and Centre for Scientific Modeling and Computation Chinese University of Hong Kong, Shatin, Hong Kong, China.
| |
Collapse
|
16
|
Qu C, Bowman JM. Assessing the Importance of the H 2-H 2O-H 2O Three-Body Interaction on the Vibrational Frequency Shift of H 2 in the sII Clathrate Hydrate and Comparison with Experiment. J Phys Chem A 2019; 123:329-335. [PMID: 30525619 DOI: 10.1021/acs.jpca.8b11675] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The vibrational frequency shift of H2 in the 512 cage of the sII clathrate hydrate with and without surrounding water molecules is reported at 0 K, using diffusion Monte Carlo calculations for the ground and first excited vibrational states of H2. Approximate 1d calculations of the frequency shift are also reported with the H2 at the equilibrium position in the clathrate hydrate. These calculations make use of full-dimensional potential energy surfaces for the H2-H2O 2-body and H2-H2O-H2O 3-body interactions. The inclusion of the 3-body interaction is shown to make roughly a 33% contribution to the frequency shift and to bring the calculated value of -40 ± 4 cm-1 to within just 3 cm-1 of the experimental value at 20 K. This level of agreement with experiment may be somewhat fortuitous; however, the importance of the 3-body interaction is firmly established by these calculations. The frequency shift reported here with 2-body interactions does not agree with a previously reported calculation using just 2-body interactions from a different ab initio potential energy surface and with a different method to obtain the frequency shift. A similar 1d calculation of the frequency shift using that potential is reported and agrees to within roughly 10% of the one previously reported. Therefore, this suggests that the difference between the present calculations and the previous one using just 2-body interactions is mainly due to differences in the potential energy surfaces.
Collapse
Affiliation(s)
- Chen Qu
- Cherry L. Emerson Center for Scientifc Computations and Department of Chemistry , Emory University , Atlanta , Georgia 30322 , United States
| | - Joel M Bowman
- Cherry L. Emerson Center for Scientifc Computations and Department of Chemistry , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
17
|
Qu C, Bowman JM. Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases? Phys Chem Chem Phys 2019; 21:3397-3413. [DOI: 10.1039/c8cp04990d] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| |
Collapse
|
18
|
Arismendi-Arrieta DJ, Valdés Á, Prosmiti R. A Systematic Protocol for Benchmarking Guest-Host Interactions by First-Principles Computations: Capturing CO 2 in Clathrate Hydrates. Chemistry 2018; 24:9353-9363. [PMID: 29600599 DOI: 10.1002/chem.201800497] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2018] [Indexed: 01/19/2023]
Abstract
Clathrate hydrates of CO2 have been proposed as potential molecular materials in tackling important environmental problems related to greenhouse gases capture and storage. Despite the increasing interest in such hydrates and their technological applications, a molecular-level understanding of their formation and properties is still far from complete. Modeling interactions is a challenging and computationally demanding task, essential to reliably determine molecular properties. First-principles calculations for the CO2 guest in all sI, sII, and sH clathrate cages were performed, and the nature of the guest-host interactions, dominated by both hydrogen-bond and van der Waals forces, was systematically investigated. Different families of density functionals, as well as pairwise CO2 @H2 O model potentials versus wavefunction-based quantum approaches were studied for CO2 clathrate-like systems. Benchmark energies for new distance-dependent datasets, consisting of potential energy curves sampling representative configurations of the systems at the repulsive, near-equilibrium, and asymptotic/long-range regions of the full-dimensional surface, were generated, and a general protocol was proposed to assess the accuracy of such conventional and modern approaches at minimum and non-minimum orientations. Our results show that dispersion interactions are important in the guest-host stabilization energies of such clathrate cages, and the encapsulation of the CO2 into guest-free clathrate cages is always energetically favorable. In addition, the orientation of CO2 inside each cage was explored, and the ability of current promising approaches to accurately describe non-covalent CO2 @H2 O guest-host interactions in sI, sII, and sH clathrates was discussed, providing information for their applicability to future multiscale computer simulations.
Collapse
Affiliation(s)
| | - Álvaro Valdés
- Departamento de Física, Universidad Nacional de Colombia, Calle 26, Cra 39, Edificio, 404, Bogotá, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
Collapse
|
19
|
Piquemal JP, Jordan KD. Preface: Special Topic: From Quantum Mechanics to Force Fields. J Chem Phys 2018; 147:161401. [PMID: 29096449 DOI: 10.1063/1.5008887] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This Special Topic issue entitled "From Quantum Mechanics to Force Fields" is dedicated to the ongoing efforts of the theoretical chemistry community to develop a new generation of accurate force fields based on data from high-level electronic structure calculations and to develop faster electronic structure methods for testing and designing force fields as well as for carrying out simulations. This issue includes a collection of 35 original research articles that illustrate recent theoretical advances in the field. It provides a timely snapshot of recent developments in the generation of approaches to enable more accurate molecular simulations of processes important in chemistry, physics, biophysics, and materials science.
Collapse
Affiliation(s)
- Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, UMR 7616 CNRS, UPMC, Sorbonne Universités, 75252 Paris Cedex 05, France
| | - Kenneth D Jordan
- Department of Chemistry, The University of Pittsburgh, 219 Parkman Avenue, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
20
|
Qu C, Yu Q, Van Hoozen BL, Bowman JM, Vargas-Hernández RA. Assessing Gaussian Process Regression and Permutationally Invariant Polynomial Approaches To Represent High-Dimensional Potential Energy Surfaces. J Chem Theory Comput 2018; 14:3381-3396. [DOI: 10.1021/acs.jctc.8b00298] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Brian L. Van Hoozen
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | | |
Collapse
|
21
|
Affiliation(s)
- Chen Qu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Qi Yu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M. Bowman
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|