Theoretical studies of host-guest interaction in gas hydrates

Ab initio calculations and atoms-in-molecules (AIM) analysis have been used to investigate the host-guest interaction in dodecahedral water cages using a variety of guest species that include monatomic (He, Ne, Ar, Kr, and Xe), diatomic (CO, H(2), N(2), O(2), and NO), triatomic (CO(2), NO(2), and O(3)), and polyatomic (CH(4) and NH(3)) molecules. Geometry optimization for the guest species, host cage, and their complexes was carried out using the second order M?ller-Plesset perturbation method with the 6-31G** basis set. Single point energy calculations using the same method but different basis sets (6-31++G**, 6-311++G**, aug-cc-pVDZ, and aug-cc-pVTZ) were carried out for the MP2/6-31G** optimized geometries. The interaction energy between the guest species and the host cage has been obtained in the complete basis set limit by basis set extrapolation.     

Molecular dynamics simulations for the growth of CH_4-CO_2 mixed hydrate

Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xCO2= 75%, xCO2= 50%, and xCO2= 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.     

Molecular dynamics simulations for the growth of CH4-CO2 mixed hydrate

Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.     

Stability and reactivity of methane clathrate hydrates: insights from density functional theory

The sI methane clathrate hydrate consists of methane gas molecules encapsulated as dodecahedron (5(12)CH(4)) and tetrakaidecahedron (5(12)6(2)CH(4)) water cages. The characterization of the stability of these cages is crucial to an understanding of the mechanism of their formation. In the present work, we perform calculations using density functional theory to calculate interaction energies, free energies, and reactivity indices of these cages. The contributions from polarization functions to interaction energies is more than diffuse functions from Pople basis sets, though both functions from the correlation-consistent basis sets contribute significantly to interaction energies. The interaction energies and free energies show that the formation of the 5(12)CH(4) cage (from the 5(12) cage) is more favored compared to the 5(12)6(2)CH(4) cage (from the 5(12)6(2) cage). The pressure-dependent study shows a spontaneous formation of the 5(12)CH(4) cage at 273 K (P ≥ 77 bar) and the 5(12)6(2)CH(4) cage (P = 100 bar). The reactivity of the 5(12)CH(4) cage is similar to that of the 5(12) cage, but the 5(12)6(2)CH(4) cage is more reactive than the 5(12)6(2) cage.     

Isomer stability of N6C6H6 cages

Recent theoretical studies have identified carbon-nitrogen cages that are potentially stable high energy density materials (HEDM). One such molecule is an N(6)C(6)H(6) cage in which a six-membered ring of nitrogen is bonded to C(3)H(3) triangles on both sides. This molecule is based on the structure of the most stable N(12) cage, with six carbon atoms substituted into the structure. In the current study, several N(6)C(6)H(6) isomers (including the previously studied cage) are examined by theoretical calculations to determine which is actually the most stable. Stability will be evaluated from two points of view: (1) thermodynamic stability of one isomer versus another and (2) kinetic stability of each isomer as determined by the energetics of bond breaking. Density functional theory (B3LYP), perturbation theory (MP2 and MP4), and coupled-cluster theory (CCSD(T)) are used in this study, along with the correlation-consistent basis sets of Dunning. Trends in thermodynamic and kinetic stability are discussed.     

Molecular-dynamics simulations of binary structure II hydrogen and tetrahydrofurane clathrates

The binary structure II hydrogen and tetrahydrofurane (THF) clathrates are studied with molecular-dynamics simulations. Simulations are done at pressures of 120 and 1.013 bars for temperatures ranging from 100 to 273 K. For the small cages of the structure II unit cell, H2 guest molecule occupancies of 0, 16 (single occupancy), and 32 (double occupancy) are considered. THF occupancies of 0-8 in the large cages are studied. For cases in which THF does not occupy all large cages in a unit cell, the remaining large cages can be occupied with sets of four H2 guest molecules. The unit-cell volumes and configurational energies are compared in the different occupancy cases. Increasing the small cage occupancy leads to an increase in the unit-cell volume and thermal-expansion coefficient. Among simulations with the same small cage occupancy, those with the large cages containing 4H2 guests have the largest volumes. The THF guest molecules have a stabilizing effect on the clathrate and the configurational energy of the unit cell decreases linearly as the THF content increases. For binary THF + H2 clathrates, the substitution of the THF molecules in the large cages with sets of 4H2 molecules increases the configurational energy. For the binary clathrates, various combinations of THF and H2 occupancies have similar configurational energies.     

