A Review of Clathrate Hydrate Nucleation - ACS Sustainable

Nov 3, 2017 - He is a recipient of the Industrial Postgraduate Program scholarship from Economic Development Board of Singapore and Lloyd's Register G...
1 downloads 16 Views 2MB Size
Subscriber access provided by READING UNIV

Perspective

A Review of Clathrate Hydrate Nucleation Maninder Khurana, Zhenyuan Yin, and Praveen Linga ACS Sustainable Chem. Eng., Just Accepted Manuscript • DOI: 10.1021/ acssuschemeng.7b03238 • Publication Date (Web): 03 Nov 2017 Downloaded from http://pubs.acs.org on November 4, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Sustainable Chemistry & Engineering is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

A Review of Clathrate Hydrate Nucleation Maninder Khurana1, Zhenyuan Yin1,2, Praveen Linga1* 1

Department of Chemical and Biomolecular Engineering, National University of Singapore, Engineering

Drive 4, Singapore 117585 2

Lloyd's Register Global Technology Centre Pte Ltd, Singapore 138522

Abstract Clathrate hydrates are crucial from the point of view of flow assurance, future energy resource as well as promising innovative and sustainable applications such as gas separation, CO2 sequestration, district and data centre cooling and seawater desalination and natural gas storage. Although proof of concept has been demonstrated, significant progress is necessary in order to achieve industrial level validation and commercialization. Most of the applications possess a common requirement of enhanced kinetics in formation and dissociation. There is a need for a broader understanding of hydrate nucleation mechanisms, cause-effect relations and investigation techniques. The stochastic nature of hydrate nucleation, confounding cause-effect relations and spatial-temporal scales have made it even more challenging to study nucleation. The use of hydrate promoters, novel reactor configurations such as porous media in a packed bed, use of nanoparticles and hydrogels necessitates us to obtain further insights about clathrate nucleation. This review provides an in-depth analysis about the characteristics of clathrate hydrate nucleation and the techniques adopted for studying nucleation from an applicationoriented perspective and enables further development of clathrate technology towards future applications. Keywords: Gas Hydrates; Nucleation; Clathrate Hydrates; Stochastic nature; Molecular simulation; Porous media. *Corresponding author: [email protected] (P. Linga)

ACS Paragon Plus Environment

1

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 89

Introduction Clathrate hydrates are crystalline compounds of water molecules acting as host and guest molecules1-2. Traditionally, natural gas hydrates have been a cause of major problems for the oil and gas industry in production lines, during drilling and in workover operations.3,4 Hammerschmidt in 1934 first discovered that the natural gas pipelines were plugged due to hydrate formation rather than ice which was originally perceived5. Hydrate formation can compromise the structural integrity of the pipelines/surface facilities and can cause disruptions in production.1 Hydrates have also been widely stated as leading source of deep-water flow assurance problem6-7. Gas hydrate is a technology enabler for a number of application in water, energy and environment research domains. Notable applications with gas hydrate as an enabler are gas separation8-15, energy storage16-22, energy transport23-27, cold energy storage28-29, CO2 sequestration30-32 and desalination applications14,

33-36

. These applications of gas hydrates

necessitated the systems to have higher throughput, lower energy consumption, and higher efficiency of separations. Hydrate formation including both, nucleation and growth is extremely crucial step of hydrate-based processes and thorough understanding of nucleation is necessary for the applications. Some applications such as desalination and gas separation requires continuous process operations and hence necessitates precise predictions and knowledge of each step. Clathrate hydrates required high pressure and low temperature for formation. The process of hydrate formation has many similarities with that of crystallization, i.e., it can be divided into a nucleation phase and a growth phase. Three common clathrate structures exist (see Table 1): structure I (cubic); structure II (cubic) and structure H (hexagonal). The structure formed is

ACS Paragon Plus Environment

2

Page 3 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

found to depend upon the size of the guest molecule.37 The pentagonal dodecahedron (512 cage) is the basic building block in the two hydrate structures. The other two types of cages are 51262 and 51264. The other two types of cages observed are attributed to the inability of the dodecahedron to tessellate a 3D space. When two 512 cages are separated by bridging water molecules, they form 51262 cages as in sI. On the other hand, sII structure hydrates are formed when 512 cages share the faces and the resulted gaps are filled by creating 51264 cages. CH4, C2H6 and CO2 form sI hydrate while propane and isobutane tend to form sII hydrate. Smaller molecules such as H2 or N2 also tend to form sII hydrate since they occupy smaller cages that are present in greater fraction in sII hydrates. The crystal structures of sI hydrate and sII hydrate were first determined in the late 1940s and early 1950s by von Stackelberg and co-workers using X-ray diffraction2. The third type of hydrate crystals (sH) are composed of three small 512 cages, two small 12-hedra 435663 cavities and one large 18-hedra 51268 cage.38 Net equation of hydrate formation can be described as:

G+n wH2O

G.nwH2O

(1)

Figure 1 shows the schematic of a typical time evolution of gas consumption for hydrate formation in a semi-batch reactor. The different stages for the process of hydrate nucleation can be observed from the schematic. Since the guest species is present in the gas phase, the process starts with the increase in gas uptake due to dissolution into the liquid phase. After the dissolution phase, the super saturation phase begins when the pressure and temperature thermodynamically favour the formation of hydrates phase but still without the appearance of critical hydrate nuclei. A critical nucleus of hydrate phase can be defined as the minimum amount of new phase that is capable of existing independently. The time interval between establishment of super saturation and the formation critical nuclei is called as the induction time.

ACS Paragon Plus Environment

3

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 89

Induction time has been confirmed experimentally in various literature. It has been shown in independent studies that supersaturation does not always guarantee hydrate formation.39,40 This critical nucleus can then act as the centre for further hydrate growth. Hydrate nucleation process is followed by the catastrophic growth process with rapid increase in the gas uptake. Finally, the gas uptake may end at different levels depending upon the mass transfer resistance and final hydrate formation as shown in Figure 1 by the dotted trajectories. While the thermodynamics of hydrate formation/dissociation is well established, there are lot of unanswered questions with respect to nucleation. Ripmeester and Alavi41 in their recent review suggested that nucleation, decomposition, and the memory effect during reformation are among the most important outstanding issues to be understood regarding clathrate hydrate science. Understanding hydrate nucleation at a molecular level is challenging due to small time/length scales of nucleation and the stochastic nature of the process. Although there have been article focussed towards reviewing hydrate nucleation, but the articles have been targeted towards particular aspects failing to draw a holistic understanding

of clathrate nucleation necessary for

application of hydrate in the fields of gas separation, energy storage and transport. Hence, this work is an attempt to document the various works, techniques and important results along with providing insights into the standing issues for clathrate nucleation.

ACS Paragon Plus Environment

4

Page 5 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Figure 1: Hydrate formation schematic represented by gas uptake in an experiment vs time showing three main phases in hydrate formation process: Dissolution phase, supersaturated phase and growth phase.

Next, we provide some basic background about nucleation and subsequently delve various aspects pertinent to nucleation. Types of nucleation and characteristics There are two types of primary nucleation: Homogeneous nucleation In homogeneous nucleation (HON), the nucleus of hydrate phase emerges directly from the parent phase. Homogeneous nucleation is generally observed in systems without any impurities and considered to be stochastic in nature, i.e., the critical nucleus is formed because of a local thermodynamic fluctuation of the system42. Before the critical radius is attained, the clusters of molecules formed may either grow or shrink as a result of thermodynamic fluctuation. At lower sub-cooling or driving force, near the equilibrium curve, homogeneous nucleation has low probability of occurrence since the critical radius can be substantially high43. For instance, nucleation can take up to thousand years in case of hydrate formation from air inclusion in ice

ACS Paragon Plus Environment

5

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 89

cores44,. Knott et al.45 estimated the homogeneous nucleation rate under realistic conditions as 10-111 nuclei cm-3 s-1 showed by performing MD simulations of methane hydrates. In practical conditions, hydrate nucleation through homogeneous nucleation is very unlikely and proceeds through heterogeneous nucleation. Heterogeneous nucleation In case of heterogeneous nucleation (HEN), the hydrate phase nucleates in contact with a third phase which can be either foreign particle or surface. When we consider the system free energy, it is more favourable to form hydrate on a two dimensional surface than a three dimensional nucleus in bulk water phase. The presence of third phase lowers the interfacial energy necessary to overcome in nucleation phenomenon. Hence, heterogeneous nucleation is more rapid than homogeneous and the critical radius for nucleation is less than homogeneous46-47. Similar effects have been reported for the case of ice nucleation as well48-49. The contact angle between the hydrate and the pre-existing surface controls the reduction in the specific superficial energy of solution-hydrate interface, which, in turn, decreases the amount of work required in the formation of the new phase50. The stabilization of hydrate nuclei in the presence of surfaces is discussed further in Section 9 on nucleation in porous media. In case of dispersed systems such as emulsions, water droplets behave as independent reactors and the nucleation can be to some extent disjoint in each droplet51,52,53. Droplet collision in shared systems or in concentrated emulsions may cause nucleation propagation between droplets and is called as secondary nucleation43. Site of Nucleation To investigate the site for hydrate nucleation, Long and Sloan54 performed a series of experiments on nucleation of natural gas and carbon dioxide in a sapphire tube. They studied the

ACS Paragon Plus Environment

6

Page 7 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

effect of precipitated amorphous silica as hydrate nucleating agent, Sodium dodecyl sulfate (SDS) as a surface inhibitor. Carbon dioxide was chosen to establish the effect of higher guest solubility as compared to natural gas. While the induction times were found to be not predictable, hydrate nucleation was always initiated at an interface, at vapour-water interface in most cases but along the sapphire tube or along the reactor walls in other cases. Higher solubility of carbon dioxide resulted into much faster hydrate growth in the solution. Various other studies in the literature have confirmed the nucleation site as the vapour-liquid interface for both methane and carbon dioxide hydrates55. The preferential hydrate nucleation at the vapour-liquid interface has also been demonstrated by means of Molecular Dynamics (MD) simulations by Moon et al.56. Apart from the lowering of the Gibbs free energy at the interface as stated before, the interface has very high concentration of the host and the guest molecules. The hydrate guest composition at the interface can be as high as 0.15 mole fraction compared to the maximum of 0.001 in aqueous phase2. The higher concentration that exists due to surface adsorption increases the probability of nucleation. Induction Time Induction time is characterised by over-saturation and metastability of the solution. As stated previously, nucleation is a stochastic process. Various studies have shown that hydrate induction times has a significant scatter for experiments performed at constant temperatures57-59. Whereas, the scatter for hydrate formation experiments at constant cooling rate show significantly less scatter2. Further, the degree of stochastic behaviour is strongly affected by the magnitude of driving force. Natarajan et al.60 showed that hydrate induction time is far more reproducible at higher pressures and subsequently formed empirical correlations for induction times as a

ACS Paragon Plus Environment

7

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 89

function of supersaturation ratio. Guo and Rodger61 studied the pre-nucleation stage of methane hydrate formation using Molecular Dynamics (MD) simulations. They studied methane solubility under metastable conditions and found that it can be increased by both reducing temperature and increasing pressure. But, lowering temperature was found more effective for promoting hydrate formation because increasing the pressure was found to reduce the number of water cages. They also calculated critical solubility of methane to be ~1.7 molecules/nm3 beyond which metastable solution form hydrate spontaneously and induction time is not observed. Similar enhanced nucleation has also been reported by Angeles and Firoozabadi62 near fluidfluid spinodal decomposition. For propane-water system, they demonstrated that propane density fluctuations near fluid-fluid spinodal decomposition produce enhanced nucleation rates and is more general than short-range attractive interactions. Nucleation pathways Hydrate nucleation mechanism has been a long-standing area of study and the opinions remain divided on the nucleation pathways. Both time-resolved experiments63-69 and simulations

56, 70-78

have been used in order to capture the pathway and obtain a deeper understanding for the nucleation process. However, due to the difficulty of obtaining direct evidence of nucleation, molecular simulations have been the preferred over experiments to study nucleation pathway. Four major conceptual theories proposed for hydrate nucleation (see Table 2). These are discussed in the subsequent section. Classical Nucleation Theory Classical nucleation theory (CNT) has been well studied for the case of crystallization. CNT was originally derived for condensation of vapour into a liquid in 1920s, it has been later applied to crystallization of supersaturated solutions by analogy79,80. CNT has been commonly used to

ACS Paragon Plus Environment

8

Page 9 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

predict the rate of nucleation and the height of free energy barrier79,

81

. CNT describes the

activation barrier to nucleation as the sum of the increase in free energy due to creation of a new interface and the decrease in the energy from creating a more stable phase. Due to its analytical simplicity, it has been widely applied not just for crystallization but also clathrates. At the end of 19th century, Gibbs developed a thermodynamic description of CNT by defining the free energy change required for cluster formation as the sum of free energy change for the phase transformation and the free energy change for the surface formation.

∆G = ∆GV + ∆Gs

(2)

The phase transformation free energy change is negative and decreases the free energy of the system. Whereas, the solid/liquid interface increase the free energy of the system and the growth of the cluster depends upon the two conflicting terms. Figure 2 provides the free energy diagram for the nucleation process.

Figure 2: Free Energy change for Hydrate nucleation process according to the Classical Nucleation Theory.

ACS Paragon Plus Environment

9

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 89

The positive surface free energy term dominates at the critical radius rc beyond which the total free energy decreases continuously and the growth becomes energetically favourable. The freeenergy barrier corresponds to the increase in the free-energy at the critical radius. The critical radius for hydrate is attained by means of random fluctuations and subsequently, the hydrate growth is spontaneous. Various studies in literature have adopted a unified approach for nucleation and growth by treating nucleation using classical theory and modelling growth phase as a chemical reaction3,82. The classical nucleation approach has two major shortcomings: 1) Macroscopic treatment of hydrate nucleus leading to substantial errors in excess free energy and critical radius estimations.83 2) No insight into the nucleation pathways and the exact structure of hydrate structures. Labile Cluster Hypothesis (LCH) To obtain insights into the nucleation pathways, Sloan and Fleyfel first introduced labile cluster hypothesis (LCH) for hydrate nucleation from ice84. Their work amongst others by MullerBongratz et al.85 and Christiansen and Sloan86. together form the foundation for LCH. According to the LCH, labile ring nucleation structures exist in pure liquid water.84,87 Upon dissolution of the hydrate former gas, critical nucleus is then formed by agglomeration of labile clusters around the guest molecules. Sloan and Fleyfel first analysed two studies reported on reproduction of induction or hydrate metastability. First study was by Barrer and Edge88 on the formation of hydrates of three inert gases namely argon, krypton, and xenon from ice. It was reported that while argon and xenon formed hydrates immediately, krypton had an induction period of about one hour. The induction phenomenon was reported to be similar to the metastability of the subcooled solutions in the absence of a seed crystal. Second study considered was of Falabella et al.89 based on determining the equilibrium and kinetic properties of hydrates of methane, ethane,

ACS Paragon Plus Environment

10

Page 11 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

ethylene, acetylene, carbon dioxide, and krypton at subatmospheric, constant pressures. Falabella et al.89 reported two distinct type of kinetic results: one with a clearly defined induction period for methane and krypton, and second for all other gases without any induction period. The observed induction period for methane and krypton was related to hydrated metastability or primary nucleation. Sloan and Fleyfel84 proposed a new induction nucleation parameter, guest to cavity size ratio. They proposed a kinetic mechanism for hydrate formation with multiple intermediary steps and modelled the system as a set of chemical reactions as shown in the Figure 3. The mechanism starts from the stable species initially as shown in Figure 3(A) and ends with the attainment of critical hydrate radius and beginning of crystal monotonic growth species Figure 3.D. Figure 3.B shows a very rapid transition between different hydrate structures if they exist. Consequentially, two different paths for construction of sI and sII hydrates were proposed. The induction time observed in methane and krypton hydrate was explained on the basis of transition in the two paths proposed. They validated their model using the experimental results of Falabella et al.89 .

ACS Paragon Plus Environment

11

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 89

Figure 3: Proposed kinetic mechanism for hydrate formation from ice by Sloan and Fleyfel84. Skovberg et al.90 were the first study to challenge the LCH. They pointed out the difference in driving force (∆Tsub ) for the CH4 hydrate nucleation experiments was 19 K whereas for C2H6 hydrate was 49 K. Natarajan et al.91 showed similar driving force difference but with a different driving force based on fugacity given below:

Driving Force =

f gV f eq

−1

(3)

Radhakrishnan and Trout71 showed, using Landau free energy surface, that labile clusters proposed by Sloan and Fleyfal84 are easily formed by an activated process only for dilute solutions. For concentrations near the CO2-H2O interface, the free energy penalty of formation of labile clusters is large and hence nucleation cannot occur through labile cluster pathway. Nucleation at interface hypothesis Long et al.92 and Kvamme et al.93 proposed a variation of the labile cluster hypothesis, that the assembly of the labile clusters takes place on the vapour side of the vapour-liquid interface. According to the hypothesis, gas molecules are transported to the interface where they get adsorbed on the aqueous surface. Subsequently, the molecules transport through surface diffusion to a suitable location where water molecules form first partial and then complete cages resulting into labile clusters. Kvamme et al.93 did not attempt to formulate a quantitative estimation of the rate of nucleation. Local Structuring Mechanism Upon providing evidence against LCH, Radhakrishnan and Trout71 proposed a local structuring (see Figure 4) mechanism for nucleation according to which, thermal fluctuation causes the local

ACS Paragon Plus Environment

12

Page 13 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

ordering of CO2 molecules. The local ordering of the guest molecule induces ordering of the host molecules and finally leads to the formation of critical nucleus. Upon attaining the critical nucleus, the water molecules rearrange to form a proper hydrate framework, leading to hydrate crystallization. Their work was based on the following assumptions: 1) Free energy barrier to nucleation remain unaltered by the limited simulation size. 2) Nucleation is governed by equilibrium thermodynamics. 3) Minimum free energy path lies within the chosen parameter space. A more fundamental difference between local structuring mechanism and labile cluster hypothesis is whether water ordering is driven by guest molecule or guest ordering is driven by water molecule.

Local ordering of guest molecules. Structing of water molecules around the ordered strucutre. Driven by thermal fluctuations and stochastic in nature.

Number of gas molecules in structure > critical nucleus; Guest-guest and host-host stable cluster similar to clathrate phase attained

Restructing of the cluster to form clathrate cages

Figure 4: Schematic representing the local structuring mechanism Rodger and co-workers73,56,94 and Walsh et al.75 have shown methane hydrate nucleation in unconstrained atomistic molecular dynamics (MD) simulations under conditions of high driving force. Guo et al.74 showed similar results and called the mechanism as cage adsorption hypothesis. The results were similar to those proposed by LCH, clusters of methane molecules

ACS Paragon Plus Environment

13

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 89

were observed to be surrounded by water molecules. The difference, however, was that the structure was amorphous rather than crystalline as proposed by LCH. To explore this further, Jacobson et al.95 studied the amorphous precursors in the nucleation of clathrate hydrates and proposed a blob formation mechanism discussed next. Blob Formation Mechanism Jacobson et al.

95

proposed the nucleation pathway for hydrophobic guest molecules via blob

formation. Jacobson et al.96 developed a coarse-grained model for molecular dynamics simulations and investigated the mechanisms of nucleation of clathrate hydrates. They introduced blob as a guest-rich amorphous precursor in the nucleation pathway of clathrates of hydrophobic guests. The blob comprises of the same polyhedral cages as crystalline hydrates, but lack their long-range order. In the blob, amorphous clathrate cage continuously forms and dissolve until a cluster of cages reaches a critical size. Figure 5 shows the mechanism of hydrate nucleation via blob formation. The critical size was proposed to be a function of the temperature and at 0.7 Tm sub-cooling it is estimated about 5 cages. After reaching the critical size, spacefilling growth is achieved through face sharing between the cages. During its lifetime, a blob can produce many hydrate crystals until it finally dissolves into the solution. The blobs themselves are not stationary but diffuse in solution. In this perspective, blobs are large analogues of labile clusters. In the literature, similar mechanisms involving amorphous metastable states have been proposed for crystallization of proteins and colloids97,98.

ACS Paragon Plus Environment

14

