Size-, Aggregation- and Oxidization-Dependent Perturbation of

C , Just Accepted Manuscript. DOI: 10.1021/acs.jpcc.9b02659. Publication Date (Web): May 1, 2019. Copyright © 2019 American Chemical Society. Cite th...
0 downloads 0 Views 2MB Size
Subscriber access provided by IDAHO STATE UNIV

C: Physical Processes in Nanomaterials and Nanostructures

Size-, Aggregation- and Oxidization-Dependent Perturbation of Methane Hydrate by Graphene Nanosheets Revealed by Molecular Dynamics Simulations Shixin Li, Rujie Lv, Yining Wu, Fang Huang, Xianren Zhang, and Tongtao Yue J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.9b02659 • Publication Date (Web): 01 May 2019 Downloaded from http://pubs.acs.org on May 1, 2019

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 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 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.

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 30 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

The Journal of Physical Chemistry

Size-, Aggregation- and Oxidization-Dependent Perturbation of Methane Hydrate by Graphene Nanosheets Revealed by Molecular Dynamics Simulations

Shixin Li†,‡, Rujie Lv‡, Yining Wu§, Fang Huang†, Xianren Zhang||, Tongtao Yue*,†,‡

†State

Key Laboratory of Heavy Oil Processing, China University of Petroleum (East

China), Qingdao, 266580, China ‡College

of Chemical Engineering, China University of Petroleum (East China),

Qingdao, 266580, China §School

of Petroleum Engineering, China University of Petroleum (East China),

Qingdao, 266580, China ||State

Key Laboratory of Organic-Inorganic Composites, Beijing University of

Chemical Technology, Beijing, 100029, China

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

ABSTRACT: Understanding and control of the methane hydrate formation are of central importance for applications ranging from natural gas exploitation to transportation. Fabricated carbon nanomaterials, owing to their outstanding physicochemical properties, are increasingly considered as additives to manipulate the hydrate formation, while little is known about the underlying molecular mechanism. Here, we investigate the methane hydrate formation in the presence of graphene nanosheets (GNs) using molecular dynamics simulation. Particular attention is placed on the effects of size, aggregation and oxidation of GNs. Individual GNs are found to play roles in a size-dependent manner, as sharp corners of GNs are preferentially anchored into cavities at the hydrate surface, exposing other segment in solvent to disturb the local hydrate structure. Once GNs form aggregates exceeding a critical size, methane molecules can be recruited to promote formation of nanobubbles, thus retarding the hydrate formation due to depletion of methane in the aqueous phase. Graphene oxide forms hydrogen bonds with water both in the aqueous phase and at the hydrate surface, thereby reducing the water activity to obstruct the hydrate growth. Our results have important consequences for regulating the methane hydrate formation, and open up new avenues for the energy application of graphene family nanomaterials.

2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30 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

The Journal of Physical Chemistry

INTRODUCTION When guest molecules encounter water under high-pressure and low-temperature, gas hydrate forms. Structures I, II and H are the three types of hydrate identified depending on the size of guest molecules and the environment. In particular, methane hydrate existing in nature primarily as the structure I is the most common type of clathrate hydrate,1, 2 and widespread in the permafrost and deep ocean regions.3, 4 The large amount storage of methane hydrate has attracted considerable efforts as it is expected as a future energy resource. Besides, the gas hydrate technology can be used in many other areas, such as water desalination,5-7 CO2 capture and sequestration,8-11 gas separation,12,

13

food industry,14-17 refrigeration,18,

19

gas storage and

transportation.20-23 However, adverse impacts exist when the gas hydrate is formed or dissolved abnormally. For example, plugging of oil and gas pipelines is a serious industrial problem caused by gas hydrate formation,24, 25 while drastic hydrate melting releases greenhouse gases that cause severe climate change.26 Therefore, understanding and control of the gas hydrate formation are crucial for both promoting applications in relevant areas and avoiding adverse impacts. An elegant solution to these problems is designing delicate inhibitors27-36 and promoters37-47 to regulate formation of gas hydrate. However, due to stochastic and dynamic natures of the gas hydrate evolution, the mode of action and efficacy of additives to influence the hydrate formation are associated with uncertainties. Even results using different equipment and testing procedures are not always repeatable. We conjecture that the inconsistency may arise from the fact that the same additives may influence the gas hydrate formation in different modes, depending on the additive concentration and aggregation, the hydrate formation stage (nucleation and growth), and the hydrate type (I, II and H) encountering the additives. Conventional thermodynamic inhibitors (THIs), such as methanol and ethylene glycol, shift the hydrate formation condition to lower temperatures and higher pressures, but at least 20 wt% of THIs is generally required to suppress the hydrate formation,48 thus increasing the cost and risks of toxicity and flammability. By 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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 30

contrast, kinetic hydrate inhibitors (KHIs) can retard the hydrate formation at concentrations below 1.0 wt%, being more economically favorable.49 KHIs delay the hydrate formation with three possible mechanisms, i.e. binding to the nucleus to prevent it from reaching the critical size, adsorbing onto the hydrate surface to inhibit the crystal growth, and structuring water molecules to restrict recrystallization.50-59 Inspired by the anti-icing activity of natural anti-freezing proteins (AFPs), which were found to anchor their methyl groups into half-cages at the hydrate surface to prevent guest molecules diffusing into cages,60 some analogues were synthesized to exhibit similar properties of hydrate inhibition.61-68 For facilitating the hydrate formation, promoters were designed, following one or more of the principles: acting as seeds to promote the hydrate nucleation,44 lowering surface tension at the gas/liquid interface to enhance the mass transfer from gas to aqueous phase,47, 69 increasing the solubility of gas in the aqueous solution,45,

46, 70

improving heat transfer efficiency,39,

44

and

altering morphology of the growing hydrate crystals.42, 46, 70 Carbon-based nanomaterials (CBNs) can be easily produced from the petroleum coke,71 and show extraordinary mechanical, electrical and thermal properties, making them desirable materials for a wide range of applications.72-77 In particular, different types of CBNs have been studied for the gas hydrate application.39,

78-82

However,

both promotion and inhibition effects of CBNs on the gas hydrate formation were reported,18,

83-87

depending on the experimental procedure, impurity, and apparatus.

Yet more essentially, the mechanisms with which CBNs influence the gas hydrate formation still remain ambiguous. Here, we use all-atom molecular dynamics (MD) simulation to investigate the methane hydrate formation in the presence of graphene nanosheets (GNs), focusing on the effects of size, aggregation and oxidation of GNs. Simulations reveal that individual GNs preferentially anchor sharp corners into cavities at the hydrate-liquid interface, with other segments exposed to disturb the local hydrate structure in a size-dependent manner. Multiple GNs form aggregates to enhance the hydrate perturbation. Once the GN or GN aggregate size exceeds a critical value, methane molecules can be recruited to form nanobubbles, thus retarding the hydrate growth due to depletion of methane in the aqueous phase. Oxidized GNs 4

ACS Paragon Plus Environment

Page 5 of 30 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

The Journal of Physical Chemistry

are found to form hydrogen bonds with water both in the aqueous phase and at the hydrate surface, thus reducing the water activity to impede the hydrate growth, being reminiscent of the anti-freezing activity of oxidized GNs discovered recently.88

