Dissociation of Methane Hydrate in Aqueous NaCl Solutions

Sep 19, 2014 - down in an early stage before bubble formation, and another is swift bubble formation ... transportation, and gas separation purposes.3...
0 downloads 0 Views 2MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Dissociation of Methane Hydrate in Aqueous NaCl Solutions Takuma Yagasaki, Masakazu Matsumoto, Yoshimichi Andoh, Susumu Okazaki, and Hideki Tanaka J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/jp507978u • Publication Date (Web): 19 Sep 2014 Downloaded from http://pubs.acs.org on September 26, 2014

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.

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

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

Dissociation of Methane Hydrate in Aqueous NaCl Solutions

Takuma Yagasaki,† Masakazu Matsumoto,† Yoshimichi Andoh,§ Susumu Okazaki,§ and Hideki Tanaka†, ‡,* †

Department of Chemistry, Faculty of Science, Okayama University, Okayama,

700-8530, Japan ‡

Research Center of New Functional Materials for Energy Production, Storage and

Transport, Okayama, 700-8530, Japan §

*

Department of Applied Chemistry, Nagoya University, Nagoya 464-8603, Japan

Email: [email protected].

Phone and Fax: +81-86-251-7769

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 Molecular dynamics simulations of the dissociation of methane hydrate in aqueous NaCl solutions are performed.

It is shown that the dissociation of the hydrate

is accelerated by the formation of methane bubbles both in NaCl solutions and in pure water.

We find two significant effects on the kinetics of the hydrate dissociation by

NaCl.

One is slowing down in an early stage before bubble formation and another is

swift bubble formation that enhances the dissociation.

These effects arise from the low

solubility of methane in NaCl solution which gives rise to a non-uniform spatial distribution of solvated methane in the aqueous phase.

We also demonstrate that

bubbles form near the hydrate interface in dense NaCl solutions, and that the hydrate dissociation proceeds inhomogeneously due to the bubbles.

Keywords Molecular dynamics, salt effects, bubble formation.

2

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28

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 Clathrate hydrate is a crystalline inclusion compound in which small molecules, such as methane and hydrogen, are held within polyhedral water cages.1

It is known

that a huge amount of natural gas is stored as clathrate hydrates in ocean floor and permafrost regions, and these naturally occurring gas hydrates are expected as a future energy resource.2-4

It has also been proposed that gas hydrates can be used for energy

storage, transportation, and gas separation purposes.3,5-11 Molecular dynamics (MD) is a powerful tool to examine microscopic features of condensed phases, and has provided much insight into static and dynamic properties of gas hydrates,

such as

thermal stability,12-23 structure

of the

interface,24-26

dissociation,27-42 formation,43-67 and molecular diffusion in the cage structure.68-74

In a

previous paper, we performed large-scale MD simulations of the dissociation of a methane hydrate cluster in pure water.41

It was shown that the concentration of

methane in the liquid phase increases with time, and this results in a decrease in the dissociation rate due to reconstruction of hydrate cages. further, bubbles of methane form. phase.

As the dissociation proceeds

Once formed, they rapidly grow in the aqueous

The dissociation rate turns to increase after the bubble formation because the

bubbles absorb surrounding methane molecules supersaturated in the solution.

We

also demonstrated that methane hydrate can be preserved as a metastable superheated solid if there are no pre-existent bubbles near the hydrate/liquid interface. It is experimentally known that properties of gas hydrates are affected by the presence of solutes in the surrounding aqueous phase.

For example, the

hydrate-liquid-gas three-phase equilibrium temperature is lowered by adding alcohols

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

into the aqueous phase.1,3

Page 4 of 28

This is due to the lower chemical potential of water in the

aqueous solution of alcohol than in pure water.

To the best of our knowledge, however,

there are only a limited number of MD simulations to examine roles of solutes other than the guest species,22,64 and thus the effects of solutes on the dissociation kinetics are still unclear at the molecular scale. In this paper, we perform MD simulations of the dissociation of methane hydrate in

aqueous

NaCl solutions.

Experimental

studies have

shown

that the

hydrate-liquid-gas three-phase equilibrium temperature, Teq, in aqueous NaCl solutions becomes lower with increasing NaCl concentrations.75,76

This fact suggests that the

dissociation of methane hydrate is also faster in solutions of a higher NaCl concentration at a given temperature.

This study demonstrates, however, that the

effects of NaCl are not so simple as what we anticipated.

We find that NaCl has both

deceleration and acceleration effects on the kinetics of hydrate dissociation.

These

effects are explained by the spatial distribution of methane molecules released from the hydrate.

We also discuss the temperature dependence of the NaCl effects.

Methods The initial structure consists of a hydrate cluster and the surrounding aqueous phase.

The hydrate cluster is a 9 × 9 × 9 unit cell replica of fully occupied structure I

methane hydrate. 0.6 mol kg-1.

Simulations are performed for two NaCl concentrations, m = 4.8 and

The concentration of m = 4.8 mol kg-1 is somewhat lower than the

experimental saturation condition at ambient temperature, m ~ 6 mol kg-1, and m = 0.6 mol kg-1 is close to the NaCl concentration of seawater. 4

ACS Paragon Plus Environment

The analyses are also

Page 5 of 28

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

performed for the m = 0.0 mol kg-1 solution, i.e., pure water.

The trajectories of the

pure water system are the same as those examined in our previous paper.41

The

number of each chemical species in the initial configuration is summarized in Table I. The sum of molecules and ions in the system is 124510 for all three concentrations. This large system size is required to treat methane bubbles that grow quite rapidly.