Page 15 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Figure 5: Schematic for blob formation mechanism as proposed by Jacobson et al.95. This two-step nucleation mechanism of amorphous nuclei to crystalline clathrate have been observed in other studies by Vatamanu et al.76, Liang et al.78,Sarupria et al.99, and He et al.100 using different methodologies. These results contrast the classical nucleation theory according to which nucleation of crystalline phase takes through a build-up of monomer already arranged with the symmetry of the crystal phase. Experimentally, amorphous tetrahydrofuran (THF) hydrate solids have been synthesized under a high pressure of 1.3 GPa101. In their later study, Jacobson et al.96 analysed the effect of guest molecule solubility on nucleation pathway. They studied the nucleation pathways for guest molecules varying over two orders of magnitude in solubility. They concluded that nucleation occurs through blob formation pathway for both hydrophobic and hyrophillic guest molecules. They concluded that that formation of a blob of hydrophobic (small) guest is a rare event and limits the rate of nucleation more as compared to more soluble guest. It was observed that blobs of hydrophobic guests are rarer and longer-lived than those for hydrophilic guest.

ACS Paragon Plus Environment

15

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 89

Figure 6: Effect of the size of guest molecule on reaction pathway for hydrate nucleation shown by Jacobson et al.96 They established different blob pathways for nucleation of different guest molecules based on the guest size. In Figure 6, three pathways for small (S), medium (M) and extra-large (XL) guest molecule can be observed. The S solute fills all the cages, M primarily fills the larger cages and XL fills exclusively large cages. They observed that formation of empty 512 cages to be a parallel competing mechanism for nucleation. All the molecular simulations discussed until now were constant temperature simulations NPT or NVT. Ripmeester, Alavi, Englezos and coworkers102,103,104 have shown by means of constant energy (NVE) simulations that the heat released due to hydrate significantly affects the mechanism and the rate of dissociation. Liang et al.105 performed NVE simulations of hydrate formation to understand similar effect for hydrate nucleation. They observed higher order of crystallinity for the case of NVE simulation but the two-step mechanism of nucleation was confirmed in their simulation as well. Zhang et al.106 performed MD simulations to form hydrates with a methane non-bubble in liquid water at 250 K and 50 MPa to study the effect of different

ACS Paragon Plus Environment

16

Page 17 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

ensembles on nucleation kinetics of methane hydrates. They reported the sequence of nucleation rate in the ensembles as NPT>NVT>NVE. They reported that the order of crystallinity was reversed, i.e. the faster the hydrate forms the lower the crystallinity is. In terms of mimicking the real phenomenon in the early stage of nucleation, NVE is most suitable because the use of thermostat and/or a barostat in NPT and NVT ensembles is artificial since the exothermic heat of hydrate formation cannot be observed105. NVE simulations however come at the cost significantly higher simulation time, ~6.5 times slower than comparable NPT or NVT simulations107. The formation of amorphous crystals in the simulation studies had one common thing though. The molecular simulations were performed at high driving force by means of very high subcooling or pressurizing. Jacobson et al. in their subsequent work108 analysed two major questions: How can crystalline hydrates arise from amorphous precursors? In addition, are the amorphous precursors stable even at lower driving force for temperatures closer to equilibrium? They showed in their study that at high driving force, the amorphous nuclei are kinetically favourable over crystalline nuclei because of lower barriers of formation. They also reported that both amorphous and crystalline nuclei lead to the formation of crystalline clathrates. An implication of this finding is that macroscopic evidence for crystalline clathrates cannot be used to rule out amorphous pre-cursors. Similar results have been reported in parallel studies74. Interestingly, they also showed that cross nucleation of sII clathrates from sI nuclei is also possible in cases where both sI and sII structures are more stable than liquid water. They were not able to completely answer the question about the stability of amorphous crystals at lower driving force.

ACS Paragon Plus Environment

17

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 89

Nucleation studies based on diffraction and spectroscopic experiments have reported formation of mixture of different crystal structures68, 109-111. Walsh et al.112 suggested that even though the presence of amorphous structure observed might be subject to the time-scales and the driving force but both experiments and simulations suggest towards post-nucleation solid-solid rearrangement. To study the transformation from amorphous crystals to hydrate crystals by means of solid-solid rearrangement, Walsh et al.112 performed 20 MD simulation of methane clathrate. Walsh et al.112 reported the 7 types of cages (Figure 7) observed in all nucleation. They identified two types of cage-cage transformations: a) insertion type and b) rotation type transformation. They concluded that templates for exclusive growth of the thermodynamically favourable crystalline phase could form almost immediately upon nucleation, different than the amorphous crystals intermediates.

Figure 7: The seven dominant cages in incipient clathrate hydrate formation as observed from MD simulations from the study of Walsh et al.112. The most common cages shaded more prominently and the number represents the number of water molecule in each cage. Grey dot (51261) represents cage with physical infeasibility.

ACS Paragon Plus Environment

18

Page 19 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Multi-pathway nucleation It is evident from the discussion that there is enough evidence to support for the intermediate amorphous precursors. The stability of amorphous nuclei as compared to crystalline nuclei, the interconversion between them are a few things that could still require more proof to be conclusively established. More recently, studies have shown evidence of the two-step hydrate nucleation and direct crystalline nuclei occurring in parallel as competing pathways for hydrate nucleation107, 113. An extremely interesting observation made by Bi et al.113 in their recent work shows that although theories vastly different than CNT have been proposed for nucleation, the free-energy curves for the nucleation in their study is very close to that predicted by CNT. On the basis of multiple non-classical pathways and yet CNT like energy curves, Bi et al.113 concluded that hydrate nucleation occurs through energetically similar yet different pathways. Moreover, it was conjectured that hydrate nucleation is entropically driven kinetic process. Future work needs to be undertaken to corroborate these conjectures but lays the foundation for coming endeavors. Nucleation rate, critical radius and work of formation. Vysniauskas and Bishnoi et al.82,114 in their seminal work reviewed the kinetics of hydrate formation and developed a semi-empirical model to correlate experimental data on methane and ethane hydrate formation. In their work, the hydrate formation modelling did not include the hydrate nucleation. The driving force for the hydrate formation was considered to be the degree of supercooling that is the difference between the equilibrium temperature and the experimental temperature. Their study was the first attempt to describe quantitatively and model the formation kinetics of gas hydrates. Similar driving force had been adopted in various other studies as well before115. Englezos et al.116 subsequently modified the experimental procedure of Vysnauskas and Bishnoi to achieve homogenous hydrate nucleation conditions. With the new obtained

ACS Paragon Plus Environment

19

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 89

experimental results, Englezos et al.116 then proposed a mechanistic model (Englezos-Bishnoi model) for hydrate growth including nucleation with only one tunable parameter r that represents the rate of reaction. The overall driving force for the process was proposed as fugacity difference between the equilibrium and the experimental conditions.

∆f = f − f eq

(4)

Englezos et al.116 also estimated the critical radius by starting from the Gibbs free energy change in formation of a new phase

∆G = ∆Gs + ∆Gv = Apσ + ∆g vV p

(5)

where σ is the surface free energy and ∆gv is the difference in phase free energies. Assuming a spherical nucleus, at critical radius d ∆G/dr=0, we obtain critical radius as following:

rc =

(−∆gv ) =

−2σ ∆gv

(6)