MATERIALS AND METHODS Atomistic models and simulation system setup. To understand the effect of GNs on the methane hydrate formation, two model systems were considered in parallel. Firstly, a cubic box of size 5 × 5 × 5 nm3 was generated (Figure 1a), in which a mixture of 256 methane molecules and 2944 water molecules was randomly placed. GNs of three sizes (Figure 1b), each consisting of 22, 58 and 238 carbon atoms, were randomly placed in the aqueous box. For ease of comparison, the number of GNs were respectively set as 11, 4 and 1 to approximately ensure the same total number of carbon atoms. Oxidized GNs were constructed by decorating carbon atoms at the GN edge with carboxyl groups. In the other system, a box of size 3.5 × 3.5 × 9.3 nm3 was constructed, being composed of a gas phase, an aqueous phase, and a hydrate phase (Figure 1c). The gas phase consists of either 271 or 188 methane molecules, corresponding to different gas concentrations. The numbers of methane and water molecules in the middle aqueous phase were 100 and 1656, respectively.47 The hydrate phase was represented by a 3 × 3 × 2 supercell of SI crystal structure to serve as a template for the directional hydrate growth. GNs were positioned in the aqueous phase 0.6 nm close to the hydrate surface (Figure 1c).

Figure 1. Atomistic models and simulation system setup. (a) The cubic box of a 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

mixture of methane and water molecules, in which GNs were randomly positioned. (b) GNs of three different sizes used in simulations. Carbon atoms at edge of the GN of 58 carbon atoms are decorated with carboxyl groups to construct the oxidized GN. (c) The other simulation system containing a preexisted hydrate phase, a middle aqueous phase, and a gas phase. The lower panel of (c) shows the initial distribution of methane molecules along the z axis. Details of MD simulations. During MD simulations, water molecules were described by the TIP4P/ice model, and methane molecules were represented by the OPLS/AA model.89, 90 The force field parameters for GNs were obtained from Tu et al.91 Short-range interactions were truncated at 1.2 nm, and long-range coulomb interactions were calculated using the particle-mesh Ewald summation method.92 The standard Lorentz-Berthelot mixing rule was used to treat cross interactions between different molecules.93 Motion equations were integrated using the leapfrog algorithm with a time step of 1.0 fs.94 After energy minimization with the steepest descent algorithm, the system was relaxed with the NVT ensemble for 30 ps. Then the system was equilibrated for hundreds of nanoseconds with the NPT ensemble. Pressure and temperature were kept constant using the Parrinello−Rahman barostat and the Nose-Hoover thermostat with time constants of 2 ps and 0.4 ps, respectively. The pressure coupling was applied along all three directions of the simulation box isotropically for both systems. We note that using thermostat and barostat can accelerate the hydrate nucleation, albeit at the cost of reducing the hydrate crystallinity. That is because a barostat retains gas-phase pressure and a thermostat removes the exothermic heat of hydrate formation.95 Considering that our simulations mainly focused on the effect of GNs on the methane hydrate formation, the NPT ensemble was used in the production run for accelerating the methane hydrate formation. For the selection of temperature, it was previously estimated that the melting temperature of the SI methane hydrate is ~ 304 K at P = 500 MPa.96 For a mixture of methane and water molecules randomly placed in a cubic simulation box, we chose a higher temperature of 275 K for reducing the hydrate nucleation time. To achieve the regular hydrate growth with less amorphous fragments, we decreased the 6

ACS Paragon Plus Environment

Page 6 of 30

Page 7 of 30 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

The Journal of Physical Chemistry

temperature to 250 K for the system containing a preexisted hydrate phase. Periodic boundary conditions were applied in all three directions. All MD simulations were performed using the open source code Gromacs 4.6.7.97 Snapshots were rendered by VMD.98 F4 parameter. The four-body structural order parameter (F4), which was defined as F4 

1 n  cos 3i , was used to evaluate degree of the hydrate formation, where n n i 1

is the total number of H2O-H2O pairs with the O-O distance smaller than 3.5 Å, and

i is the torsion angle formed by the oxygen atoms and two outermost hydrogen atoms in the i-th H2O-H2O pair. In theory, the value of F4 parameter varies from -0.4 to 0.7, with specific values of -0.4, -0.04 and 0.7 representing ice, liquid and hydrate phases, respectively.99

Hydrate cage identification. The hydrate cage structure was identified based on the connection of water molecules and the topology of rings they formed.100 The whole procedure of the hydrate cage identification is described in detail as follows: (i) Two water molecules were defined to be connected if the O-O distance was less than 3.5 Å. (ii) Search for all water molecules to identify the non-repetitive triplets of successive vertices. (iii) List all possible five-membered and six-membered rings formed by connected water molecules using the code by Matsumoto et al.24 (iv) Search for five-membered and six-membered rings having five or six other rings connected on each edge. This structure was defined as the half-cage. (v) Other structures with at least four inter-connected rings were referred to as amorphous fragments. (vi) If edges of two half-cages coincide with each other, a full hydrate cage was identified.

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

RESULTS Methane hydrate formation in the absence of GNs. We started our simulations by examining the methane hydrate formation in the absence of GNs. To that end, a mixture of 256 methane and 2944 water molecules was randomly placed in a cubic box of size 5 × 5 × 5 nm3 (Figure 1a). Note that the solvated methane concentration can influence the nucleation time with a direct effect on the methane growth rate.101-103 As predicted by Lauricella et al. using the nucleation theory, the hydrate can form in months when the methane concentration is approximately 0.01.102 Even in the experimental range between 0.015 and 0.0175, the predicted nucleation time reaches 1-60 seconds, which is still not accessible by current MD simulations. Therefore, in our simulations, the methane concentration was increased to 0.08 to decrease the nucleation time for reduction of the simulation time. Shown in Figure 2a is the time sequence of typical snapshots depicting how the methane hydrate was formed from a mixture of methane and water. Under conditions of T = 275 K and P = 500 bar,104 a number of fragments first appeared via forming hydrogen bonds between water molecules. Such amorphous structures were unstable and gradually rearranged to decrease the system free energy (Figure 2b). At t = 15 ns, two half-cages were formed, from which the hydrate started to nucleate and grow into a larger structure. We monitored numbers of water molecules participating in the formation of both amorphous fragments and regular cages, and found that the whole process can be divided into two stages (Figure 2c). The first stage was featured with formation of fragments, as reflected by striking increases of both the F4 parameter and the number of water molecules in the fragment (Figure 2c, d). Once stable half-cages appeared to initiate the hydrate formation, the number of water molecules belonging to the fragment stopped increasing, and the increase of the F4 parameter became less apparent, suggesting that the growth of hydrate was accomplished via transformation from fragmental to cage-like structures. At the end of the simulation, a large hydrate structure was formed, in which most methane molecules were encapsulated, as reflected by the distribution of methane molecules along the x direction (Figure 2e). 8

ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30 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

The Journal of Physical Chemistry

Figure 2. Methane hydrate formation in the absence of GN. (a) Time sequence of typical snapshots, with amorphous clathrate fragments consisting of at least four five membered and six membered rings sharing edges displayed in blue, and half cages in red. (b) Time evolution of the system potential energy. (c) Time evolutions of the number of water molecules participating in the formation of fragmental and cage-like structures. (d) Time evolution of the average F4 parameter for all water molecules. (e) Distributions of the methane concentration along the x direction at different time points. The system temperature and pressure are 275 K and 500 bar, respectively. The total number of methane and water molecules are 256 and 2944, respectively. 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