Table I.

The number of each chemical species in the initial configuration. Hydrate

Liquid phase

m / mol kg-1

CH4

H2O

Na+

Cl-

H2 O

0.0

5832

33534

0

0

85144

0.6

5832

33534

900

900

83344

4.8

5832

33534

6273

6273

72598

In accordance with the previous study,41 we employ the TIP4P/2005 model for water,77 and the OPLS united model for methane molecules.78 parameters for Na+ and Cl- are taken from ref. 79.

The MD simulations are carried out

using the MODYLAS package with a time step of 2 fs.80 MPa in all simulations. directions.

The potential

The pressure is kept at 0.1

Periodic boundary conditions are applied in all three

Long-range Coulomb interactions are treated with the fast multipole

method.80 The equilibration procedure is the same as that of our previous paper.41

First,

the aqueous phase is equilibrated at 300 K with fixing all the degrees of freedom of the hydrate cluster.

The simulation is continued for 2 ns at 220 K without any constrains

to relax the structure of the hydrate cluster.

Then, the temperature is suddenly changed 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

to a higher temperature.

Page 6 of 28

We set this instant as the time origin, t = 0, and analyze the

following dissociation process.

Methane molecules are classified into three types,

methane in the hydrate cluster, solvated in the aqueous phase, and in the gas phase at each time step based on a procedure given in the previous paper.41

The F3 parameter is

used for the classification.27

Results and Discussion In Figure 1, we present the number of methane molecules in the hydrate cluster, Ng, at T = 292 K.

The hydrate-liquid-gas three-phase equilibrium temperature, Teq, in

aqueous NaCl solutions is lower for higher NaCl concentrations.75,76

If the

dissociation rate is simply proportional to the degree of superheating, T – Teq, Ng would be smaller for solutions of higher concentrations at any moment. demonstrates that this simple expectation is not true.

Figure 1 clearly

It is seen that Ng is larger for 10

ns < t < 90 ns for the 0.6 mol kg-1 solution than for the 0.0 mol kg-1 solution. indicates that the presence of NaCl slows down the hydrate dissociation.

This

The slowing

down is also seen in the 4.8 mol kg-1 solution.

Figure 1. T = 292 K.

Time evolution of the number of methane molecules in the hydrate cluster at Black, green, and red curves are the results of the 0.0, 0.6, and 4.8 mol kg-1 6

ACS Paragon Plus Environment

Page 7 of 28

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

solutions, respectively.

It is well known that the formation of methane hydrate is enhanced by pressure. This is because the pressurization increases the possibility of finding methane molecules near the hydrate surface in water.

In our simulations, the dissociation of the hydrate

cluster gradually increases the number of methane molecules at the surface.

This

results in a similar effect to the pressurization, i.e., the enhancement of the construction of hydrate cages.

Since the temperature is higher than Teq, this is seen as the decrease

in the dissociation rate.

The slowing down of the dissociation in the 0.0 mol kg-1

solution for t < 80 ns shown in Figure 1 arises from this mechanism.

The dissociation

rate turns to increase at t ~ 80 ns because two bubbles form at t = 73 and 84 ns. similar behavior is observed in the NaCl solutions.

A

An increase in the dissociation rate

is seen at ~70 ns in the 0.6 mol kg-1 solution and at ~20 ns in the 4.8 mol kg-1 solution. This result suggests that bubble formation occurs earlier in denser NaCl solutions. NaCl has two effects on the kinetics of the hydrate dissociation: slowing down at an early stage before bubble formation and a rapid bubble formation that enhances the dissociation.

Hereafter, we focus on the difference between the dissociation processes

in the 4.8 and 0.0 mol kg-1 solutions because both effects are more significant in the denser NaCl solution.

Figure 2 presents snapshots of the dissociation process in the

0.0 mol kg-1 solution at T = 292 K. initial condition.

The hydrate cluster is cubic at t = 0 owing to the

The hydrate cluster changes its shape from a cube to a sphere (Figure

2b) because acute parts of the hydrate cluster dissociate faster than planar parts due to the Gibbs-Thomson effect81 (note that the contribution of the Gibbs-Thomson effect to 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

Page 8 of 28

the dissociation rate of methane hydrate was examined in detail in our previous paper41). The spherical shape of the cluster is maintained during the dissociation for t > ~15 ns. The concentration of methane in the aqueous phase increases as the dissociation proceeds until the formation of bubbles at t = 73 ns and 84 ns.

The number of solvated

methane molecules in the aqueous phase begins to decrease after the formation of bubbles because the bubbles absorb the surrounding methane molecules.

Figure 2.

Several snapshots along the hydrate dissociation process in the 0.0 mol kg-1

solution at T = 292 K.

Gray, blue, and red particles represent methane molecules in the

hydrate cluster, solvated in the aqueous phase, and in bubbles, respectively.

Water

molecules are not shown.

Figure 3 shows snapshots of the hydrate dissociation in the 4.8 mol kg-1 solution at T = 292 K.

There are several significant differences in dissociation behavior between

the 0.0 and 4.8 mol kg-1 solutions.

One is the location of bubbles: All bubbles form in 8

ACS Paragon Plus Environment

Page 9 of 28

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 vicinity of the hydrate-liquid interface in the dense NaCl solution. elapsed time before bubble formation.

Another is the

As expected from Figure 1, the bubbles form

much earlier in the 4.8 mol kg-1 solution than in the 0.0 mol kg-1 solution. the shape of the hydrate cluster.