Cage occupancies in the high pressure structure H methane hydrate: a neutron diffraction study

A neutron diffraction study was performed on the CD(4) : D(2)O structure H clathrate hydrate to refine its CD(4) fractional cage occupancies. Samples of ice VII and hexagonal (sH) methane hydrate were produced in a Paris-Edinburgh press and in situ neutron diffraction data collected. The data were analyzed with the Rietveld method and yielded average cage occupancies of 3.1 CD(4) molecules in the large 20-hedron (5(12)6(8)) cages of the hydrate unit cell. Each of the pentagonal dodecahedron (5(12)) and 12-hedron (4(3)5(6)6(3)) cages in the sH unit cell are occupied with on average 0.89 and 0.90 CD(4) molecules, respectively. This experiment avoided the co-formation of Ice VI and sH hydrate, this mixture is more difficult to analyze due to the proclivity of ice VI to form highly textured crystals, and overlapping Bragg peaks of the two phases. These results provide essential information for the refinement of intermolecular potential parameters for the water-methane hydrophobic interaction in clathrate hydrates and related dense structures.     

Molecular dynamics study of structure H clathrate hydrates of methane and large guest molecules

Methane storage in structure H (sH) clathrate hydrates is attractive due to the relatively higher stability of sH as compared to structure I methane hydrate. The additional stability is gained without losing a significant amount of gas storage density as happens in the case of structure II (sII) methane clathrate. Our previous work has showed that the selection of a specific large molecule guest substance (LMGS) as the sH hydrate former is critical in obtaining the optimum conditions for crystallization kinetics, hydrate stability, and methane content. In this work, molecular dynamics simulations are employed to provide further insight regarding the dependence of methane occupancy on the type of the LMGS and pressure. Moreover, the preference of methane molecules to occupy the small (5(12)) or medium (4(3)5(6)6(3)) cages and the minimum cage occupancy required to maintain sH clathrate mechanical stability are examined. We found that thermodynamically, methane occupancy depends on pressure but not on the nature of the LMGS. The experimentally observed differences in methane occupancy for different LMGS may be attributed to the differences in crystallization kinetics and/or the nonequilibrium conditions during the formation. It is also predicted that full methane occupancies in both small and medium clathrate cages are preferred at higher pressures but these cages are not fully occupied at lower pressures. It was found that both small and medium cages are equally favored for occupancy by methane guests and at the same methane content, the system suffers a free energy penalty if only one type of cage is occupied. The simulations confirm the instability of the hydrate when the small and medium cages are empty. Hydrate decomposition was observed when less than 40% of the small and medium cages are occupied.     

Self-assembly of discrete M6L8 coordination cages based on a conformationally flexible tripodal phosphoric triamide ligand

Reactions of a tripodal ligand, N,N',N″-tris(3-pyridinyl)phosphoric triamide (TPPA), and a series of transition-metal ions result in the assembly of five discrete M(6)L(8) coordination cages [M(6)(TPPA)(8)(H(2)O)(12)](ClO(4))(12)·57H(2)O [M = Ni(2+) (1), Co(2+) (2), Zn(2+) (3), Cd(2+) (4)] and [Pd(6)(TPPA)(8)]Cl(12)·22H(2)O (5). X-ray structural analyses reveal that the cages have large internal cavities and flexible windows. The flexible ligand TPPA adopts the syn conformation in cages 1-4, but it transforms to the anti conformation in cage 5. Because of the conformational transformation, the sizes of the windows and the volume of the internal cavity of cage 5 are increased. (1)H NMR and electrospray mass spectrometric studies show that cage 5 maintains its structural integrity in solution. Additionally, compounds 3 and 4 exhibit strong blue fluorescent emissions, which are 1 order of magnitude higher than that of the free ligand.     