For the system containing a hydrate phase at bottom of the simulation box, we decreased the temperature to 250 K to probe the directional growth of the methane hydrate, skipping the first step of hydrate nucleation. Note that the lower temperature can effectively inhibit formation of amorphous fragments for better accomplishing regular hydrate growth.104 As expected, we observed a continuous growth of the hydrate along the z direction (Figure S1a), as reflected by gradual increases of both the F4 parameter and the cage number (Figure S1b and c). Different from the case under a higher temperature of 275 K, where amorphous fragments preferentially appeared in the early stage, both the presence of a hydrate phase and the lower temperature led to a directional growth of regular hydrates along the templet. By contrast, less fragments formed during the simulation (Figure S1d). As methane molecules in the aqueous phase were depleted, the rate of hydrate growth was decreased and mainly dependent on the rate of transferring methane molecules from the gas to the aqueous phase.47 If we decreased the number of methane molecules in the gas phase from 271 to 188, water molecules in the aqueous phase were found to penetrate into the gas phase to open a channel, which resulted in a curved gas-water interface. It has been demonstrated by Walsh et al. that the concentration of methane in water increased with the interfacial curvature due to the effective pressure increase in the methane phase as governed by the Young-Laplace equation.69 As expected, the appearance of a water channel facilitated the gas-water mixing to promote the hydrate formation (Figure S2a), as quantified by the F4 parameter (Figure S2b), the cage number (Figure S2c), and the system potential energy (Figure S2d). Given ubiquity of the gas-water interfacial curvature, such effect is thought to play roles frequently in real systems. For example, the porous media with small pore size was found to shorten the induction time of hydrate formation, arising from the enhanced water-methane mass transfer efficiency.105

Minor effect of small GNs on the methane hydrate formation. Figure 3a depicts the methane hydrate formation in the presence of 11 small GNs, each 10

ACS Paragon Plus Environment

Page 10 of 30

Page 11 of 30 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

The Journal of Physical Chemistry

consisting of 22 carbon atoms. Both the F4 parameter and the hydrate cage number evolved resembling those in the absence of GNs (Figure 3b and c), suggesting that the rate of methane hydrate formation is barely affected by small GNs. The lower potential energy of the system was ascribed to extra interactions involving GNs (Figure 3d). Initially, all GNs located randomly in the aqueous box. Similar with the case in the absence of GNs, amorphous fragments first appeared, but distant from and repelling GNs into the diminishing aqueous phase. As the simulation proceeded, more hydrate structures formed to further decrease the volume of the aqueous phase accommodating insoluble GNs. To decrease unfavorable interactions with surrounding water, GNs formed aggregates in the aqueous phase. Both the methane hydrate formation and the GN aggregation caused decreases of GN contacts with water and methane molecules (Figure 3e).

Figure 3. The methane hydrate formation in the presence of 11 small GNs, each consisting of 22 carbon atoms. (a) Time sequence of the typical snapshot. (b) Time evolutions of the average F4 parameter in the absence and presence of GNs. (c) Time evolutions of the cage number. (d) Time evolutions of the system potential energy. (e) Time evolutions of the number of methane and water molecules in contact with GNs. Amorphous fragments are shown in blue, half cages in red, and GNs in black.

Uncertainties exist in above simulations due to the stochastic nature of methane hydrate nucleation and the heterogeneity of hydrate growth. We next skipped the step of hydrate nucleation and directly positioned a small GN of 22 carbon atoms 0.6 nm above the surface of a preexisted hydrate phase at bottom of the simulation box. 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

Shown in Figure 4a is the time sequence of typical snapshots depicting how the directional growth of methane hydrate was affected by a small GN. We observed an early hydrate growth along the z axis, except the region around the GN. Under the lower concentration of methane molecules in the gas phase, a water channel was rapidly formed, through which a number of water molecules penetrated into the gas phase to promote the gas-water mixing. As the simulation proceeded, we observed coalescence of the hydrate growing respectively from top and bottom the box. Such rapidly growing structure was found to wrap the GN, though the local hydrate structure was disturbed by the small GN (Figure 4b). Apart from that, the lower value of F4 parameter at the late simulation period compared to that in the absence of GN was mainly ascribed to the hydrate formation at top of the box (Figure 4c), which surrounded the gas phase to retard diffusion of methane molecules into water to sustain further growth of the hydrate. With such occasional event being eliminated, it can be concluded that rather small GNs play minor roles in the methane hydrate formation.

Figure 4. The methane hydrate growth in the presence of a small GN of 22 carbon atoms initially positioned 0.6 nm above the preexisted hydrate surface. (a) Time sequence of the typical snapshot. (b) The enlarged structure depicting the perturbed local hydrate structure by the GN. (c) Time evolutions of the average F4 parameter for water in systems with and without the GN.

Enhanced perturbation of methane hydrate by larger GNs. Next, a larger GN of 58 carbon atoms was positioned above the methane hydrate surface. Interestingly, the methane hydrate stopped growing upon encountering the GN (Figure 5a). Moreover, a half cage structure was identified, in which a sharp corner of the GN was 12

ACS Paragon Plus Environment

Page 12 of 30

Page 13 of 30 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

The Journal of Physical Chemistry

stably anchored, exposing other segment in the surrounding aqueous phase to disturb formation of hydrate around the GN (Figure 5b). Nearly no change in the GN position along the z direction was observed in the whole simulation time (Figure 5c). Again, at the lower concentration of methane in the gas phase, water molecules were found to penetrate into the gas phase through opening a water channel, thus promoting the gas-water mixing due to effect of the Laplace pressure. Nevertheless, no cage-like structures were formed near the GN in the rest of the simulation period, as further reflected by evolution of the hydrate cage number (Figure 5d), manifesting that the methane hydrate perturbation can be enhanced by increasing the GN size.

Figure 5. The enhanced perturbation of methane hydrate by a larger GN of 58 carbon atoms. (a) Time sequence of the typical snapshot. (b) Local enlarged structure showing stable anchorage of a sharp corner of the GN into a cavity. (c) Time evolution of the GN position along the z axis. (d) Time evolutions of the system potential energy in the absence and presence of GNs.

As four GNs of 232 carbon atoms in total were randomly positioned in the 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

aqueous box of 256 methane and 2944 water molecules, interestingly, the methane hydrate formation was strongly inhibited by GNs, but with a different mechanism. In the early stage, amorphous fragments first appeared distant from GNs (Figure 6a), resulting in an increase of the F4 parameter (Figure 6b). Concurrently, GNs formed an aggregate in the aqueous phase to reduce unfavorable contacts with water. Also due to the hydrophobicity, nearby methane molecules were adsorbed by the GN aggregate (Figure 6c), with the methane concentration in the aqueous phase accordingly decreased to retard the hydrate formation (Figure S3). It has been proved by earlier simulations that a concentration lower than 0.03 in molar fraction of methane solvated in water would make the hydrate nucleation a rare event in the time accessible by standard MD simulations.106 As the fraction was still larger than 0.03, we observed a gradual increase of the F4 parameter in the early 100 ns, informing a continuous growth of the methane hydrate, albeit with a lower rate compared to that in the absence of GNs (Figure 6b). Shortly after t = 124 ns, four GNs were rearranged to further decrease the interaction energy between GNs (Figure S4). Upon forming a layered aggregate, a larger number of methane molecules were instantly recruited and adsorbed by the GN aggregate to form a nanobubble (Figure 6c). Accordingly, the preformed hydrate structure was rapidly melted, with only several amorphous fragments remained in the aqueous box.

Figure 6. The methane hydrate formation in the presence of four GNs, each consisting of 58 carbon atoms. (a) Time sequence of the typical snapshot. (b) Time evolutions of the F4 parameter for systems in the absence and presence of GNs. (c) 14

ACS Paragon Plus Environment

Page 14 of 30

Page 15 of 30 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

The Journal of Physical Chemistry

Time evolution of the number of methane molecules inside the bubble around the GNs.