The other is

As shown in Figures 3d and 3e, the hydrate cluster is

nonspherical in the dense NaCl solution.

Figure 3.

Snapshots of the hydrate dissociation in 4.8 mol kg-1 solution at T = 292 K.

Gray, blue, and red particles represent methane molecules in the hydrate cluster, solvated in the aqueous phase, and in bubbles, respectively. Water molecules and ions are not shown.

To analyze the spatial distribution of each chemical species, we consider a cylindrical coordinate system whose origin is set to the position of G0, where G0 is the guest molecule that is the closest to the center of mass (COM) of the hydrate cluster at t = 0.

The number density profile along the cylinder axis is calculated with a grid size 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

of 0.5 Å.

Page 10 of 28

The number of molecules located within 15 Å from the axis in each grid is

counted to construct the density profile (i.e., we consider a cylinder of a radius of 15 Å). Figures 4a-4d present the number density profiles in a cylinder in the 0.0 mol kg-1 solution.

We refer to this cylinder as Cylinder-0.

As shown in Figure 5a, the

direction of Cylinder-0 is so chosen as to include a given bubble. Before the formation of the bubble, the cylinder axis passes through a point where the bubble emerges at t = 73 ns.

The cylinder axis passes through the COM of the bubble when t > 73 ns.

Figure 4a shows that the hydrate-liquid interface resides in c = 80 Å, i.e., 80 Å away from G0, at t = 1 ns.

Since the initial structure is prepared without caring for open

cages at the interface, a certain amount of guest molecules is quickly released from the hydrate.

This gives rise to a peak in the distribution of the solvated methane near the

interface.

As the dissociation proceeds, the interface shifts inward and the number of

solvated methane molecules increases.

Figure 4c shows that the density profile of

solvated methane is almost flat everywhere in the liquid phase.

This indicates that the

timescale of the diffusion of methane in the aqueous phase is faster than that of the release of methane from the hydrate.

The concentration of methane in the aqueous

phase increases with time, and a bubble forms at t = 73 ns.

Because the distribution of

the solvated methane is uniform, there is no particular place where bubbles would form more easily.

In this case, the bubble forms at c ~ 75 Å, about 40 Å away from the

hydrate/liquid interface, due to the concentration fluctuation. supersaturation is roughly 0.002 Å-3.

The limit of

Note that similar results are obtained for another

bubble in the 0.0 mol kg-1 solution because the dissociation proceeds isotropically.

10

ACS Paragon Plus Environment

Page 11 of 28

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

The number density profiles along the cylinder axis.

Black, red, and green

curves are the number density profiles of methane molecules in the hydrate cluster, aqueous phase, and bubbles, respectively.

Blue curves show the number density of the

Na+ ion. Curves for Cl- are not shown because they are quite similar to the profiles of Na+.

A moving average is taken with a window of 1 ns to reduce noisy fluctuations.

Figure 5.

(a) Cylinder-0 at t = 80 ns.

The blue arrow is the cylinder axis that passes 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

through the COM of a bubble.

Page 12 of 28

(b) Cylinder-4.8a (green) and Cylinder-4.8b (blue) at t

= 20 ns.

The number density profiles in a cylinder in the 4.8 mol kg-1 solution are shown in Figures 4e-4g.

This cylinder, referred to as Cylinder-4.8a, does not include bubbles

(Figure 5b), and can be straightforwardly compared with Cylinder-0 in which no bubble exists at t = 1, 7, and 20 ns.

Note that the results do not depend on the angle of the

cylinder when it does not include any bubble.

At t = 1 ns, the density profiles of

methane molecules in the 4.8 mol kg-1 solution (Figures 4e) are similar to those in the 0.0 mol kg-1 system (Figure 4a).

Water and methane molecules are released into the

aqueous phase as the dissociation proceeds. hydrate interface (Figures 4f and 4g).

This results in an ion-poor region near the

In pure water, as mentioned above, the timescale

of the diffusion of methane in the aqueous phase is faster than that of the release of methane from the hydrate.

The distribution of methane in solution therefore becomes

flat at t = 20 ns (Figure 4c).

In contrast, the solvated methane molecules in the NaCl

solution prefer to stay near the interface, and do not enter into the ion-rich region as shown in Figure 4g.

This non-uniform distribution is caused by the low solubility of

methane in aqueous NaCl solutions. The lower solubility cause by increasing NaCl concentration is examined for the present model interactions by a simple particle insertion method.82

MD simulations

are performed for aqueous solutions containing 1,000 molecules (ions) using the GROMACS package.83,84 The excess chemical potentials are calculated by trial insertion of methane for the configurations generated by MD simulations. 12

ACS Paragon Plus Environment

As shown

Page 13 of 28

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

in Figure 6, the calculated excess chemical potential is higher for solutions of higher NaCl concentrations.

This result is consistent with experimental observations.85,86

Salt changes the structure of water in a similar way as pressurization.87

The increase in

NaCl concentration reduces the number of large voids in the solution which can accommodate methane molecules.

As a result, the excess chemical potential increases

as the NaCl concentration increases.

Figure 6.

Excess chemical potential of methane in aqueous NaCl solutions at 300 K.

The presence of solvated methane near the hydrate/liquid interface facilitates the reconstruction of hydrate cages.41

This fact, together with the comparison between

Figures 4c and 4g, suggests that the hydrate dissociation rate is slower in Cylinder-4.8a than in Cylinder-0.

