1
|
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]
|
2
|
Michalis VK, Economou IG, Stubos AK, Tsimpanogiannis IN. Phase equilibria molecular simulations of hydrogen hydrates via the direct phase coexistence approach. J Chem Phys 2022; 157:154501. [PMID: 36272800 DOI: 10.1063/5.0108738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report the three-phase (hydrate-liquid water-vapor) equilibrium conditions of the hydrogen-water binary system calculated with molecular dynamics simulations via the direct phase coexistence approach. A significant improvement of ∼10.5 K is obtained in the current study, over earlier simulation attempts, by using a combination of modifications related to the hydrogen model that include (i) hydrogen Lennard-Jones parameters that are a function of temperature and (ii) the water-guest energy interaction parameters optimized further by using the Lorentz-Berthelot combining rules, based on an improved description of the solubility of hydrogen in water.
Collapse
Affiliation(s)
| | - Ioannis G Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Athanasios K Stubos
- Environmental Research Laboratory, National Center for Scientific Research "Demokritos," 15310 Aghia Paraskevi Attikis, Greece
| | - Ioannis N Tsimpanogiannis
- Chemical Process and Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH), 57001 Thermi-Thessaloniki, Greece
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
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
|
8
|
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
|
9
|
Shilpi V, Kaur SP, Ramachandran CN. Density functional studies of fused dodecahedral and irregular-dodecahedral water cages. RSC Adv 2015. [DOI: 10.1039/c5ra13268a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
The fused cages of dodecahedral and irregular-dodecahedral water cages with the maximum number of t1d hydrogen bonds were studied using the dispersion corrected density functional method.
Collapse
Affiliation(s)
- V. Shilpi
- Department of Chemistry
- Indian Institute of Technology Roorkee
- Roorkee
- India-247667
| | - Surinder Pal Kaur
- Department of Chemistry
- Indian Institute of Technology Roorkee
- Roorkee
- India-247667
| | - C. N. Ramachandran
- Department of Chemistry
- Indian Institute of Technology Roorkee
- Roorkee
- India-247667
| |
Collapse
|
10
|
Moustafa SG, Schultz AJ, Kofke DA. Effects of Finite Size and Proton Disorder on Lattice-Dynamics Estimates of the Free Energy of Clathrate Hydrates. Ind Eng Chem Res 2014. [DOI: 10.1021/ie504008h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sabry G. Moustafa
- Department
of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| | - Andrew J. Schultz
- Department
of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| | - David A. Kofke
- Department
of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| |
Collapse
|
11
|
Papadimitriou NI, Tsimpanogiannis IN, Economou IG, Stubos AK. Influence of combining rules on the cavity occupancy of clathrate hydrates by Monte Carlo simulations. Mol Phys 2014. [DOI: 10.1080/00268976.2014.902136] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
12
|
|
13
|
Albertí M, Pirani F, Laganà A. Carbon Dioxide Clathrate Hydrates: Selective Role of Intermolecular Interactions and Action of the SDS Catalyst. J Phys Chem A 2013; 117:6991-7000. [DOI: 10.1021/jp3126158] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- M. Albertí
- IQTCUB, Departament
de Química
Física, Universitat de Barcelona, Barcelona, Spain
| | - F. Pirani
- Dipartimento
di Chimica, Università di Perugia, Perugia, Italy
| | - A. Laganà
- Dipartimento
di Chimica, Università di Perugia, Perugia, Italy
| |
Collapse
|
14
|
Tsimpanogiannis IN, Papadimitriou NI, Stubos AK. On the limitation of the van der Waals-Platteeuw-based thermodynamic models for hydrates with multiple occupancy of cavities. Mol Phys 2012. [DOI: 10.1080/00268976.2012.666278] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
15
|
Albertí M, Costantini A, Laganá A, Pirani F. Are Micelles Needed to Form Methane Hydrates in Sodium Dodecyl Sulfate Solutions? J Phys Chem B 2012; 116:4220-7. [DOI: 10.1021/jp301124z] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- M. Albertí
- IQTCUB, Departament
de Química
Física, Universitat de Barcelona, Barcelona, Spain
| | - A. Costantini
- Dipartimento
di Chimica, Università di Perugia, Perugia, Italy
| | - A. Laganá
- Dipartimento
di Chimica, Università di Perugia, Perugia, Italy
| | - F. Pirani
- Dipartimento
di Chimica, Università di Perugia, Perugia, Italy
| |
Collapse
|
16
|
Srivastava HK, Sastry GN. Viability of Clathrate Hydrates as CO2 Capturing Agents: A Theoretical Study. J Phys Chem A 2011; 115:7633-7. [DOI: 10.1021/jp203599g] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Hemant Kumar Srivastava
- Molecular Modelling Group, Indian Institute of Chemical Technology, Tarnaka, Hyderabad 500 607, India
| | - G. Narahari Sastry
- Molecular Modelling Group, Indian Institute of Chemical Technology, Tarnaka, Hyderabad 500 607, India
| |
Collapse
|
17
|
Zhang H, Singer SJ. Analysis of the subcritical carbon dioxide-water interface. J Phys Chem A 2011; 115:6285-96. [PMID: 21438534 DOI: 10.1021/jp200190k] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We follow the evolution of the H(2)O/CO(2) interface at 300 K from the low pressure limit to near-critical pressures in molecular dynamics simulations using the SPC water and EPM2 carbon dioxide models. The intrinsic structure of the interface is elucidated by accumulating density profiles relative to the fluctuating capillary wave surface. Our main finding is that a carbon dioxide film of increasing density and thickness grows in two stages at the interface while the structure of the water surface barely changes. At low density, the entire film density profile grows linearly with the bulk CO(2) density. This regime continues up to a bulk CO(2) density of roughly 0.00095 Å(-3). At pressures above this point, we observe a distinct second peak in the CO(2) density, along with a tail of excess density that decays exponentially with distance from the interface. The decay length of the exponential tail diverges with increasing CO(2) pressure according to an inverse power law decay. Over the entire range of pressures, the CO(2) film had no detectable effect on the orientational order of the water surface. As expected, when the film of excess CO(2) at the interface grows, we find that the surface tension drops with increasing pressure. This is in qualitative accord with existing measurements, although the rate at which the surface tension falls with increasing pressure according to the SPC and EPM2 models is too small, indicating that the surface excess of CO(2) is underestimated by these models.
Collapse
Affiliation(s)
- Hui Zhang
- Department of Chemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | |
Collapse
|
18
|
Lenz A, Ojamäe L. Structures of the I-, II- and H-Methane Clathrates and the Ice−Methane Clathrate Phase Transition from Quantum-Chemical Modeling with Force-Field Thermal Corrections. J Phys Chem A 2011; 115:6169-76. [DOI: 10.1021/jp111328v] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Annika Lenz
- Department of Chemistry, IFM, Linköping University, SE-58183 Linköping, Sweden
| | - Lars Ojamäe
- Department of Chemistry, IFM, Linköping University, SE-58183 Linköping, Sweden
| |
Collapse
|
19
|
Brodskaya EN, Sizov VV. Molecular simulation of gas hydrate nanoclusters in water shell: Structure and phase transitions. COLLOID JOURNAL 2009. [DOI: 10.1134/s1061933x09050020] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
20
|
Hartke B. Methane-water clusters under pressure: Are clathrate cages optimal clusters? J Chem Phys 2009; 130:024905. [DOI: 10.1063/1.3058479] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
21
|
Papadimitriou N, Tsimpanogiannis I, Papaioannou A, Stubos A. Monte Carlo study of sII and sH argon hydrates with multiple occupancy of cages. MOLECULAR SIMULATION 2008. [DOI: 10.1080/08927020802101734] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- N.I. Papadimitriou
- a Environmental Research Laboratory , National Center for Scientific Research ‘Demokritos’ , Agia Paraskevi, 15310, Greece
| | - I.N. Tsimpanogiannis
- a Environmental Research Laboratory , National Center for Scientific Research ‘Demokritos’ , Agia Paraskevi, 15310, Greece
| | - A.Th. Papaioannou
- b School of Chemical Engineering, National Technical University of Athens , Zografou, 15780, Greece
| | - A.K. Stubos
- a Environmental Research Laboratory , National Center for Scientific Research ‘Demokritos’ , Agia Paraskevi, 15310, Greece
| |
Collapse
|
22
|
Whitman CA, Mysyk R, White MA. Investigation of factors affecting crystallization of cyclopentane clathrate hydrate. J Chem Phys 2008; 129:174502. [DOI: 10.1063/1.3005379] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Alavi S, Dornan P, Woo TK. Determination of NMR lineshape anisotropy of guest molecules within inclusion complexes from molecular dynamics simulations. Chemphyschem 2008; 9:911-9. [PMID: 18386265 DOI: 10.1002/cphc.200700805] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Nonspherical cages in inclusion compounds can result in non-uniform motion of guest species in these cages and anisotropic lineshapes in NMR spectra of the guest. Herein, we develop a methodology to calculate lineshape anisotropy of guest species in cages based on molecular dynamics simulations of the inclusion compound. The methodology is valid for guest atoms with spin 1/2 nuclei and does not depend on the temperature and type of inclusion compound or guest species studied. As an example, the nonspherical shape of the structure I (sI) clathrate hydrate large cages leads to preferential alignment of linear CO(2) molecules in directions parallel to the two hexagonal faces of the cages. The angular distribution of the CO(2) guests in terms of a polar angle theta and azimuth angle phi and small amplitude vibrational motions in the large cage are characterized by molecular dynamics simulations at different temperatures in the stability range of the CO(2) sI clathrate. The experimental (13)C NMR lineshapes of CO(2) guests in the large cages show a reversal of the skew between the low temperature (77 K) and the high temperature (238 K) limits of the stability of the clathrate. We determine the angular distributions of the guests in the cages by classical MD simulations of the sI clathrate and calculate the (13)C NMR lineshapes over a range of temperatures. Good agreement between experimental lineshapes and calculated lineshapes is obtained. No assumptions regarding the nature of the guest motions in the cages are required.
Collapse
Affiliation(s)
- Saman Alavi
- Centre for Catalysis Research and Innovation, Department of Chemistry, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada.
| | | | | |
Collapse
|
24
|
Alavi S, Ripmeester JA, Klug DD. Molecular dynamics study of the stability of methane structure H clathrate hydrates. J Chem Phys 2007; 126:124708. [PMID: 17411153 DOI: 10.1063/1.2710261] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Molecular dynamics simulations are used to study the stability of structure H (sH) methane clathrate hydrates in a 3 x 3 x 3 sH unit cell replica. Simulations are performed at experimental conditions of 300 K and 2 GPa for three methane intermolecular potentials. The five small cages of the sH unit cell are assigned methane guest occupancies of one and large cage guest occupancies of one to five are considered. Radial distribution functions, unit cell volumes, and configurational energies are studied as a function of large cage CH(4) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies. Large cage occupancy of five is the most stable configuration for a Lennard-Jones united-atom potential and the Tse-Klein-McDonald potential parametrized for condensed methane phases and two for the most stable configuation for the Murad and Gubbins potential.
Collapse
Affiliation(s)
- Saman Alavi
- Department of Chemistry, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada.
| | | | | |
Collapse
|