1
|
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
|
2
|
Yanes-Rodríguez R, Prosmiti R. Confining He Atoms in Diverse Ice-Phases: Examining the Stability of He Hydrate Crystals through DFT Approaches. Molecules 2023; 28:7893. [PMID: 38067621 PMCID: PMC10708247 DOI: 10.3390/molecules28237893] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2023] [Revised: 11/24/2023] [Accepted: 11/28/2023] [Indexed: 10/16/2024] Open
Abstract
In the realm of solid water hydrostructures, helium atoms have a tendency to occupy the interstitial spaces formed within the crystal lattice of ice structures. The primary objective of this study is to examine the stability of various ice crystals when influenced by the presence of He atoms. Presenting a first attempt at a detailed computational description of the whole energy components (guest-water, water-water, guest-guest) in the complete crystal unit cells contributes to enhancing the knowledge available about these relatively unexplored helium-water systems, which could potentially benefit future experiments. For this purpose, two different ice structures were considered: the previously established He@ice II system, and the predicted (but currently nonexistent) He@ice XVII system. One of the main features of these He-filled structures is the stability conferred by the weak van der Waals dispersion forces that occur between the host lattice and the guest atoms, in addition to the hydrogen bonds established among the water molecules. Hence, it is crucial to accurately describe these interactions. Therefore, the first part of this research is devoted examining the performance and accuracy of various semi-local and non-local DFT/DFT-D functionals, in comparison with previous experimental and/or high-level computational data. Once the best-performing DFT functional has been identified, the stability of these empty and He-filled structures, including different number of He atoms within the lattices, is analysed in terms of their structural (lattice deformation), mechanical (pressure compression effects) and energetic properties (binding and saturation energies). In this manner, the potential formation of these structures under zero temperature and pressure conditions can be evaluated, while their maximum storage capacity is also determined. The obtained results reveal that, despite the weak underlying interactions, the He encapsulation has a rather notable effect on both lattice parameters and energetics, and therefore, the guest-host interactions are far from being negligible. Besides, both ice crystals are predicted to remain stable when filled with He atoms, with ice XVII exhibiting a higher capacity for accommodating a larger number of guest atoms within its interstitial spaces.
Collapse
Affiliation(s)
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain;
| |
Collapse
|
3
|
Yanes-Rodríguez R, Prosmiti R. Computational investigations of stable multiple-cage-occupancy He clathrate-like hydrostructures. Phys Chem Chem Phys 2023. [PMID: 37314248 DOI: 10.1039/d3cp00603d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
One of the several possibilities offered by the interesting clathrate hydrates is the opportunity to encapsulate several atoms or molecules, in such a way that more efficient storage materials could be explored or new molecules that otherwise do not exist could be created. These types of applications are receiving growing attention from technologists and chemists, given the future positive implications that they entail. In this context, we investigated the multiple cage occupancy of helium clathrate hydrates, to establish stable novel hydrate structures or ones similar to those predicted previously by experimental and theoretical studies. To this purpose, we analyzed the feasibility of including an increased number of He atoms inside the small (D) and large (H) cages of the sII structure through first-principles properly assessed density functional approaches. On the one hand, we have computed energetic and structural properties, in which we examined the guest-host and guest-guest interactions in both individual and two-adjacent clathrate-like sII cages by means of binding and evaporation energies. On the other hand, we have carried out a thermodynamical analysis on the stability of such He-containing hydrostructures in terms of changes in enthalpy, ΔH, Gibbs free energy, ΔG, and entropy, ΔS, during their formation process at various temperature and pressure values. In this way, we have been able to make a comparison with experiments, reaffirming the ability of computational DFT approaches to describe such weak guest-host interactions. In principle, the most stable structure involves the encapsulation of one and four He atoms inside the D and H sII cages, respectively; however, more He atoms could be entrapped under lower temperature and/or higher pressure thermodynamic conditions. We foresee such accurate computational quantum chemistry approaches contributing to the current emerging machine-learning model development.
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
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
4
|
Li M, Chen B, Li K, Song Y, Yang M. Stability and structure of multiply occupied sII CO2 clathrate hydrates: a possibility for carbon capturing. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/31/2023]
|
5
|
Cabrera-Ramírez A, Prosmiti R. Modeling of Structure H Carbon Dioxide Clathrate Hydrates: Guest-Lattice Energies, Crystal Structure, and Pressure Dependencies. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:14832-14842. [PMID: 36110497 PMCID: PMC9465682 DOI: 10.1021/acs.jpcc.2c04140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 07/26/2022] [Indexed: 06/15/2023]
Abstract
We performed first-principles computations to investigate the complex interplay of molecular interaction energies in determining the lattice structure and stability of CO2@sH clathrate hydrates. Density functional theory computations using periodic boundary conditions were employed to characterize energetics and the key structural properties of the sH clathrate crystal under pressure, such as equilibrium lattice volume and bulk modulus. 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. Structural relaxations of the fully CO2-filled and empty sH unit cells yield crystal structure and lattice energies, while their compressibility parameters were derived by including the pressure dependencies. The present quantum chemistry computations suggest anisotropy in the compressibility of the sH clathrate hydrates, with the crystal being less compressible along the a-axis direction than along the c-axis one, in distinction from nearly isotropic sI and sII structures. The detailed results presented here give insight into the complex nature of the underlying guest-host interactions, checking earlier assumptions, providing critical tests, and improving estimates. Such entries may eventually lead to better predictions of thermodynamic properties and formation conditions, with a direct impact on emerging hydrate-based technologies.
Collapse
Affiliation(s)
- 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, 28049, Madrid, Spain
| | - Rita Prosmiti
- Institute
of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
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: 0.7] [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
|
Yanes-Rodríguez R, Prosmiti R. Assessment of DFT approaches in noble gas clathrate-like clusters: stability and thermodynamics. Phys Chem Chem Phys 2021; 24:1475-1485. [PMID: 34935011 DOI: 10.1039/d1cp04935f] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have assessed the performance and accuracy of different wavefunction-based electronic structure methods, such as DFMP2 and domain-based local pair-natural orbital (DLPNO-CCSD(T)), as well as a variety of density functional theory (DFT) approaches on He@(H2O)N cage systems. We have selected representative clathrate-like structures corresponding to the building blocks present in each of the sI, sII and sH natural gas clathrate hydrates, and we have carefully studied the interaction between a He atom with each of their individual cages. We reported well-converged DFMP2 and DLPNO-CCSD(T) reference data, together with interaction and cohesive energies of four different density functionals (two GGA, revPBE and PW86PBE, and two hybrids, B3LYP and PBE0), including diverse dispersion correction schemes (D3(0), D3(BJ), D4 and XDM) for both He-filled and empty clathrate-like cages. After the analysis of the results, we came to the conclusion that the PW86PBE functional, with both XDM and D4 corrections, and the PBE0-D4 functional present reasonably adequate approaches to describe the guest-host noncovalent interactions that take place in such He clathrate hydrates. Taking into account that the He@sII is the only helium clathrate that scientists have been able to synthesize recently, we have performed a thermodynamic study on the individual 512 and 51264 cages present in the sII crystal. We determined the change in enthalpy, ΔH, and in Gibbs free energy, ΔG, at various temperatures and pressures, and we found out that in the range of experimental conditions the reactions associated with the encapsulation of the He atom inside the cages are exothermic and spontaneous. Finally, we highlighted the importance of an accurate description of the interaction in He@water mixtures, as a crucial component in construction of reliable data-driven models.
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
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
8
|
Carrillo-Bohórquez O, Valdés Á, Prosmiti R. Encapsulation of a Water Molecule inside C 60 Fullerene: The Impact of Confinement on Quantum Features. J Chem Theory Comput 2021; 17:5839-5848. [PMID: 34420292 PMCID: PMC8444341 DOI: 10.1021/acs.jctc.1c00662] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Indexed: 11/30/2022]
Abstract
We introduce an efficient quantum fully coupled computational scheme within the multiconfiguration time-dependent Hartree (MCTDH) approach to handle the otherwise extremely costly computations of translational-rotational-vibrational states and energies of light-molecule endofullenes. Quantum calculations on energy levels are reported for a water molecule inside C60 fullerene by means of such a systematic approach that includes all nine degrees of freedom of H2O@C60 and does not consider restrictions above them. The potential energy operator is represented as a sum of natural potentials employing the n-mode expansion, along with the exact kinetic energy operator, by introducing a set of Radau internal coordinates for the H2O molecule. On the basis of the present rigorous computations, various aspects of the quantized intermolecular dynamics upon confinement of H2O@C60 are discussed, such as the rotational energy level splitting and the significant frequency shifts of the encapsulated water molecule vibrations. The impact of water encapsulation on quantum features is explored, and insights into the nature of the underlying forces are provided, highlighting the importance of a reliable first-principles description of the guest-host interactions.
Collapse
Affiliation(s)
- Orlando Carrillo-Bohórquez
- Departamento
de Física, Universidad Nacional
de Colombia, Calle 26, Cra 39, 404 Edificio, Bogotá, Colombia
- Institute
of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, 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
|