In Figure 7, we show the position of the hydrate/liquid interface

against time for these cylinders.

The position of the interface is defined as the most

distant grid from G0 in which the number density of the methane molecules in the hydrate is larger than 0.005 Å-3.

The dissociation rate is indeed slower in

Cylinder-4.8a than in Cylinder-0 for t > 10 ns.

The dissociation rates are the almost

same for t < 10 ns because both cylinders include an edge of the cubic hydrate cluster and the dissociation kinetics is dominated by the Gibbs-Thomson effect in this period. 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

Page 14 of 28

The present result demonstrates that the slowing down effect caused by NaCl arises from the non-uniform distribution of the solvated methane.

Figure 7.

Position of the hydrate interface against time.

Data for the 4.8 mol kg-1

solution is truncated at t = 41 ns because G0 dissolves into the aqueous phase at this time.

A bubble forms when the local concentration of solvated methane exceeds the limit of supersaturation.

Because the solvated methane molecules can stay only in a

small area near the interface, the bubble forms at an early stage of the dissociation in the dense NaCl solution.

In Figures 4h-4j, we present the number density profile of a

cylinder in the 4.8 mol kg-1 solution that includes a bubble. Cylinder-4.8b and shown in Figure 5b.

This cylinder is named

Similar results are obtained for other bubbles.

The local concentration exceeds the limit of supersaturation, 0.002 Å-3 as early as t = 7 ns and a bubble forms near the hydrate interface. bubbles form at t = 7, 8, 13, and 21 ns.

In the 4.8 mol kg-1 solution, four

Acceleration of the dissociation found at t ~ 20

ns in Figure 1 is caused by these bubbles.

Note that bubble formation is a stochastic

phenomenon, and thus it is not surprising that no bubble occurs in Cylinder-4.8a even though the local concentration exceeds the limit of supersaturation. 14

ACS Paragon Plus Environment

It should also be

Page 15 of 28

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

noted that the probability of bubble formation would be the same anywhere on the hydrate/liquid interface if the initial hydrate cluster is spherical.

In this study, however,

the hydrate cluster is initially cubic, and each bubble forms near a vertex of the cube which dissociates faster than planar parts.

In macroscopic systems, there must be

asperity on the hydrate interface where the dissociation proceeds faster than other regions.

The present simulation suggests that methane bubbles tend to form near such

places in dense NaCl solutions. The acceleration of the hydrate dissociation due to bubbles is explained from the chemical potential of methane in the hydrate (µH), liquid (µL), and gas phases (µG). Before bubble formation, the driving force for the dissociation can be given as F0 = C(µH – µL) + Fw, where C is a constant and Fw is the driving force arising from the difference between the chemical potentials of water in the hydrate and liquid phases. After bubble formation, the driving force can reach F’ = C(µH – µG) + Fw.

Because of

the nature of the metastable state, the chemical potential of methane in water is higher than that in the gas phase, µL > µG, and hence F’ > F0. is accelerated after the bubble formation.

Thus, the hydrate dissociation

The driving force is maximized when a

bubble is in contact with the hydrate interface.

This is evidently seen in Figure 8.

The configurations presented in Figure 8 are the same as those shown in Figure 3d and 3e but are drawn in a different way to show the shape of hydrate cluster more clearly. Figure 8 demonstrates that the cluster interface near each bubble is concave inward in the 4.8 mol kg-1 solution.

Figure 7 also shows that the dissociation is much faster in

Cylinder-4.8b than in Cylinder-4.8a.

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

Figure 8.

Page 16 of 28

Shape of the hydrate cluster in the 4.8 mol kg-1 solution at (a) t = 30 ns and

(b) t = 40 ns.

Large particles are the methane molecules in the hydrate cluster that are

colored based on the distance from the COM of the cluster.

Green methane molecules

indicate the sphere surface of the equivalent volume with the cluster, and red and blue ones indicate the 20% elevated and concaved surface regions, respectively.

Small

white particles are methane molecules in bubbles.

The condition of bubble formation depends on the system size and the ratio of the number of methane molecules to that of water, R.88,89

As shown above, bubbles form

in the bulk region of the aqueous phase in the 0.0 mol kg-1 solution.

If the simulation

is performed with much smaller R, i.e., larger aqueous region, the dissociation kinetics of this system would change drastically because the system cannot reach the limit of supersaturation and bubbles never form.

In contrast, the dissociation behavior of the

4.8 mol kg-1 solution is expected to be insensitive to the ratio R because the bubbles form in the vicinity of the hydrate interface. still a size effect.

Even when R is kept constant, there is

It is known that the free energy barrier of bubble formation

decreases with increasing system size.88

In the 0.0 mol kg-1 solution, the first bubble

forms at t = 73 ns, at which about 70 % of methane molecules dissolved into the

16

ACS Paragon Plus Environment

Page 17 of 28

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

aqueous phase.

If both the aqueous part and the initial hydrate crystal of this system

are equally enlarged, the bubble formation would occur on the earlier stage in the dissociation process at which less than 70 % of methane molecules dissolved into water (note that the elapsed time before bubble formation might be longer than 73 ns because it takes longer time for methane molecule in larger systems to diffuse throughout the system).

This size effect is also expected to be small in the 4.8 mol kg-1 system for the

same reason. This study assumes that the hydrate phase is well distant from the gas phase, which is also a stable phase at the present thermodynamic condition.

This assumption

would be relevant to real bulk systems, in which bubble formation dominates the dissociation kinetics as shown above.

However, there is another possible condition