The cages, dynamics, and structuring of incipient methane clathrate hydrates

Interest in describing clathrate hydrate formation mechanisms spans multiple fields of science and technical applications. Here, we report findings from multiple molecular dynamics simulations of spontaneous methane clathrate hydrate nucleation and growth from fully demixed and disordered two-phase fluid systems of methane and water. Across a range of thermodynamic conditions and simulation geometries and sizes, a set of seven cage types comprises approximately 95% of all cages formed in the nucleated solids. This set includes the ubiquitous 5(12) cage, the 5(12)6(n) subset (where n ranges from 2-4), and the 4(1)5(10)6(n) subset (where n also ranges from 2-4). Transformations among these cages occur via water pair insertions/removals and rotations, and may elucidate the mechanisms of solid-solid structural rearrangements observed experimentally. Some consistency is observed in the relative abundance of cages among all nucleation trajectories. 5(12) cages are always among the two most abundant cage types in the nucleated solids and are usually the most abundant cage type. In all simulations, the 5(12)6(n) cages outnumber their 4(1)5(10)6(n) counterparts with the same number of water molecules. Within these consistent features, some stochasticity is observed in certain cage ratios and in the long-range ordering of the nucleated solids. Even when comparing simulations performed at the same conditions, some trajectories yield swaths of multiple adjacent sI unit cells and long-range order over 5 nm, while others yield only isolated sI unit cells and little long-range order. The nucleated solids containing long-range order have higher 5(12)6(2)/5(12) and 5(12)6(3)/4(1)5(10)6(2) cage ratios when compared to systems that nucleate with little long-range order. The formation of multiple adjacent unit cells of sI hydrate at high driving forces suggests an alternative or addition to the prevailing hydrate nucleation hypotheses which involve formation through amorphous intermediates.     

NMR shielding constants for hydrogen guest molecules in structure II clathrates

Proton NMR shielding constants and chemical shifts for hydrogen guests in small and large cages of structure II clathrates are calculated using density-functional theory and the gauge-invariant atomic-orbital method. Shielding constants are calculated at the B3LYP level with the 6-311++G(d,p) basis set. The calculated chemical shifts are corrected with a linear regression to reproduce the experimental chemical shifts of a set of standard molecules. The calculated chemical shifts of single hydrogen molecules in the small and large structure II cages are 4.94 and 4.84 ppm, respectively, which show that within the error range of the method the H2 guest molecules in the small and large cages cannot be distinguished. Chemical shifts are also calculated for double occupancy of the hydrogen guests in small cages, and double, triple, and quadruple occupancy in large cages. Multiple occupancy changes the chemical shift of the hydrogen guests by approximately 0.2 ppm. The relative effects of other guest molecules and the cage on the chemical shift are studied for the cages with multiple occupancies.     

Selective formation of a self-assembling homo or hetero cavitand cage via metal coordination based on thermodynamic or kinetic control

The selective formation of a homo or hetero cavitand cage composed of two molecules of tetra(4-pyridyl)-cavitand (1), tetrakis(4-cyanophenyl)-cavitand (2), or tetrakis(4-pyridylethynyl)-cavitand (3), and four molecules of Pd(dppp)(OTf)(2) (4) or Pt(dppp)(OTf)(2) (5) has been studied. A 1:1:4 mixture of 1 with more steric restriction, 2 with less coordination ability, and 4 or 5 specifically self-assembled into a hetero cavitand cage 6 or 7, respectively. In contrast, a 1:1:4 mixture of 2, 3, and 4 in CDCl(3) at room temperature assembled into the most labile homo cyanophenyl cavitand cage 8 and the most stable homo pyridylethynyl cavitand cage 9 in a 1:1 ratio. Upon heating at 50 degrees C, the thermodynamic equilibrium was shifted to a 1:1:1 mixture of 8, 9, and a hetero cavitand cage 10. When 1 equiv of 3 was added to 8 at room temperature, 8, 9, and 10 were formed initially in a 1:1:3 ratio and finally shifted to a 1:1:1 ratio. In the Pt-system, upon addition of 1 equiv of 3 to homo cyanophenyl cavitand cage 11 in CDCl(3) at room temperature, the ratio of hetero to homo cavitand cage (13/12) initially attained was 8.7 and remained above 5.6 at room temperature. Upon heating at 50 degrees C, 13 was finally converted to 11 and 12. Thus, the selectivity for the self-assembly of the homo or hetero cavitand cage is controlled by the balance between kinetic and thermodynamic stabilities of cages based on a combination of factors such as coordination ability and steric demand of the cavitands.     