fb, j n v ( P − P∞)  RT  2 )+ w w ∑1 θ j ln(  vh  f∞, j RT 

(7)

where vh and vw are the molar volumes of hydrate and water respectively; fb,j and f∞,j are the bulk phase and equilibrium fugacities, nw is the number of water molecules per gas molecule. Using the above expression, Englezos et al.116 calculated the critical radius of methane to be 30-170 Å. Skovberg et al.90 expressed the driving force for nucleation as the difference in the chemical potential of water in hydrate phase and water in water phase. They showed that induction times

ACS Paragon Plus Environment

20

Page 21 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

for methane, ethane and, mixed methane-ethane hydrate formation varies exponentially with the driving force. ∆ µ water = µ water liq − µ water hyd

(8)

Bishnoi et al.47 followed up the study of Englezos et al.116 and provided a unified description of the kinetics of hydrate nucleation, growth and decomposition. They modelled the induction time for Methane, Ethane and Carbon dioxide as a function of driving force expressed as fugacity difference between experimental conditions and equilibrium. Natarajan et al.91 adopted similar approach in modelling the induction times.

tind

 fV  = K  g − 1  f   eq 

−m

(9)

Where K is a constant and m is the index and both are species dependent. The work of Natarajan et al.91 was limited to modelling pure component nucleation. This was in light of the experimental observation from Flabella et al. that methane-ethane mixtures do not show induction time whereas pure methane hydrate have an non-zero induction time. Kaschiev et al.117 considered the crystallization of hydrates analogous to precipitation of salts. The driving force for the formation of new phase was considered the difference in chemical potential between the existing phase and the new phase42.

∆µ = µ gs + nw µ w − µ h

(10)

They obtained a general expression for the driving force as

∆µ = kT ln[γ ( P, T , C )υW C ] + µ gs* ( P, T ) + nW ( P, T ) µW ( P, T ) − µh ( P, T )

ACS Paragon Plus Environment

(11)

21

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 89

Further expressions were derived for the cases when the solution is in chemical equilibrium with the gas phase and depending on isothermal or isobaric regime of super-saturation in the solution. Upon establishing driving force based on chemical potential, Kashchiev and Firoozabadi estimated the size of hydrate nucleus and the work of nucleus formation as a function of super saturation in their subsequent work.118 Based on classical theory of nucleation, Work W (J) to form a hydrate cluster for 3 dimensional nucleation of one component condensed phase to hydrate nucleation at a given super saturation ∆µ (J) is given by42:

W (n) = −n∆µ + cν h2/3σ ef n 2/3

(12)

In the formula, ∆ µ is a known function of the pressure P or the temperature T taken from their previous study118. The number, n* of building units constituting the hydrate nucleus and nucleation work (the value of W at n=n*) can be estimated from the above equation. At n=n*, dW/dn=0 and results into the following general formula for critical size and the work done as: n* = 8c3vh2σ ef3 / 27 ∆µ 3 W * = 4c3vh2σ ef3 / 27 ∆µ 2

(13) (14)

Based on the estimations of work done, Kashchiev and Firoozabadi then established a general expression for the nucleation rate of one-component gas hydrate as: −4 c3 vH2 σ ef3 ∆g 2 kT 27 kT ( ∆ g)

J nuc = Ae e

ACS Paragon Plus Environment

(15)

22

Page 23 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Where c is the shape factor, σ ef is the effective superficial energy for the hydrate-solution interface. The kinetic factor, A (m-3 s-1) accounts for the attachment mechanism oy hydrate building units. The kinetic factor (A) contains the information about the kind of nucleation and the mechanism through which the hydrate building blocks attach to the hydrate nucleus. Typically, A is orders of magnitude lower for HEN than HON. With appropriate value of the driving force ∆µ calculated from their previous study117, the above equation provides the nucleation rate for different cases. Kashchiev and Firoozabadi117 performed an analysis on the effect of pressure and temperature on the nucleation rate of the methane hydrate under isothermal (273.2K) and isobaric (19.4 MPa) supersaturation respectively. A rapid increase in the nucleation rate with increase in the supersaturation was reported. They also reported that for smaller values of supersaturation, homogenous nucleation was order of magnitude smaller than heterogeneous nucleation. While at higher supersaturation, homogenous nucleation becomes dominant. Due to the rapid reduction of nucleation rate with reducing super saturation, a minimum supersaturation was identified below which the nucleation is virtually arrested.

Multicomponent nucleation: Anklam and Firoozabadi model for mixture predictions The model is based on the difference in chemical potentials as the driving force. As in their previous study117, hydrate formation was considered as an aqueous phase reaction for multicomponent system:

∑ n G +n i

i

w

H2O

n1G1n2G 2 ...nng G ng .n w H 2O

(16)

i

The driving force per unit cell is given by:

ACS Paragon Plus Environment

23

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 89

∆µ = ∑ ni ( P, T , x) µgi ( P,T ) + nw ( P,T )µsw ( P, T , x) − µH ( P, T ,z) i

In the equation, µgi is the chemical potential of gas species i in the solution, potential of water in the solution,

µ H is

(17)

µ sw is

the chemical

the chemical potential of hydrate. The composition and

the type of hydrate was shown as a function of z where

zi =

ni

∑n

. j

j

They first derived the expression for driving force for isothermal operation, with the assumptions of gas-saturated liquid, negligible compressibility of coefficient for liquid and hydrate phases and most importantly, composition of hydrate phase fixed as the equilibrium composition. The expression was effectively the same expression as given by Christiansen and Sloan87:  fi G (T , Pop , y)  ∆g = ∑ ni (T , Peq , y )kT ln  G  i =1  fi (T , Peq , y)  NG

(18)

Anklam and Firoozabadi119 challenged the assumption of constant hydrate phase composition even for super-saturation especially for the conditions where the conditions are far away from the gas-liquid-hydrate equilibrium. In their subsequent work, Anklam and Firoozabadi incorporated varying fractional cage occupancy to account for hydrate phase composition different than at equilibrium. They derived expressions for the critical nucleus size and the composition by taking the partial derivatives of the work of formation with respect to the number of molecules of individual components in the new phase. The generic expression derived for estimating the work done and the critical radius was: 2σ Vβ j R*

− ∆µ j = 0

ACS Paragon Plus Environment

(19)

24

Page 25 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

The final expression for an isothermal supersaturation process,

 1 − ∑k θ jk (T , Pop , y)   ∆g = nw ( Pop − Peq )(vw − vHw ) − kT ∑ v j ln   j 1 − ∑k θ jk (T , Peq , y)   

(20)

Where θ jk is the fraction of species j in cavity k and is given by Langmuir adsorption theory as

θ jk =

C jk f gj

(21)

1 + ∑ C jk f gj j

The expression derived for the critical radius in the study was as below:

rc =

2σ nwvHw ∆g

(22)

Anklam and Firoozabadi119 analysed the difference in driving force estimation by the two approaches, fixed hydrate phase composition vs. varying hydrate phase composition. The difference was found to be negligible for the case of pure methane gas and at lower pressure for the case of 75 % methane and 25 % ethane mixture (sI hydrate). The difference at higher pressure was found to be even larger for the case of 75% methane, 20% ethane, 5 % propane mixture (sII hydrate). Larger driving force in the case of mixtures leads to faster nucleation rates, lower work required for hydrate phase formation, smaller critical radius and reduced induction times. However, the results were not confirmed by experiments and required further inspection. In another study, Ma et al.120 studied hydrate formation kinetics of methane+ethylene+THF+H2O by experiments and modelling. They adopted two-step hydrate formation model by Chen-Guo121 to calculate the chemical potential based driving force for a mixture. First step is the formation of basic hydrate cage in which the basic cavity (larger cavity) is filled. The second step comprises

ACS Paragon Plus Environment

25

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 89

of adsorption of the gas molecules into smaller empty cages. Based on this model, they proposed the driving force as: ∆µ f0 = λ1 ln(1 − ∑θ j ) + λ2 ∑ xi* ln i RT fi j i

xi* =

fi

f i [1 − ∑θ j ]λ1 / λ2 0

(23)

(24)

j

It is evident from Equation 23 that filling of small cavities by small guest molecules can increase the driving force. In THF+CH4 or THF+H2 clathrates, THF fills the basic or the larger cavity first and subsequently the smaller molecule either H2 or CH4 occupies the smaller cages. This is highly relevant for the application of storage and transport of CH4 as solidified natural gas (SNG) and hydrogen. In these applications, THF acts as a promoter and results in formation of binary clathrates with THF occupying the larger cages. Ma et al.120 also evaluated the difference between fixed vs. variable fractional filling of the cavity. Their results corroborated those previously by Anklam and Firoozabadi. Lee et al.22 studied the formation of H2-THF binary hydrates and reported that they exhibit “tuning-effect” and the wt% of H2 can be controlled. Lee et al.22 observed a maximum storage capacity of 4.0 wt % for 0.15 mol% of THF concentration.. This has very strong implication for H2 storage application as the capacity and composition of H2 can be controlled as shown by Lee et al.22. This control was proposed due to partial occupancy of large cages by THF and hence increasing the sites available for H2 storage which can occupy the large cages quadruply or the smaller cages double. Strobel et al.122 showed that the amount of H2 stored exhibits Langmuir behaviour and saturates asymptotically to 1 wt % H2 with increase in pressure. 1 wt % H2 corresponds to one H2 molecule each occupying the smaller cages (512) and one THF molecule

ACS Paragon Plus Environment

26

Page 27 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

occupying the large cage (51264). The results of Strobel et al. contradicted the study of Lee et al. who were able to attain as high as 4 wt % H2 loading. Several other studies also have failed to attain the 4 wt% loading and have observed a maximum of 1 wt%123-125. Song et al.126 studied the growth of binary clathrates using molecular simulations. Their results showed that the stabilization energy of binary clathrate with larger cages occupied by XL molecule and smaller cages by SS molecule was larger than the sum of stabilization energies of a clathrate with only larger cage occupied by XL molecule and smaller cage occupied by SS. This is somewhat different than pure-component nucleation mechanism of Jacobson et al.96 as discussed in previous section. The additional stability was attributed to a decrease in lowfrequency modes of water network when all cages are occupied as stated by Rodger et al127. In their study, XL molecule was similar to the size of THF while SS molecule was corresponding to H2. Another major finding from their result was that the composition of clathrates at higher driving force is not necessarily the same as of the most stable phase. This should be kept into consideration while estimating the fractional cage occupancy. Their result showed that THF can act as both thermodynamic and kinetic promoter for H2-THF clathrates. However, there lies a trade-off between the promoting effect of THF and the storage capacity. For the application of H2 storage, the storage capacity is crucial. While THF facilitates formation of hydrates at milder conditions and increases the rate of formation, it reduces the H2 uptake capacity by occupying the larger cages. Song et al.126 showed that the fractional occupancy can be tuned by varying the growth temperature. Veluswamy et al.128 in their in-depth review of H2 storage in clathrate hydrates suggested that combined use of NMR and Raman spectroscopy simultaneously to analyse H2 inside the hydrate phase might help resolve the ambiguity of “tuning-effect” further.

ACS Paragon Plus Environment

27

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 89

Presence of binary components or higher number of hydrate formers can not only affect the cage occupancy but also the structure of hydrate formed. For example, addition of 1 mol% propane to methane can change from type I to type II hydrates129. Since natural gas at pipeline conditions mostly has small amounts of propane, hydrate plugs in pipelines usually have sII structure130. Abay and Svartaas131 studied multi-component gas hydrate nucleation for the effect of cooling rate and gas composition. They used two synthetic natural gas mixtures, for studying the secondary nucleation rate at equilibrium (secondary nucleation) and the stochasticity of nucleation. They showed that that systems with different gas composition respond differently to change in cooling rates. The implications of the results involve the uncertainty in systems employing varying amount of additives in the form of KHIs (discussed in next section) in the system. The studies discussed above treat the interface boundary as distinct whereas the actual interface for the hydrate nuclei constitutes a diffuse region. The sharp interface assumed by the CNT based studies is valid for the case where the interface is significantly smaller than the critical radius of the nuclei. Within the diffuse interface region, the properties vary continuously between the values of the individual phases. There have been many studies reporting the diffuse nature of the nuclei132,133,134,135. The diffuse nature of the interface can be incorporated within the modelling studies by means of phase field theories as shown by Kvamme et al.136,137, Granasy et al.138, Svandal et al.139. Next, we briefly discuss the phase field theory applied by kvamme et al.136.

Phase field theory based model: Kvamme Model Kvamme et al.136 developed a generic hydrate nucleation model based on phase field theory for describing the nucleation of CO2 hydrate in aqueous solutions. The phase field theory developed

ACS Paragon Plus Environment

28

Page 29 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

in the study was shown to be considerably more accurate than the sharp-interface droplet model of the classical nucleation theory. Based on the thermodynamic and interfacial properties, it was shown that the size of the critical fluctuations (nuclei) is comparable to the interface thickness, implying that the droplet model should be rather inaccurate. They proposed the local state of matter to be characterized by two fields (i) a structural order parameter, m, called as the phase field. This field describes the transition between the disordered liquid and ordered crystalline structures. (ii) A conserved field, χ, which maybe the coarse-grained density, ρ, or the solute concentration, c. Subsequently, they formulated Helmholtz free energy as a function of these fields. The resultant formula upon analysis was observed to have two minima. And the critical fluctuation (hydrate nucleus), being at is found at the extreme of the free energy function. The final expression for steady state provided by their model was:

J SS = J 0e −W

*

/ kT

(25)

There are two key parameters necessary for accurate prediction of the nucleation rate, superficial energy and the thickness of the interfacial region. It was proposed that the two can be obtained using molecular simulations for a given system.

Critical radius estimation using molecular simulations Radhakrishnan and Trout71 in their study of CO2 hydrate nucleation using Monte-Carlo (MC) simulations and non-Boltzmann sampling proposed an alternative approach for the estimation of critical radius. They implanted nuclei of two sizes 9.6 Å and 19.3 Å and analysed the trajectory of the two implants. Out of the two, the smaller crystal dissolved with time whereas the larger crystal grew and the entire system was observed to be converted to hydrates. Hence, it was concluded that the critical radius lies between the two values. They further calculated Landau-

ACS Paragon Plus Environment

29

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 89

free energy hypersurface as a function of a set of order parameters and plotted the first-order distribution functions for implants of size 9.6 Å, 14.5 Å and 19.3 Å. Upon observing the first order distribution, it was observed that the minima in the free energy shifts towards values of clathrate phase, as the system appears more clathrate like. It was concluded that the critical radius should lie between 9.6-14.5 Å. This was significantly smaller than ~32 Å predicted using CNT in an independent study by Larson and Garside140. Gibbs-Thomsom equation for spherical particles108:

Tm ( R ) = Tmbulk −

R *(T ) =

T bulk γν 2 K GT where KGT = m ∆H m R

−2γ 2γ ≈ ∆µ (T ) ρl ( ρ∆S m (Tmbulk − T ))

(26)

(27)

Where γ is the liquid-solid surface tension, ∆H is the bulk enthalpy of melting and v is the molar volume. The expression correlates melting temperature and the radius of the solid phase. At critical radius, the value of Tm equates with the bulk temperature T. Using the above approach, Jacobson et al.108 calculated the critical radii for both amorphous and crystalline nuclei in their study as discussed previously in the nucleation pathway section. Barnes et al.141 used the order parameter (OP) MCG-1 for estimating the critical nucleus size and nucleation rates. They estimated the critical nucleus size in terms of the OP MCG-1 by following the trajectory of each nucleus. Upon the observation, they proposed a MCG-1 value of 16 as the critical nuclei. They evaluated the quality of the proposed MCG value by performing a pB histogram test and showed that it is binomially shaped peaking near 0.5. Hence concluding that the proposed value is a good estimate.

ACS Paragon Plus Environment

30

Page 31 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Yuhara et al.142 analysed the results of Barnes et al.141 and calculated the nucleation rates and critical nucleus size of methane hydrates by implementing the work originally applied to the domain of vapour-liquid nucleation143. In attempt to reduce the computational load of direct numerical simulations performed by Barnes et al.141, Yuhara et al.142 used mean first-passage time (MFPT)144 and survival probability (SP)145 methods for estimating the nucleation rates and the critical nucleus size of methane hydrates. They obtained similar nucleation rates while higher critical nucleus but since MFPT and SP only require the simulation trajectories they are computationally less intensive.

Effect of Hydrate Inhibitors The discussion on inhibitors, its mechanism and properties discussed in the current work are restricted to the transport of the hydrocarbons in the oil and gas industry. Various methods have been adopted to curb the formation of hydrates such as raising the temperature/heating, lowering the pressure, removal of water. However, these options are not cost effective and the quest for lowering the cost has driven towards new techniques. To prevent the formation of hydrates in pipelines, thermodynamic hydrate inhibitors (THI) such as methanol, glycols etc. have been used traditionally. This type of inhibitor acts by shifting the hydrate formation phase boundary away from the operating conditions of the process in question. However, the use of THI was not feasible due to economic and environmental reasons. As per the current industry trends in oil and gas production industry, concentrations as high as 60 vol% methanol could be required to effective hydrate control. Koh146 in her review on natural gas hydrates estimated that for a small two-well satellite field producing 5.66x106 m3 gas per day at 379 K and 33 MPa, the amount of methanol required would be 20 tonnes per day ($5 million per year). This clearly demonstrate the economic infeasibility of using THIs.

ACS Paragon Plus Environment

31

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 89

As an alternative to high volume and subsequently high cost THI’s, low dosage hydrate inhibitors (LDHIs) have been developed to reduce hydrate formation. These include major two types: Kinetic hydrate inhibitors (KHIs) and Anti-agglomerates (AAs). Both classes of additives are added at low concentrations, typically around 0.1-1 wt%. LDHIs do not lower the three phase equilibrium but rather kinetically inhibit the formation of hydrates and crystal growth so that hydrate formation does not occur in the transport and operations146-147. It has been suggested that LDHIs act by adsorbing onto the hydrate surface and hinder hydrate growth.148,149,150 The efficiency of LDHIs, subsequently, depends upon the adsorption affinity. Anderson et al.149 identified the two molecular characteristics that lead to strongly binding inhibitors as charge distribution on the edge of the inhibitor and congruence of the size of the inhibitor with respect to the available space at the hydrate-surface binding site. KHIs are water-soluble polymers or copolymers such as PVP, poly (N-vinylpyrrolidonem Nvinylcaprolactam, dimethylaminoethyl acrylate etc. Efficient KHIs are commonly reported to have both amide groups and hydrophobic parts. Kelland et al.59 in their study reported that the major issues with application of KIs are: Performance, overall cost, environmental impact and compatibility. Conventionally, it has been proposed that the adsorption of KHIs is mainly due to the hydrogen bonds between the amide groups and water molecule on the hydrate surface.151,147 Carver et al.152 demonstrated the adsorption caused as a result of hydrogen bonding in hydrate/gas interface by means of Monte carlo simulations. Yagasaki et al.153 subsequently suggested that this adsorption effect would be affected for the case of hydrate/water interface since KHIs are water-soluble polymers. To further understand the adsorption of KHIs on the hydrate cages, Yagasaki et al.153 studied the role of hydrogen bonding and the effect of hydratewater interface by means of MD simulations. They concluded that contrary to prior suppositions,

ACS Paragon Plus Environment

32

Page 33 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

hydrogen bonding does not affect the affinity of KHIs adsorption on hydrate cages. It was reported that the adsorption affinities of KHIs is a resultant of entropic stabilization arising from the presence of cavities at the surface of the hydrates. Amide groups present on the KHIs do not enhance the adsorption affinity but enhance the solubility of the KHIs and hence are essential. Although certain KHIs have been shown to act through surface adsorption mechanism,154, there have been studies that show the mechanism might be molecule dependent rather than general.72 While molecules like PVCap19149, PDMAEMA155 have been shown to adsorb on the hydrate surface, Moon et al.72 showed that this is not the case for PVP molecule in their study based on MD simulation of methane hydrate. It was observed that PVP molecule remained at least 5-10 Å away from the surface of hydrate crystal. However, Moon et al.72 explained the results in terms of increased interfacial energy due to PVP. An increase in the hydrate/liquid interfacial energy results into increase in critical radius and decrease the stability of the particles below critical radius and thereby increasing the induction time of the system and delaying the nucleation. It was argued that the hydrate-water interface has been observed to be diffuse at a molecular level with the structure of water changing to a distance of 10-15 Å.55 Hence, this would place the PVP molecule in the middle of the hydrate-water interface and hence increasing the interfacial surface energy. Authors suggested that, this alternate pathway of PVP inhibition might also be responsible for synergistic effects upon combination of PVP and other KHIs such as PVPCap156. The results corroborated those by Anderson et al.149 which showed an absence of free energy driving force (∆Gads = 0.4 ±3.9 kcal/mol) for the case of PVP, as the authors quoted. However, it can be clearly seen that the result cited from the work of Anderson et al.149 was of Polyethylene Oxide (PEO) molecule mistook for PVP molecule. In the study, PVP was shown to have a binding energy of -20.6 kcal/mol which should be considered as physisorption. Even though not

ACS Paragon Plus Environment

33

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 89

corroborated by the results of Anderson et al.149 as claimed in their work, results of Moon et al.72 and the proposed hypothesis should be investigated further for veracity. Certain organisms that can survive at low temperatures do so by means of producing anti-freeze proteins (AFP) which adsorb to microscopic ice-crystals and prevent their further growth. These AFPs have also been used as KHIs for hydrate formation and are called as “green inhibitors”157. AFPs are generally produced by purifying from the organism that produces it such as bacteria, yeast, algae etc. While AFPs have been found as effective as other commercial KHIs such as PVP157-164, large-scale production, non-reproducible performance remain some of the major issues towards their implementation. Walker et al.165 provided an insightful review on the application of AFPs as gas hydrate inhibitors. For the purpose of evaluating KHIs and obtaining an understanding of synthesizing better KHIs, various experimental and simulation techniques have been employed in the literature. Most common of them being rocking cell150 and high-pressure cell or auto-claves166-167 or both in combination have been used168-169. Techniques such as nuclear magnetic resonance (NMR) microscopy160 and NMR combined techniques such as in-situ powder X-ray diffraction170, Highperformance differential scanning calorimetry53 (HP-DSC) and ultrasonic testing techniques171 have all been used to study the impact of KHIs on the hydrate nucleation and growth kinetics. KHIs have also been evaluated on a larger setup by means of high-pressure flow loops by Talaghat et al.172. High pressure automated lag time apparatus (HP-ALTA) is another technique that have been recently adopted for the purpose of studying hydrate formation. The main advantage of HP-ALTA is the ability to perform large number of experiments in a relatively short time as compared to the conventional techniques such as rocking cell or auto-claves. HPALTA provides a setup for measuring sufficient data in the hydrate formation system and allows

ACS Paragon Plus Environment

34

Page 35 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

representation in terms of probability distribution function173 and hence reduces the effects of intrinsic stochasticity of the process. Maeda et al.174-175 adopted the original version of HPALTA from the work of Heneghan,Heymat and Wilson176-178. HP-ALTA allows large number (>100) of nucleation and growth events to be recorded for a given sample under controlled pressure-temperature conditions. HP-ALTA cools a given sample at a specified rate until it detects formation of hydrates from a sudden reduction in the transmitted light through the sample. The formation temperature, Tf is is recorded and then the sample is heated to 10-15 K above equilibrium dissociation temperature and maintained for some time. After this, the sample cooling is started again and the system is repeated. Tf is recorded for each run and can be used for further analysis and representation by using a probability density function. The traditional HP-ALTA was modified for studying hydrate formation in which gas- liquid phase is separate and there is either a requirement for enhancing the mass transfer by mixing or focus the detection system at the gas-water interface. They modified the HP-ALTA for being able to operate in two configurations: Bulk-transmittance configuration that detects hydrate formation in bulk in conjunction with stirring and Interfacial-transmittance configuration that detects gas-hydrate surface. There are two major shortcomings in HP-ALTA, relative smaller size and quiescent setup. The size of HP-ALTA is smaller than other conventional techniques such as autoclave or a flow loop. As the sample size increases, the probability of nucleation occurring also increases. This could result into having a wider variation in onset of nucleation. Hence, HP-ALTA might lead to overestimation of the stochasticity. This is paradoxical since the intention of using HP-ALTA is to predict the stochasticity in the beginning. The use of linear cooling ramp causes the solution to be under-saturated and since the system is quiescent, the mass transfer to from the gas phase

ACS Paragon Plus Environment

35

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 89

represents a resistance179. We also need to keep in consideration that these setups might still be order of magnitude smaller than the industrial or actual implementation. It is peremptory to understand the effect of scale-up on stochasticity. Anti-agglomerants (AAs) on the other hand allow hydrates to form but the prevent them from agglomerating and accumulating into large masses. Some of the best performing AAs operate at higher subcooling than KHIs. There are two mechanisms by which AA operates depending on their physical structure. Quaternary AA’s180,181 are designed with hyrdophillic head group and hydrophobic tail group. The hydrophilic group at the quaternary centre binds to the hydrate particles whereas the long hydrophobic tail prevents the hydrate from continuing to grow on the surface. The hydrophobic tail also makes the surface more attractive to hydrocarbon phase and hence easily disperses the particles in the hydrocarbon phase. The seconds class of AAs form special water-in-oil emulsion.182 Hydrate formation is subsequently restricted to within the water phase and the end product is a slurry of oil with dispersed water-hydrate phase. In general, the upper water limit for the AAs to be effective is approximately 50% otherwise the slurry becomes too viscous for transporting.147 While KHIs are effective in delaying hydrate formation, a major downside of KHIs is the possibility of significantly enhanced hydrate growth (catastrophic growth) after hydrates nucleate183-184. Recently, Sharifi and Englezos185 introduced hydrate catastrophic index (HCI) to quantify the onset of the enhanced hydrate growth phenomenon based on laboratory data and the type of experimental conducted. Ohno et al.186 in their study of AFPs with hydrocarbon hydrates also showed that nucleation and growth inhibition to be independent processes. In principal, materials for controlling hydrate plugging may operate either by inhibiting nucleation or growth or both. The difference in the two modes is suggestive towards two separate adsorption

ACS Paragon Plus Environment

36

Page 37 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

characteristics of KHIs. One, adsorption with pre-nucleation clusters and foreign particles to “poison” the nucleation sites in case of heterogeneous nucleation. Two, adsorption on hydrate crystal surface. Hence, caution should be taken with the application of KHIs to limit pipeline plugging. A multi-scale approach is necessary for a complete evaluation of KHIs. Walker et al.165 outlined such a methodology involving the use of apparatuses ranging from XRD, Raman and NMR to Stirred Autoclaves. Apart from THIs and KHIs, applying localised high temperature2,

130

, external electric and

electro-magnetic fields have also been tried to disrupt the already-formed hydrates. Makogon130 reported that electric field intensity of 107 V/m would have to be employed to have any significant effect on equilibrium. English and MacElroy187 performed non-equilibrium MD simulation to evaluate how external electro-magnetic field disrupts the formation of methane hydrate nano-crystals They concluded that shifting dipolar alignment inside the crystals weakens the structure. The effect is intensity and frequency dependent and required at least 0.01-0.05 V/Å for tangible results.

Effect of Hydrate promoters Surfactants are one class of molecules that have been used in the literature as hydrate promoters because of their diverse agglomeration characteristics. Surfactants are organic compounds that lower the interfacial tension in a liquid-gas system by adsorbing at the interface. They are amphiphilic molecules i.e. have both hydrophobic (non-polar tail) and hydrophilic (polar head) groups. Three major types of surfactants are: Anionic, Cationic and Nonionic surfactants are classified based on the charge that the hydrophilic head possesses, negative, positive and neutral respectively. The presence of dual centres in surfactants results into diverse surface impacts in gas-liquid, solid-liquid and hydrate-liquid interfaces. Recently, Kumar et al.188 published a

ACS Paragon Plus Environment

37

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 89

review on the role of surfactants in promoting gas hydrate formation. In their work, Kumar et al.188 compiled a list of various molecules that have been used as surfactants for promoting gas hydrates. The focus of the current work has been restricted on the phenomenon of promoter mechanism and its effects pertaining to nucleation. We recommend other works for a further elaborate account on surfactants. Kalogerakis et al.189 experimentally investigated the effect of four surfactants on the kinetics of methane hydrates. In their study, they included anionic and non-ionic surfactants. The liquid side mass transfer coefficient was found to reduce for both cases, more so for anionic surfactant. Anionic surfactants enhanced the rate of hydrate formation more as compared to non-ionic surfactants. Karaaslan et al.190 also studied the performance of cationic, anionic and non-ionic surfactants in varying concentration for natural gas hydrate formation. It was reported that the hydrate kinetics was enhanced with anionic surfactants under all concentrations. Whereas for cationic surfactants, enhanced kinetics was observed at lower concentrations and contrasting effect (hydrate inhibition) at higher concentrations. Non-ionic surfactants had the least affect out of them all. Kumar et al.191 studied the effect of three different porous media and the effect of nonionic surfactant Tween-80 (T-80), cationic dodecyltrimethylammonium chloride (DTACl) and anionic Sodium Dodecyl Sulphate (SDS) on CO2 hydrates at constant pressure (3.55 MPa) and temperature (274 K). SDS was found to be most effective in reducing the induction time and enhancing the rate of hydrate formation. Before we discuss further about the mechanism of action for surfactants, let us first review the solubility and micelles formation in surfactants.

Krafft point and Critical micelle concentration Zhong and Rogers192 studied the effect of Sodium dodecyl sulfate (SDS) as a surfactant in methane/water and natural gas/water system. It was reported that SDS can increase the hydrate

ACS Paragon Plus Environment

38

Page 39 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

formation rates upto 700 times in quiescent systems. Zhong and Rogers calculated the critical micellar concentration (CMC) of SDS molecule to be 242 ppm at hydrate forming conditions. Critical micelle concentration is defined as the minimum concentration that is required at a given temperature to form micelle. They reported that above CMC, hydrate initiated subsurface around the micelle-solubilized hydrocarbon gas. The significantly enhanced hydrate growth was attributed to micelle solution with the micelles providing more nucleation sites apart from interfacial region which was observed as the initiation for hydrate growth quiescent systems. These developing hydrate particles migrated subsequently due to buoyancy to the water-gas interface and water-wet cell walls. Ramaswamy et al.193 in their study also concluded that SDS significantly reduces the induction time above CMC. Ganji et al. 18 reported that for cationic and non-ionic surfactants, much higher than 1000 ppm is required to have the desired promotor effect. However, the solubility of many surfactants can be significantly lower than it and hence the desired effect is not attainable.194

Figure 8: Effect of Temperature on surfactant concentration demonstrating the critical micelle concentration195.

ACS Paragon Plus Environment

39

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 89

To have a better understanding of the surfactant system, let us look at the different phases present for a surfactant/ water system (Figure 8). TK as shown in the figure is the Krafft point of a given ionic surfactant below which the solute does not form micelle rather precipitates as a hydrated solid196. The hypothesis of Zhong and Rogers was challenged by various studies later 194,195,197

. Watanabe et al.195 studied the effect of SDS on hydrate formation in a quiescent system

using HFC-32. Apart from other reasons in argument to the micelle hypothesis proposed by Zhong and Rogers, the main argument provided was that the Kraft point of SDS was significantly higher than the hydrate forming conditions adopted in their work (~279K). Hence, micelle formation was not possible and subsequently, the results should be considered cautiously. Di Profio et al.194 later corroborated the absence of micelle formation in the hydrate forming conditions when using SDS, Sodium oleate (SO) amongst other anionic surfactants. Watanabe et al.195 also challenged the CMC calculated in the study of Zhong and Rogers and Sun et al.198 since it was an order of magnitude lower than reported at atmospheric conditions in absence of methane/ethane. An important consideration generally forgotten is the impact that surfactants can have on the foam formation during hydrate dissociation. Higher concentrations of surfactants would expectedly increase the formation of foam during dissociation and hinder the process performance. Hence, it is important that surfactant has the desired effect on hydrate formation with lower concentration.

Effect of promoters on surface tension The hydrate formation for immiscible liquids is found to be highly dependent on the surface properties of the phase boundary. Surfactants reduce the interfacial surface tension of the gasliquid systems and increase the interfacial area of the system199,200,201. Surfactants have also been

ACS Paragon Plus Environment

40

Page 41 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

found to reduce the liquid side mass transfer coefficients in those studies. There is a divided opinion on cause of this reduction. Some studies report that surfactants induce a local modification in slip velocity at the interface which is responsible for reduced mass transfer.202,203 Some studies have reported that surfactants create a hydrodynamic change at the interface and an addition of a new resistance at the boundary layer film due to reduced local diffusion200, 204. And the third reasoning proposed is that by reducing the surface tension at the interface, promoters reduce the interfacial renewal rate and hence have lower mass transfer rate.205 The structure of the surfactant has been found to play a key role in its effect on surface tension. Surfactants with large hydrophobic and hydrophilic heads show lower interfacial tension than similar surfactants with smaller heads.206-207

Contrasting effects of SDS Across the various studies reported, Sodium dodecyl sulfate (SDS) has been shown to perform better than other anionic surfactants in terms of both rate of hydrate formation and the total hydrate formation in quiescent systems.208,209,210,191, 211,212 It was claimed by Zhang et al.94 that the effect of SDS is restricted to light hydrocarbons. For CO2 hydrate, SDS was found to have had no promoting effect.94 To understand the absence of promoting effect of SDS, Zhang et al.213 later studied the competitive adsorption between SDS and carbonate on THF hydrates. The reasoning provided was that the mechanism of SDS promotion is through surface adsorption of DS- on the hydrate crystals. Since at higher concentration of CO2, the Carbonate ions are sufficient to provide the adsorption promoted hydrate formation, the promotion effect of SDS is not observed for CO2 hydrates. SDS was also reported to have no promoting effect for methane-propane clathrates on the bubble surface214.

ACS Paragon Plus Environment

41

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 89

Contrary to the above studies, there also have been various studies across literature in which the effect of SDS has been clearly observed even for CO2 hydrates.215-217 Hence the attribution of the absence of SDS promotion effect to CO2 seems very likely to be incorrect. Ho et al.218 in their study on CO2 hydrates in an unstirred reactor showed that SDS in presence of Cyclopentane (CP) has no promotor effect in their system. Since CP has a density of ~750 Kg/m3, it forms a layer above water in quiescent system. The layer of CP above water inhibits the effect of SDS as a promoter and hence, indirectly indicates the importance of effect of surfactant on lowering the surface tension. It was also shown by Zhang et al.219 in their study that SDS does not act as a promoter in the case of heterogeneous nucleation for CH4+THF system. It was further conjectured that the enhancement effect of SDS might be eclipsed due to the higher nucleation rates present in heterogeneous systems. It may be argued that this acts as a proof that promoters function by adsorption on the surface of hydrate crystals. Lo et al.220 qualitatively investigated the adsorption of DS- on cyclopentane (CP) hydrates and proposed that the orientation of DS- on CP hydrates changes from “lie-down” to “stand-up” mode as the concentration of SDS increases. In the “stand-up” configuration, the head groups attach themselves to the hydrate surface while the tail orients towards the liquid phase forming a hydrophobic microdomains on the hydrate surface219,220. Hydrate formers such as methane can then solubilize in the hydrophobic microdomains, increasing their concentration at the interface and hence enhance the growth rate. To quantify the adsorption effects of surfactants, Lo et al.221 studied the adsorption of SDS and Dodecyl-trimethylammonium bromide (DTAB), a cationic surfactant on CP hydrates. CP hydrates were chosen because of their stability under ambient pressure at less than 280K suitable for the analytical techniques adopted.

ACS Paragon Plus Environment

42

Page 43 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Based on the abundant literature, it can be clearly seen that SDS has been proposed to be the best available promoter in certain cases while has been found to have no impact in some. The effect of SDS, apart from its concentration, can be seen to depend on the guest molecule (and hence hydrate structure) and reactor configuration as well. Dependence of reactor configuration and the guest molecule clearly indicates towards competing effects of surfactant on hydrate nucleation and growth. Surfactant act as a promotor of hydrate nucleation and growth by reducing the surface tension, increase the solubility of guest molecule, and adsorption on hydrate cages. It acts as an inhibitor forming a layer at the gas liquid interface and increase the mass transfer resistance. The pathway mechanism for surfactant can be seen as a resultant of the competing effects and hence is system dependent. Figure 9 shows the various impacts of surfactants and the system configuration impacting the overall effect of surfactant. Apart from the nucleation rates and the effects discussed above, surfactants have also been observed to affect the morphology of the gas hydrates formed195,189, 222,223,224,225. In current work we have not covered the morphology aspects of hydrate growth and readers are requested to refer to the cited articles for further readings.

ACS Paragon Plus Environment

43

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 89

Adsorption on Hydrate cages to promote hydrate growth

Reactor configuration and media

Effect of Surfactant on Hydrate Nucleation and Growth

Adsorption at gas-liquid interface to reduce surface tension

Morphological changes in hydrates

Figure 9: Various parameters to be considered while understanding the effect of surfactants on hydrate nucleation and growth.

Memory effect and hypothesis Memory effect is termed as the phenomenon in which hydrate nucleation is faster in hydratemelt or ice-melt systems than from fresh water systems. This concept originated from empirical observations such as that of Bishnoi and co-workers226,82 and Takeya et al.227 in which the induction times for methane hydrates was found to drastically reduce in the system with prior hydrate formation. This effect was also empirically observed in systems with prior ice-melt apart from hydrate-melt.227,228 Sefidroodi et al.229 examined the hydrate of cyclopentane from dissociated water and showed that a small amount of dissociated water was enough to induce the memory effect comparable to 100% dissociated water. The general consensus amongst researchers in the field is that if the system is heated to sufficiently high temperatures or long enough, the memory effect will be destroyed. A further complicating attribute of nucleation

ACS Paragon Plus Environment

44

Page 45 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

observed is the intrinsic stochasticity.227,145 Stochasticity nature of the memory effect reported from empirical studies implies that two systems can have different memory effects even while having the same thermal history. These attributes have made memory effect as a mysterious phenomenon but in the pursuit of explaining it, two major hypotheses have been put forward and tested. These are the residual cage structure hypothesis and higher local concentration in the solution. There have three primary methods to capture the nature of memory effect in the literature: (i) measuring of physical properties to detect the residual cage structures existing in the system, (ii) Spectroscopic methods for direct microscopic measurement of the nature of the system, and (iii) use of molecular simulation studies to obtain pathway insights. These will be discussed with the supporting theories in the next section.

Residual hydrate cage structures Upon these empirical observations, it was hypothesized that memory effect originates from residual clusters of water molecules after hydrate/ice dissociation. Surviving residual hydratelike structures provide the precursor crystals for hydrates for nucleation once its thermodynamically favourable to form hydrates again. Makogon et al.130 in their work provided data in support of this concept showing that hydrates do not completely decompose on dissociation but leave a partial structure that enables hydrate formation in the subsequent cycles. There have been various experimental techniques that have been adopted to capture the change in the liquid water post hydrate formation that would corroborate the residual cage hypothesis. Sloan et al.230 estimated effective kinematic viscosity of the liquid water phase by measuring the time taken by a stainless steel ball to travel one inch in the liquid. They reported an increase in the kinematic viscosity of liquid water phase post hydrate formation. Kato et al.231 performed simultaneous surface tension and kinematic viscosity measurements using a surface light

ACS Paragon Plus Environment

45

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 89

scattering method. It was observed that the two properties increase unusually elevated just after the hydrate dissociation. In contrast to the above stated supporting studies, Bylov and Rasmussen232 reported no detectable difference in the refractive index of the liquid water before hydrate formation and just after dissociation at a water surface in contact with methane or natural gas. Following these studies, Ohmura et al.233 measured the interfacial tension at the surface a sessile drop of a dense hydrate-forming liquid immersed in a liquid water pool. The memory effect was experimentally validated in their system without any change in the interfacial tension. Hence it was concluded that the memory effect did not arise from the molecular structuring in the liquid exerting a considerable effect on the surface tension. Apart from measuring physical properties, molecular simulations are the second means that have been used for studying the memory effect. In one such study, Baez et al.70 used MD simulations to study crystal growth and dissolution of natural gas hydrates. They concluded the persistence of dodecahedron cage (512 cage) preferentially over other cages as the reason for the memory effect. Subsequently, Yasuoka et al.234 showed that all dissociation is agnostic to the cage structure by mean of MD simulations of methane hydrate. However, it was conjectured in the study that memory effect is due to the residual structures of hydrate cages. The third means of understanding the memory effect has been spectroscopic techniques in an attempt for direct microscopic measurement. In one such study, Gao et al.235 used NMR to monitor the hydrate formation and decomposition in THF-heavy water mixtures. They studied the spin-lattice relaxation time, T1 to monitor the dynamics of THF during formation and decomposition. A marked increased in T1 was shown to be an evidence for the residual structures in the solution.

ACS Paragon Plus Environment

46

Page 47 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Higher local concentration of guest molecules An opposing hypothesis to residual hydrate cage structure that has been proposed states that memory effect is due to the presence of higher concentration of guest molecule and retarded diffusion upon melting. Although less in number as compared to the studies supporting residual cage hypothesis, there are various studies which support the hypothesis. It has been admitted that the higher local concentration of guest molecules would not explain the presence of memory effect in ice-melt systems. Rodger et al.236 studied memory effect phenomenon in methane hydrates by using long timescale MD simulations of a methane hydrate/methane gas interface. The thermodynamic conditions (15-20 oC above the equilibrium temperature) were chosen to ensure gentle melting of the hydrates and enabling the capture of memory effect, if any. They observed a significant increase of crystal structures in the melt, both hydrate-like and ice-like. The structures were predominantly ice-like. However, it was reported that there is no evidence of structured water molecules clusters and hence the authors concluded that memory effect is not related to longlived metastable hydrate precursors. As a corollary, the authors concluded that memory effect is attributed to higher methane concentration and retarded diffusion which they observed in their studies. Chapman et al.237 in their unpublished study based on NMR spectroscopy oh tetrahydrofuran also concluded the higher concentration and retarded diffusion is more related to the memory effect than the structural residuals. Buchanan et al.238 used neutron diffraction with H/D isotopes in combination with gas consumption measurements. They concluded that there is no significant structural difference in water before the hydrate formation and after hydrate dissociation and there is some dissolved residual gas post melting. Overall, the evidence

ACS Paragon Plus Environment

47

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 48 of 89

supporting the hypothesis remains thin and circumstantial and hence should be considered cautiously.

Surface Imprint Hypothesis Zeng et al.157,

239

studied the effect of Antifreeze proteins on the nucleation, growth and the

memory effect of THF hydrates. They reported that AFPs were able to eliminate the memory effect. They first estimated the homogenous nucleation temperature of THF hydrate to be 237 K using a DSC. Since the temperatures considered in the study were far above 237K, the authors concluded that they are operating in the regime of heterogeneous nucleation. They did not observe any memory effect for homogeneous nucleation in THF system. They argued against residual structure hypothesis since that should be an intrinsic property and should have been observed in their system. They provided an alternative hypothesis for memory effect in the form of surface imprint hypothesis. The hypothesis states that growth of hydrate might alter the surface properties of the impurity. In heterogeneous nucleation, these altered sites promote further nucleation resulting into the memory effect. Contact with the aqueous solution will allow the surface to return to its original state and limits the memory effect to a certain temperature and timescale. An important point to note is that their results on the effect of AFPs cannot distinguish whether they are hindering nucleation or growth. Effect of AFPs might be restricted to growth and might be mistaken for nucleation.

Factors affecting memory effect As discussed before, memory effect has been found to be mostly restricted to conditions close to equilibrium curve and the increasing the system temperature for a longer duration has been observed to remove the memory of the system240. Takeya et al.227 reported that nucleation rate was an order of magnitude of higher for CO2 hydrates at 1.9 K of superheating than at a higher

ACS Paragon Plus Environment

48

Page 49 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

temperature. The memory effect was completely lost at 298 K (25 K superheating). Ohmura et al. measured the induction time distributions of HCFC-141b hydrate from dissociated water at different superheating (0.5, 1, 1.5 K)145. They reported that the memory effect had started to wane at 1.5 K superheating. They also reported that there remains a degree of stochasticity to memory effect as different samples with same thermal history behave differently. Sefidroodi et al.229 investigated the strength and source of memory effect for cyclopentane hydrates. They studied the effect of melting temperatures and melting duration on the onset temperature of hydrate re-formation. The effect of memory effect is not limited to increasing the onset temperature of hydrate re-formation but also reduces the stochasticity of hydrate nucleation. The spread in the onset temperature reduces as the melting temperature increases. Their results showed that memory effect reduces with increasing melting temperature and duration. Sefidroodi et al.229 also showed that memory effect is transferable between solutions by showing its presence when they spiked fresh solution with only 3 ml and 1 ml of hydrate-melted water. They showed that the average onset temperature of the solutions was identical for the spiked systems and completely hydrate-melted solution. This provides an extremely strong evidence for residual cage hypothesis. Memory effect has been found to be affected by various additives as well. THF has been found to not affect the memory effect, whereas Antifreeze proteins have been found to remove the memory effect157. May et al.241 in their study quantitatively compared the effect of various KHIs and memory effect using HP-ALTA. Performance of six structurally different KHIs was evaluated along with the memory effect in the systems. Interestingly, it was empirically shown that memory effect is found to vary with the structure of KHIs given the same thermal history. Hence, it was conjectured that the result corroborates the residual structure hypothesis which

ACS Paragon Plus Environment

49

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 50 of 89

would imply that residual structures interact differently with the KHIs depending upon the structure of the KHI. While there have been significant studies demonstrating the presence of memory effect, there are also studies which have shown no memory effect at all232, 242-244. These studies either involved the superheating temperature significantly higher243 or involved guest molecules such as THF with high solubility so that supersaturation hypothesis cannot be refuted242. Recently, Sowa and Maeda studied the memory effect in model natural gas hydrate systems179 using HP-ALTA. They used 4 identical boat-shaped custom made glass sample cell (“boats”). Their results varied across each boat. They argued that the walls of the glass sample have effect on the memory effect. They also critically analysed the three hypothesis proposed for memory effect. They were not able to conclusive prove and showed that none of the hypothesis. In future, maybe a combination of hypothesis can be proposed to explain the elusive memory effect completely.

Molecular Simulation techniques and Order parameters Due to the complex nature of hydrate nucleation, different molecular simulations techniques have been developed to obtain qualitative and quantitative insights into the process. The molecular simulation domain has evolved rapidly from gaining insights on trajectories through direct numerical simulation under high driving force to obtaining statistical insights under more realistic conditions. These improvements has been achieved by means of advanced sampling methods such as forward flux sampling

245

, well-tempered metadynamics246, equilibrium path

sampling141, and generalized replica exchange method247. English et al.248 in their review provided insightful summary on the various molecular simulation techniques employed for hydrates. In current work, we have not covered in such details considering the broad scope of the study and the readers are recommended to find the molecular simulations in depth in the work by

ACS Paragon Plus Environment

50

Page 51 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

English et al.248 amongst other works. Next we discuss about the order parameters which form a key component of the molecular simulation studies but one that has been generally overlooked by literature review studies. Order parameters (OPs) are guidance parameters used in molecular simulation studies for estimating the progress of hydrate formation and other key properties. In general, order parameters are a quantitative measure of the degree of order in the system. The order parameters are system dependent, and are usually determined based on the symmetry of the ordered phase and how it differs from the disordered phase. Hydrate nucleation and growth is more complex and may involve multiple intermediary species. A broad spectrum of order parameters have been used in the study of gas hydrates such as detecting solvent structure of water to classify as being ice, liquid, or hydrate phase73,249 or to classify hydrate cages250 etc. OPs for hydrate nucleation have been based on guest species251,252,253, host molecules (water)245,252, or both141 (Table 3). Another classification of OPs is whether it considers a system structure metric or Global metric250, or it considers a local cluster structuring metric. First we discuss about three OPs based on water structuring F3, F4φ, F4t. These OPs were specifically designed to distinguish between different types of water structuring. They are derived from a two-level description of water molecule, a central molecule at location A is used to identify set of water molecules (Bi) in its solvation shell. F3 is a three body OP designed for probing angle deviation from those of tetrahedral hydrogen bonding network. The value of F3 is zero for tetrahedral structures and larger than zero otherwise. F4φ and F4t were originally designed as a test for clathrate-like structure by probing into torsional angles within a hydrogen bonding network and maybe defined in terms of quartet of oxygen atoms. Analogous parameter can be defined in terms of dimer of water molecule rather than through all water tetramers by

ACS Paragon Plus Environment

51

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 52 of 89

assuming that those O-H bonds not involved in the dimer and is computationally advantageous. F4φ is 0.7 for clathrate like structure; 0.3 for ice and close to 0 for liquid water55. Moon et al.72 adopted the three OPs for evaluating the effect of additives on nucleation and growth by means of direct simulations. They used instantaneous values of OPs to assign the water molecules to local phase and determine the extent of hydrate formation. They showed that the method is 95 % accurate in estimating the phase of molecule for those within the stable hydrate and ice phases. Radhakrishnan and Trout71 in their study showed that the Steinhardt bond-orientational OP254 was rendered useless for the case of crystalline clathrate and aqueous CO2 solution. They used tetrahedral OP as the leading OP for their analysis. As shown by Radhakrishnan and Trout, the combined use of the above OPs allows for the identification of regions of solid clathrates. However, these OP cannot distinguish between amorphous and crystalline or different crystalline clathrate regions. For the purpose of comprehensive pursuit of hydrate nucleation, Jacobson et at.253 used a portfolio of OPs. As a first step, they employed largest cluster of solvent-separated guests (LCSSG) OP for tracking identifying the densification of guest molecules. LCSSG tracks the formation of blob in the system but is insensitive towards the ordering of SSG and cannot distinguish between blob/amorphous crystals/crystalline clathrates. A cage OP was subsequently used to identify 5126n polyhedral water cages. The cage OP, however, does not distinguish between crystalline and amorphous nucleus255. Finally, for distinguishing between amorphous and crystalline and between different crystalline phases, they used vertex order parameter. In clathrate crystals, each water molecule is hydrogen bonded four other water molecules and constitutes the vertex shared by four polyhedral cages. The vertex order parameter differentiates between the different clathrate structures on the basis of types of cages in the vertices. The vertex

ACS Paragon Plus Environment

52

Page 53 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

OP developed in their study was based solely on water structuring. They provided an argument against guest-based OP that although not any more difficult to implement, it is not advisable due to two reasons. Firstly, the dual cages formed by guest molecule is inherently larger than the water cages, it is not suitable for identifying small or nascent crystallites. Secondly, not all cages in clathrates are necessarily filled with guests and hence reduces the symmetry of the guest ordering. Moreover, since it requires solely guest coordinates, it is computationally inexpensive to implement. Chakraborty et al.252 also utilized the computational advantage only solely guest coordinates in their work based on application of Voronoi tessellation analysis of clathrate hydrates252. Voronoi tessellation has been applied to other domains of research such as Leonard-Jones solids256, glassforming liquids257, phospholipid membranes258 and more, to study local structuring but is new for hydrate MD simulation studies. The simulation space is first tessellated (divided) based on guest molecule positions and results into Voronoi polyhedra (VP). The topological and metric properties of VP are further used to provide local ordering analysis and is used as OP. Chakrobarty et al.252 showed that the real advantage of using VP lies in computation speed (two order of magnitude faster than cage count code) in detecting ordered environment. However, VP lacks in differentiating between cage structures or detecting empty or doubly occupied cages. Ideally, VP could be combined with cage differentiating OP such as vertex OP developed by Jacobson et al.253 to provide an optimum combination of computational speed and desired cage detection. All the OPs discussed above have been single component OP that are based on either guest or host structuring. Motivated from the work of Walsh et al.75 showing adsorbed guest molecule on the water cages, Barnes et al.141 were the first to propose a two-component OP and quantified

ACS Paragon Plus Environment

53

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 54 of 89

nucleation and growth. Barnes et al.141 proposed an order parameter combining the use of both guest and solvent ordering, Mutually Coordinated Guest (MCG). MCG quantifies the appearance and connectivity of molecular clusters composed of guests separated by water clusters. Based on the set parameters, the unit MCG monomers are identified and the monomers adjacent to other monomers are classified as a cluster. The guests must be within distance Rgcut of each other, and the water molecules must be within distance Rwcut of each guest. A given minimum number of water molecules, Nw, must be within a specified angle φ of the guest-guest vector as shown in Figure 10. Barnes et al.141 compared the performance of MCG with other order parameters such as F4ϴ global metric, 512 water cage count250,259,112, potential energy trace, and the LCSSG OP proposed by Jacobson et al.253 . It was shown that MCG and global metric have higher resolution as compared to cage count to detect events such as nucleation, and critical nucleus size. MCG OP is structure neutral and detects blobs and crystalline nuclei.

Figure 10: Mutually Coordinated guest order parameter along with the tuneable model parameters253. Apart from their work on development of OPs for clathrate nucleation, Barnes et al.260 also conducted review for the advances in molecular simulation advancement. They identified key advanced sampling methods, developing optimal OPs, better understanding of gas-water interface as some of the major issues that remains to be tackled in the coming times. MD

ACS Paragon Plus Environment

54

Page 55 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

simulations of clathrate hydrates nucleation till date have been mostly been limited to conditions not identical with real experimental conditions in terms of driving force, time scales and homogenous nucleation. The use of advanced sampling methods along with a suitable order parameter and computations involving free energy landscape formulation and subsequent analysis form the barebones of the MD simulation studies. To move towards more realistic simulations, it requires an understanding and improvements on each component apart from the above listed shortcomings. As stated previously, under realistic conditions, homogeneous nucleation is rate is extremely low and nucleation is almost certainly to proceed through heterogeneous nucleation. Hence, the future efforts should be focussed towards understanding of heterogeneous nucleation, effect of various surfaces on hydrate nucleation and the mechanism of stabilization of hydrate nuclei in presence of surfaces. Some part of this is discussed in the next section on nucleation in porous media.

Nucleation in the presence of porous media As clathrate technology has moved towards the potential applications of gas separation and energy storage and transport, there is a pressing need for improving the mass transfer kinetics and heat transfer properties of the reactor configuration. Traditional stirred tank reactors face the problems of agglomeration of hydrate crystals at the surface and reduce the mass transfer coefficients19. Hydrate promoters are one of the agents that have been employed to improve the hydrate formation rate and increase process throughput. Along with promoters, various novel reactor configuration has been tried that might enable higher surface area for hydrate nucleation, better mass transfer properties and higher hydrate conversions. As a result, water saturated porous media such as in packed bed has gained preference over other reactor configuration261-264.

ACS Paragon Plus Environment

55

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 56 of 89

Further, applications such as CO2 sequestration involve replacing CH4 with CO2 in the reservoirs necessitate a better understanding hydrate formation in porous media. Hydrate formation in presence of solid surfaces or porous media has been reported to be much faster than bulk solution because they provide nucleation sites for hydrate formation265. An interesting phenomenon observed for the case of hydrate formation in porous media is the presence of multiple nucleation sites that have been experimentally observed. Apart from nucleation being heterogeneous, spatial heterogeneity has been observed in hydrate nucleation in the case of porous media264, 266. This is somewhat in line with nucleation in dispersed systems where each droplet can act as independent reactor and due to the stochasticity of nucleation, nucleation does not occur simultaneously in different locations. Linga et al.264 monitored the temperature inside a packed bed by means of multiple thermocouples. Figure 11 shown below the gas uptake measurement curves along with six thermocouple data for the hydrate formation. Hydrate nucleation and catastrophic growth process was detected by the heat release due to the process shown by the peaks in the recorded temperature in Figure 11A. The peak in the temperature was asynchronous for the thermocouples demonstrating the spatial heterogeneity in hydrate nucleation and nucleation in one part of the bed does not ensure or spread into other regions. Similar results have been shown in other studies as well. Linga et al.264 attributed higher overall hydrate formation in the experiment and hence higher throughput as the advantage of multiple nucleation.

ACS Paragon Plus Environment

56

Page 57 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Figure 11: Gas uptake measurement curve along with thermocouple data provided by Linga et al.264

Linga and Clarke267 also reported that the new porous media materials employed need to be tested within the scope of industrial application. Scaling up and further testing is fundamental towards the development of industrially viable hydrate based processes. Experimental research on hydrate formation in porous media encounters difficulty because of temporal and spatial limitations of monitoring techniques268 to provide detailed nucleation mechanism and a molecular level understanding. To obtain molecular level insights to hydrate formation in porous media, there have been various MD simulation studies to understand hydrate nucleation and growth on solid surfaces 269-270. Bai et al.269, 271 performed MD simulations of CO2 hydrate formation from hydroxylated silica surfaces towards the application of CO2 sequestration. They reported that nucleation tends to occur on the silica layer due to the stabilization effect of the silica layer. They concluded that solid surfaces with different

ACS Paragon Plus Environment

57

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 58 of 89

hydrophilic characters will affect the nucleation pathways in different way and hence further studies need to be conducted. They proposed a three step nucleation mechanism for hydrate growth in presence of solid surfaces269. First on a nanosecond timescale, an ice-like layer is formed closest to the substrate. On the microsecond scale, CO2 hydrate motif layer is observed that acts as nucleation seeds. While the intermediate step is formed by intermediate structure that transforms the ice-like layer and the final motif. Cygan et al.272 investigated the behaviour of CH4 hydrate in a clay interlayer. They reported that the CH4 hydrate structure in clays was different than in aqueous solution and bulk CH4 hydrate. Following their work, Yan et al.270 performed MD simulations of CH4 hydrate nucleation and growth process in a system containing bulk solution layer and clay layer. As per the study, CH4 first diffuses first from the bulk to the clay surface that promotes the formation of semi-cages and subsequently, hydrate grows along with stacking fault in the bulk solution region. CH4 molecules diffuse into the clay nanopore to form the ‘interlayer’ hydrate. The molecular diffusion of CH4 molecule into the nanopore is slow due to steric hindrance provided by hydrate crystals blocking the entrance of the pore. He et al.273 studied the effect of substrates on CH4 hydrate formation in Microsecond MD simulations. They investigated CH4 hydrate formation from gas/water two-phase system between hydrophilic silica surface and hydrophobic graphite surface. CH4 is preferentially adsorbed to the graphite surface while water bonds with silanol groups in silica. A cylindrical nano-bubble is formed that leads to higher aqueous CH4 concentration and hydrate formation. The studies clearly have demonstrated that hydrate nucleation and growth are more complex in porous media and require deeper understanding of the phenomenon. There also have been efforts to perform macroscopic numerical simulations consisting of multiphase flows in porous media and modelling of hydrate formation274-277. These works either use

ACS Paragon Plus Environment

58

Page 59 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

equilibrium reaction models or kinetic models based on intrinsic rate kinetics such as by Kim and Bishnoi47. These works focus on the macro-level parameters and can provide as useful tools for evaluating the process level performance. However, they do use empirical parameters for the simulations based on the physical setup such as characteristics of porous media or hydrate formation rates amongst others and require prior knowledge of the systems.

Future Perspectives Further development of experimental techniques There have been various spectroscopic studies adopted for studying the hydrate nucleation process and obtain direct evidence in support of theories or exploratory studies. However, as the application of hydrates have expanded towards the use of porous media and other reactor configurations, experimental techniques need to be adapted towards these applications. As discussed for the case of studying KHIs and in multi-component nucleation for studying the “tuning effect”, multi-scale approach needs to be adopted ranging from ranging from XRD, NMR, Raman Spectroscopy to Stirred Autoclaves. Further, in-situ experimental techniques for detection and study in porous media needs to be further developed.

Understanding the spatial-temporal stochastic nature of hydrate nucleation in porous media. The stochastic nature of hydrate nucleation, induction time and memory effect makes it challenging to design a process Application such as gas separation, desalination that requires continuous process face even bigger challenge for hydrate application. Experimental setups such as HP-ALTA have been adopted to statistically represent induction time in bulk solutions. The limitation of HP-ALTA to small size makes it challenging to use for quantifying temporal stochasticity desirably for industrial applications. It is very important to quantify the effect of

ACS Paragon Plus Environment

59

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 60 of 89

scale-up on the stochasticity. In porous media, this along with spatial heterogeneity adds the complexity further. Efforts need to be undertaken to understand the phenomenon of spatial and temporal stochastic nature of hydrate nucleation in porous media in order to develop hydrate based processes towards industrial applications.

Molecular Simulation improvement for simulating closer to experimental conditions. As discussed in the article, most of the molecular simulations due to computational limitations simulate systems with artificially high driving forces or lacking potential models. The nucleation pathways and results obtained are reliant on the conditions and can lead to artificial pathways non-existent in real world. Hence, the results should be considered carefully along with an effort to bring the MD simulations closer to the real conditions. More emphasis should be placed on the study of heterogeneous nucleation since under realistic conditions nucleation is most-likely to occur through this pathway. As stated by English et al.248, the massively parallel implementation of MD simulations on supercomputers and the increasing use of GPU278 for parallel implementation of MD are techniques that will play important roles.

Establishing the effects of hydrate promoters and porous media on the process. The effect of hydrate promoters is clearly multi-dimensional and the overall effect varies with the setups. The porous media also plays a key role in the hydrate based process and requires to be optimized for the application. In order to develop hydrate based technology, better understanding of these variables not simply in isolation but confounding effects needs to be clearly established.

Multi-scale simulation approach to hydrate based processes. There have been molecular simulation studies focused at obtaining molecular level understanding, experimental studies in search of validation at lab-scale setups and process or

ACS Paragon Plus Environment

60

Page 61 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

macro-level simulation studies to provide process performance indicator. Developing hydrate based technology towards industrial application is truly a multi-disciplinary effort and requires a common thread to bind the different domains of studies. There is a need for integration of molecular simulation studies of realistic hydrate processes to design experimental setups and seek validation.

Taking on the challenges towards industrial implementation Hydrate processes are extremely complex due to all the reasons listed in the current article and more. Apart from technical challenges such as optimizing the conditions, materials and additives, there is a requirement for overcoming the stochastic nature of the process by means of either statistical modelling or reducing it within bounds that do not alter the process performance significantly. Applications such as gas separation or desalination or energy storage and transport require a semi-batch or continuous process with precise control over the process and little scope for stochasticity. Scaling-up is a significant part of the challenge towards industrial implementation of any technology and generally a reduction in performance is observed upon scaling-up. Veluswamy et al.279 provided a multi-scale experimental validation of rapid methane hydrate formation for developing a cost-effective large scale energy storage system. They studied methane uptake in presence of THF-water system at three scales: small scale (2ml), medium scale (53 ml) and large scale (220 ml). Their result showed a reduction in induction time with increasing reactor volume. Also, they showed that the error bars spread for the induction time reduced which represents that the stochasticity of induction time reduced with increasing scales. They were able to attain this scale-up without significant loss in methane uptake. Figure 12 shows the results taken from their study. Their study shows that scaling-up in certain cases might turn out to be in favour of hydrate based technology which might be counter-intuitive at

ACS Paragon Plus Environment

61

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 62 of 89

first. The use of THF, however, has strong impacts on environment and hence has been questioned. So apart from scaling-up, the environmental aspects of the system need to be evaluated. Along with the use of novel reactor configuration and hydrate promoter, novel techniques such as electronucleation have also been proposed recently280-281. Carpenter and Bahadur showed that for pure THF hydrates can be reduced to an order of few minutes by applying 100 V for a few minutes pre-nucleation. They performed experiments at -5 oC and in a 5 ml system. In a follow up study, Shahriari281 et al. showed that using Al foam as anode reduced the induction time further to only tens of seconds. They proposed that the effect is due to two mechanisms, bubble formation and metal compound formation at anodes, both of which promote hydrate formation. Interestingly, the effect of increasing voltage on the induction times is observed to have similar effect of increasing reactor size shown by Veluswamy et al.279. The induction times reduce and also the stochasticity. There is a need for developing such novel techniques and demonstrate the effect of scale-up aspects.

ACS Paragon Plus Environment

62

Page 63 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Figure 12: Positive Effect of process scale-up on nucleation time (or induction time) in methane storage system experimentally validated by Veluswamy et al.279 Figure 13 shows the schematic representing the facets of hydrate nucleation discussed in the current article and clathrate process level performance parameters. The translation is not simple but effective implementation can only be brought about by keeping the broader picture into perspective.

ACS Paragon Plus Environment

63

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 64 of 89

Figure 13: Schematic representing facets of nucleation and its translation into overall clathrate process performance.

Conclusion From being traditionally considered a source of trouble in flow assurance and pipelines, clathrates are now seen as a potential for various applications such as gas separation, energy storage and transport, desalination and CO2 sequestration. Hydrate nucleation is a key step in the process and hence requires detailed understanding. While much efforts have been undertaken towards the nucleation pathways still remain up for debate and it seems likely that nucleation occurs through multiple competing pathways with varying degree of each depending upon the studies. Molecular simulation studies have provided significant insights to the nucleation pathways, nucleation rates and related effects and will continue to do so. While the effect of hydrate inhibitors on nucleation is more aligned, the effect of hydrate promoters is observed to vary significantly. Hydrate promoters are found to have promoting effects on hydrate nucleation in form of reducing the surface tension, increase the solubility of

ACS Paragon Plus Environment

64

Page 65 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

guest molecule, and adsorption on hydrate cages. They also act as an inhibitor forming a layer at the gas liquid interface and increase the mass transfer resistance. The overall effect is an outcome of the two contrasting effects. Memory effect has been experimentally validated but the presence of additives and the reactor media alters the memory effect. These effects have been studied in isolation but clearly cannot be reliable when other parameters change as well. The outlook for clathrate has shifted focus towards new applications and experimental proof of concepts have been conducted. Newer reactor configurations have been adopted towards better mass/heat transfer properties and faster kinetics. Significant progress is still required in terms of our understanding of clathrate processes, scaling up of processes and process development for these applications.

Acknowledgements The work was funded in part under the Energy Innovation Research Programme (EIRP, Award No. NRF2015EWTEIRP002-002), administrated by the Energy Market Authority (EMA), Singapore. We also acknowledge Lloyd’s Register Global Technology Center Singapore for the funding support.

ACS Paragon Plus Environment

65

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 66 of 89

Tables Table 1: Structural properties of most common clathrate hydrates.

Property

sI

sII

sH

Lattice type

Primitive cubic

Face centered cubic

Hexagonal a=1.21,

a=1.20

a=1.70

Unit cell parameters

c=1.01 3 (512) (S)

Unit Cell structure

2 (512) (S)

16 (512) (S)

2 (435663) (S)

6 (51262) (L)

8 (51262) (L)

1 (51268) (L) 0.391 (S)

Average cavity

0.395 (S)

0.391 (S)

0.406 (S)

radius

0.433 (L)

0.473 (L)

0.571 (L)

46

136

34

Number of water molecules per unit cell

ACS Paragon Plus Environment

66

Page 67 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Table 2 Compilation of nucleation pathways reported in the literature.

Studies Pathway

Comments

Nucleation type

supporting/Adopting Classical nucleation Theory

Mostly adopted from other established

3

Sloan et al. , Bishnoi domains

47

such

as

crystallization

and

et al. , Englezos et

Homogeneous applied analogously to the case of hydrate

al.116

formation. 84

Labile Cluster

Sloan and Fleyfel ,

Hypothesis

Christiansen and

(LCH)

Sloan86

It has been shown in literature that labile clusters

are

not

thermodynamically Homogeneous

favoured and the pathway is not correct for most cases of nucleation. Variation of LCH where labile clusters are

Nucleation at

Long et al.92 and

Interface

93

Kvamme et al.

formed at the interface and nucleation Homogeneous occurs at interface prior to diffusing into

Mechanism bulk phase. Thermal fluctuations are the cause of initial appearance of water cages and Local Structuring

Radhakrishnan and

results into structuring of guest molecules.

Trout71

Conceptual differs from LCH by stating

Homogeneous

mechanism that

guest

structuring

follows

water

structuring. Blob

Jacobson et al.95,

First hypothesis to propose amorphous

hypothesis: 2

Vatamanu et al.76,

precursors to clathrate crystals in a two-

Homogeneous

ACS Paragon Plus Environment

67

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Step formation

Lauricella et al.246,

step mechanism. A few other studies have

He et al.100

supported the occurrence of amorphous

Page 68 of 89

precursors prior to amorphous-crystalline transition. Reports of direct crystalline cages like Multi-pathway

Bi et al.113

nucleation

CNT along with two-step Blob hypothesis Homogeneous suggesting these mechanisms might be occurring in parallel. The nucleation occurs along the threephase line rather than vapour- liquid

2 Step formation Mechanism

271

Bai et al.

interface in a CO2/H2O/Silica system.

Heterogeneous

Silica surface acts as stabilizer for the CO2 amorphous crystals through hydrogen bonding and act as nucleation sites.

ACS Paragon Plus Environment

68

Page 69 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Table 3 Compilation of the order parameters used in the molecular simulation studies for establishing the nucleation pathways.

OP

Study

Scale

Comment

Based on host ordering F3i (3 body parameter) is used to identify Rodger et al.236, F3i, F4φi, F4ti

water molecules in a perfect tetrahedral

Moon56, Hawtin et al.

Global Ordering

69

network while F4φi, F4ti (4 body parameters) are used to distinguish between ice and hydrate crystals.

Face-saturated incomplete cage

Guo et al.

73

Global

Identifies all possible face-saturated

Ordering

complete and incomplete cages

analysis (FSICA) Cage Identifies 5 member and 6 member rings identification algorithm based

Matsumo et al.207,

Local

first and then dodecahedron cages (512

Jacobson et al. 216

ordering

cages) and tetrahedron cages (51262 cages)

on Ring referred as cups. perception Vertex OP based on Cage

Jacobson et al.253

Local/Global

Identification of sI, sII and amorphous

ordering

types lattice upon identification of Cage

identification

Based on guest ordering

ACS Paragon Plus Environment

69

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 70 of 89

Vatamanu et al.76 proposed this as critical Local Methane

72

Local

OP in early stages of methane nucleation

ordering

and explanation for the presence of

Vatmanu et al. density

memory effects. Calculates the Coordination number of

Guest coordintation

Jacobson et al. 253

Local guests molecules upon identifiction of ordering LCSSG

(GC) OP

Based on both guest and host ordering Bond Orientational OP,

Steinhardt et al.254,

Local/Global

Sensitive to the degree of global

Wolde et al.282

Ordering

orientational order in the system.

Q6 First identify solvent separated guests Largest cluster of solvent-separated

Jacobson et al.253

Local

(SSG), Then use clustering algorithm for

Ordering

connecting adjacent SSG and finding the

guests (LCSSG) largest cluster. Compatible with defective cages;

Mutually coordinated guest

Barnes et al.141

Local Combines the use of both guest and host Ordering

(MCG)

molecules Sensitive to both perfect crystals and

Voronoi

Chakraborty et

Local

amorphous crystals, computationally less

Tessellation

al.209

Ordering

intensive since uses only guest coordinates.

ACS Paragon Plus Environment

70

Page 71 of 89

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

71

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 72 of 89

References 1. Englezos, P., Clathrate hydrates. Ind. Eng. Chem. Res. 1993, 32 (7), 1251-1274. 2. Sloan Jr, E. D.; Koh, C., Clathrate hydrates of natural gases. CRC press: 2007. 3. Sloan, E. D., Fundamental principles and applications of natural gas hydrates. Nature 2003, 426 (6964), 353-363. 4. Jamaluddin, A.; Kalogerakis, N.; Bishnoi, P., Hydrate plugging problems in undersea natural gas pipelines under shutdown conditions. J. Pet. Sci. Eng. 1991, 5 (4), 323-335. 5. Hammerschmidt, E., Formation of gas hydrates in natural gas transmission lines. Ind. Eng. Chem. 1934, 26 (8), 851-855. 6. Sloan, E. D.; Koh, C. A.; Sum, A., Natural gas hydrates in flow assurance. Gulf Professional Publishing: 2010. 7. Sloan, E. D., A changing hydrate paradigm—from apprehension to avoidance to risk management. Fluid Phase Equilib. 2005, 228, 67-74. 8. Englezos, P.; Lee, J. D., Gas hydrates: A cleaner source of energy and opportunity for innovative technologies. Korean J. Chem. Eng. 2005, 22 (5), 671-681. 9. Linga, P.; Kumar, R.; Englezos, P., The clathrate hydrate process for post and pre-combustion capture of carbon dioxide. J. Hazard. Mater. 2007, 149 (3), 625-629. 10. Kang, S.-P.; Lee, H., Recovery of CO2 from flue gas using gas hydrate: thermodynamic verification through phase equilibrium measurements. Environ. Sci. Technol. 2000, 34 (20), 4397-4400. 11. Babu, P.; Linga, P.; Kumar, R.; Englezos, P., A review of the hydrate based gas separation (HBGS) process for carbon dioxide pre-combustion capture. Energy 2015, 85, 261-279. 12. Xu, C.-G.; Li, X.-S., Research progress of hydrate-based CO 2 separation and capture from gas mixtures. RSC Advances 2014, 4 (35), 18301-18316. 13. Li, X.-S.; Xu, C.-G.; Chen, Z.-Y.; Wu, H.-J., Tetra-n-butyl ammonium bromide semi-clathrate hydrate process for post-combustion capture of carbon dioxide in the presence of dodecyl trimethyl ammonium chloride. Energy 2010, 35 (9), 3902-3908. 14. Park, S.; Lee, S.; Lee, Y.; Lee, Y.; Seo, Y., Hydrate-based pre-combustion capture of carbon dioxide in the presence of a thermodynamic promoter and porous silica gels. INT. J. GREENHOUSE GAS CONTROL 2013, 14, 193-199. 15. Tajima, H.; Yamasaki, A.; Kiyono, F., Energy consumption estimation for greenhouse gas separation processes by clathrate hydrate formation. Energy 2004, 29 (11), 1713-1729. 16. Gudmundsson, J. S.; Parlaktuna, M.; Khokhar, A., Storage of natural gas as frozen hydrate. SPE PROD FACIL 1994, 9 (01), 69-73. 17. Khokhar, A.; Gudmundsson, J.; Sloan, E., Gas storage in structure H hydrates. Fluid Phase Equilib. 1998, 150, 383-392. 18. Ganji, H.; Manteghian, M.; Omidkhah, M.; Mofrad, H. R., Effect of different surfactants on methane hydrate formation rate, stability and storage capacity. Fuel 2007, 86 (3), 434-441. 19. Mori, Y. H., Recent advances in hydrate-based technologies for natural gas storage—a review. J. Chem. Ind. Eng.(China) 2003, 54 (1), 1-17. 20. Chong, Z. R.; Yang, S. H. B.; Babu, P.; Linga, P.; Li, X.-S., Review of natural gas hydrates as an energy resource: Prospects and challenges. APPL. ENERG. 2016, 162, 1633-1652. 21. Collett, T. S., Energy resource potential of natural gas hydrates. AAPG bulletin 2002, 86 (11), 1971-1992. 22. Lee, H.; Lee, J.-w.; Park, J.; Seo, Y.-T.; Zeng, H.; Moudrakovski, I. L.; Ratcliffe, C. I.; Ripmeester, J. A., Tuning clathrate hydrates for hydrogen storage. Nature 2005, 434 (7034), 743-746. 23. Gudmundsson, J. S., Method for production of gas hydrates for transportation and storage. Google Patents: 1996. 24. Kanda, H. In Economic study on natural gas transportation with natural gas hydrate (NGH) pellets, 23rd world gas conference, Amsterdam, 2006.

ACS Paragon Plus Environment

72

Page 73 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

25. Rehder, G.; Eckl, R.; Elfgen, M.; Falenty, A.; Hamann, R.; Kähler, N.; Kuhs, W. F.; Osterkamp, H.; Windmeier, C., Methane hydrate pellet transport using the self-preservation effect: a techno-economic analysis. Energies 2012, 5 (7), 2499-2523. 26. Kim, N.-J.; Lee, J. H.; Cho, Y. S.; Chun, W., Formation enhancement of methane hydrate for natural gas transport and storage. Energy 2010, 35 (6), 2717-2722. 27. Rossi, F.; Filipponi, M.; Castellani, B., Investigation on a novel reactor for gas hydrate production. APPL. ENERG. 2012, 99, 167-172. 28. Bi, Y.; Guo, T.; Zhu, T.; Zhang, L.; Chen, L., Influences of additives on the gas hydrate cool storage process in a new gas hydrate cool storage system. Energy Convers. Manage. 2006, 47 (18), 29742982. 29. Xie, Y.; Li, G.; Liu, D.; Liu, N.; Qi, Y.; Liang, D.; Guo, K.; Fan, S., Experimental study on a small scale of gas hydrate cold storage apparatus. APPL. ENERG. 2010, 87 (11), 3340-3346. 30. Kvamme, B.; Graue, A.; Buanes, T.; Kuznetsova, T.; Ersland, G., Storage of CO 2 in natural gas hydrate reservoirs and the effect of hydrate as an extra sealing in cold aquifers. INT. J. GREENHOUSE GAS CONTROL 2007, 1 (2), 236-246. 31. Hirai, S.; Tabe, Y.; Kuwano, K.; Ogawa, K.; Okazaki, K., MRI measurement of hydrate growth and an application to advanced CO2 sequestration technology. Ann. N.Y. Acad. Sci. 2000, 912 (1), 246253. 32. Koide, H.; Takahashi, M.; Shindo, Y.; Tazaki, Y.; Iijima, M.; Ito, K.; Kimura, N.; Omata, K., Hydrate formation in sediments in the sub-seabed disposal of CO2. Energy 1997, 22 (2-3), 279-283. 33. Kang, K. C.; Linga, P.; Park, K.-n.; Choi, S.-J.; Lee, J. D., Seawater desalination by gas hydrate process and removal characteristics of dissolved ions (Na+, K+, Mg 2+, Ca 2+, B 3+, Cl−, SO 4 2−). Desalination 2014, 353, 84-90. 34. Babu, P.; Kumar, R.; Linga, P., Unusual behavior of propane as a co-guest during hydrate formation in silica sand: Potential application to seawater desalination and carbon dioxide capture. Chem. Eng. Sci. 2014, 117, 342-351. 35. Javanmardi, J.; Moshfeghian, M., Energy consumption and economic evaluation of water desalination by hydrate phenomenon. Appl. Therm. Eng. 2003, 23 (7), 845-857. 36. Knox, W. G.; Hess, M.; Jones, G.; Smith, H., The hydrate process. Chem. Eng. Prog 1961, 57 (2), 66-71. 37. Strobel, T. A.; Hester, K. C.; Koh, C. A.; Sum, A. K.; Sloan, E. D., Properties of the clathrates of hydrogen and developments in their applicability for hydrogen storage. Chem. Phys. Lett. 2009, 478 (4), 97-109. 38. Ripmeester, J. A.; Tse, J. S.; Ratcliffe, C. I.; Powell, B. M., A new clathrate hydrate structure. Nature 1987, 325 (6100), 135-136 DOI: 10.1038/325135a0. 39. Barrer, R.; Ruzicka, D., Non-stoichiometric clathrate compounds of water. Part 4.—Kinetics of formation of clathrate phases. Trans. Faraday Soc. 1962, 58, 2262-2271. 40. Maini, B. B.; Bishnoi, P., Experimental investigation of hydrate formation behaviour of a natural gas bubble in a simulated deep sea environment. Chem. Eng. Sci. 1981, 36 (1), 183-189. 41. Ripmeester, J. A.; Alavi, S., Some current challenges in clathrate hydrate science: Nucleation, decomposition and the memory effect. Curr. Opin. Solid State Mater. Sci. 2016, 20 (6), 344-351. 42. Kashchiev, D., Nucleation. Butterworth-Heinemann: 2000. 43. Davies, S. R.; Hester, K. C.; Lachance, J. W.; Koh, C. A.; Sloan, E. D., Studies of hydrate nucleation with high pressure differential scanning calorimetry. Chem. Eng. Sci. 2009, 64 (2), 370-375. 44. Ohno, H.; Lipenkov, V. Y.; Hondoh, T., Air bubble to clathrate hydrate transformation in polar ice sheets: a reconsideration based on the new data from Dome Fuji ice core. Geophys. Res. Lett. 2004, 31 (21). 45. Knott, B. C.; Molinero, V.; Doherty, M. F.; Peters, B., Homogeneous nucleation of methane hydrates: Unrealistic under realistic conditions. J. Am. Chem. Soc. 2012, 134 (48), 19544-19547. 46. Jensen, L.; Thomsen, K.; von Solms, N., Propane hydrate nucleation: Experimental investigation and correlation. Chem. Eng. Sci. 2008, 63 (12), 3069-3080.

ACS Paragon Plus Environment

73

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 74 of 89

47. Bishnoi, P. R.; Natarajan, V., Formation and decomposition of gas hydrates. Fluid Phase Equilib. 1996, 117 (1-2), 168-177. 48. Wilson, P.; Heneghan, A.; Haymet, A., Ice nucleation in nature: supercooling point (SCP) measurements and the role of heterogeneous nucleation. Cryobiology 2003, 46 (1), 88-98. 49. Salamatin, A. N.; Hondoh, T.; Uchida, T.; Lipenkov, V. Y., Post-nucleation conversion of an air bubble to clathrate air–hydrate crystal in ice. J. Cryst. Growth. 1998, 193 (1), 197-218. 50. Kashchiev, D.; Firoozabadi, A., Induction time in crystallization of gas hydrates. J. Cryst. Growth. 2003, 250 (3), 499-515. 51. Dalmazzone, D.; Hamed, N.; Dalmazzone, C.; Rousseau, L., Application of high pressure DSC to the kinetics of formation of methane hydrate inwater-in-oil emulsion. J. Therm. Anal. Calorim 2006, 85 (2), 361-368. 52. Clausse, D.; Gomez, F.; Dalmazzone, C.; Noik, C., A method for the characterization of emulsions, thermogranulometry: Application to water-in-crude oil emulsion. J. Colloid Interface Sci. 2005, 287 (2), 694-703. 53. Lachance, J. W. Investigation of gas hydrates using differential scanning calorimetry with waterin-oil emulsions. Colorado School of Mines Golden, CO, 2008. 54. Long, J.; Sloan, E., Hydrates in the ocean and evidence for the location of hydrate formation. Int. J. Thermophys. 1996, 17 (1), 1-13. 55. Rodger, P.; Forester, T.; Smith, W., Simulations of the methane hydrate/methane gas interface near hydrate forming conditions conditions. Fluid Phase Equilib. 1996, 116 (1), 326-332. 56. Moon, C.; Taylor, P. C.; Rodger, P. M., Molecular dynamics study of gas hydrate formation. J. Am. Chem. Soc. 2003, 125 (16), 4706-4707. 57. Nerheim, A. R.; Svartaas, T. M.; Samuelsen, E. J. In Laser light scattering studies of gas hydrate formation kinetics, The Fourth International Offshore and Polar Engineering Conference, International Society of Offshore and Polar Engineers: 1994. 58. Bansal, N.; Drummond, C., Comment on “Kinetic study on the hexacelsian-celsian phase transformation”. J. Mater. Sci. Lett. 1994, 13 (6), 423-424. 59. Kelland, M. A.; Svartaas, T. M.; Øvsthus, J.; Namba, T., A new class of kinetic hydrate inhibitor. Ann. N.Y. Acad. Sci. 2000, 912 (1), 281-293. 60. Venkateswaran, N., Thermodynamics and nucleation kinetics of gas hydrates. Chemical and Petroleum Engineering, University of Calgary: 1993. 61. Guo, G.-J.; Rodger, P. M., Solubility of aqueous methane under metastable conditions: Implications for gas hydrate nucleation. J. Phys. Chem. B 2013, 117 (21), 6498-6504. 62. Jimenez-Angeles, F.; Firoozabadi, A., Enhanced hydrate nucleation near the limit of stability. J. Phys. Chem. C 2015, 119 (16), 8798-8804. 63. Koh, C. A.; Savidge, J. L.; Tang, C. C., Time-resolved in-situ experiments on the crystallization of natural gas hydrates. J. Phys. Chem. 1996, 100 (16), 6412-6414. 64. Subramanian, S.; Sloan, E., Molecular measurements of methane hydrate formation. Fluid Phase Equilib. 1999, 158, 813-820. 65. Koh, C. A.; Wisbey, R. P.; Wu, X.; Westacott, R. E.; Soper, A. K., Water ordering around methane during hydrate formation. J. Chem. Phys 2000, 113 (15), 6390-6397. 66. Moudrakovski, I. L.; Sanchez, A. A.; Ratcliffe, C. I.; Ripmeester, J. A., Nucleation and growth of hydrates on ice surfaces: new insights from 129Xe NMR experiments with hyperpolarized xenon. J. Phys. Chem. B 2001, 105 (49), 12338-12347. 67. Uchida, T.; Takeya, S.; Wilson, L.; Tulk, C.; Ripmeester, J.; Nagao, J.; Ebinuma, T.; Narita, H., Measurements of physical properties of gas hydrates and in situ observations of formation and decomposition processes via Raman spectroscopy and X-ray diffraction. Can. J. Phys. 2003, 81 (1-2), 351-357. 68. Murshed, M. M.; Kuhs, W. F., Kinetic studies of methane–ethane mixed gas hydrates by neutron diffraction and Raman spectroscopy. J. Phys. Chem. B 2009, 113 (15), 5172-5180.

ACS Paragon Plus Environment

74

Page 75 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

69. Ohno, H.; Strobel, T. A.; Dec, S. F.; Sloan, J., E Dendy; Koh, C. A., Raman studies of methaneв€’ ethane hydrate metastability. J. Phys. Chem. A 2009, 113 (9), 1711-1716. 70. Baez, L. A.; Clancy, P., Computer simulation of the crystal growth and dissolution of natural gas hydratesa. Ann. N.Y. Acad. Sci. 1994, 715 (1), 177-186. 71. Radhakrishnan, R.; Trout, B. L., A new approach for studying nucleation phenomena using molecular simulations: application to CO 2 hydrate clathrates. J. Chem. Phys 2002, 117 (4), 1786-1796. 72. Moon, C.; Hawtin, R.; Rodger, P. M., Nucleation and control of clathrate hydrates: insights from simulation. Faraday Discuss. 2007, 136, 367-382. 73. Hawtin, R. W.; Quigley, D.; Rodger, P. M., Gas hydrate nucleation and cage formation at a water/methane interface. PCCP 2008, 10 (32), 4853-4864. 74. Guo, G.-J.; Li, M.; Zhang, Y.-G.; Wu, C.-H., Why can water cages adsorb aqueous methane? A potential of mean force calculation on hydrate nucleation mechanisms. PCCP 2009, 11 (44), 1042710437. 75. Walsh, M. R.; Koh, C. A.; Sloan, E. D.; Sum, A. K.; Wu, D. T., Microsecond simulations of spontaneous methane hydrate nucleation and growth. Science 2009, 326 (5956), 1095-1098. 76. Vatamanu, J.; Kusalik, P. G., Observation of two-step nucleation in methane hydrates. PCCP 2010, 12 (45), 15065-15072. 77. Guo, G.-J.; Zhang, Y.-G.; Liu, C.-J.; Li, K.-H., Using the face-saturated incomplete cage analysis to quantify the cage compositions and cage linking structures of amorphous phase hydrates. PCCP 2011, 13 (25), 12048-12057. 78. Liang, S.; Kusalik, P. G., Exploring nucleation of H 2 S hydrates. Chem. Sci 2011, 2 (7), 12861292. 79. Volmer, M.; Weber, Α., Keimbildung in übersättigten Gebilden. Zeitschrift für physikalische Chemie 1926, 119 (1), 277-301. 80. Becker, R.; Döring, W., Kinetische behandlung der keimbildung in übersättigten dämpfen. Annalen der Physik 1935, 416 (8), 719-752. 81. Kelton, K., Crystal nucleation in liquids and glasses. Solid state physics 1991, 45, 75-177. 82. Vysniauskas, A.; Bishnoi, P., A kinetic study of methane hydrate formation. Chem. Eng. Sci. 1983, 38 (7), 1061-1072. 83. Thompson, S.; Gubbins, K.; Walton, J.; Chantry, R.; Rowlinson, J., A molecular dynamics study of liquid drops. J. Chem. Phys 1984, 81 (1), 530-542. 84. Sloan, E.; Fleyfel, F., A molecular mechanism for gas hydrate nucleation from ice. AlChE J. 1991, 37 (9), 1281-1292. 85. Muller-Bongartz, B.; Wildeman, T.; Sloan Jr, R. In A Hypothesis For Hydrate Nucleation Phonemena, The Second International Offshore and Polar Engineering Conference, International Society of Offshore and Polar Engineers: 1992. 86. Christiansen, R. L.; Sloan, E. D., Mechanisms and kinetics of hydrate formation. Ann. N.Y. Acad. Sci. 1994, 715 (1), 283-305. 87. Christiansen, R. L.; Sloan, E. D. J., A compact model for hydrate formation. Gas Processors Association, Tulsa, OK (United States): 1995; p Medium: X; Size: pp. 15-22. 88. Barrer, R. M., Edge, A.V.J.,, Gas hydrates containing argon krypton and xenon: kinetics, and energetics of formation and equilibria. Proceedings of the Royal Society of London 1967, A 300, 1-24. 89. Falabella, B. J., A study of natural gas hydrates. 1975. 90. Skovborg, P.; Ng, H.; Rasmussen, P.; Mohn, U., Measurement of induction times for the formation of methane and ethane gas hydrates. Chem. Eng. Sci. 1993, 48 (3), 445-453. 91. Natarajan, V.; Bishnoi, P.; Kalogerakis, N., Induction phenomena in gas hydrate nucleation. Chem. Eng. Sci. 1994, 49 (13), 2075-2087. 92. Long, J., Gas hydrate formation mechanism and kinetic inhibition. Colerado School of Mines: 1994. 93. Kvamme, B., A new theory for the kinetics of hydrate formation. . Proceedings of 2nd International Conference on Natural Gas Hydrates 1996, Monfort, J.P. (Ed.) p. 139.

ACS Paragon Plus Environment

75

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 76 of 89

94. Zhang, J.; Lo, C.; Somasundaran, P.; Lu, S.; Couzis, A.; Lee, J., Adsorption of sodium dodecyl sulfate at THF hydrate/liquid interface. J. Phys. Chem. C 2008, 112 (32), 12381-12385. 95. Jacobson, L. C.; Hujo, W.; Molinero, V., Amorphous precursors in the nucleation of clathrate hydrates. J. Am. Chem. Soc. 2010, 132 (33), 11806-11811. 96. Jacobson, L. C.; Hujo, W.; Molinero, V., Nucleation pathways of clathrate hydrates: effect of guest size and solubility. J. Phys. Chem. B 2010, 114 (43), 13796-13807. 97. Erdemir, D.; Lee, A. Y.; Myerson, A. S., Nucleation of crystals from solution: classical and twostep models. Acc. Chem. Res. 2009, 42 (5), 621-629. 98. Whitelam, S., Control of pathways and yields of protein crystallization through the interplay of nonspecific and specific attractions. Phys. Rev. Lett. 2010, 105 (8), 088102. 99. Sarupria, S.; Debenedetti, P. G., Homogeneous nucleation of methane hydrate in microsecond molecular dynamics simulations. J. Phys. Chem. Lett. 2012, 3 (20), 2942-2947. 100. He, Z.; Linga, P.; Jiang, J., What are the key factors governing the nucleation of CO2 hydrate? PCCP 2017, 19 (24), 15657-15661 DOI: 10.1039/C7CP01350G. 101. Suzuki, Y., Evidence of pressure-induced amorphization of tetrahydrofuran clathrate hydrate. Phys. Rev. B 2004, 70 (17), 172108. 102. Alavi, S.; Ripmeester, J., Nonequilibrium adiabatic molecular dynamics simulations of methane clathrate hydrate decomposition. J. Chem. Phys 2010, 132 (14), 144703. 103. Bagherzadeh, S. A.; Englezos, P.; Alavi, S.; Ripmeester, J. A., Molecular modeling of the dissociation of methane hydrate in contact with a silica surface. J. Phys. Chem. B 2012, 116 (10), 31883197. 104. Bagherzadeh, S. A.; Englezos, P.; Alavi, S.; Ripmeester, J. A., Molecular simulation of nonequilibrium methane hydrate decomposition process. J. Chem. Thermodyn. 2012, 44 (1), 13-19. 105. Liang, S.; Kusalik, P. G., Nucleation of gas hydrates within constant energy systems. J. Phys. Chem. B 2013, 117 (5), 1403-1410. 106. Zhang, Z.; Liu, C.-J.; Walsh, M. R.; Guo, G.-J., Effects of ensembles on methane hydrate nucleation kinetics. PCCP 2016, 18 (23), 15602-15608. 107. Zhang, Z.; Walsh, M. R.; Guo, G.-J., Microcanonical molecular simulations of methane hydrate nucleation and growth: evidence that direct nucleation to sI hydrate is among the multiple nucleation pathways. PCCP 2015, 17 (14), 8870-8876. 108. Jacobson, L. C.; Molinero, V., Can amorphous nuclei grow crystalline clathrates? The size and crystallinity of critical clathrate nuclei. J. Am. Chem. Soc. 2011, 133 (16), 6458-6463. 109. Staykova, D. K.; Kuhs, W. F.; Salamatin, A. N.; Hansen, T., Formation of porous gas hydrates from ice powders: Diffraction experiments and multistage model. J. Phys. Chem. B 2003, 107 (37), 10299-10311. 110. Klapproth, A.; Goreshnik, E.; Staykova, D.; Klein, H.; Kuhs, W. F., Structural studies of gas hydrates. Can. J. Phys. 2003, 81 (1-2), 503-518. 111. Schicks, J. M.; Ripmeester, J. A., The coexistence of two different methane hydrate phases under moderate pressure and temperature conditions: Kinetic versus thermodynamic products. Angew. Chem. Int. Ed. 2004, 43 (25), 3310-3313. 112. Walsh, M. R.; Rainey, J. D.; Lafond, P. G.; Park, D.-H.; Beckham, G. T.; Jones, M. D.; Lee, K.H.; Koh, C. A.; Sloan, E. D.; Wu, D. T., The cages, dynamics, and structuring of incipient methane clathrate hydrates. PCCP 2011, 13 (44), 19951-19959. 113. Bi, Y.; Porras, A.; Li, T., Free energy landscape and molecular pathways of gas hydrate nucleation. J. Chem. Phys 2016, 145 (21), 211909. 114. Vysniauskas, A.; Bishnoi, P., Kinetics of ethane hydrate formation. Chem. Eng. Sci. 1985, 40 (2), 299-303. 115. Makogon, Iпё U. F., Hydrates of natural gas. PennWell Books Tulsa, Oklahoma: 1981. 116. Englezos, P.; Kalogerakis, N.; Dholabhai, P.; Bishnoi, P., Kinetics of formation of methane and ethane gas hydrates. Chem. Eng. Sci. 1987, 42 (11), 2647-2658.

ACS Paragon Plus Environment

76

Page 77 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

117. Kashchiev, D.; Firoozabadi, A., Driving force for crystallization of gas hydrates. J. Cryst. Growth. 2002, 241 (1), 220-230. 118. Kashchiev, D.; Firoozabadi, A., Nucleation of gas hydrates. J. Cryst. Growth. 2002, 243 (3), 476489. 119. Anklam, M. R.; Firoozabadi, A., Driving force and composition for multicomponent gas hydrate nucleation from supersaturated aqueous solutions. J. Chem. Phys 2004, 121 (23), 11867-11875. 120. Ma, Q.-L.; Chen, G.-J.; Zhang, L.-W., Experimental and modeling study on gas hydrate formation kinetics of (methane+ ethylene+ tetrahydrofuran+ H2O). J. Chem. Eng. Data 2009, 54 (9), 2474-2478. 121. Chen, G.-J.; Guo, T.-M., A new approach to gas hydrate modelling. Chem. Eng. J. 1998, 71 (2), 145-151. 122. Strobel, T. A.; Taylor, C. J.; Hester, K. C.; Dec, S. F.; Koh, C. A.; Miller, K. T.; Sloan, E., Molecular hydrogen storage in binary THFв€’ H2 clathrate hydrates. J. Phys. Chem. B 2006, 110 (34), 17121-17125. 123. Anderson, R.; Chapoy, A.; Tohidi, B., Phase relations and binary clathrate hydrate formation in the system H2в€’ THFв€’ H2O. Langmuir 2007, 23 (6), 3440-3444. 124. Ogata, K.; Hashimoto, S.; Sugahara, T.; Moritoki, M.; Sato, H.; Ohgaki, K., Storage capacity of hydrogen in tetrahydrofuran hydrate. Chem. Eng. Sci. 2008, 63 (23), 5714-5718. 125. Struzhkin, V. V.; Militzer, B.; Mao, W. L.; Mao, H.-k.; Hemley, R. J., Hydrogen storage in molecular clathrates. Chem. Rev. 2007, 107 (10), 4133-4151. 126. Song, B.; Nguyen, A. H.; Molinero, V., Can guest occupancy in binary Clathrate hydrates be tuned through control of the growth temperature? J. Phys. Chem. C 2014, 118 (40), 23022-23031. 127. Rodger, P. M., Stability of gas hydrates. J. Phys. Chem. 1990, 94 (15), 6080-6089 DOI: 10.1021/j100378a082. 128. Veluswamy, H. P.; Kumar, R.; Linga, P., Hydrogen storage in clathrate hydrates: current state of the art and future directions. APPL. ENERG. 2014, 122, 112-132. 129. Susilo, R.; Alavi, S.; Ripmeester, J.; Englezos, P., Tuning methane content in gas hydrates via thermodynamic modeling and molecular dynamics simulation. Fluid Phase Equilib. 2008, 263 (1), 6-17. 130. Makogon, Iпё U. F., Hydrates of hydrocarbons. Pennwell Books: 1997. 131. Abay, H. K.; Svartaas, T. M., Multicomponent gas hydrate nucleation: the effect of the cooling rate and composition. Energy Fuels 2010, 25 (1), 42-51 DOI: 10.1021/ef1011879. 132. Mandell, M.; McTague, J.; Rahman, A., Crystal nucleation in a three‐dimensional Lennard‐Jones system. II. Nucleation kinetics for 256 and 500 particles. J. Chem. Phys 1977, 66 (7), 3070-3075. 133. Harrowell, P.; Oxtoby, D. W., A molecular theory of crystal nucleation from the melt. J. Chem. Phys 1984, 80 (4), 1639-1646. 134. Honeycutt, J. D.; Andersen, H. C., Small system size artifacts in the molecular dynamics simulation of homogeneous crystal nucleation in supercooled atomic liquids. J. Phys. Chem. 1986, 90 (8), 1585-1589. 135. Bagdassarian, C. K.; Oxtoby, D. W., Crystal nucleation and growth from the undercooled liquid: A nonclassical piecewise parabolic free‐energy model. J. Chem. Phys 1994, 100 (3), 2139-2148. 136. Kvamme, B.; Graue, A.; Aspenes, E.; Kuznetsova, T.; GrГЎnГЎsy, L.; TГіth, G.; Pusztai, T.; Tegze, G., Kinetics of solid hydrate formation by carbon dioxide: Phase field theory of hydrate nucleation and magnetic resonance imaging. PCCP 2004, 6 (9), 2327-2334. 137. Kvamme, B.; Qasim, M.; Baig, K.; KivelГ¤, P.-H.; Bauman, J., Hydrate phase transition kinetics from Phase Field Theory with implicit hydrodynamics and heat transport. INT. J. GREENHOUSE GAS CONTROL 2014, 29, 263-278. 138. GrГЎnГЎsy, L.; Pusztai, T.; TГіth, G.; Jurek, Z.; Conti, M.; Kvamme, B., Phase field theory of crystal nucleation in hard sphere liquid. J. Chem. Phys 2003, 119 (19), 10376-10382. 139. Svandal, A.; Kvamme, B.; GrГ nГ sy, LГ s.; Pusztai, TamГ s.; Buanes, T.; Hove, J., The phasefield theory applied to CO 2 and CH 4 hydrate. J. Cryst. Growth. 2006, 287 (2), 486-490.

ACS Paragon Plus Environment

77

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 78 of 89

140. Larson, M.; Garside, J., Solute clustering in supersaturated solutions. Chem. Eng. Sci. 1986, 41 (5), 1285-1289. 141. Barnes, B. C.; Beckham, G. T.; Wu, D. T.; Sum, A. K., Two-component order parameter for quantifying clathrate hydrate nucleation and growth. J. Chem. Phys 2014, 140 (16), 164506. 142. Yuhara, D.; Barnes, B. C.; Suh, D.; Knott, B. C.; Beckham, G. T.; Yasuoka, K.; Wu, D. T.; Sum, A. K., Nucleation rate analysis of methane hydrate from molecular dynamics simulations. Faraday Discuss. 2015, 179, 463-474. 143. Yasuoka, K.; Matsumoto, M., Molecular dynamics of homogeneous nucleation in the vapor phase. I. Lennard-Jones fluid. J. Chem. Phys 1998, 109 (19), 8451-8462. 144. Wedekind, J.; Strey, R.; Reguera, D., New method to analyze simulations of activated processes. J. Chem. Phys 2007, 126 (13), 134103. 145. Ohmura, R.; Ogawa, M.; Yasuoka, K.; Mori, Y. H., Statistical study of clathrate-hydrate nucleation in a water/hydrochlorofluorocarbon system: Search for the nature of the “memory effect”. J. Phys. Chem. B 2003, 107 (22), 5289-5293. 146. Koh, C. A., Towards a fundamental understanding of natural gas hydrates. Chem. Soc. Rev. 2002, 31 (3), 157-167. 147. Kelland, M. A., History of the development of low dosage hydrate inhibitors. Energy & Fuels 2006, 20 (3), 825-847. 148. Kelland, M. A.; Mønig, K.; Iversen, J. E.; Lekvam, K., Feasibility study for the use of kinetic hydrate inhibitors in deep-water drilling fluids. Energy & Fuels 2008, 22 (4), 2405-2410. 149. Anderson, B. J.; Tester, J. W.; Borghi, G. P.; Trout, B. L., Properties of inhibitors of methane hydrate formation via molecular dynamics simulations. J. Am. Chem. Soc. 2005, 127 (50), 17852-17862. 150. Lederhos, J.; Long, J.; Sum, A.; Christiansen, R.; Sloan, E., Effective kinetic inhibitors for natural gas hydrates. Chem. Eng. Sci. 1996, 51 (8), 1221-1229. 151. Perrin, A.; Musa, O. M.; Steed, J. W., The chemistry of low dosage clathrate hydrate inhibitors. Chem. Soc. Rev. 2013, 42 (5), 1996-2015. 152. Michael, G.; MarkГЎRodger, P., Inhibition of crystal growth in methane hydrate. J. Chem. Soc., Faraday Trans. 1995, 91 (19), 3449-3460. 153. Yagasaki, T.; Matsumoto, M.; Tanaka, H., Adsorption mechanism of inhibitor and guest molecules on the surface of gas hydrates. J. Am. Chem. Soc. 2015, 137 (37), 12079-12085. 154. Carver, T. J.; Drew, M. G.; Rodger, P., Configuration‐biased Monte Carlo simulations of poly (vinylpyrrolidone) at a gas hydrate crystal surface. Ann. N.Y. Acad. Sci. 2000, 912 (1), 658-668. 155. Hawtin, R. W.; Rodger, P. M., Polydispersity in oligomeric low dosage gas hydrate inhibitors. J. Mater. Chem. 2006, 16 (20), 1934-1942. 156. Larsen, R.; Knight, C. A.; Rider, K. T.; Sloan, E. D., Melt growth and inhibition of ethylene oxide clathrate hydrate. J. Cryst. Growth. 1999, 204 (3), 376-381. 157. Zeng, H.; Wilson, L. D.; Walker, V. K.; Ripmeester, J. A., Effect of antifreeze proteins on the nucleation, growth, and the memory effect during tetrahydrofuran clathrate hydrate formation. J. Am. Chem. Soc. 2006, 128 (9), 2844-2850. 158. Wathen, B.; Kuiper, M.; Walker, V.; Jia, Z., New simulation model of multicomponent crystal growth and inhibition. Chem. Eur. J. 2004, 10 (7), 1598-1605. 159. Daraboina, N.; Linga, P.; Ripmeester, J.; Walker, V. K.; Englezos, P., Natural gas hydrate formation and decomposition in the presence of kinetic inhibitors. 2. Stirred reactor experiments. Energy and Fuels 2011, 25 (10), 4384-4391 DOI: 10.1021/ef200813v. 160. Daraboina, N.; Moudrakovski, I. L.; Ripmeester, J. A.; Walker, V. K.; Englezos, P., Assessing the performance of commercial and biological gas hydrate inhibitors using nuclear magnetic resonance microscopy and a stirred autoclave. Fuel 2013, 105, 630-635 DOI: 10.1016/j.fuel.2012.10.007. 161. Daraboina, N.; Ripmeester, J.; Walker, V. K.; Englezos, P., Natural gas hydrate formation and decomposition in the presence of kinetic inhibitors. 3. Structural and compositional changes. Energy and Fuels 2011, 25 (10), 4398-4404 DOI: 10.1021/ef200814z.

ACS Paragon Plus Environment

78

Page 79 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

162. Daraboina, N.; Ripmeester, J.; Walker, V. K.; Englezos, P., Natural gas hydrate formation and decomposition in the presence of kinetic inhibitors. 1. High pressure calorimetry. Energy and Fuels 2011, 25 (10), 4392-4397 DOI: 10.1021/ef200812m. 163. Perfeldt, C. M.; Chua, P. C.; Daraboina, N.; Friis, D.; Kristiansen, E.; Ramløv, H.; Woodley, J. M.; Kelland, M. A.; Von Solms, N., Inhibition of gas hydrate nucleation and growth: Efficacy of an antifreeze protein from the longhorn beetle rhagium mordax. Energy and Fuels 2014, 28 (6), 3666-3672 DOI: 10.1021/ef500349w. 164. Daraboina, N.; Malmos Perfeldt, C.; von Solms, N., Testing antifreeze protein from the longhorn beetle Rhagium mordax as a kinetic gas hydrate inhibitor using a high-pressure micro differential scanning calorimeter. Can. J. Chem 2015, 93 (9), 1025-1030. 165. Walker, V. K.; Zeng, H.; Ohno, H.; Daraboina, N.; Sharifi, H.; Bagherzadeh, S. A.; Alavi, S.; Englezos, P., Antifreeze proteins as gas hydrate inhibitors. Can. J. Chem 2015, 93 (8), 839-849. 166. Svartaas, T.; Kelland, M.; Dybvik, L., Experiments related to the performance of gas hydrate kinetic inhibitors. Ann. N.Y. Acad. Sci. 2000, 912 (1), 744-752. 167. Chua, P. C.; Kelland, M. A., Tetra (iso-hexyl) ammonium Bromideо—ё The Most Powerful Quaternary Ammonium-Based Tetrahydrofuran Crystal Growth Inhibitor and Synergist with Polyvinylcaprolactam Kinetic Gas Hydrate Inhibitor. Energy Fuels 2012, 26 (2), 1160-1168. 168. Daraboina, N.; Linga, P., Experimental investigation of the effect of poly-N-vinyl pyrrolidone (PVP) on methane/propane clathrates using a new contact mode. Chem. Eng. Sci. 2013, 93, 387-394. 169. Daraboina, N.; Linga, P.; Ripmeester, J.; Walker, V. K.; Englezos, P., Natural gas hydrate formation and decomposition in the presence of kinetic inhibitors. 2. Stirred reactor experiments. Energy Fuels 2011, 25 (10), 4384-4391. 170. Ohno, H.; Moudrakovski, I.; Gordienko, R.; Ripmeester, J.; Walker, V. K., Structures of hydrocarbon hydrates during formation with and without inhibitors. J. Phys. Chem. A 2012, 116 (5), 1337-1343. 171. Yang, J.; Tohidi, B., Characterization of inhibition mechanisms of kinetic hydrate inhibitors using ultrasonic test technique. Chem. Eng. Sci. 2011, 66 (3), 278-283. 172. Talaghat, M.; Esmaeilzadeh, F.; Fathikaljahi, J., Experimental and theoretical investigation of simple gas hydrate formation with or without presence of kinetic inhibitors in a flow mini-loop apparatus. Fluid Phase Equilib. 2009, 279 (1), 28-40. 173. Wu, R.; Kozielski, K. A.; Hartley, P. G.; May, E. F.; Boxall, J.; Maeda, N., Probability distributions of gas hydrate formation. AlChE J. 2013, 59 (7), 2640-2646. 174. Maeda, N.; Wells, D.; Becker, N. C.; Hartley, P. G.; Wilson, P. W.; Haymet, A. D.; Kozielski, K. A., Development of a high pressure automated lag time apparatus for experimental studies and statistical analyses of nucleation and growth of gas hydrates. Rev. Sci. Instrum. 2011, 82 (6), 065109. 175. Maeda, N.; Wells, D.; Hartley, P. G.; Kozielski, K. A., Statistical analysis of supercooling in fuel gas hydrate systems. Energy Fuels 2012, 26 (3), 1820-1827. 176. Heneghan, A.; Haymet, A., Liquid-to-crystal nucleation: A new generation lag-time apparatus. J. Chem. Phys 2002, 117 (11), 5319-5327. 177. Heneghan, A.; Wilson, P.; Haymet, A., Heterogeneous nucleation of supercooled water, and the effect of an added catalyst. Proceedings of the National Academy of Sciences 2002, 99 (15), 9631-9634. 178. Heneghan, A.; Wilson, P.; Wang, G.; Haymet, A., Liquid-to-crystal nucleation: Automated lagtime apparatus to study supercooled liquids. J. Chem. Phys 2001, 115 (16), 7599-7608. 179. Sowa, B.; Maeda, N., Statistical study of the memory effect in model natural gas hydrate systems. J. Phys. Chem. A 2015, 119 (44), 10784-10790. 180. Ulfert Cornelis, K. R. R. Method for inhibiting the plugging of conduits by gas hydrates. 1996. 181. Urdahl, O. L., A. In IIR Conference on Waxes, Hydrates and Asphaltenes, Aberdeen, Scotland, Aberdeen, Scotland, 1998. 182. Bourgmayer, P. S., A; Behar., E. In 1989, 4th Multiphase Flow Conference, BHR Group. 183. Lee, J. D.; Englezos, P., Unusual kinetic inhibitor effects on gas hydrate formation. Chem. Eng. Sci. 2006, 61 (5), 1368-1376.

ACS Paragon Plus Environment

79

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 80 of 89

184. Cha, M.; Shin, K.; Seo, Y.; Shin, J.-Y.; Kang, S.-P., Catastrophic growth of gas hydrates in the presence of kinetic hydrate inhibitors. J. Phys. Chem. A 2013, 117 (51), 13988-13995. 185. Sharifi, H.; Englezos, P., Accelerated hydrate crystal growth in the presence of low dosage additives known as kinetic hydrate inhibitors. J. Chem. Eng. Data 2015, 60 (2), 336-342. 186. Ohno, H.; Susilo, R.; Gordienko, R.; Ripmeester, J.; Walker, V. K., Interaction of antifreeze proteins with hydrocarbon hydrates. Chem. Eur. J. 2010, 16 (34), 10409-10417. 187. English, N. J.; MacElroy, J., Theoretical studies of the kinetics of methane hydrate crystallization in external electromagnetic fields. J. Chem. Phys 2004, 120 (21), 10247-10256. 188. Kumar, A.; Bhattacharjee, G.; Kulkarni, B.; Kumar, R., Role of surfactants in promoting gas hydrate formation. Ind. Eng. Chem. Res. 2015, 54 (49), 12217-12232. 189. Kalogerakis, N.; Jamaluddin, A.; Dholabhai, P.; Bishnoi, P. In Effect of surfactants on hydrate formation kinetics, SPE international symposium on oilfield chemistry, Society of Petroleum Engineers: 1993. 190. Karaaslan, U. u.; Parlaktuna, M., Surfactants as hydrate promoters? Energy Fuels 2000, 14 (5), 1103-1107. 191. Kumar, A.; Sakpal, T.; Linga, P.; Kumar, R., Influence of contact medium and surfactants on carbon dioxide clathrate hydrate kinetics. Fuel 2013, 105, 664-671. 192. Zhong, Y.; Rogers, R., Surfactant effects on gas hydrate formation. Chem. Eng. Sci. 2000, 55 (19), 4175-4187. 193. Ramaswamy, D.; Sharma, M. M. In The Effect of Surfactants on the Kinetics of Hydrate Formation, SPE International Symposium on Oilfield Chemistry, Society of Petroleum Engineers: 2011. 194. Di Profio, P.; Arca, S.; Germani, R.; Savelli, G., Surfactant promoting effects on clathrate hydrate formation: Are micelles really involved? Chem. Eng. Sci. 2005, 60 (15), 4141-4145. 195. Watanabe, K.; Imai, S.; Mori, Y. H., Surfactant effects on hydrate formation in an unstirred gas/liquid system: An experimental study using HFC-32 and sodium dodecyl sulfate. Chem. Eng. Sci. 2005, 60 (17), 4846-4857. 196. Kunieda, H.; Shinoda, K., Krafft points, critical micelle concentrations, surface tension, and solubilizing power of aqueous solutions of fluorinated surfactants. J. Phys. Chem. 1976, 80 (22), 24682470. 197. Zhang, J.; Lee, S.; Lee, J. W., Does SDS micellize under methane hydrate-forming conditions below the normal Krafft point? J. Colloid Interface Sci. 2007, 315 (1), 313-318. 198. Sun, C.-Y.; Chen, G.-J.; Yang, L.-Y., Interfacial tension of methane+ water with surfactant near the hydrate formation conditions. J. Chem. Eng. Data 2004, 49 (4), 1023-1025. 199. GГіmez-DГ-az, D.; Navaza, J.; Sanjurjo, B., Mass-Transfer Enhancement or Reduction by Surfactant Presence at a Gasв€’ Liquid Interface. Ind. Eng. Chem. Res. 2009, 48 (5), 2671-2677. 200. Sardeing, R.; Painmanakul, P.; HГ©brard, G., Effect of surfactants on liquid-side mass transfer coefficients in gas–liquid systems: a first step to modeling. Chem. Eng. Sci. 2006, 61 (19), 6249-6260. 201. Hebrard, G.; Zeng, J.; Loubiere, K., Effect of surfactants on liquid side mass transfer coefficients: a new insight. Chem. Eng. J. 2009, 148 (1), 132-138. 202. Cuenot, B.; Magnaudet, J.; Spennato, B., The effects of slightly soluble surfactants on the flow around a spherical bubble. J. Fluid Mech. 1997, 339, 25-53. 203. Alves, S.; Orvalho, S.; Vasconcelos, J., Effect of bubble contamination on rise velocity and mass transfer. Chem. Eng. Sci. 2005, 60 (1), 1-9. 204. Painmanakul, P.; LoubiГЁre, K.; HГ©brard, G.; Mietton-Peuchot, M.; Roustan, M., Effect of surfactants on liquid-side mass transfer coefficients. Chem. Eng. Sci. 2005, 60 (22), 6480-6491. 205. Rosso, D.; Huo, D. L.; Stenstrom, M. K., Effects of interfacial surfactant contamination on bubble gas transfer. Chem. Eng. Sci. 2006, 61 (16), 5500-5514. 206. Barakat, Y.; Fortney, L. N.; Schechter, R. S.; Wade, W. H.; Yiv, S. H.; Graciaa, A., Criteria for structuring surfactants to maximize solubilization of oil and water: II. Alkyl benzene sodium sulfonates. J. Colloid Interface Sci. 1983, 92 (2), 561-574.

ACS Paragon Plus Environment

80

Page 81 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

207. Kunieda, H.; Shinoda, K. o. z. o., Correlation between critical solution phenomena and ultralow interfacial tensions in a surfactant/water/oil system. Bull. Chem. Soc. Jpn. 1982, 55 (6), 1777-1781. 208. Lin, W.; Chen, G.-J.; Sun, C.-Y.; Guo, X.-Q.; Wu, Z.-K.; Liang, M.-Y.; Chen, L.-T.; Yang, L.-Y., Effect of surfactant on the formation and dissociation kinetic behavior of methane hydrate. Chem. Eng. Sci. 2004, 59 (21), 4449-4455. 209. Link, D. D.; Ladner, E. P.; Elsen, H. A.; Taylor, C. E., Formation and dissociation studies for optimizing the uptake of methane by methane hydrates. Fluid Phase Equilib. 2003, 211 (1), 1-10. 210. Xie, Y.; Guo, K.; Liang, D.; Fan, S.; Gu, J., Steady gas hydrate growth along vertical heat transfer tube without stirring. Chem. Eng. Sci. 2005, 60 (3), 777-786. 211. Li, J.; Liang, D.; Guo, K.; Wang, R., The influence of additives and metal rods on the nucleation and growth of gas hydrates. J. Colloid Interface Sci. 2005, 283 (1), 223-230. 212. Zhang, J.; Lee, S.; Lee, J. W., Kinetics of methane hydrate formation from SDS solution. Ind. Eng. Chem. Res. 2007, 46 (19), 6353-6359. 213. Zhang, J.; Lo, C.; Somasundaran, P.; Lee, J., Competitive adsorption between SDS and carbonate on tetrahydrofuran hydrates. J. Colloid Interface Sci. 2010, 341 (2), 286-288. 214. Lee, S. Y.; Kim, H. C.; Lee, J. D., Morphology study of methane–propane clathrate hydrates on the bubble surface in the presence of SDS or PVCap. J. Cryst. Growth. 2014, 402, 249-259. 215. Kang, S.-P.; Lee, J.-W., Kinetic behaviors of CO 2 hydrates in porous media and effect of kinetic promoter on the formation kinetics. Chem. Eng. Sci. 2010, 65 (5), 1840-1845. 216. Torré, J.-P.; Ricaurte, M.; Dicharry, C.; Broseta, D., CO 2 enclathration in the presence of water-soluble hydrate promoters: hydrate phase equilibria and kinetic studies in quiescent conditions. Chem. Eng. Sci. 2012, 82, 1-13. 217. Dicharry, C.; Duchateau, C.; Asbaï, H.; Broseta, D.; Torré, J.-P., Carbon dioxide gas hydrate crystallization in porous silica gel particles partially saturated with a surfactant solution. Chem. Eng. Sci. 2013, 98, 88-97. 218. Ho, L. C.; Babu, P.; Kumar, R.; Linga, P., HBGS (hydrate based gas separation) process for carbon dioxide capture employing an unstirred reactor with cyclopentane. Energy 2013, 63, 252-259. 219. Zhang, J.; Lee, J. W., Effect of sodium dodecyl sulfate on the supercooling point of ice and clathrate hydrates. Energy Fuels 2009, 23 (6), 3045-3047. 220. Lo, C.; Zhang, J.; Somasundaran, P.; Lu, S.; Couzis, A.; Lee, J., Adsorption of surfactants on two different hydrates. Langmuir 2008, 24 (22), 12723-12726. 221. Lo, C.; Zhang, J.; Couzis, A.; Somasundaran, P.; Lee, J., Adsorption of cationic and anionic surfactants on cyclopentane hydrates. J. Phys. Chem. C 2010, 114 (31), 13385-13389. 222. Gayet, P.; Dicharry, C.; Marion, G.; Graciaa, A.; Lachaise, J.; Nesterov, A., Experimental determination of methane hydrate dissociation curve up to 55MPa by using a small amount of surfactant as hydrate promoter. Chem. Eng. Sci. 2005, 60 (21), 5751-5758. 223. Yoslim, J.; Linga, P.; Englezos, P., Enhanced growth of methane–propane clathrate hydrate crystals with sodium dodecyl sulfate, sodium tetradecyl sulfate, and sodium hexadecyl sulfate surfactants. J. Cryst. Growth. 2010, 313 (1), 68-80. 224. Veluswamy, H. P.; Chen, J. Y.; Linga, P., Surfactant effect on the kinetics of mixed hydrogen/propane hydrate formation for hydrogen storage as clathrates. Chem. Eng. Sci. 2015, 126, 488499. 225. Lim, Y.-A.; Babu, P.; Kumar, R.; Linga, P., Morphology of carbon dioxide–hydrogen–cyclopentane hydrates with or without sodium dodecyl sulfate. Cryst. Growth Des. 2013, 13 (5), 2047-2059. 226. Parent, J.; Bishnoi, P., Investigations into the nucleation behaviour of methane gas hydrates. Chem. Eng. Commun. 1996, 144 (1), 51-64. 227. Takeya, S.; Hori, A.; Hondoh, T.; Uchida, T., Freezing-memory effect of water on nucleation of CO2 hydrate crystals. J. Phys. Chem. B 2000, 104 (17), 4164-4168.

ACS Paragon Plus Environment

81

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 82 of 89

228. Hwang, M.; Wright, D.; Kapur, A.; Holder, G., An experimental study of crystallization and crystal growth of methane hydrates from melting ice. J. Incl. Phenom. Macrocycl. Chem. 1990, 8 (1), 103-116. 229. Sefidroodi, H.; Abrahamsen, E.; Kelland, M. A., Investigation into the strength and source of the memory effect for cyclopentane hydrate. Chem. Eng. Sci. 2013, 87, 133-140. 230. Sloan, E. D.; Subramanian, S.; Matthews, P.; Lederhos, J.; Khokhar, A., Quantifying hydrate formation and kinetic inhibition. Ind. Eng. Chem. Res. 1998, 37 (8), 3124-3132. 231. KATO, T., Y. YOKOTA & Y. NAGASAKA. In Observation of the dynamic variations of interfacial tension and kinematic viscosity of Xe–water by the surface lightscattering method., Proc. 5th Asian Thermophysical Properties Conf., M.S. Kim & S.T. Ro, E., Ed. 1998. 232. Bylov, M.; Rasmussen, P., Experimental determination of refractive index of gas hydrates. Chem. Eng. Sci. 1997, 52 (19), 3295-3301. 233. Ohmura, R.; Shigetomi, T.; Mori, Y., Mechanical Properties of Water/Hydrate‐Former Phase Boundaries and Phase‐Separating Hydrate Films. Ann. N.Y. Acad. Sci. 2000, 912 (1), 958-966. 234. Yasuoka, K.; Murakoshi, S., Molecular dynamics simulation of dissociation process for methane hydrate. Ann. N.Y. Acad. Sci. 2000, 912 (1), 678-684. 235. S. Gao, W. C., and W. House Proceedings of the Fifth International Conference on Gas Hydrates: Volume 1 Kinetics and Transport Phenomena, Trondheim, Norway, 13–15 June 2005 Trondheim, Norway, 2005. 236. Rodger, P., Methane hydrate: melting and memory. Ann. N.Y. Acad. Sci. 2000, 912 (1), 474-482. 237. W.G. Chapman, S. G., M. Yarrison, K. Song and W. House In Proceedings of the Tenth International Conference on PPEPPD, , Snowbird, UT,, Engineering Conferences International, NY: Snowbird, UT,, 2004. 238. Buchanan, P.; Soper, A. K.; Thompson, H.; Westacott, R. E.; Creek, J. L.; Hobson, G.; Koh, C. A., Search for memory effects in methane hydrate: structure of water before hydrate formation and after hydrate decomposition. J. Chem. Phys 2005, 123 (16), 164507. 239. Zeng, H.; Moudrakovski, I. L.; Ripmeester, J. A.; Walker, V. K., Effect of antifreeze protein on nucleation, growth and memory of gas hydrates. AlChE J. 2006, 52 (9), 3304-3309. 240. Duchateau, C.; Peytavy, J.-L.; GlГ©nat, P.; Pou, T.-E.; Hidalgo, M.; Dicharry, C., Laboratory evaluation of kinetic hydrate inhibitors: a procedure for enhancing the repeatability of test results. Energy Fuels 2009, 23 (2), 962-966. 241. May, E. F.; Wu, R.; Kelland, M. A.; Aman, Z. M.; Kozielski, K. A.; Hartley, P. G.; Maeda, N., Quantitative kinetic inhibitor comparisons and memory effect measurements from hydrate formation probability distributions. Chem. Eng. Sci. 2014, 107, 1-12. 242. Wilson, P.; Haymet, A., Hydrate formation and re-formation in nucleating THF/water mixtures show no evidence to support a “memory” effect. Chem. Eng. J. 2010, 161 (1), 146-150. 243. FandiГ±o, O.; Ruffine, L., Methane hydrate nucleation and growth from the bulk phase: Further insights into their mechanisms. Fuel 2014, 117, 442-449. 244. Becker, N.; Kozielski, K.; Haymet, A.; Hartley, P.; Wilson, P. In Nucleation of clathrates from supercooled THF/water mixtures shows that no memory effect exists, 6th International Conference on Gas Hydrates, 2008. 245. Bi, Y.; Li, T., Probing methane hydrate nucleation through the forward flux sampling method. J. Phys. Chem. B 2014, 118 (47), 13324-13332. 246. Lauricella, M.; Meloni, S.; English, N. J.; Peters, B.; Ciccotti, G., Methane clathrate hydrate nucleation mechanism by advanced molecular simulations. J. Phys. Chem. C 2014, 118 (40), 2284722857. 247. MaЕ‚olepsza, E.; Keyes, T., Pathways through equilibrated states with coexisting phases for gas hydrate formation. J. Phys. Chem. B 2015, 119 (52), 15857-15865. 248. English, N. J.; MacElroy, J., Perspectives on molecular simulation of clathrate hydrates: Progress, prospects and challenges. Chem. Eng. Sci. 2015, 121, 133-156.

ACS Paragon Plus Environment

82

Page 83 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

249. Fidler, J.; Rodger, P., Solvation structure around aqueous alcohols. J. Phys. Chem. B 1999, 103 (36), 7695-7703. 250. Matsumoto, M.; Baba, A.; Ohmine, I., Topological building blocks of hydrogen bond network in water. J. Chem. Phys 2007, 127 (13), 134504. 251. Matsumoto, M., Four-body cooperativity in hydrophonic association of methane. J. Phys. Chem. Lett. 2010, 1 (10), 1552-1556. 252. Chakraborty, S. N.; Grzelak, E. M.; Barnes, B. C.; Wu, D. T.; Sum, A. K., Voronoi tessellation analysis of clathrate hydrates. J. Phys. Chem. C 2012, 116 (37), 20040-20046. 253. Jacobson, L. C.; Matsumoto, M.; Molinero, V., Order parameters for the multistep crystallization of clathrate hydrates. J. Chem. Phys 2011, 135 (7), 074501. 254. Steinhardt, P. J.; Nelson, D. R.; Ronchetti, M., Bond-orientational order in liquids and glasses. Phys. Rev. B 1983, 28 (2), 784. 255. Dellago, C.; Bolhuis, P. G., Transition Path Sampling and Other Advanced Simulation Techniques for Rare Events. In Advanced Computer Simulation Approaches for Soft Matter Sciences III, Holm, C.; Kremer, K., Eds. Springer Berlin Heidelberg: Berlin, Heidelberg, 2009; pp 167-233. 256. Nosé, S.; Yonezawa, F., Isothermal–isobaric computer simulations of melting and crystallization of a Lennard‐Jones system. J. Chem. Phys 1986, 84 (3), 1803-1814. 257. Starr, F. W.; Sastry, S.; Douglas, J. F.; Glotzer, S. C., What do we learn from the local geometry of glass-forming liquids? Phys. Rev. Lett. 2002, 89 (12), 125501. 258. Alinchenko, M. G.; Voloshin, V. P.; Medvedev, N. N.; Mezei, M.; Pártay, L. v.; Jedlovszky, P., Effect of cholesterol on the properties of phospholipid membranes. 4. Interatomic voids. J. Phys. Chem. B 2005, 109 (34), 16490-16502. 259. Jacobson, L. C.; Hujo, W.; Molinero, V., Thermodynamic stability and growth of guest-free clathrate hydrates: a low-density crystal phase of water. J. Phys. Chem. B 2009, 113 (30), 10298-10307. 260. Barnes, B. C.; Sum, A. K., Advances in molecular simulations of clathrate hydrates. Curr Opin Chem Eng. 2013, 2 (2), 184-190. 261. Babu, P.; Kumar, R.; Linga, P., Pre-combustion capture of carbon dioxide in a fixed bed reactor using the clathrate hydrate process. Energy 2013, 50, 364-373. 262. Babu, P.; Kumar, R.; Linga, P., Medium pressure hydrate based gas separation (HBGS) process for pre-combustion capture of carbon dioxide employing a novel fixed bed reactor. INT. J. GREENHOUSE GAS CONTROL 2013, 17, 206-214. 263. Linga, P.; Daraboina, N.; Ripmeester, J. A.; Englezos, P., Enhanced rate of gas hydrate formation in a fixed bed column filled with sand compared to a stirred vessel. Chem. Eng. Sci. 2012, 68 (1), 617623. 264. Linga, P.; Haligva, C.; Nam, S. C.; Ripmeester, J. A.; Englezos, P., Gas hydrate formation in a variable volume bed of silica sand particles. Energy Fuels 2009, 23 (11), 5496-5507. 265. Cha, S.; Ouar, H.; Wildeman, T.; Sloan, E., A third-surface effect on hydrate formation. J. Phys. Chem. 1988, 92 (23), 6492-6494. 266. Kneafsey, T. J.; Tomutsa, L.; Moridis, G. J.; Seol, Y.; Freifeld, B. M.; Taylor, C. E.; Gupta, A., Methane hydrate formation and dissociation in a partially saturated core-scale sand sample. J. Pet. Sci. Eng. 2007, 56 (1), 108-126. 267. Linga, P.; Clarke, M., A Review of Reactor Designs and Materials Employed for Increasing the Rate of Gas Hydrate Formation. Energy Fuels 2016, 31 (1), 1-13. 268. Shen, Y. R.; Ostroverkhov, V., Sum-frequency vibrational spectroscopy on water interfaces: polar orientation of water molecules at interfaces. Chem. Rev. 2006, 106 (4), 1140-1154. 269. Bai, D.; Chen, G.; Zhang, X.; Wang, W., Microsecond molecular dynamics simulations of the kinetic pathways of gas hydrate formation from solid surfaces. Langmuir 2011, 27 (10), 5961-5967. 270. Yan, K.-F.; Li, X.-S.; Chen, Z.-Y.; Xia, Z.-M.; Xu, C.-G.; Zhang, Z., Molecular Dynamics Simulation of the Crystal Nucleation and Growth Behavior of Methane Hydrate in the Presence of the Surface and Nanopores of Porous Sediment. Langmuir 2016, 32 (31), 7975-7984.

ACS Paragon Plus Environment

83

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 84 of 89

271. Bai, D.; Chen, G.; Zhang, X.; Wang, W., Nucleation of the CO2 hydrate from three-phase contact lines. Langmuir 2012, 28 (20), 7730-7736. 272. Cygan, R. T.; Guggenheim, S.; Koster van Groos, A. F., Molecular models for the intercalation of methane hydrate complexes in montmorillonite clay. J. Phys. Chem. B 2004, 108 (39), 15141-15149. 273. He, Z.; Linga, P.; Jiang, J., CH4 Hydrate Formation between Silica and Graphite Surfaces: Insights from Microsecond Molecular Dynamics Simulations. Langmuir 2017, 33 (43), 11956-11967 DOI: 10.1021/acs.langmuir.7b02711. 274. Kowalsky, M. B.; Moridis, G. J., Comparison of kinetic and equilibrium reaction models in simulating gas hydrate behavior in porous media. Energy Convers. Manage. 2007, 48 (6), 1850-1863. 275. Li, B.; Li, X.-S.; Li, G., Kinetic studies of methane hydrate formation in porous media based on experiments in a pilot-scale hydrate simulator and a new model. Chem. Eng. Sci. 2014, 105, 220-230. 276. Gamwo, I. K.; Liu, Y., Mathematical modeling and numerical simulation of methane production in a hydrate reservoir. Ind. Eng. Chem. Res. 2010, 49 (11), 5231-5245. 277. Nazridoust, K.; Ahmadi, G., Computational modeling of methane hydrate dissociation in a sandstone core. Chem. Eng. Sci. 2007, 62 (22), 6155-6177. 278. Varini, N.; English, N. J.; Trott, C. R., Molecular dynamics simulations of clathrate hydrates on specialised hardware platforms. Energies 2012, 5 (9), 3526-3533. 279. Veluswamy, H. P.; Wong, A. J. H.; Babu, P.; Kumar, R.; Kulprathipanja, S.; Rangsunvigit, P.; Linga, P., Rapid methane hydrate formation to develop a cost effective large scale energy storage system. Chem. Eng. J. 2016, 290, 161-173. 280. Carpenter, K.; Bahadur, V., Electronucleation for Rapid and Controlled Formation of Hydrates. J. Phys. Chem. Lett. 2016, 7 (13), 2465-2469 DOI: 10.1021/acs.jpclett.6b01166. 281. Shahriari, A.; Acharya, P. V.; Carpenter, K.; Bahadur, V., Metal-Foam-Based Ultrafast Electronucleation of Hydrates at Low Voltages. Langmuir 2017, 33 (23), 5652-5656 DOI: 10.1021/acs.langmuir.7b00913. 282. Wolde, P.-R., Simulation of homogeneous crystal nucleation close to coexistence. Faraday Discuss. 1996, 104, 93-110.

ACS Paragon Plus Environment

84

Page 85 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

TOC/Abstract Graphic 3.33 inches (8.47 cm) wide by 1.875 inches (4.76 cm) deep

Non-technical Synopsis (20 words) A comprehensive review of clathrate hydrate nucleation for the development of sustainable technologies for water, energy and environment is presented

ACS Paragon Plus Environment

85

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 86 of 89

Biography of Authors Maninder Khurana

Dr. Maninder Khurana received his Ph.D. in Chemical Engineering from the National University of Singapore in 2016 under the supervision of Professor Farooq Shamsuzzaman on the topic of integrated optimization of vacuum swing adsorption for post-combustion carbon capture. His Ph.D. research involved adsorbent evalution, first principles modelling, developing in-house process simulation tools and economic analysis of CO2 capture. He is currently working as a Post-doctoral Research Fellow with Professor Praveen Linga and works on gas hydrate kinetic modelling, thermodynamic analysis and numerical simulations for hydrate processes. His current research interests are process systems engineering, developing simulation platforms for novel technologies with application towards CO2 capture, energy storage and clean energy.

ACS Paragon Plus Environment

86

Page 87 of 89

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Sustainable Chemistry & Engineering

Zhenyuan Yin:

Mr. Zhenyuan Yin is currently a Ph.D candidate in the Department of Chemical and Biomolecular Engineering at the National University of Singapore (NUS), Singapore. He is a recipient of the Industrial Postgraduate Program scholarship from Economic Development Board of Singapore and Lloyd’s Register Global Technology Centre. He holds a 1st Class honors degree in Chemical and Biomolecular Engineering from NUS. The primary focus of his research is numerical modelling and reservoir modelling of the kinetic behaviour of methane hydrate formation and dissocaition in porous medium. He is the finalist for the 2017 IChemE Singapore Young Industrialist Award.

ACS Paragon Plus Environment

87

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 88 of 89

Praveen Linga:

Professor Praveen Linga is an associate professor in the Department of Chemical and Biomolecular Engineering at the National University of Singapore (NUS). He is also the co-lead for natural gas research in the centre for energy research & technology (CERT) at the Faculty of Engineering, NUS. His research interests are in the areas of clathrate (gas) hydrates, storage and transport of fuels, carbon dioxide capture, storage & utilization (CCS & U), seawater desalination and recovery of energy. His research group at NUS particularly focuses on enhancing the kinetics of hydrate formation for several applications of interest by developing novel reactor designs, experimental methods and techniques. Up to date, he has published more than 75 research articles and delivered about 50 keynote/invited talks and seminars. He has won numerous local and international awards including the 2017 NUS Young Researcher Award, 2017 NUS Engineering Young Researcher Award, 2017 Energies Young Investigator Award and the 2017 Donald W. Davidson Award for outstanding contributions to gas hydrate research.

ACS Paragon Plus Environment

88

Page 89 of 89

ACS Sustainable Chemistry & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

89