where the gas phase initially exists near the hydrate interface.37

This could be

observed, for example, in dissociation of gas hydrates in porous media, in which the dissociation rate is determined by the timescale of transfer of methane molecules from the vicinity of the hydrate to the gas phase.

In between these two extreme cases, both

the bubble formation and the mass transfer would affect the dissociation kinetics. We finally discuss the temperature dependence of the dissociation rate.

Figure 9

presents the number of methane molecules in the hydrate at T = 296, 308, and 324 K. It is shown that the effects of NaCl decrease with increasing temperature.

The rate of

cage decomposition grows exponentially with temperature because it follows the Arrhenius law with an activation energy that is related to the mechanical stability of the hydrate.41

Other effects on the dissociation rate become negligible at higher

temperatures.

Thus, the dissociation kinetics does not depend on the NaCl

concentration at the highest temperature, T = 324 K.

At T = 296K, only the slowing

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

down effect is observed in the 0.6 mol kg-1 solution.

In contrast, the deceleration

effect is small while the acceleration effect is evident in the 4.8 mol kg-1 solution at T = 308 K.

There is no simple explanation why only one of the two effects of NaCl is

dominant under these conditions.

The kinetics of hydrate dissociation involves many

factors, including the activation energy of cage breaking, the rate of cage reconstruction, free energy barrier for bubble formation, and diffusivity of methane in the aqueous phase.

A model accounting for these factors, as well as the dependence on temperature

and the NaCl concentration, and knowledge as to how much they contribute to the non-equilibrium dissociation process is required for a quantitative understanding of the effects of NaCl.

18

ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28

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

Time evolution of the number of methane molecules in the hydrate cluster at

T = 296, 308, and 324 K.

Black, green, and red curves are the results of the 0.0, 0.6,

and 4.8 mol kg-1 solutions, respectively.

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

Page 20 of 28

Conclusions We have investigated the dissociation of methane hydrate in aqueous NaCl solutions using MD simulations.

The dissociation rate decreases with time, and then it

turns to increase by the formation of the methane bubbles both in pure water and in NaCl solutions.

It is found that two effects of NaCl on the dissociation kinetics.

One

is the slowing down which is significant before bubble formation, and another is the acceleration due to rapid bubble formation.

These effects arise from the low solubility

of methane in NaCl solutions which results in a locally high concentration of solvated methane near the hydrate interface. We also find that the dissociation occurs rather inhomogeneously in dense NaCl solutions because the methane bubbles, which enhance the dissociation, form in proximity to the interface. Other salts, such as CaCl2 and K2SO4, also decrease the solubility of methane.85 These salts would have similar effects on the kinetics of hydrate dissociation. interesting to consider the effects of other types of solute.

It is

Amphiphilic solutes may

exhibit a different mechanism because they can stabilize bubbles in aqueous solutions. The presence of solutes should affect the formation mechanism as well as the dissociation rate.

Controlling formation/dissociation kinetics using solutes might

become a useful technique in the industrial use of gas hydrates in the future.

Acknowledgments The present work was supported by a Grant-in-Aid by JSPS and by HPCI Strategic Programs for Innovative Research (SPIRE) and Computational Materials Science 20

ACS Paragon Plus Environment

Page 21 of 28

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

Initiative (CMSI), MEXT, Japan.

Calculations were performed on the K computer at

the RIKEN Advanced Institute for Computational Science (Project ID: hp140216).

References (1) Sloan, E. D.; Koh, C. A.: Clathrate Hydrates of Natural Gases; CPC Press: Boca Raton, 2008. (2) Kvenvolden, K. A. Gas hydrate and humans. Ann. N.Y. Acad. Sci. 2000, 2000 912, 17-22. (3) Sloan, E. D. Fundamental principles and applications of natural gas hydrates. Nature 2003, 2003 426, 353-359. (4) Boswell, R. Resource potential of methane hydrate coming into focus. J. Pet. Sci. Eng. 2007, 2007 56, 9-13. (5) Struzhkin, V. V.; Militzer, B.; Mao, W. L.; Mao, H. K.; Hemley, R. J. Hydrogen storage in molecular clathrates. Chem. Rev. 2007, 2007 107, 4133-4151. (6) Mao, W. L.; Mao, H. K. Hydrogen storage in molecular compounds. Proc. Natl. Acad. Sci.

U.S.A. 2004, 2004 101, 708-710. (7) 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, 2004 306, 469-471. (8) Lee, H.; Lee, J. W.; Kim, D. Y.; Park, J.; Seo, Y. T.; Zeng, H.; Moudrakovski, I. L.; Ratcliffe, C. I.; Ripmeester, J. A. Tuning clathrate hydrates for hydrogen storage. Nature 2005, 2005 434, 743-746. (9) Linga, P.; Kumar, R.; Englezos, P. The clathrate hydrate process for post and pre-combustion capture of carbon dioxide. J. Hazard. Mater. 2007, 2007 149, 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, 2000 34, 4397-4400. (11) Seo, Y.-T.; Moudrakovski, I. L.; Ripmeester, J. A.; Lee, J.-w.; Lee, H. Efficient Recovery of CO2 from Flue Gas by Clathrate Hydrate Formation in Porous Silica Gels. Environ. Sci.

Technol. 2005, 2005 39, 2315-2319. (12) Tanaka, H.; Kiyohara, K. On the Thermodynamic Stability of Clathrate Hydrate .1. J.

Chem. Phys. 1993, 1993 98, 4098-4109. 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