How much carbon dioxide can be stored in the structure H clathrate hydrates?: a molecular dynamics study

The stability of structure H (sH) carbon dioxide clathrate hydrates at three temperature-pressure conditions are determined by molecular dynamics simulations on a 3x3x3 sH unit cell replica. Simulations are performed at 100 K at ambient pressure, 273 K at 100 bars and also 300 K and 5.0 kbars. The small and medium cages of the sH unit cell are occupied by a single carbon dioxide guest and large cage guest occupancies of 1-5 are considered. Radial distribution functions are given for guests in the large cages and unit cell volumes and configurational energies are studied as a function of large cage CO(2) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at three temperature/pressure conditions stated above. At the low temperature, large cage occupancy of 5 is the most stable while at the higher temperature, the occupancy of 3 is the most favored. Calculations are also performed to show that the CO(2) sH clathrate is more stable than the methane clathrate analog. Implications on CO(2) sequestration by clathrate formation are discussed.     

Phase equilibrium measurements and crystallographic analyses on structure-H type gas hydrate formed from the CH4-CO2-neohexane-water system

Phase equilibrium conditions and the crystallographic properties of structure-H type gas hydrates containing various amounts of methane (CH4), carbon dioxide (CO2), neohexane (2,2-dimethylbutane; NH), and liquid water were investigated. When the CH4 concentration was as high as approximately 70%, the phase equilibrium pressure of the structure-H hydrate, which included NH, was about 1 MPa lower at a given temperature than that of the structure-I hydrate with the same composition (except for a lack of NH). However, as the CO2 concentration increased, the pressure difference between the structures became smaller and, at CO2 concentrations below 50%, the phase equilibrium line for the structure-H hydrate crossed that for the structure I. This cross point occurred at a lower temperature at higher CO2 concentration. Extrapolating this relation between the cross point and the CO2 concentration to 100% CO2 suggests that the cross-point temperature would be far below 273.2 K. It is then difficult to form structure-H hydrates in the CO2-NH-liquid water system. To examine the structure, guest composition, and formation process of structure-H hydrates at various CH4-CO2 compositions, we used the methods of Raman spectroscopy, X-ray diffraction, and gas chromatography. Raman spectroscopic analyses indicated that the CH4 molecules were found to occupy both 5(12) and 4(3)5(6)6(3) cages, but they preferably occupied only the 5(12) cages. On the other hand, the CO2 molecules appeared to be trapped only in the 4(3)5(6)6(3) cages. Thus, the CO2 molecules aided the formation of structure-H hydrates even though they reduced the stability of that structure. This encaged condition of guest molecules was also compared with the theoretical calculations. In the batch-type reactor, this process may cause the fractionation of the remaining vapor composition in the opposite sense as that for CH4-CO2 hydrate (structure-I), and thus may result in an alternating formation of structure-H hydrates and structure-I in the same batch-type reactor.     

p-Xylene-selective metal-organic frameworks: a case of topology-directed selectivity

Para-disubstituted alkylaromatics such as p-xylene are preferentially adsorbed from an isomer mixture on three isostructural metal-organic frameworks: MIL-125(Ti) ([Ti(8)O(8)(OH)(4)(BDC)(6)]), MIL-125(Ti)-NH(2) ([Ti(8)O(8)(OH)(4)(BDC-NH(2))(6)]), and CAU-1(Al)-NH(2) ([Al(8)(OH)(4)(OCH(3))(8)(BDC-NH(2))(6)]) (BDC = 1,4-benzenedicarboxylate). Their unique structure contains octahedral cages, which can separate molecules on the basis of differences in packing and interaction with the pore walls, as well as smaller tetrahedral cages, which are capable of separating molecules by molecular sieving. These experimental data are in line with predictions by molecular simulations. Additional adsorption and microcalorimetric experiments provide insight in the complementary role of the two cage types in providing the para selectivity.     

