1
|
Zhao Z, Filip MA, Thom AJW. Rapidly convergent quantum Monte Carlo using a Chebyshev projector. Faraday Discuss 2024. [PMID: 39083360 DOI: 10.1039/d4fd00035h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2024]
Abstract
The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.
Collapse
Affiliation(s)
- Zijun Zhao
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK.
| | - Maria-Andreea Filip
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK.
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK.
| |
Collapse
|
2
|
Pham KN, Modrzejewski M, Klimeš J. Contributions beyond direct random-phase approximation in the binding energy of solid ethane, ethylene, and acetylene. J Chem Phys 2024; 160:224101. [PMID: 38856055 DOI: 10.1063/5.0207090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Accepted: 05/22/2024] [Indexed: 06/11/2024] Open
Abstract
The random-phase approximation (RPA) includes a subset of higher than second-order correlation-energy contributions, but stays in the same complexity class as the second-order Møller-Plesset perturbation theory (MP2) in both Gaussian-orbital and plane-wave codes. This makes RPA a promising ab initio electronic structure approach for the binding energies of molecular crystals. Still, some issues stand out in practical applications of RPA. Notably, compact clusters of nonpolar molecules are poorly described, and the interaction energies strongly depend on the reference single-determinant state. Using the many-body expansion of the binding energy of a crystal, we investigate those issues and the effect of beyond-RPA corrections. We find the beneficial effect of quartic-scaling exchange and non-ring coupled-cluster doubles corrections. The nonadditive interactions in compact trimers of molecules are improved by using the self-consistent Hartree-Fock orbitals instead of the usual Kohn-Sham states, but this kind of orbital input also leads to underestimated dimer energies. Overall, a substantial improvement over the RPA with a renormalized singles approach is possible at a modest quartic-scaling cost, which encourages further research into additional RPA corrections.
Collapse
Affiliation(s)
- Khanh Ngoc Pham
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
3
|
Gómez S, Flórez E, Acelas N, Cappelli C, Hadad C, Restrepo A. Encapsulation of charged halogens by the 5 12 water cage. Phys Chem Chem Phys 2024; 26:15426-15436. [PMID: 38747303 DOI: 10.1039/d4cp01340a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | - Elizabeth Flórez
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| | - Nancy Acelas
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | - Cacier Hadad
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| |
Collapse
|
4
|
Yanes-Rodríguez R, Prosmiti R. Analysing the stability of He-filled hydrates: how many He atoms fit in the sII crystal? Phys Chem Chem Phys 2024; 26:2519-2528. [PMID: 38170811 DOI: 10.1039/d3cp05410a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Clathrate hydrates have the ability to encapsulate atoms and molecules within their cavities, and thus they could be potentially large storage capacity materials. The present work studies the multiple cage occupancy effects in the recently discovered He@sII crystal. On the basis of previous theoretical and experimental findings, the stability of He(1/1)@sII, He(1/4)@sII and He(2/4)@sII crystals was analysed in terms of structural, mechanical and energetic properties. For this purpose, first-principles DFT/DFT-D computations were performed by using both semi-local and non-local functionals, not only to elucidate which configuration is the most energetically favoured, but also to scrutinize the relevance of the long-range dispersion interactions in these kinds of compounds. We have encountered that dispersion interactions play a fundamental role in describing the underlying interactions, and different tendencies were observed depending on the choice of the functional. We found that PW86PBE-XDM shows the best performance, while the non-local functionals tested here were not able to correctly account for them. The present results reveal that the most stable configuration is the one presenting singly occupied D cages and tetrahedrally occupied H cages (He(1/4)@sII) in line with the experimental observation.
Collapse
Affiliation(s)
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
5
|
Cieśliński D, Tucholska AM, Modrzejewski M. Post-Kohn-Sham Random-Phase Approximation and Correction Terms in the Expectation-Value Coupled-Cluster Formulation. J Chem Theory Comput 2023; 19:6619-6631. [PMID: 37774375 PMCID: PMC10569055 DOI: 10.1021/acs.jctc.3c00496] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Indexed: 10/01/2023]
Abstract
Using expectation-value coupled-cluster theory and many-body perturbation theory (MBPT), we formulate a series of corrections to the post-Kohn-Sham (post-KS) random-phase approximation (RPA) energy. The beyond-RPA terms are of two types: those accounting for the non-Hartree-Fock reference and those introducing the coupled-cluster doubles non-ring contractions. The contributions of the former type, introduced via the semicanonical orbital basis, drastically reduce the binding strength in noncovalent systems. The good accuracy is recovered by the attractive third-order doubles correction referred to as Ec2g. The existing RPA approaches based on KS orbitals neglect most of the proposed corrections but can perform well thanks to error cancellation. The proposed method accounts for every contribution in the state-of-the-art renormalized second-order perturbation theory (rPT2) approach but adds additional terms which initially contribute in the third order of MBPT. The cost of energy evaluation scales as noniterative O ( N 4 ) in the implementation with low-rank tensor decomposition. The numerical tests of the proposed approach demonstrate accurate results for noncovalent dimers of polar molecules and for the challenging many-body noncovalent cluster of CH4···(H2O)20.
Collapse
Affiliation(s)
- Dominik Cieśliński
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | | | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
6
|
Yanes-Rodríguez R, Cabrera-Ramírez A, Prosmiti R. Delving into guest-free and He-filled sI and sII clathrate hydrates: a first-principles computational study. Phys Chem Chem Phys 2022; 24:13119-13129. [PMID: 35587105 DOI: 10.1039/d2cp00701k] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dynamics of the formation of a specific clathrate hydrate as well as its thermodynamic transitions depend on the interactions between the trapped molecules and the host water lattice. The molecular-level understanding of the different underlying processes benefits not only the description of the properties of the system, but also allows the development of multiple technological applications such as gas storage, gas separation, energy transport, etc. In this work we investigate the stability of periodic crystalline structures, such as He@sI and He@sII clathrate hydrates by first-principles computations. We consider such host water networks interacting with a guest He atom using selected density functional theory approaches, in order to explore the effects on the encapsulation of a light atom in the sI/sII crystals, by deriving all energy components (guest-water, water-water, guest-guest). Structural properties and energies were first computed by structural relaxations of the He-filled and empty sI/sII unit cells, yielding lattice and compressibility parameters comparable to experimental and theoretical values available for those hydrates. According to the results obtained, the He enclathration in the sI/sII unit cells is a stabilizing process, and both He@sI and He@sII clathrates, considering single cage occupancy, are predicted to be stable whatever the XDM or D4 dispersion correction applied. Our results further reveal that despite the weak underlying interactions the He encapsulation has a rather notable effect on both lattice parameters and energetics, with the He@sII being the most energetically favorable in accord with recent experimental observations.
Collapse
Affiliation(s)
- Raquel Yanes-Rodríguez
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain. .,Doctoral Programme in Theoretical Chemistry and Computational Modelling, Doctoral School, Universidad Autónoma de Madrid, Spain
| | - Adriana Cabrera-Ramírez
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain. .,Doctoral Programme in Theoretical Chemistry and Computational Modelling, Doctoral School, Universidad Autónoma de Madrid, Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
7
|
Cabrera-Ramírez A, Arismendi-Arrieta DJ, Valdés Á, Prosmiti R. Exploring CO 2 @sI Clathrate Hydrates as CO 2 Storage Agents by Computational Density Functional Approaches. Chemphyschem 2021; 22:359-369. [PMID: 33368985 DOI: 10.1002/cphc.202001035] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Indexed: 12/21/2022]
Abstract
The formation of specific clathrate hydrates and their transformation at given thermodynamic conditions depends on the interactions between the guest molecule/s and the host water lattice. Understanding their structural stability is essential to control structure-property relations involved in different technological applications. Thus, the energetic aspects relative to CO2 @sI clathrate hydrate are investigated through the computation of the underlying interactions, dominated by hydrogen bonds and van der Waals forces, from first-principles electronic structure approaches. The stability of the CO2 @sI clathrate is evaluated by combining bottom-up and top-down approaches. Guest-free and CO2 guest-filled aperiodic cages, up to the gradually CO2 occupation of the entire sI periodic unit cells were considered. Saturation, cohesive and binding energies for the systems are determined by employing a variety of density functionals and their performance is assessed. The dispersion corrections on the non-covalent interactions are found to be important in the stabilization of the CO2 @sI energies, with the encapsulation of the CO2 into guest-free/empty cage/lattice being always an energetically favorable process for most of the functionals studied. The PW86PBE functional with XDM or D3(BJ) dispersion corrections predicts a lattice constant in accord to the experimental values available, and simultaneously provides a reliable description for the guest-host interactions in the periodic CO2 @sI crystal, as well as the energetics of its progressive single cage occupancy process. It has been found that the preferential orientation of the single CO2 in the large sI crystal cages has a stabilizing effect on the hydrate, concluding that the CO2 @sI structure is favored either by considering the individual building block cages or the complete sI unit cell crystal. Such benchmark and methodology cross-check studies benefit new data-driven model research by providing high-quality training information, with new insights that indicate the underlying factors governing their structure-driven stability, and triggering further investigations for controlling the stabilization of these promising long-term CO2 storage materials.
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 Medellín, A. A., 3840, Medellíın, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
Collapse
|
8
|
Modrzejewski M, Yourdkhani S, Śmiga S, Klimeš J. Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage. J Chem Theory Comput 2021; 17:804-817. [PMID: 33445879 DOI: 10.1021/acs.jctc.0c00966] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.
Collapse
Affiliation(s)
- Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland.,Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Sirous Yourdkhani
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziądzka 5, 87-100 Toruń, Poland
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
9
|
Cabrera-Ramírez A, Yanes-Rodríguez R, Prosmiti R. Computational density-functional approaches on finite-size and guest-lattice effects in CO 2@sII clathrate hydrate. J Chem Phys 2021; 154:044301. [PMID: 33514100 DOI: 10.1063/5.0039323] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We performed first-principles computations to investigate guest-host/host-host effects on the encapsulation of the CO2 molecule in sII clathrate hydrates from finite-size clusters up to periodic 3D crystal lattice systems. Structural and energetic properties were first computed for the individual and first-neighbors clathrate-like sII cages, where highly accurate ab initio quantum chemical methods are available nowadays, allowing in this way the assessment of the density functional (DFT) theoretical approaches employed. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. On this basis, structural relaxations of the CO2-filled and empty sII unit cells yield lattice and compressibility parameters comparable to experimental and previous theoretical values available for sII hydrates. According to these data, the CO2 enclathration in the sII clathrate cages is a stabilizing process, either by considering both guest-host and host-host interactions in the complete unit cell or only the guest-water energies for the individual clathrate-like sII cages. CO2@sII clathrates are predicted to be stable whatever the dispersion correction applied and in the case of single cage occupancy are found to be more stable than the CO2@sI structures. Our results reveal that DFT approaches could provide a good reasonable description of the underlying interactions, enabling the investigation of formation and transformation processes as a function of temperature and pressure.
Collapse
Affiliation(s)
| | | | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
10
|
Cladek BR, Everett SM, McDonnell MT, Tucker MG, Keffer DJ, Rawn CJ. Local structure and distortions of mixed methane-carbon dioxide hydrates. Commun Chem 2021; 4:6. [PMID: 36697523 PMCID: PMC9814247 DOI: 10.1038/s42004-020-00441-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Accepted: 12/07/2020] [Indexed: 01/28/2023] Open
Abstract
A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.
Collapse
Affiliation(s)
- Bernadette R. Cladek
- grid.411461.70000 0001 2315 1184Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996-2100 USA
| | - S. Michelle Everett
- grid.135519.a0000 0004 0446 2659Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6475 USA
| | - Marshall T. McDonnell
- grid.135519.a0000 0004 0446 2659Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6475 USA
| | - Matthew G. Tucker
- grid.135519.a0000 0004 0446 2659Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6475 USA
| | - David J. Keffer
- grid.411461.70000 0001 2315 1184Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996-2100 USA
| | - Claudia J. Rawn
- grid.411461.70000 0001 2315 1184Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996-2100 USA
| |
Collapse
|
11
|
Bier I, O'Connor D, Hsieh YT, Wen W, Hiszpanski AM, Han TYJ, Marom N. Crystal structure prediction of energetic materials and a twisted arene with Genarris and GAtor. CrystEngComm 2021. [DOI: 10.1039/d1ce00745a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
A molecular crystal structure prediction workflow, based on the random structure generator, Genarris, and the genetic algorithm (GA), GAtor, is successfully applied to two energetic materials and a chiral arene.
Collapse
Affiliation(s)
- Imanuel Bier
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Dana O'Connor
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Yun-Ting Hsieh
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Wen Wen
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Anna M. Hiszpanski
- Materials Science Division, Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
| | - T. Yong-Jin Han
- Materials Science Division, Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
- Department of Physics, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| |
Collapse
|
12
|
Tsimpanogiannis IN, Michalis VK, Economou IG. Novel methodology for the calculation of the enthalpy of enclathration of methane hydrates using molecular dynamics simulations. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1711976] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH) Thermi-Thessaloniki, Greece
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Vasileios K. Michalis
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
13
|
Wang Z, Yang L, Yang K, Deng R, Lin S. Effects of cage occupancy in hydrate by first-principles calculation and modification of the van der Waals-Platteeuw hypothesis. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:095402. [PMID: 31689692 DOI: 10.1088/1361-648x/ab5481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The changes of methane hydrate lattice with the decrease of cage occupancy were calculated by first-principles methods. The calculation results show that the decrease of the cages occupancy in sII and sH hydrate does not lead to large deformation in the lattice. Even if all the methane molecules are removed so that the hydrates have become new types of ice, the sII and sH lattices remain stable. The same conclusion is also true when the occupancy of the small cages in sI hydrate is reduced. However, the sI hydrate lattice will deform and almost collapse as the large cage occupancy decreases. These calculation results suggest that sI hydrate cannot exist with empty cages. Since the van der Waals-Platteeuw theory is based on the assumption that the stability of host lattice is independent of the occupancy of guest molecule, it would be applicable to sII and sH lattices, but not to sI hydrates. We propose a modification to the van der Waals-Platteeuw hypothesis so that the theory seems more reasonable.
Collapse
Affiliation(s)
- Zhao Wang
- School of Science, Hainan University, Haikou 570228, People's Republic of China
| | | | | | | | | |
Collapse
|
14
|
Molecular Insights into Cage Occupancy of Hydrogen Hydrate: A Computational Study. Processes (Basel) 2019. [DOI: 10.3390/pr7100699] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Density functional theory calculations and molecular dynamics simulations were performed to investigate the hydrogen storage capacity in the sII hydrate. Calculation results show that the optimum hydrogen storage capacity is ~5.6 wt%, with the double occupancy in the small cage and quintuple occupancy in the large cage. Molecular dynamics simulations indicate that these multiple occupied hydrogen hydrates can occur at mild conditions, and their stability will be further enhanced by increasing the pressure or decreasing the temperature. Our work highlights that the hydrate is a promising material for storing hydrogen.
Collapse
|
15
|
Pétuya C, Martin-Gondre L, Aurel P, Damay F, Desmedt A. Unraveling the metastability of the SI and SII carbon monoxide hydrate with a combined DFT-neutron diffraction investigation. J Chem Phys 2019; 150:184705. [PMID: 31091912 DOI: 10.1063/1.5093202] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Clathrate hydrates are crystalline compounds consisting of water molecules forming cages (so-called "host") inside of which "guest" molecules are encapsulated depending on the thermodynamic conditions of formation (systems stable at low temperature and high pressure). These icelike systems are naturally abundant on Earth and are generally expected to exist on icy celestial bodies. Carbon monoxide hydrate might be considered an important component of the carbon cycle in the solar system since CO gas is one of the predominant forms of carbon. Intriguing fundamental properties have also been reported: the CO hydrate initially forms in the sI structure (kinetically favored) and transforms into the sII structure (thermodynamically stable). Understanding and predicting the gas hydrate structural stability then become essential. The aim of this work is, thereby, to study the structural and energetic properties of the CO hydrate using density functional theory (DFT) calculations together with neutron diffraction measurements. In addition to the comparison of DFT-derived structural properties with those from experimental neutron diffraction, the originality of this work lies in the DFT-derived energy calculations performed on a complete unit cell (sI and sII) and not only by considering guest molecules confined in an isolated water cage (as usually performed for extracting the binding energies). Interestingly, an excellent agreement (within less than 1% error) is found between the measured and DFT-derived unit cell parameters by considering the Perdew-Burke-Ernzerhof (denoted PBE) functional. Moreover, a strategy is proposed for evaluating the hydrate structural stability on the basis of potential energy analysis of the total nonbonding energies (i.e., binding energy and water substructure nonbonding energy). It is found that the sII structure is the thermodynamically stable hydrate phase. In addition, increasing the CO content in the large cages has a stabilizing effect on the sII structure, while it destabilizes the sI structure. Such findings are in agreement with the recent experimental results evidencing the structural metastability of the CO hydrate.
Collapse
Affiliation(s)
- Claire Pétuya
- Institut des Sciences Moléculaires, UMR5255 CNRS-Université de Bordeaux, Talence, France
| | | | - Philippe Aurel
- Institut des Sciences Moléculaires, UMR5255 CNRS-Université de Bordeaux, Talence, France
| | - Françoise Damay
- Laboratoire Léon Brillouin UMR 12 CEA-CNRS, Gif-sur-Yvette, France
| | - Arnaud Desmedt
- Institut des Sciences Moléculaires, UMR5255 CNRS-Université de Bordeaux, Talence, France
| |
Collapse
|
16
|
Tetrahydrofuran (THF)-Mediated Structure of THF·(H2O)n=1–10: A Computational Study on the Formation of the THF Hydrate. CRYSTALS 2019. [DOI: 10.3390/cryst9020073] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Tetrahydrofuran (THF) is well known as a former and a promoter of clathrate hydrates, but the molecular mechanism for the formation of these compounds is not yet well understood. We performed ab initio calculations and ab initio molecular dynamics simulations to investigate the formation, structure, and stability of THF·(H2O)n=1–10 and its significance to the formation of the THF hydrate. Weak hydrogen bonds were found between THF and water molecules, and THF could promote water molecules from the planar pentagonal or hexagonal ring. As a promoter, THF could increase the binding ability of the CH4, CO2, or H2 molecule onto a water face, but could also enhance the adsorption of other THF molecules, causing an enrichment effect.
Collapse
|
17
|
Al-Hamdani YS, Tkatchenko A. Understanding non-covalent interactions in larger molecular complexes from first principles. J Chem Phys 2019; 150:010901. [PMID: 30621423 PMCID: PMC6910608 DOI: 10.1063/1.5075487] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Accepted: 12/05/2018] [Indexed: 01/02/2023] Open
Abstract
Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
18
|
|
19
|
Košata J, Merkl P, Teeratchanan P, Hermann A. Stability of Hydrogen Hydrates from Second-Order Møller-Plesset Perturbation Theory. J Phys Chem Lett 2018; 9:5624-5629. [PMID: 30179509 DOI: 10.1021/acs.jpclett.8b02274] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The formation of gas hydrates and clathrates critically depends on the interaction between the host water network and the guest gas species. Density functional calculations can struggle to quantitatively capture these dispersion-type interactions. Here, we report wave function-based calculations on hydrogen hydrates that combine periodic Hartree-Fock with a localized treatment of electronic correlation. We show that local second-order Møller-Plesset perturbation theory (LMP2) reproduces the stability of the different filled-ice-like hydrates in excellent agreement with experimental data. In contrast to various dispersion-corrected density functional theory implementations, LMP2 correctly identifies the pressures needed to stabilize the C0, C1, and C2 hydrates and does not find a spurious region of stability for an ice-Ih-based dihydrate. Our results suggest that LMP2 or similar approaches can provide quantitative insights into the mechanisms of formation and eventual decomposition of molecular host-guest compounds.
Collapse
Affiliation(s)
- Jan Košata
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| | - Padryk Merkl
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| | - Pattanasak Teeratchanan
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| | - Andreas Hermann
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| |
Collapse
|
20
|
Tsimpanogiannis IN, Economou IG. Monte Carlo simulation studies of clathrate hydrates: A review. J Supercrit Fluids 2018. [DOI: 10.1016/j.supflu.2017.12.017] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
21
|
Hapka M, Rajchel Ł, Modrzejewski M, Schäffer R, Chałasiński G, Szczęśniak MM. The nature of three-body interactions in DFT: Exchange and polarization effects. J Chem Phys 2017; 147:084106. [PMID: 28863509 DOI: 10.1063/1.4986291] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a physically motivated decomposition of density functional theory (DFT) 3-body nonadditive interaction energies into the exchange and density-deformation (polarization) components. The exchange component represents the effect of the Pauli exclusion in the wave function of the trimer and is found to be challenging for density functional approximations (DFAs). The remaining density-deformation nonadditivity is less dependent upon the DFAs. Numerical demonstration is carried out for rare gas atom trimers, Ar2-HX (X = F, Cl) complexes, and small hydrogen-bonded and van der Waals molecular systems. None of the tested semilocal, hybrid, and range-separated DFAs properly accounts for the nonadditive exchange in dispersion-bonded trimers. By contrast, for hydrogen-bonded systems, range-separated DFAs achieve a qualitative agreement to within 20% of the reference exchange energy. A reliable performance for all systems is obtained only when the monomers interact through the Hartree-Fock potential in the dispersion-free Pauli blockade scheme. Additionally, we identify the nonadditive second-order exchange-dispersion energy as an important but overlooked contribution in force-field-like dispersion corrections. Our results suggest that range-separated functionals do not include this component, although semilocal and global hybrid DFAs appear to imitate it in the short range.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Łukasz Rajchel
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Rainer Schäffer
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | | |
Collapse
|
22
|
Al-Hamdani YS, Rossi M, Alfè D, Tsatsoulis T, Ramberger B, Brandenburg JG, Zen A, Kresse G, Grüneis A, Tkatchenko A, Michaelides A. Properties of the water to boron nitride interaction: From zero to two dimensions with benchmark accuracy. J Chem Phys 2017; 147:044710. [DOI: 10.1063/1.4985878] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yasmine S. Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Mariana Rossi
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Theodoros Tsatsoulis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Wien, Austria
| | - Jan Gerit Brandenburg
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AH, United Kingdom
| | - Andrea Zen
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Wien, Austria
| | - Andreas Grüneis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
23
|
Dubecký M. Noncovalent Interactions by Fixed-Node Diffusion Monte Carlo: Convergence of Nodes and Energy Differences vs Gaussian Basis-Set Size. J Chem Theory Comput 2017; 13:3626-3635. [DOI: 10.1021/acs.jctc.7b00537] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701
03 Ostrava, Czech Republic
- ATRI, Faculty of Materials
Science and Technology, Slovak University of Technology, Paulínska
16, 917 24 Trnava, Slovakia
| |
Collapse
|
24
|
|
25
|
Zen A, Roch LM, Cox SJ, Hu XL, Sorella S, Alfè D, Michaelides A. Toward Accurate Adsorption Energetics on Clay Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2016; 120:26402-26413. [PMID: 27917256 PMCID: PMC5126707 DOI: 10.1021/acs.jpcc.6b09559] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 10/31/2016] [Indexed: 05/31/2023]
Abstract
Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.
Collapse
Affiliation(s)
- Andrea Zen
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Loïc M. Roch
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Xiao Liang Hu
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Sandro Sorella
- SISSA−International
School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- INFM
Democritos National Simulation Center, 34151 Trieste, Italy
| | - Dario Alfè
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United
Kingdom
| | - Angelos Michaelides
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
26
|
Alavi S, Ohmura R. Understanding decomposition and encapsulation energies of structure I and II clathrate hydrates. J Chem Phys 2016; 145:154708. [DOI: 10.1063/1.4964673] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Azadi S, Cohen RE. Low-pressure phase diagram of crystalline benzene from quantum Monte Carlo. J Chem Phys 2016. [DOI: 10.1063/1.4960434] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Azadi S, Cohen RE. Chemical accuracy from quantum Monte Carlo for the benzene dimer. J Chem Phys 2016; 143:104301. [PMID: 26374029 DOI: 10.1063/1.4930137] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
Collapse
Affiliation(s)
- Sam Azadi
- Department of Earth Science and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - R E Cohen
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom and Extreme Materials Initiative, Geophysical Laboratory, Carnegie Institution of Washington, Washington, D.C. 20015, USA
| |
Collapse
|
29
|
Zeng Q, Li J, Huang H, Wang X, Yang M. Polarization response of clathrate hydrates capsulated with guest molecules. J Chem Phys 2016; 144:204308. [PMID: 27250307 DOI: 10.1063/1.4952417] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Clathrate hydrates are characterized by their water cages encapsulating various guest atoms or molecules. The polarization effect of these guest-cage complexes was studied with combined density functional theory and finite-field calculations. An addition rule was noted for these systems whose total polarizability is approximately equal to the polarizability sum of the guest and the cage. However, their distributional polarizability computed with Hirshfeld partitioning scheme indicates that the guest-cage interaction has considerable influence on their polarization response. The polarization of encapsulated guest is reduced while the polarization of water cage is enhanced. The counteraction of these two opposite effects leads to the almost unchanged total polarizability. Further analysis reveals that the reduced polarizability of encapsulated guest results from the shielding effect of water cage against the external field and the enhanced polarizability of water cage from the enhanced bonding of hydrogen bonds among water molecules. Although the charge transfer through the hydrogen bonds is rather small in the water cage, the polarization response of clathrate hydrates is sensitive to the changes of hydrogen bonding strength. The guest encapsulation strengthens the hydrogen bonding network and leads to enhanced polarizability.
Collapse
Affiliation(s)
- Qun Zeng
- Institute of Chemical Materials, China Academy of Engineering Physics (CAEP), Mianyang 621900, China
| | - Jinshan Li
- Institute of Chemical Materials, China Academy of Engineering Physics (CAEP), Mianyang 621900, China
| | - Hui Huang
- Institute of Chemical Materials, China Academy of Engineering Physics (CAEP), Mianyang 621900, China
| | - Xinqin Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| | - Mingli Yang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| |
Collapse
|
30
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
31
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
32
|
Gillan MJ, Alfè D, Manby FR. Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations. J Chem Phys 2015; 143:102812. [PMID: 26374005 DOI: 10.1063/1.4926444] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.
Collapse
Affiliation(s)
- M J Gillan
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - D Alfè
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - F R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
33
|
Teeratchanan P, Hermann A. Computational phase diagrams of noble gas hydrates under pressure. J Chem Phys 2015; 143:154507. [DOI: 10.1063/1.4933371] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Pattanasak Teeratchanan
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| | - Andreas Hermann
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| |
Collapse
|
34
|
Azadi S, Foulkes WMC. Systematic study of finite-size effects in quantum Monte Carlo calculations of real metallic systems. J Chem Phys 2015; 143:102807. [DOI: 10.1063/1.4922619] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Sam Azadi
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - W. M. C. Foulkes
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| |
Collapse
|
35
|
Quigley D, Alfè D, Slater B. Communication: On the stability of ice 0, ice i, and I(h). J Chem Phys 2015; 141:161102. [PMID: 25362263 DOI: 10.1063/1.4900772] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I(h). The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.
Collapse
Affiliation(s)
- D Quigley
- Department of Physics and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - D Alfè
- Department of Earth Sciences and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - B Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
36
|
Gillan MJ. Many-body exchange-overlap interactions in rare gases and water. J Chem Phys 2014; 141:224106. [DOI: 10.1063/1.4903240] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|