(13) Tanaka, H.; Kiyohara, K. The Thermodynamic Stability of Clathrate Hydrate .2. Simultaneous Occupation of Larger and Smaller Cages. J. Chem. Phys. 1993, 1993 98, 8110-8118. (14) Koyama, Y.; Tanaka, H.; Koga, K. On the thermodynamic stability and structural transition of clathrate hydrates. J. Chem. Phys. 2005, 2005 122, 074503. (15) Katsumasa, K.; Koga, K.; Tanaka, H. On the thermodynamic stability of hydrogen clathrate hydrates. J. Chem. Phys. 2007, 2007 127, 044509. (16) Nakayama, T.; Koga, K.; Tanaka, H. Augmented stability of hydrogen clathrate 2009 131, 214506. hydrates by weakly polar molecules. J. Chem. Phys. 2009, (17) Matsumoto, M.; Tanaka, H. On the Structure Selectivity of Clathrate Hydrates. J. Phys.

Chem. B 2011, 2011 115, 8257-8265. (18) Alavi, S.; Ripmeester, J. A.; Klug, D. D. Molecular-dynamics study of structure II hydrogen clathrates. J. Chem. Phys. 2005, 2005 123, 024507. (19) Alavi, S.; Ripmeester, J. A.; Klug, D. D. Molecular-dynamics simulations of binary structure II hydrogen and tetrahydrofurane clathrates. J. Chem. Phys. 2006, 2006 124, 014704. (20) Conde, M. M.; Vega, C.; Tribello, G. A.; Slater, B. The phase diagram of water at negative pressures: Virtual ices. J. Chem. Phys. 2009, 2009 131, 034510. (21) Jensen, L.; Thomsen, K.; von Solms, N.; Wierzchowski, S.; Walsh, M. R.; Koh, C. A.; Sloan, E. D.; Wu, D. T.; Sum, A. K. Calculation of Liquid Water−Hydrate−Methane Vapor Phase Equilibria from Molecular Simulations. J. Phys. Chem. B 2010, 2010 114, 5775-5782. (22) Qi, Y. X.; Wu, W. D.; Liu, Y. F.; Xie, Y. M.; Chen, X. The influence of NaCl ions on hydrate structure and thermodynamic equilibrium conditions of gas hydrates. Fluid Phase

Equilib. 2012, 2012 325, 6-10. (23) 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, 2005

127, 17852-17862. (24) Rodger, P. M.; Forester, T. R.; Smith, W. Simulations of the methane hydrate/methane gas interface near hydrate forming conditions conditions. Fluid Phase Equilib. 1996, 1996 116, 326-332. (25) Chihaia, V.; Adams, S.; Kuhs, W. F. Molecular dynamics simulations of properties of a (0 0 1) methane clathrate hydrate surface. Chem. Phys. 2005, 2005 317, 208-225. (26) Mastny, E. A.; Miller, C. A.; de Pablo, J. J. The effect of the water/methane interface on methane hydrate cages: The potential of mean force and cage lifetimes. J. Chem. Phys. 2008, 2008

129, 034701. (27) Baez, L. A.; Clancy, P. Computer-Simulation of the Crystal-Growth and Dissolution of Natural-Gas Hydrates. Ann. N.Y. Acad. Sci. 1994, 1994 715, 177-186. (28) Yasuoka, K.; Murakoshi, S. Molecular dynamics simulation of dissociation process for

22

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28

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

methane hydrate. Ann. N.Y. Acad. Sci. 2000, 2000 912, 678-684. (29) Tse, J. S.; Klug, D. D. Formation and decomposition mechanisms for clathrate hydrates.

J. Supramol. Chem. 2002, 2002 2, 467-472. (30) English, N. J.; Johnson, J. K.; Taylor, C. E. Molecular-dynamics simulations of methane hydrate dissociation. J. Chem. Phys. 2005, 2005 123, 244503. (31) English, N. J.; Phelan, G. M. Molecular dynamics study of thermal-driven methane hydrate dissociation. J. Chem. Phys. 2009, 2009 131, 074704. (32) Myshakin, E. M.; Jiang, H.; Warzinski, R. P.; Jordan, K. D. Molecular Dynamics Simulations of Methane Hydrate Decomposition. J. Phys. Chem. A 2009, 2009 113, 1913-1921. (33) Conde, M. M.; Vega, C. Determining the three-phase coexistence line in methane hydrates using computer simulations. J. Chem. Phys. 2010, 2010 133, 064507. (34) Alavi, S.; Ripmeester, J. A. Nonequilibrium adiabatic molecular dynamics simulations of methane clathrate hydrate decomposition. J. Chem. Phys. 2010, 2010 132, 144703. (35) 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, 2012

116, 3188-3197. (36) Bagherzadeh, S. A.; Englezos, P.; Alavi, S.; Ripmeester, J. A. Molecular simulation of non-equilibrium methane hydrate decomposition process. J. Chem. Thermodyn. 2012, 2012 44, 13-19. (37) Bagherzadeh, S. A.; Alavi, S.; Ripmeester, J. A.; Egnglezos, P. Evolution of methane during gas hydrate dissociation. Fluid Phase Equilib. 2013, 2013 358, 114-120. (38) Sarupria, S.; Debenedetti, P. G. Molecular Dynamics Study of Carbon Dioxide Hydrate Dissociation. J. Phys. Chem. A 2011, 2011 115, 6102-6111. (39) Smirnov, G. S.; Stegailov, V. V. Melting and superheating of sI methane hydrate: Molecular dynamics study. J. Chem. Phys. 2012, 2012 136, 044523. (40) Liu, Y.; Zhao, J. J.; Xu, J. C. Dissociation mechanism of carbon dioxide hydrate by molecular dynamic simulation and ab initio calculation. Comput. Theor. Chem. 2012, 2012 991, 165-173. (41) Yagasaki, T.; Matsumoto, M.; Andoh, Y.; Okazaki, S.; Tanaka, H. Effect of Bubble Formation on the Dissociation of Methane Hydrate in Water: A Molecular Dynamics Study.