Betaine-induced assembly of neutral infinite columns and chains of linked silver(I) polyhedra with embedded acetylenediide

Ten polymeric silver(I) double salts containing embedded acetylenediide: [(Ag2C2)2(AgCF3CO2)9(L1)3] (1), [(Ag2C2)2(AgCF3CO2)10(L2)3]H2O (2), [(Ag2C2)(AgCF3CO2)4(L3)(H2O)]0.75 H2O (3), [(Ag2C2)(1.5)(AgCF3CO2)7(L4)2] (4), [(Ag2C2)(AgCF3CO2)7(L5)2(H2O)] (5), [(Ag2C2) (AgC2F5CO2)7(L1)3(H2O)] (6), [(Ag2C2)(AgCF3CO2)7(L1)3(H2O)]2 H2O (7), [(Ag2C2)(AgC2F5CO2)6(L3)2] (8), [(Ag2C2)2(AgC2F5CO2)12(L4)2(H2O)4]H2O (9), and [(Ag2C2)(AgCF3CO2)6(L3)2(H2O)]H2O (10) have been isolated by varying the types of betaines, the perfluorocarboxylate ligands employed, and the reaction conditions. Single-crystal X-ray analysis has shown that 1-4 all have a columnar structure composed of fused silver(I) double cages, with C2(2-) species embedded in its stem and an exterior coat comprising anionic and zwitterionic carboxylates. For 5 and 6, single silver(I) cages are linked into a beaded chain through both types of carboxylate ligands. In 7, two different coordination modes of L1 connect the silver(I) polyhedra into a chain. For 8, the mu(2)-O,O' coordination mode of L3 connects the silver(I) double cages into a chain. Compound 9 exhibits a two-dimensional architecture generated from the cross-linkage of double cages by C2F5CO2-, L4, and [Ag2(C2F5CO2)2] units. Similar to 9, 10 is also a two-dimensional structure, which is formed by connecting the chains of linked double cages through [Ag2(CF3CO2)2] bridging.     

Heterogeneous crystal growth of methane hydrate on its sII [001] crystallographic face

This paper presents a systematic molecular simulation study of the heterogeneous crystal growth of methane hydrate sII from supersaturated aqueous methane solutions. The growth of sII hydrate on the [001] crystallographic face is achieved through utilization of a recently proposed methodology, and rates of crystal growth of 1 A/ns were sustained for the molecular models and specific conditions employed in this work. Characteristics of the crystals grown as well as properties and structure of the interface are examined. Water cages with a 5(12)6(3) arrangement, which are improper to both sI and sII structures, are identified during the heterogeneous growth of sII methane hydrate. We show that the growth of a [001] face of sII hydrate can produce an sI crystalline structure, confirming that cross-nucleation of methane hydrate structures is possible. Defects consisting of two methane molecules trapped in large 5(12)6(4) cages and water molecules trapped in small and large cages are observed, where in one instance we have found a large 5(12)6(4) cage containing three water molecules.     

Argentophilicity and solvent-induced structural diversity in double salts of silver acetylide with silver perfluoroalkyl carboxylates