Note that smaller GNs also aggregated but did not exhibit inhibition effect on the methane hydrate formation (Figure 3), indicating that the methane adsorbing capability of GNs contributing to the hydrate suppression is dependent on the GN aggregate size and stability. Such property can be reflected in the evolution of methane hydrate in the presence of four larger GNs (Figure 6a), showing competitive adsorption of methane by the hydrate structure and the GN aggregate. In particular, before t = 25 ns, four GNs first formed a disordered aggregate that recruited surrounding methane molecules to form a nanobubble which, however, was unstable and rapidly dissolved for the hydrate growth in the following 100 ns (Figure 6b, c). Upon forming an ordered GN aggregate via rearrangement (Figure S3), a sudden decrease of the F4 parameter and an increase of the nanobubble size were observed (Figure 6b, c), suggesting that the preformed hydrate structure was rapidly melted due to the enhanced capability of methane recruitment by the ordered GN aggregate. Note that the aggregation of GNs can generate a local confined space, which may ease the condition for the methane hydrate growth as revealed by earlier experiments.107 However, due to the small size and mobility of GNs in our simulations, such subtle confinement can be ignored when analyzing the effect of GNs on the methane hydrate formation. We asked whether the formed structure of nanobubbles encapsulating GN aggregates is stable for sufficiently long time to generate a long-term inhibition effect. Shown in Figure 7a is the transient distribution of methane molecules at t = 182 ns as a function of the distance to the center of aggregated GNs. Combining the typical snapshot (Figure 7b), no methane molecules located at the layer gap, and the major peak reveals the dense adsorption of methane molecules on the aggregated GN surface. We tracked positions of these molecules and surprisingly found that they diffused freely inside the nanobubble, rather than constantly adsorbing onto the GN surface (Figure 7c). To check stability of the nanobubble, we extended the above 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

simulation to 500 ns and no apparent change of either the F4 parameter and the number of methane molecules around the GN aggregate was observed in the following 250 ns, until the hydrate structure reappeared at t = 450 ns around the nanobubble, as reflected by an increase of the F4 parameter and a decrease of the bubble size (Figure 6b and c). Besides the freely diffusing methane molecules inside the nanobubble making the local methane concentration higher at the air-water interface (Figure 7), the π-π stacking of GNs became less stable inside the hydrophobic nanobubble environment than that in the aqueous phase. A further rearrangement of GNs was observed (Figure 6a), which in turn reduced the stabilizing effect on the nanobubble. As a consequence, more methane molecules were found to diffuse from the nanobubble into the aqueous phase to restart the hydrate formation (Figure S3a).

Figure 7. Characterization of methane molecules inside the bubble around a GN aggregate. (a) The number of methane molecules around GNs as a function of the distance to the GN aggregate center along the layer normal direction. (b) Transient structure at t = 182 ns showing methane molecules around the aggregate. Methane molecules adsorbed on the first layer are shown in red for tracking their movement. (c) Transient structure at t = 183 ns, showing that methane molecules diffused freely inside the nanobubble.

As we further increased the GN size to that of 238 carbon atoms, nearly no stable cage-like hydrate structures formed in the finite simulation period (Figure 8a). Compared to the former case of four GNs of 232 carbon atoms in total, the single large GN of 238 carbon atoms absorbed a larger number of methane molecules more 16

ACS Paragon Plus Environment

Page 16 of 30

Page 17 of 30 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

The Journal of Physical Chemistry

rapidly (Figure 8b), as further reflected by a continuous decrease of the molar fraction of methane solvated in water (Figure S3b). Due to the decreased methane concentration in water, the formation of hydrate was retarded from start of the simulation, as no increase of the F4 parameter was observed during the simulation (Figure 8c). The molar fraction of methane solvated in water was found to fluctuate around the value of 0.03, below which the hydrate nucleation was predicted to rarely occur in the time accessible by standard MD simulations.106 Therefore, no methane hydrate was formed in the finite simulation time (Figure 8). Such inhibition effect of larger GNs on the methane hydrate formation was further verified by another simulation, where a large GN was positioned close to a preformed hydrate surface (Figure S5). Simulation results showed that, before onset of the methane hydrate growth, the GN rapidly absorbed a large number of methane molecules in the aqueous phase (Figure S5a), as reflected by a peak of the methane concentration at around 2 nm at t = 300 ns (Figure S5b). Then, the GN carrying a methane nanobubble moved randomly in the aqueous phase. Upon getting contact with the nearby gas phase, the nanobubble coalesced to increase the layer thickness of the gas phase (Figure S5b). Simultaneously, we observed a sudden decrease of the F4 parameter indicating perturbation of the hydrate structure (Figure S5c). In the rest of the simulation time, nearly no methane molecules were found to diffuse back into the aqueous phase to support the hydrate growth, manifesting the size-dependent stabilizing effect of GNs on the nanobubble.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

Figure 8. Inhibition of methane hydrate formation by a larger GN of 238 carbon atoms. (a) Typical snapshots showing that the larger GN rapidly adsorb a large number of methane molecules to form a stable nanobubble, accordingly decreasing the gas concentration in the aqueous phase to retard the hydrate formation. (b) Time evolution of the number of methane molecules around the GN. (c) Time evolutions of the F4 parameter showing the inhibition effect of large GNs on the hydrate formation.

Above simulations convincingly showed that the effect of GNs on the methane hydrate formation is strongly influenced by the size of GNs or GN aggregates. Besides the size dependent perturbation of local hydrate structures by individual GNs (Figures 4 and 5), multiple GNs form aggregates exhibiting an ability of recruiting methane molecules to form nanobubbles that reduce the methane concentration in the aqueous phase to retard the hydrate formation. We showed the detailed arrangement of methane molecules around GNs of different sizes, and found that nearly no methane molecules were adsorbed by a rather small GN of 22 carbon atoms (Figure 9a). By contrast, an average of 12 methane molecules was found to adsorb on the GN of 58 carbon atoms (Figure 9b). As we increased the GN size to that of 238 carbon atoms, the average number of methane molecules closely adsorbed onto each GN surface was increased to 34 (Figure 9c). Despite the finding that methane molecules diffused freely inside the nanobubble rather than constantly attaching to the GN surface (Figure 7), we have calculated the distribution of distance between neighboring methane molecules both adsorbed on the GN surface and in a nanobubble. Interestingly, the inter-methane distance on the largest GN surface was found to range from 0.36 nm to 0.42 nm (Figure 9d). By contrast, the inter-methane distance was distributed wider in the nanobubble phase (Figure S6), suggesting stabilization of larger GNs on the nanobubble via forming ordered and regular distribution of methane molecules. Note that the lateral size of the small GN is only 0.3 nm, smaller than the shortest distance measured between neighboring methane molecules either adsorbing on the GN surface or in the nanobubble, thus explaining why rather small GNs, albeit forming aggregates in the aqueous phase, showed no inhibition effect on 18

ACS Paragon Plus Environment

Page 18 of 30

Page 19 of 30 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

The Journal of Physical Chemistry

the methane hydrate formation.

Figure 9. The size effect of GNs on adsorbing methane molecules. (a) Methane molecules around smaller GNs, each consisting of 22 carbon atoms. (b) Methane molecules in close contact with larger GNs of 58 carbon atoms. (c) Methane molecules adsorbed on the surface of a larger GN of 238 carbon atoms. (d) Distribution of the distance between adjacent methane molecules adsorbed on the GNs averaged from (b) and (c).