J. Phys. Chem. B 2014, 2014 118, 1900-1906. (42) Uddin, M.; Coombe, D. Kinetics of CH and CO Hydrate Dissociation and Gas Bubble Evolution via MD Simulation. J. Phys. Chem. A 2014, 2014 118, 1971-1988. (43) Radhakrishnan, R.; Trout, B. L. A new approach for studying nucleation phenomena using molecular simulations: Application to CO2 hydrate clathrates. J. Chem. Phys. 2002, 2002

117, 1786-1796. 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

(44) Moon, C.; Taylor, P. C.; Rodger, P. M. Molecular Dynamics Study of Gas Hydrate Formation. J. Am. Chem. Soc. 2003, 2003 125, 4706-4707. (45) Hawtin, R. W.; Quigley, D.; Rodger, P. M. Gas hydrate nucleation and cage formation at a water/methane interface. Phys. Chem. Chem. Phys. 2008, 2008 10, 4853-4864. (46) English, N. J.; MacElroy, J. M. D. Theoretical studies of the kinetics of methane hydrate crystallization in external electromagnetic fields. J. Chem. Phys. 2004, 2004 120, 10247-10256. (47) English, N. J.; Lauricella, M.; Meloni, S. Massively parallel molecular dynamics simulation of formation of clathrate-hydrate precursors at planar water-methane interfaces: Insights into heterogeneous nucleation. J. Chem. Phys. 2014, 2014 140, 204714. (48) Vatamanu, J.; Kusalik, P. G. Unusual Crystalline and Polycrystalline Structures in Methane Hydrates. J. Am. Chem. Soc. 2006, 2006 128, 15588-15589. (49) Vatamanu, J.; Kusalik, P. G. Molecular insights into the heterogeneous crystal growth of sI methane hydrate. J. Phys. Chem. B 2006, 2006 110, 15896-15904. (50) Vatamanu, J.; Kusalik, P. G. Heterogeneous Crystal Growth of Methane Hydrate on Its sII [001] Crystallographic Face. J. Phys. Chem. B 2008, 2008 112, 2399-2404. (51) Vatamanu, J.; Kusalik, P. G. Observation of two-step nucleation in methane hydrates.

Phys. Chem. Chem. Phys. 2010, 2010 12, 15065-15072. (52) Liang, S.; Kusalik, P. G. Exploring nucleation of H2S hydrates. Chemical Science 2011, 2011

2, 1286. (53) Pirzadeh, P.; Kusalik, P. G. Molecular insights into clathrate hydrate nucleation at an ice-solution interface. J. Am. Chem. Soc. 2013, 2013 135, 7278-7287. (54) Nada, H. Growth mechanism of a gas clathrate hydrate from a dilute aqueous gas solution: A molecular dynamics simulation of a three-phase system. J. Phys. Chem. B 2006, 2006

110, 16526-16534. (55) Nada, H. Anisotropy in Growth Kinetics of Tetrahydrofuran Clathrate Hydrate: A Molecular Dynamics Study. J. Phys. Chem. B 2009, 2009 113, 4790-4798. (56) 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, 2011

115, 21241-21248. (57) 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, 2009 326, 1095-1098. (58) Jacobson, L. C.; Hujo, W.; Molinero, V. Amorphous Precursors in the Nucleation of Clathrate Hydrates. J. Am. Chem. Soc. 2010, 2010 132, 11806-11811. (59) Jacobson, L. C.; Hujo, W.; Molinero, V. Nucleation Pathways of Clathrate Hydrates: Effect of Guest Size and Solubility. J. Phys. Chem. B 2010, 2010 114, 13796-13807.

24

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28

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

(60) Jacobson, L. C.; Matsumoto, M.; Molinero, V. Order parameters for the multistep crystallization of clathrate hydrates. J. Chem. Phys. 2011, 2011 135, 074501. (61) Nguyen, A. H.; Jacobson, L. C.; Molinero, V. Structure of the Clathrate/Solution Interface and Mechanism of Cross-Nucleation of Clathrate Hydrates. J. Phys. Chem. C 2012, 2012

116, 19828-19838. (62) Knott, B. C.; Molinero, V.; Doherty, M. F.; Peters, B. Homogeneous nucleation of methane hydrates: unrealistic under realistic conditions. J. Am. Chem. Soc. 2012, 2012 134, 19544-19547. (63) Tung, Y. T.; Chen, L. J.; Chen, Y. P.; Lin, S. T. The Growth of Structure I Methane Hydrate from Molecular Dynamics Simulations. J. Phys. Chem. B 2010, 2010 114, 10804-10813. (64) Tung, Y.-T.; Chen, L.-J.; Chen, Y.-P.; Lin, S.-T. Molecular Dynamics Study on the Growth of Structure I Methane Hydrate in Aqueous Solution of Sodium Chloride. J. Phys.