A series of novel double salts of silver(I) were isolated by dissolving Ag(2)C(2) in a concentrated aqueous solution of R(F)CO(2)Ag (R(F) = CF(3), C(2)F(5)) and AgBF(4). Different ancillary solvento ligands such as H(2)O, CH(3)CN, and C(2)H(5)CN were found to affect the crystallization process that led to the assembly of various silver(I) cages with embedded C(2)(2-) ions. 2Ag(2)C(2) x 12CF(3)CO(2)Ag x 5H(2)O (1) consists of two independent C(2)@Ag(7) cages, each having the shape of a basket with a square base. Ag(2)C(2) x 6CF(3)CO(2)Ag x 3CH(3)CN (2) contains a zigzag chain of edge-sharing triangulated dodecahedra, and 4Ag(2)C(2) x 23CF(3)CO(2)Ag x 7C(2)H(5)CN x 2.5H(2)O (3) features an unusual double-walled silver column constructed from the fusion of four different kinds of irregular polyhedra. Ag(2)C(2) x 10C(2)F(5)CO(2)Ag x 9.5H(2)O (4), Ag(2)C(2) x 9C(2)F(5)CO(2)Ag x 3CH(3)CN x H(2)O (5), and Ag(2)C(2) x 6C(2)F(5)CO(2)Ag x 2C(2)H(5)CN (6) all contain an edge-sharing double cage with each single cage in the shape of a square antiprism, a capped square antiprism, and a triangulated dodecahedron, respectively.     

New 13-vertex metallacarborane sandwich compounds; synthetic and structural studies

Reduction of 1,12-closo-C2B10H12 followed by reaction with the appropriate metal halide and metathesis with either [K(18-crown-6)]Br or [BTMA]Cl ([BTMA] = [C6H5CH2N(CH3)3]+) affords isolable salts of the supraicosahedral metallacarborane sandwich anions [4,4-M-(1,10-closo-C2B10H12)2]n- in moderate to good yield. Compounds prepared are [BTMA][4,4-Co-(1,10-closo-C2B10H12)2] ( 1), [K(18-crown-6)][4,4-Co-(1,10-closo-C2B10H12)2] ( 2), [K(18-crown-6)]2[4,4-Ni-(1,10-closo-C2B10H12)2] ( 3), [K(18-crown-6)]2[4,4-Fe-(1,10-closo-C2B10H12)2] ( 4), [BTMA]2[4,4-Fe-(1,10-closo-C2B10H12)2] ( 5) and [K(18-crown-6)]2[4,4-Ti-(1,10-closo-C2B10H12)2] ( 6). Oxidation of the iron(II) species 4 and 5 with FeCl3 in THF generates the iron(III) analogues [K(18-crown-6)][4,4-Fe-(1,10-closo-C2B10H12)2] ( 7) and [BTMA][4,4-Fe-(1,10-closo-C2B10H12)2] ( 8), respectively. All diamagnetic compounds were characterised spectroscopically and the structures of 1, 3, 4, 6, 7 and 8 were established by single crystal X-ray diffraction. All anions have the anticipated cluster structures with two docosahedral 13-vertex cages joined at the central metal atom (the common degree-six vertex 4). Carbon atoms occupy the degree-four vertex 1 and the degree-five vertex 10. 11B NMR spectroscopy suggests the anions have, on the NMR timescale, C2h symmetry in solution at room temperature, consistent with free rotation, or at least substantial libration, of cage units about the long molecular axis. In the solid state the relative conformations of the two cages may be rationalised by simple bonding arguments, the single exception being the conformation of 4, in which both cages are subject to directional B-H...K+ interactions to the [K(18-crown-6)]+ counterion. The salts 3, 6 and 7 also show B-H...K+ interactions but involving one cage only.     

Molecular dynamics study of carbon dioxide hydrate dissociation

We present results from a molecular dynamics study of the dissociation behavior of carbon dioxide (CO(2)) hydrates. We explore the effects of hydrate occupancy and temperature on the rate of hydrate dissociation. We quantify the rate of dissociation by tracking CO(2) release into the liquid water phase as well as the velocity of the hydrate-liquid water interface. Our results show that the rate of dissociation is dependent on the fractional occupancy of each cage type and cannot be described simply in terms of overall hydrate occupancy. Specifically, we find that hydrates with similar overall occupancy differ in their dissociation behavior depending on whether the small or large cages are empty. In addition, individual cages behave differently depending on their surrounding environment. For the same overall occupancy, filled small and large cages dissociate faster in the presence of empty large cages than when empty small cages are present. Therefore, hydrate dissociation is a collective phenomenon that cannot be described by focusing solely on individual cage behavior.