Graphene oxide. Although aggregation of GNs has been proved to retard the methane hydrate formation by extracting methane molecules to form nanobubbles, in real systems, such behavior may induce the GN deposition that reduces the inhibition effect and generates adverse impacts, thus requiring caution to control the GN dosage to prevent large scale GNP aggregation. Aggregation of GNs can be reduced by oxidation. Shown in Figure 10 is the methane hydrate evolution in the presence of a GN of 58 carbon atoms with all edge atoms decorated with carboxyl groups (Figure 1b). Before encountering the oxidized GN, the hydrate was found to grow along the normal direction (Figure 10a), as reflected by an increase of the F4 parameter before t = 100 ns (Figure 10b). However, a number of hydrogen bonds were formed between carboxyl groups of the GN and water both in the aqueous phase and at the hydrate surface (Figure 10c). We infer that both reduction of the water activity in solution and 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

occupancy of hydrate-liquid interface by GNs prevent water forming new hydrogen bonds with the hydrate surface to support the hydrate growth. In the rest of the simulation time, nearly no further growth of the hydrate was observed, suggesting that the hydrate structure can be perturbed by oxidized GNs.

Figure 10. The effect of an oxidized GN on the methane hydrate growth. (a) Typical snapshots. (b) Time evolutions of the F4 parameter for water in systems with and without oxidized GNs. (c) Locally enlarged configuration depicting hydrogen bonds formed between carboxyl groups of the GN and water molecules both in the aqueous phase and at the hydrate surface.

Finally, four oxidized GNs were positioned randomly in the aqueous box of 256 methane and 2944 water molecules. As expected, no aggregation of oxidized GNs formed in the finite simulation period (Figure S7). Since individual oxidized GNs reduced activity of surrounding water molecules, methane hydrate was found to form distant from GNs. After hydrate nucleation, GNs form hydrogen bonds with water at the hydrate surface to retard growth of hydrate only along specific directions. Different from the case with a preexisted hydrate phase that was regulated to grow along only the surface normal direction, the hydrate growth in the aqueous phase is 20

ACS Paragon Plus Environment

Page 20 of 30

Page 21 of 30 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

The Journal of Physical Chemistry

more heterogeneous, thus reducing the inhibition effect of oxidized GNs on the hydrate formation. Besides, no methane molecules were recruited to form nanobubbles that were found to retard the hydrate formation in the case of pristine GNs.

DISCUSSION Advances have been made in understanding methane hydrate formation in bulk phase,108 while the effect of nanoparticles, such as those present under conditions relevant to industry and environment or acting as additives to manipulate the hydrate formation, remains ambiguous. Combining neutron scattering experiments and MD simulations, it was found by Cox et al. that the formation of methane hydrate is insensitive to addition of nanoparticles. They explained that methane and water have different chemical natures that make it unlikely for the two species simultaneously attracted by nanoparticles so as to promote their mixing.78 Such explanation, although straightforward, should apply to large particles as they used planar solid surfaces in simulations to find the hydrate nucleation away from the particle surface. By contrast, the GNs considered in this study have sizes comparable to that of the hydrate cage. We found that sharp corners of GNs can be anchored into cavities at the hydrate surface, exposing other segment in the aqueous phase to disturb the local hydrate structure (Figure 5). For oxidized GNs, they were found to form hydrogen bonds with water both in the aqueous phase and at the hydrate surface, thus reducing the water activity to impede the hydrate growth (Figure 9), being reminiscent of the recently discovered anti-freezing activity.88 Earlier experimental results showed that the methane hydrate formation can be promoted by a variety of nanoparticles, with such promotion effect reflected as both the decreased induction time and the increased rate of gas consumption.39, 44, 79 NPs were expected to provide seeds for heterogeneous hydrate nucleation,83 different from simulation results showing that the methane hydrate nucleation preferentially occurred distant from nanoparticles.78 Besides, no further decrease of the induction time was observed by increasing hydrophilicity of nanoparticles, though the dispersity 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

can be enhanced,109 suggesting that other mechanisms may exist for NPs affecting the hydrate nucleation. During the hydrate growth, the rate of gas consumption was generally measured to reflect the effect of nanoparticles on the hydrate formation. However, the increase of gas consumption does not strictly reflect larger amount of hydrate. Our simulations suggested that GNs can form aggregates via π-π stacking interactions, and once the size of GNs or GN aggregates exceeds a critical value, surprisingly, a large number of methane molecules dissolved in the aqueous phase are recruited by GNs to form nanobubbles (Figure 6). Such behavior can also be relevant to an increase of the gas consumption, while the formation of gas hydrate is actually retarded due to the lower local gas concentration. Since the methane hydrate formation is regarded as an exothermic event occurring preferentially at the gas-liquid interface, any agent being able to enhance heat and mass transfer at the interface should enhance the rate of hydrate formation. Once forming aggregates, the heat transfer enhancement by nanoparticles can be lowered.110 In this regard, increasing the hydrophilicity of nanoparticles may enhance the thermal conductivity by preventing their aggregation. However, simulations presented in this work were conducted with the NPT ensemble, thus failing to measure the heat transfer influenced of GNs.

CONCLUSIONS In summary, we have demonstrated, using atomistic MD simulations, that the methane hydrate formation is influenced by GNs, depending on the GN size, aggregation and oxidation. Rather small GNs play minor roles in the hydrate formation as they can be readily wrapped by surrounding hydrate structures. For larger GNs, only sharp corners are preferentially anchored into open cavities at the hydrate surface, exposing other segment in the aqueous phase to disturb the local hydrate structure. Due to hydrophobicity, GNs form aggregate in the aqueous phase. Once the size of GNs or GN aggregates exceeds a critical value, methane molecules can be recruited by GNs to form nanobubbles, thus decreasing the local gas concentration in the aqueous phase to retard the hydrate formation. Aggregation of 22

ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30 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

The Journal of Physical Chemistry

GNs can be effectively reduced by oxidation, which influences the methane hydrate formation via forming hydrogen bonds with water both in the aqueous phase and at the hydrate surface. These findings help us to understand the molecular mechanism with which the gas hydrate formation is influenced by carbon-based nanomaterials, thus shedding lights on the design of novel additives to regulate the gas hydrate formation.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website. Additional simulation results, including methane hydrate formation in the absence of GN under different methane concentrations in the gas phase, solvated methane mole fractions, ordered GN aggregation via rearrangement and rapid formation of nanobubbles, directional methane hydrate growth in the presence of a large GN of 238 carbon atoms, distribution of inter-methane distance in a nanobubble, and heterogeneous methane hydrate formation in the presence of multiple oxidized GNs.

ACKNOWLEDGMENTS We are grateful for financial support from the National Natural Science Foundation of China (no. 31871012), and the Natural Science Foundation of Shandong Province (ZR2018MC004). This work was also partly supported by the Fundamental Research Funds for the Central Universities (no. 19CX07002A). Simulations were performed at the National Supercomputing Center in Shenzhen.

AUTHOR INFORMATION Corresponding author *E-mail: [email protected] Notes 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

The authors declare no competing financial interest.