Chem. B 2012, 2012 116, 14115-14125. (65) Sarupria, S.; Debenedetti, P. G. Homogeneous Nucleation of Methane Hydrate in Microsecond Molecular Dynamics Simulations. J. Phys. Chem. Lett. 2012, 2012 3, 2942-2947. (66) Jiménez-Ángeles, F.; Firoozabadi, A. Nucleation of Methane Hydrates at Moderate Subcooling by Molecular Dynamics Simulations. J. Phys. Chem. C 2014, 2014 118, 11310-11318. (67) Bai, J.; Zeng, X. C. Polymorphism and polyamorphism in bilayer water confined to slit nanopore under high pressure. Proceedings of the National Academy of Sciences 2012, 2012 109, 21240-21245. (68) Alavi, S.; Ripmeester, J. A. Hydrogen-Gas Migration through Clathrate Hydrate Cages.

Angew. Chem. 2007, 2007 119, 6214-6217. (69) Frankcombe, T. J.; Kroes, G.-J. Molecular Dynamics Simulations of Type-sII Hydrogen Clathrate Hydrate Close to Equilibrium Conditions. J. Phys. Chem. C 2007, 2007 111, 13044-13052. (70) Peters, B.; Zimmermann, N. E. R.; Beckham, G. T.; Tester, J. W.; Trout, B. L. Path Sampling Calculation of Methane Diffusivity in Natural Gas Hydrates from a Water-Vacancy Assisted Mechanism. J. Am. Chem. Soc. 2008, 2008 130, 17342-17350. (71) Liang, S.; Kusalik, P. G. The Mobility of Water Molecules through Gas Hydrates. J. Am.

Chem. Soc. 2011, 2011 133, 1870-1876. (72) Gorman, P. D.; English, N. J.; MacElroy, J. M. Dynamical and energetic properties of hydrogen and hydrogen-tetrahydrofuran clathrate hydrates. Phys. Chem. Chem. Phys. 2011, 2011

13, 19780-19787. (73) Gorman, P. D.; English, N. J.; MacElroy, J. M. Dynamical cage behaviour and hydrogen migration in hydrogen and hydrogen-tetrahydrofuran clathrate hydrates. J. Chem. Phys. 2012, 2012 136, 044506.

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

(74) Cao, H.; English, N. J.; MacElroy, J. M. D. Diffusive hydrogen inter-cage migration in hydrogen and hydrogen-tetrahydrofuran clathrate hydrates. J. Chem. Phys. 2013, 2013 138, 094507. (75) Maekawa, T.; Itoh, S.; Sakata, S.; Igari, S.; Imai, N. Pressure and temperature conditions for methane hydrate dissociation in sodium chloride solutions. Geochem. J. 1995, 1995

29, 325-329. (76) Jager, M. D.; Sloan, E. D. The effect of pressure on methane hydration in pure water 2001 185, 89-99. and sodium chloride solutions. Fluid Phase Equilib. 2001, (77) Abascal, J. L. F.; Vega, C. A general purpose model for the condensed phases of water: TIP4P/2005. J. Chem. Phys. 2005, 2005 123, 234505. (78) Jorgensen, W. L.; Madura, J. D.; Swenson, C. J. Optimized intermolecular potential functions for liquid hydrocarbons. J. Am. Chem. Soc. 1984, 1984 106, 6638-6646. (79) Smith, D. E.; Dang, L. X. Computer-Simulations of NaCl Association in Polarizable Water. J. Chem. Phys. 1994, 1994 100, 3757-3766. (80) Andoh, Y.; Yoshii, N.; Fujimoto, K.; Mizutani, K.; Kojima, H.; Yamada, A.; Okazaki, S.; Kawaguchi, K.; Nagao, H.; Iwahashi, K. et al., M. MODYLAS: A Highly Parallelized General-Purpose Molecular Dynamics Simulation Program for Large-Scale Systems with Long-Range Forces Calculated by Fast Multipole Method (FMM) and Highly Scalable Fine-Grained New Parallel Processing Algorithms. J. Chem. Theory Comput. 2013, 2013 9, 3201-3209. (81) Johnson, C. A. Generalization of the Gibbs-Thomson equation. Surf Sci. 1965, 1965 3, 429-444. (82) Widom, B. Some Topics in the Theory of Fluids. J. Chem. Phys. 1963, 1963 39, 2808-2812. (83) 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, 2005 26, 1701-1718. (84) Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 2008, 2008

4, 435-447. (85) Duan, Z.; Mao, S. A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar. Geochim. Cosmochim. Acta 2006, 2006 70, 3369-3386. (86) Duan, Z.; Møller, N.; Greenberg, J.; Weare, J. H. The prediction of methane solubility in natural waters to high ionic strength from 0 to 250°C and from 0 to 1600 bar. Geochim.

Cosmochim. Acta 1992, 1992 56, 1451-1460. (87) Mancinelli, R.; Botti, A.; Bruni, F.; Ricci, M. A.; Soper, A. K. Perturbation of water structure due to monovalent ions in solution. Phys. Chem. Chem. Phys. 2007, 2007 9, 2959-2967.

26

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

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

(88) Wedekind, J.; Reguera, D.; Strey, R. Finite-size effects in simulations of nucleation. J.

Chem. Phys. 2006, 2006 125, 214505. (89) Weijs, J. H.; Seddon, J. R. T.; Lohse, D. Diffusive Shielding Stabilizes Bulk Nanobubble Clusters. Chemphyschem 2012, 2012 13, 2197-2204.

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

TOC image

28

ACS Paragon Plus Environment

Page 28 of 28