REFERENCES 1. Makogon, I. F.; Makogon, I. U. r. F.; Makogon, Y. F. Hydrates of hydrocarbons. Pennwell Books: 1997. 2. Sloan, E. D.; Koh, C. A. Clathrate hydrates of natural gases. CRC press: 2007. 3. Koh, C. A.; Sloan, E. D. Natural gas hydrates: Recent advances and challenges in energy and environmental applications. Aiche. J. 2007, 53, 1636-1643. 4. Ríos-Mercado, R. Z.; Borraz-Sánchez, C. Optimization problems in natural gas transportation systems: A state-of-the-art review. Appl. Energ. 2015, 147, 536-555. 5. Javanmardi, J.; Moshfeghian, M. Energy consumption and economic evaluation of water desalination by hydrate phenomenon. Appl. Therm. Eng. 2003, 23, 845-857. 6. Park, K.-n.; Hong, S. Y.; Lee, J. W.; Kang, K. C.; Lee, Y. C.; Ha, M.-G.; Lee, J. D. A new apparatus for seawater desalination by gas hydrate process and removal characteristics of dissolved minerals (Na+, Mg2+, Ca2+, K+, B3+). Desalination 2011, 274, 91-96. 7. Babu, P.; Nambiar, A.; He, T.; Karimi, I. A.; Lee, J. D.; Englezos, P.; Linga, P. A review of clathrate hydrate based desalination to strengthen energy–water nexus. ACS Sustain. Chem. Eng. 2018, 6, 8093-8107. 8. Babu, P.; Linga, P.; Kumar, R.; Englezos, P. A review of the hydrate based gas separation (HBGS) process forcarbon dioxide pre-combustion capture. Energy 2015, 85, 261-279. 9. Dashti, H.; Yew, L. Z.; Lou, X. Recent advances in gas hydrate-based CO2 capture. J. Nat. Gas. Sci. Eng. 2015, 23, 195-207. 10. Fan, S.; Li, S.; Wang, J.; Lang, X.; Wang, Y. Efficient capture of CO2 from simulated flue gas by formation of TBAB or TBAF semiclathrate hydrates. Energ. Fuel 2009, 23, 4202-4208. 11. Linga, P.; Kumar, R.; Englezos, P. The clathrate hydrate process for post and pre-combustion capture of carbon dioxide. J. Hazard. Mater. 2007, 149, 625-629. 12. Seo, Y.; Tajima, H.; Yamasaki, A.; Takeya, S.; Ebinuma, T.; Kiyono, F. A new method for separating HFC-134a from gas mixtures using clathrate hydrate formation. Environ. Sci. Technol. 2004, 38, 4635. 13. Zhong, D.; Englezos, P. Methane separation from coal mine methane gas by tetra-n-butyl ammonium bromide semiclathrate hydrate formation. Energ. Fuel 2012, 26, 2098-2106. 14. Takeya, S.; Nakano, K.; Thammawong, M.; Umeda, H.; Yoneyama, A.; Takeda, T.; Hyodo, K.; Matsuo, S. CO2 processing and hydration of fruit and vegetable tissues by clathrate hydrate formation. Food Chem. 2016, 205, 122-128. 15. Zemelman, V. B.; Kleiner, F.; Kuchman, M. J. Gasified ice product and process having improved storage stability. Google Patents: 1982. 16. Sato, T.; Takeya, S.; Nagashima, H. D.; Ohmura, R. Preservation of carbon dioxide clathrate hydrate coexisting with sucrose under domestic freezer conditions. J. Food. Eng. 2014, 120, 69-74. 17. Englezos, P. The freeze concentration process and its applications. Asia-Pac. J. Chem. Eng. 1994, 2, 3-15. 18. Shi, X. J.; Zhang, P. A comparative study of different methods for the generation of tetra-n-butyl ammonium bromide clathrate hydrate slurry in a cold storage air-conditioning 24

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30 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

The Journal of Physical Chemistry

system. Appl. Energ. 2013, 112, 1393-1402. 19. Zhang, P.; Shi, X. J.; Ma, Z. W. Solid fraction determination in cold storage by tetra-n-butyl ammonium bromide clathrate hydrate slurry. Int. J. Refrig. 2013, 36, 809-819. 20. Florusse, L. J.; Peters, C. J.; Schoonman, J.; Hester, K. C.; Koh, C. A.; Dec, S. F.; Marsh, K. N.; Sloan, E. D. Stable low-pressure hydrogen clusters stored in a binary clathrate hydrate. Science 2004, 306, 469-471. 21. 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. 22. Veluswamy, H. P.; Kumar, S.; Kumar, R.; Rangsunvigit, P.; Linga, P. Enhanced clathrate hydrate formation kinetics at near ambient temperatures and moderate pressures: Application to natural gas storage. Fuel 2016, 182, 907-919. 23. 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. 24. Matsumoto, M.; Baba, A.; Ohmine, I. Topological building blocks of hydrogen bond network in water. J. Chem. Phys. 2007, 127, 134504. 25. Chen, J.; Yan, K.-L.; Chen, G.-J.; Sun, C.-Y.; Liu, B.; Ren, N.; Shen, D.-J.; Niu, M.; Lv, Y.-N.; Li, N. Insights into the formation mechanism of hydrate plugging in pipelines. Chem. Eng. Sci. 2015, 122, 284-290. 26. Rodhe, H. A comparison of the contribution of various gases to the greenhouse effect. Science 1990, 248, 1217-1219. 27. Zheng, Z. Molecular dynamics simulations on the inhibition of methane hydrates. Digital Repository 2010. 28. Nakarit, C.; Kelland, M. A.; Liu, D.; Chen, Y. X. Cationic kinetic hydrate inhibitors and the effect on performance of incorporating cationic monomers into N-vinyl lactam copolymers. Chem. Eng. Sci. 2013, 102, 424-431. 29. Villano, L. D.; Kelland, M. A. An investigation into the kinetic hydrate inhibitor properties of two imidazolium-based ionic liquids on Structure II gas hydrate. Chem. Eng. Sci. 2010, 65, 5366-5372. 30. Xiao, C.; Adidharma, H. Dual function inhibitors for methane hydrate. Chem. Eng. Sci. 2009, 64, 1522-1527. 31. Creek, J. L. Efficient hydrate plug prevention. Energ. Fuel 2012, 26, 4112-4116. 32. Tariq, M.; Rooney, D.; Othman, E.; Aparicio, S.; Atilhan, M.; Khraisheh, M. Gas hydrate inhibition: A review of the role of ionic liquids. Ind. Eng. Chem. Res. 2016, 53, 17855-17868. 33. Huo, Z.; Freer, E.; Lamar, M.; Sannigrahi, B.; Knauss, D. M.; Jr, E. D. S. Hydrate plug prevention by anti-agglomeration. Chem. Eng. Sci. 2001, 56, 4979-4991. 34. Reyes, F. T.; Kelland, M. A. Investigation of the kinetic hydrate inhibitor performance of a series of copolymers of N-vinyl azacyclooctanone on structure II gas hydrate. Energ. Fuel 2013, 27, 1314-1320. 35. Nazari, K.; Moradi, M. R.; Ahmadi, A. N. Kinetic modeling of methane hydrate formation in the presence of low ‐ dosage water ‐ soluble ionic liquids. Chem. Eng. Technol. 2013, 36, 1915-1923. 36. Tohidi, B.; Anderson, R.; Chapoy, A.; Yang, J.; Burgass, R. W. Do we have new solutions to the old problem of gas hydrates? Energ. Fuel 2012, 26, 4053–4058. 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

37. Karaaslan, U.; Parlaktuna, M. Surfactants as hydrate promoters? Energ. Fuel 2000, 14, 1103-1107. 38. Bhattacharjee, G.; Choudhary, N.; Kumar, A.; Chakrabarty, S.; Kumar, R. Effect of the amino acid l-histidine on methane hydrate growth kinetics. J. Nat. Gas. Sci. Eng. 2016, 35, 1453-1462. 39. Hosseini, M.; Ghozatloo, A.; Shariaty-Niassar, M. Effect of CVD graphene on hydrate formation of natural gas. J. Nanostruct. Chem. 2015, 5, 219-226. 40. Kumar, A.; Bhattacharjee, G.; Kulkarni, B. D.; Kumar, R. Role of surfactants in promoting gas hydrate formation. Ind. Eng. Chem. Res. 2015, 54, 12217-12232. 41. Liu, Y.; Chen, B.; Chen, Y.; Zhang, S.; Guo, W.; Cai, Y.; Tan, B.; Wang, W. Methane storage in a hydrated form as promoted by Leucines for possible application to natural gas transportation and storage. Energy Technol- Ger. 2015, 3, 815-819. 42. Rogers, R.; Zhang, G.; Dearman, J.; Woods, C. Investigations into surfactant/gas hydrate relationship. J. Petrol. Sci. Eng. 2007, 56, 82-88. 43. Veluswamy, H. P.; Hong, Q. W.; Linga, P. Morphology study of methane hydrate formation and dissociation in the presence of amino acid. Cryst. Growth Des. 2016, 16, 5932-5945. 44. Wang, F.; Meng, H. L.; Guo, G.; Luo, S. J.; Guo, R. B. Methane hydrate formation promoted by −SO3–-coated graphene oxide nanosheets. ACS Sustain. Chem. Eng. 2017, 5, 6597-6604. 45. 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, 68-80. 46. Zhong, Y.; Rogers, R. E. Surfactant effects on gas hydrate formation. Chem. Eng. Sci. 2000, 55, 4175-4187. 47. Ji, H.; Chen, D.; Wu, G. Molecular mechanisms for cyclodextrin-promoted methane hydrate formation in water. J. Phys. Chem. C 2017, 121, 20967-20975. 48. 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, 17852-62. 49. Kelland, M. A. History of the development of low dosage hydrate inhibitors. Energ. Fuel 2006, 20, 825-847. 50. 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, 12079-85. 51. Yagasaki, T.; Matsumoto, M.; Tanaka, H. Adsorption of kinetic hydrate inhibitors on growing surfaces: A molecular dynamics study. J. Phys. Chem. B 2018, 122, 3396. 52. Kamal, M. S.; Hussein, I. A.; Sultan, A. S.; Solms, N. V. Application of various water soluble polymers in gas hydrate inhibition. Renew. Sust. Energ. Rev. 2016, 60, 206-225. 53. Abrahamsen, E.; Heyns, I. M.; Solms, N. V.; Pfukwa, R.; Klumperman, B.; Kelland, M. A. First study of poly(3-methylene-2-pyrrolidone) as a kinetic hydrate inhibitor. Energ. Fuel 2017, 31, 13572-13577. 54. Ke, W.; Kelland, M. A. Kinetic hydrate inhibitor studies for gas hydrate systems – A review of experimental equipment and test methods. Energ. Fuel 2016, 30, 10015-10028. 55. Bellucci, M. A.; Walsh, M. R.; Trout, B. L. Molecular dynamics analysis of anti-agglomerant surface adsorption in natural gas hydrates. J. Phys. Chem. C 2018, 125, 2673-2683. 56. Li, Z.; Jiang, F.; Qin, H.; Liu, B.; Sun, C.; Chen, G. Molecular dynamics method to simulate the process of hydrate growth in the presence/absence of KHIs. Chem. Eng. Sci. 2017, 164, 26

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30 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

The Journal of Physical Chemistry

307-312. 57. Kvamme, B.; Kuznetsova, T.; Aasoldsen, K. Molecular dynamics simulations for selection of kinetic hydrate inhibitors. J. Mol. Graph. Model. 2005, 23, 524. 58. Sun, T.; Davies, P. L.; Walker, V. K. Structural basis for the inhibition of gas hydrates by α-helical antifreeze proteins. Biophys. J. 2015, 109, 1698-705. 59. Xu, J.; Li, L.; Liu, J.; Wang, X.; Yan, Y.; Zhang, J. The molecular mechanism of the inhibition effects of PVCaps on the growth of sI hydrate: an unstable adsorption mechanism. Phys. Chem. Chem. Phys. 2018, 20, 8326. 60. Bagherzadeh, S. A.; Alavi, S.; Ripmeester, J. A.; Englezos, P. Why ice-binding type I antifreeze protein acts as a gas hydrate crystal inhibitor. Phys. Chem. Chem. Phys. 2015, 17, 9984. 61. 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, 150203143345006. 62. Cruz‐Torres, A.; Romero‐Martínez, A.; Galano, A. Computational study on the antifreeze glycoproteins as inhibitors of clathrate‐hydrate formation. Chemphyschem 2010, 9, 1630-1635. 63. Myran, D.; Middleton, A.; Choi, J.; Gordienko, R.; Ohno, H.; Ripmeester, J. A.; Walker, V. K. Genetically-engineered mutant antifreeze proteins provide insight into hydrate inhibition. 2011. 64. Sharifi, H.; Walker, V. K.; Ripmeester, J.; Englezos, P. Inhibition activity of antifreeze proteins with natural gas hydrates in saline and the light crude oil mimic, heptane. Energ. Fuel 2014, 28, 3712–3717. 65. Perfeldt, C. M.; Chua, P. C.; Daraboina, N.; Friis, D.; Kristiansen, E.; Ramløv, H.; Woodley, J. M.; Kelland, M. A.; Solms, N. V. Inhibition of gas hydrate nucleation and growth: Efficacy of an antifreeze protein from the longhorn beetle rhagium mordax. Energ. Fuel 2014, 28, 3666-3672. 66. Sharifi, H.; Walker, V. K.; Ripmeester, J.; Englezos, P. Insights into the behavior of biological clathrate hydrate inhibitors in aqueous saline solutions. Cryst. Growth Des 2014, 14, 2923–2930. 67. Ohno, H.; Susilo, R.; Gordienko, R.; Ripmeester, J.; Walker, V. K. Interaction of antifreeze proteins with hydrocarbon hydrates. Chem-Eur. J. 2010, 16, 10409-10417. 68. Gordienko, R.; Ohno, H.; Singh, V. K.; Jia, Z.; Ripmeester, J. A.; Walker, V. K. Towards a green hydrate inhibitor: imaging antifreeze proteins on clathrates. Plos One 2010, 5, e8953. 69. Walsh, M. R.; Beckham, G. T.; Koh, C. A.; Sloan, E. D.; Wu, D. T.; Sum, A. K. Methane hydrate nucleation rates from molecular dynamics simulations: effects of aqueous methane concentration, interfacial curvature, and system size. J. Phys. Chem. C 2011, 115, 21241-21248. 70. Kalogerakis, N.; Jamaluddin, A. K. M.; Dholabhai, P. D.; Bishnoi, P. R. Effect of surfactants on hydrate formation kinetics. Spe International Symposium on Oilfield Chemistry 1993. 71. Wu, M.; Wang, Y.; Wu, W.; Hu, C.; Wang, X.; Zheng, J.; Li, Z.; Jiang, B.; Qiu, J. Preparation of functionalized water-soluble photoluminescent carbon quantum dots from petroleum coke. Carbon 2014, 78, 480-489. 72. Allen, B. L.; Kichambare, P. D.; Star, A. Carbon nanotube field-effect-transistor-based biosensors . Adv. Mater. 2010, 19, 1439-1451. 73. Allen-King, R. M.; Grathwohl, P.; Ball, W. P. New modeling paradigms for the sorption of hydrophobic organic chemicals to heterogeneous carbonaceous matter in soils, sediments, and rocks. Adv. Water Resour. 2002, 25, 985-1016. 74. Cai, Y.; Jiang, G.; Liu, J.; Zhou, Q. Multiwalled carbon nanotubes as a solid-phase extraction 27

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

adsorbent for the determination of bisphenol A, 4-n-nonylphenol, and 4-tert-octylphenol. Anal. Chem. 2003, 75, 2517-21. 75. Consales, M.; Crescitelli, A.; Campopiano, S.; Cutolo, A.; Penza, M.; Aversa, P.; Giordano, M.; Cusano, A. Chemical detection in water by single-walled carbon nanotubes-based optical fiber sensors. IEEE. Sens. J. 2007, 7, 1004-1005. 76. Majumder, M.; Chopra, N.; Hinds, B. J. Effect of tip functionalization on transport through vertically oriented carbon nanotube membranes. J. Am. Chem. Soc. 2005, 127, 9062. 77. Mauter, M. S.; Elimelech, M. Environmental applications of carbon-based nanomaterials. Environ. Sci. Technol. 2008, 42, 5843. 78. Cox, S. J.; Taylor, D. J.; Youngs, T. G.; Soper, A. K.; Totton, T. S.; Chapman, R. G.; Arjmandi, M.; Hodges, M. G.; Skipper, N. T.; Michaelides, A. Formation of methane hydrate in the presence of natural and synthetic nanoparticles. J. Am. Chem. Soc. 2018, 140, 3277-3284. 79. Ghozatloo, A.; Hosseini, M.; Shariaty-Niassar, M. Improvement and enhancement of natural gas hydrate formation process by Hummers' graphene. J. Nat. Gas Sci. Eng. 2015, 27, 1229-1233. 80. Nashed, O.; Partoon, B.; Lal, B.; Sabil, K. M.; Shariff, A. M. Review the impact of nanoparticles on the thermodynamics and kinetics of gas hydrate formation. J. Nat. Gas Sci. Eng. 2018, 55, 452-465. 81. Paulista Neto, A. J.; Fileti, E. E. Impact of edge groups on the hydration and aggregation properties of graphene oxide. J. Phys. Chem. B 2018, 122, 2578. 82. Rezaei, E.; Manteghian, M.; Tamaddondar, M. Kinetic study of ethylene hydrate formation in presence of graphene oxide and sodium dodecyl sulfate. J. Petrol. Sci. Eng. 2016, 147, 857-863. 83. Park, S. S.; Lee, S. B.; Kim, N. J. Effect of multi-walled carbon nanotubes on methane hydrate formation. J. Ind. Eng. Chem. 2010, 16, 551-555. 84. Kim, D.; Kim, D. W.; Lim, H. K.; Jeon, J.; Kim, H.; Jung, H. T.; Lee, H. Inhibited phase behavior of gas hydrates in graphene oxide: influences of surface and geometric constraints. Phys. Chem. Chem. Phys. 2014, 16, 22717-22722. 85. Song, Y.; Wang, F.; Liu, G. Q.; Luo, S. J.; Guo, R. B. Promotion effect of carbon nanotubes-doped SDS on methane hydrate formation. Energ. Fuel 2017, 31, 1850-1857. 86. Zhou, S.; Jiang, K.; Zhao, Y.; Chi, Y.; Wang, S.; Zhang, G. Experimental investigation of CO2 hydrate formation in the water containing graphite nanoparticles and tetra-n-butyl ammonium bromide. J. Chem. Eng. Data 2018, 63, 389-394. 87. Zhou, S. D.; Yu, Y. S.; Zhao, M. M.; Wang, S. L.; Zhang, G. Z. Effect of graphite nanoparticles on promoting CO2 hydrate formation. Energ. Fuel 2014, 28, 4694-4698. 88. Geng, H.; Liu, X.; Shi, G.; Bai, G.; Ma, J.; Chen, J.; Wu, Z.; Song, Y.; Fang, H.; Wang, J. Graphene oxide restricts growth and recrystallization of ice crystals. Angew. Chem. Int. Edit. 2017, 56, 997-1001. 89. Jorgensen, W. L.; Madura, J. D.; Swenson, C. J. Optimized intermolecular potential functions for liquid hydrocarbons. J. Am. Chem. Soc. 1984, 106:22, 6638-6646. 90. Abascal, J. L. F.; Sanz, E.; Fernández, R. G.; Vega, C. A potential model for the study of ices and amorphous water: TIP4P/Ice. J. Chem. Phys. 2005, 122, 234511. 91. Tu, Y.; Lv, M.; Xiu, P.; Huynh, T.; Zhang, M.; Castelli, M.; Liu, Z.; Huang, Q.; Fan, C.; Fang, H. Destructive extraction of phospholipids from Escherichia coli membranes by graphene nanosheets. Nat. Nanotechnol. 2013, 8, 594. 92. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N⋅ log (N) method for Ewald 28

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30 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

The Journal of Physical Chemistry

sums in large systems. J. Chem. Phys. 1993, 98, 10089-10092. 93. Allen, M.; Tildesley, D. J. Computer simulation of liquids. Clarendon. Oxford: 1987. 94. Gunsteren, W. F. V.; Berendsen, H. J. C. A Leap-frog algorithm for stochastic dynamics. Mol. Simulat. 1988, 1, 173-185. 95. Zhang, Z.; Liu, C.-J.; Walsh, M. R.; Guo, G.-J. Effects of ensembles on methane hydrate nucleation kinetics. Phys. Chem. Chem. Phys. 2016, 18, 15602-15608. 96. Conde, M.; Vega, C. Determining the three-phase coexistence line in methane hydrates using computer simulations. J. Chem. Phys. 2010, 133, 064507. 97. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. C. GROMACS: fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701. 98. Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph. 1996, 14, 33-38. 99. Moon, C.; Hawtin, R. W.; Rodger, P. M. Nucleation and control of clathrate hydrates: insights from simulation. Faraday Discuss. 2007, 136, 367-382. 100. 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, 10298-10307. 101. Guo, G.-J.; Rodger, P. M. Solubility of aqueous methane under metastable conditions: implications for gas hydrate nucleation. J. Phys. Chem. B 2013, 117, 6498-6504. 102. Lauricella, M.; Ciccotti, G.; English, N. J.; Peters, B.; Meloni, S. Mechanisms and nucleation rate of methane hydrate by dynamical nonequilibrium molecular dynamics. J. Phys. Chem. C 2017, 121, 24223-24234. 103. 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. 104. Jiménez-Ángeles, F.; Firoozabadi, A. Nucleation of methane hydrates at moderate subcooling by molecular dynamics simulations. J. Phys. Chem. C 2014, 118, 11310-11318. 105. Liu, Z.; Pan, Z.; Zhang, Z.; Liu, P.; Shang, L.; Li, B. Effect of porous media and sodium dodecyl sulphate complex system on methane hydrate formation. Energ. Fuel 2018, 32, 5736-5749. 106. 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, 1095-1098. 107. Casco, M. E.; Silvestre-Albero, J.; Ramírez-Cuesta, A. J.; Rey, F.; Jordá, J. L.; Bansode, A.; Urakawa, A.; Peral, I.; Martínez-Escandell, M.; Kaneko, K. Methane hydrate formation in confined nanospace can surpass nature. Nat. Commun. 2015, 6, 6432. 108. Warrier, P.; Khan, M. N.; Srivastava, V.; Maupin, C. M.; Koh, C. A. Overview: Nucleation of clathrate hydrates. J. Chem. Phys. 2016, 145, 211705. 109. Park, S. S.; An, E. J.; Lee, S. B.; Chun, W. G.; Kim, N. J. Characteristics of methane hydrate formation in carbon nanofluids. J. Ind. Eng. Chem. 2012, 18, 443-448. 110. Prasher, R.; Phelan, P. E.; Bhattacharya, P. Effect of aggregation kinetics on the thermal conductivity of nanoscale colloidal solutions (nanofluid). Nano Lett. 2006, 6, 1529.

29

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

TOC Graphic

30

ACS Paragon Plus Environment

Page 30 of 30