Why the Gas Uptake Behavior of Dry-Salt Water is Vastly Different

37 mins ago - By means of simulations, it was realized that the orientation of water and the hydrogen-bond network structure near the silica surface a...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Winnipeg Library

C: Surfaces, Interfaces, Porous Materials, and Catalysis

Why the Gas Uptake Behavior of Dry-Salt Water is Vastly Different above 279 K? A Dynamics-Controlled Process and Can be Intensified by Cooling Stimulation Method Jingpeng Hou, Wei Zhou, Xinrui Wang, and Dongsheng Bai J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b07612 • Publication Date (Web): 09 Nov 2018 Downloaded from http://pubs.acs.org on November 23, 2018

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

1

Why the Gas Uptake Behavior of Dry-Salt Water is Vastly Different

2

above 279 K? A Dynamics-Controlled Process and Can be Intensified

3

by Cooling stimulation Method

4

Jingpeng Houa, Wei Zhoua, Xinrui Wanga, Dongsheng Baia*

5

a

6

Industry, Beijing Technology and Business University, Beijing 100048, PR China

7

* [email protected]

Department of Chemistry, School of Science / Key Laboratory of Cosmetic, China National Light

8 9

Abstract

10

The effect of temperature on CO2 uptake behavior of dry-salt water was investigated by

11

kinetics measurements and molecular dynamics simulations. The gas uptake capacity was found

12

to be dramatically decreased at 279 K by direct cooling adsorption procedure. By means of

13

simulations, it was realized that the orientation of water and the hydrogen-bond network structure

14

near the silica surface are sensitive to temperature. With the temperature decreased, the water layer

15

becomes local-structured by forming more hydrogen bonds, facilitating the capture of CO2 by

16

forming hydrate precursor easily. The change of the structure is irreversible once the CO2 uptake

17

process occurred, hence the gas uptake can carry on continuously even if the temperature rises. A

18

cooling stimulation procedure was found to efficiently trigger gas uptake process at temperature

19

higher than 279 K. By using this method, gas adsorption can be induced due to the energy barrier

20

of nucleation is reduced. This technology can save a certain amount of energy used to keep the low

21

temperature. Therefore, it might be a potential industrial method for efficient gas storage.

22

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

23

1. Introduction

24

The efficient gas storage and transport, such as energy gas storage, become more and more

25

important in recent years since the consumption of traditional fossil energies is greatly accelerated

26

by the industrialization process. Porous materials, such as activated carbon, silica, zeolite, metal

27

organic frameworks (MOF), porous polymer networks (PPN), can be introduced to enhance the

28

adsorption of gases.1-4 Recent years have witnessed an ever-growing interest in gas/water interfaces

29

stabilized by particles because it contributes to new concepts and materials such as dry water.5 Dry

30

water (DW) is a water-in-air inverse foam produced by mixing water with hydrophobic silica

31

nanoparticles.6,7 Because it contains more than 95 wt% water, DW looks like a powder but exhibits

32

a free-flowing property. In comparison to bulk water and liquid marbles,8 DW has a much higher

33

surface-to-volume ratio. As such, the finely dispersed water droplets lead to greatly enhanced

34

kinetics of adsorption in a gaseous system,9,10 resulting in increased rate of methane hydrate

35

formation in comparison to bulk water.11-13 On the basis of efficient adsorptivity, thus, DW can be

36

considered as a material for gas capture, storage and transport, such as CH4, CO2 and Kr.7,11,14 More

37

importantly, compared with porous materials such as MOF, covalent organic frameworks (COF)

38

and PPN, DW has the advantages of cheap raw materials (only water and silica), simple preparation

39

process, complete desorption of gas after destruction of its structure, and low cost on re-preparation.

40

Therefore, it is considered to be the most likely gas storage material for industrial use.

41

At present, improving both the adsorption capacity and the formation rate of gas hydrate are

42

still key issues for practical application of DW,15,16 and many methods have been intensively

43

investigated to increase the formation rate of hydrate. Such as the use of vigorous mixing devices,

44

the addition of promoting materials, and increasing the gas-water contact area by using finely

45

ground ice particles.17,18 In addition, a mixed colloidal system made of hydrogel particles and DW

46

particles has been reported with a high adsorption capacity.19 Recently, a dry-K2CO3-containing

47

water prepared by surface modified mesoporous silica particles can enhance the CO2 capture in both

48

adsorption rate and capacity.20 Moreover, cyclodextrin was also found to be a kinetic promoter for

49

gas hydrate formation in both simulation21 and our previous experimental study.22 We should note

50

that although progressive advances made so far, gas storage by these methods is still hard to be

51

applied in industrial fields.

52

In energy industry field, the practical application of DW is constrained to both raw materials

ACS Paragon Plus Environment

Page 2 of 21

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

53

supply and energy saving.23 Since dry water is a free-flowing powder composed of water droplets,

54

numerous water is needed as raw material. If DWs can be prepared by mixing seawater (or non-

55

purified water containing electrolytes similar to seawater) with hydrophobic silica nanoparticles and

56

show an uptake capacity similar to that prepared by pure water, lots of cost on water purification

57

will be saved on the industrial scale. Furthermore, the effect of temperature on CH4 and CO2

58

adsorption kinetics has been reported,11,22 and only the temperature range from 273 K to 277 K was

59

considered to be optimal, which means a lot of energy will be consumed to keep the low temperature

60

during the whole gas adsorption period.24 If DWs have an acceptable uptake capacity when it works

61

at higher temperature, a large amount of energy consumption will be reduced in industrial field. As

62

far as we know, almost no related study concerning with these problems has been reported.

63

In this work, experiments have been carried out to investigate the CO2 uptake kinetics of DW

64

under different conditions. Considering that seawater in which the total amount of Na+ and Cl- is

65

~3wt%25 are abundant in earth, we prepared DWs by using sodium chloride solution (3 wt%) instead

66

of pure water. The DWs prepared with NaCl solutions are named dry-salt water (DSW). We found

67

that the gas uptake capacity is dramatically decreased to 0 when temperature is higher than 279 K.

68

Molecular dynamics (MD) simulations were performed to illustrate the microscopic mechanism.

69

Moreover, based on the results of temperature-variation simulations, a cooling stimulation method

70

was suggested for gas uptake by DSW, so that the gas adsorption can be easily induced at

71

temperature higher than 279 K, saving a certain amount of energy. It provides a new insight for the

72

industrial application of DW in gas storage.

73 74

2. Experimental Section

75

2.1. Materials and Synthesis

76

SIPERNAT D10, hydrophobic silica particles with the static contact angle of 119°, were

77

supplied by EVONIK (Germany). Pure carbon dioxide gas (99.99%) was supplied by Beijing

78

Yanglilai Chemical Gas Co., Ltd. (China), and sodium chloride (AR, 99.5%) was purchased from

79

Sinopharm Chemical Reagent Co., Ltd. (China).

80

The sodium chloride solution with concentration of 3 wt% was prepared firstly, and DSW was

81

prepared by mixing 5 g of hydrophobic D10 silica powders with 95 g of sodium chloride solution

82

at ~19000 rpm for 90 s in a domestic blender.

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

83

2.2. Kinetics Measurement

84

The prepared DSW of 25.0 g was loaded into a miniature high-pressure reactor (SLM-200, 300

85

mL, China), which is connected to a circulator bath (DTY-10B, China) for temperature controlling.

86

When DSW sample was cooled to a specified temperature, the reactor was pressurized to 3.5 MPa

87

with pure CO2 and sealed. Gas uptake kinetics in DSW was studied by observing the time evolution

88

of gas pressure. The temperature and pressure were recorded every 15 minutes until the adsorption

89

equilibrium is reached. To provide a qualitative assessment of the kinetics, we measured the CO2

90

uptake kinetics in DSW under different temperatures.

91 92

3. Model and Simulation Method

93

To understand the microscopic mechanism of gas uptake, we performed molecular dynamics

94

(MD) simulations by using LAMMPS.26 As is shown in Figure 1, the initial configuration contains

95

a NaCl solution/silica/carbon dioxide three-phase system, which was placed into a simulation box

96

with 20 × 8 3 × 16 nm3. In our simulations, we represent a DSW particle by a layer of

97

hydrophobic silica spheres with a certain amount of water and NaCl molecules, because DSW

98

particle is far greater than the scale of the classic MD simulations. The silica layers were prepared

99

by packing the amorphous silica spheres with diameter of 4 nm to a face-centered cubic motif

100

(Figure 1), whereas the dangling bonds on the surfaces of each silica sphere were saturated by –CH3

101

groups, showing hydrophobic property. In the initial configuration, 62800 water molecules and 600

102

Na+ and Cl- ions were added to the central part of the box to keep the 3 wt% salt solution has a

103

density of 0.997 g/cm3, and 2450 CO2 molecules were also added to the left side of the silica layer

104

and the right side of the salt solution to reach the pressure of 5 MPa (slightly higher than the

105

experiments). In this way, we can investigate gas uptake process at two kinds of interfaces at the

106

same time in same conditions: DSW interface and gas-liquid interface. We note that the direct gas-

107

liquid phase interface at the right side is just designed to compare the difference between DSW

108

interface and free liquid water interface. For a DSW particle, CO2 can only enter the aqueous water

109

core through silica pores, and there is no direct water-gas interface in perfect structure of DSW.

ACS Paragon Plus Environment

Page 4 of 21

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

110 111

Figure 1. Initial configuration of the simulation system at the projection on (a) xy plane and (b) yz plane. The

112

silica, water and CO2 molecules are represented by the stick models, NaCl ions are represented by purple and

113

green spheres, whereas hydrogen bonds are denoted by blue-dashed lines. Note that in panel (b), the water, CO2

114

and NaCl molecules are not shown for clarity.

115

In our simulations, the TIP4P/2005 water model27 in which the rigidity was restricted with

116

SHAKE algorithm28 was used. CO2 molecules were represented by the EPM2 model.29 The

117

hydroxylated silica model30 was adopted for the inner core of the silica spheres, whereas the –CH3

118

groups on the surface were represented by a single-point model.31 The potential parameters were

119

taken from ref. 32 to represent Na+ and Cl-. The unlike parameters of Lennard-Jones interactions

120

were obtained by the Lorentz-Berthelot mixing rule. The short-range interactions were cut at 12 Å,

121

whereas the long-range Coulomb interactions were evaluated using the PPPM algorithm.33 Periodic

122

boundary conditions were imposed in all three Cartesian directions.

123

A whole simulation run includes two steps: An NVT relaxation of 1 ns was performed at a

124

specified temperature of 260 K, 265 K, 270 K, 275 K, and 280 K, respectively, to eliminate the

125

effect of the initial configuration; then, a simulation at the same temperature with 500 ns was

126

performed with a time step of 2 fs to study the gas uptake mechanism. The temperature was

127

maintained by using the Nosé-Hoover algorithm,34-36 with a relaxation parameter of 0.2 ps. During

128

the simulations, the position of the silica spheres was fixed.

129 130

4. Results and Discussion

131

Like DW, the prepared DSW shows a good dispersibility, stability and flowability. Since there

132

is a certain amount of salt in the water core, the structure and the hydrogen-bond network of the

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

133

inner water will be different from pure water in low temperature, which may change the mechanism

134

and kinetics of gas uptake.

135

4.1. Uptake Kinetics of DSW at Constant Temperature

136

Many studies have shown that the gas uptake capacity is affected by temperature.11,22,37-39 Take

137

methane as an example, owing to the hydrate formation, the optimal temperature is about 273~277

138

K, in which it shows the largest uptake capacity.11 Similar to methane, carbon dioxide is also an

139

insoluble gas that can form hydrate too. We investigated the temperature effect for CO2 uptake in

140

DSW, and the gas uptake kinetics is shown in Figure 2. In this paper, the gas uptake capacity of

141

DSW is defined volumetrically, i.e. the volume of adsorbed gas in standard condition divided by

142

the volume of 25.0 g DSW we used (40 cm3), V/V. Considering the deviation from ideality, the

143

molar number of adsorbed gas was calculated first by using the Peng-Robinson equation of state

144

(EOS) according to the current experimental pressure and temperature,40 and then the volume of gas

145

in standard condition was converted by ideal gas EOS. We note that the compressibility factor of

146

CO2 in our experimental conditions (273~288 K, 3.0~3.5 MPa) is between 0.68 and 0.79 calculated

147

by Peng-Robinson EOS. One can see clearly that the optimal temperature (with the maximum

148

uptake capacity) is ~277 K in our experimental system, consistent with that in pure DW.22 The rapid

149

uptake at the initial stage at low temperature such as 273 K might be attributed to the hydrate

150

formation.

151 152

Figure 2. CO2 uptake kinetics in DSW at different temperatures. The initial pressure of CO2 in the reactor is 3.5

153

MPa.

ACS Paragon Plus Environment

Page 6 of 21

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

154

Surprisingly, when temperature is higher than 279 K, without warning, the gas uptake capacity

155

of DSW is dramatically decreased to 0. And we found that at 279 K, experiments under the same

156

conditions show completely different results: in most cases (about 75%) there is no adsorption in a

157

long time, and in a small part of cases (25%) there is a strong adsorption (~25 V/V), as marked in

158

the green band in Figure 2. We speculate that the CO2 adsorption in DSW is a dynamics-controlled

159

process at the moderate temperature of 279 K. The energy barrier restricts the adsorption process,

160

and the induction time for gas uptake probably increased obviously due to the uncertainty of

161

nucleation. By the aid of MD simulations, we explored the source of energy barrier in details.

162

4.2. Nucleation Mechanism for Gas Uptake

163

To illustrate the microscopic mechanism of gas uptake near the DSW interface, we performed

164

a series of MD simulations at 260 K, 265 K, 270 K, 275 K, and 280 K, respectively. We note that

165

the temperature we used does not match the experimental temperature, because the force field of

166

TIP4P/2005 represents a lower melting point of 252.1 K for water molecules.27 In addition, one

167

should note that during CO2 uptake process, there are two aspects that might need to be taken into

168

account to quantify the total adsorbed CO2: the uptake of CO2 by water based on the Henry constant,

169

and the acid-base reaction of carbon dioxide in water. Classical MD simulation cannot account for

170

the second one. Since the ionization equilibrium constants of carbonate and bicarbonate are so small

171

in the experimental conditions (~10-7 and ~10-11, respectively), the amount of ions formed during

172

CO2 uptake can be negligible, and it is acceptable to study the microscopic mechanism of CO2

173

uptake by using classical MD.

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

174 175

Figure 3. Density profile of CO2 at the end of the simulation time.

176

The density profiles of CO2 at the end of 500 ns are shown in Figure 3. One can see clearly

177

that at lower temperature, more CO2 molecules entered into the deep inside of the DSW. However,

178

at 275 K and 280 K, a large amount of CO2 is accumulated only near the silica particles to decrease

179

its surface free energy due to its hydrophobic nature, and is unable to enter the inner core of DSW.

180

We define the area of x range from -5 nm to -1 nm as region A (as is marked in Figure 3), and the

181

gas uptake dynamics in region A (i.e. the CO2 molecules get into region A from the left gas phase

182

by passing through the gap between silica spheres) is shown in Figure 4. At the initial stage, the

183

adsorption rate of each system is very similar. After ~50 ns, the low temperature systems can still

184

adsorb CO2 molecules, while the adsorption rate of the high temperature systems decreased

185

obviously until the equilibrium is reached at ~200 ns.

ACS Paragon Plus Environment

Page 8 of 21

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

186 187

Figure 4. Time evolution of the number of CO2 adsorbed in region A.

188

We think that the structure of water layers near the silica interface will be changed at different

189

temperatures, which can affect the mass transfer of CO2. Therefore, we calculated the orientation

190

distribution of water molecules in the first hydration layer of silica particles (i.e. water molecules

191

within 3.7 Å from the surface of silica, obtained from RDF curve shown in Figure 5). The orientation

192

of water molecule is defined as the angle between the water dipole and the connection between

193

centroid of silica and oxygen atom of the water, and the results are shown in Figure 5(a). At 275 K

194

and 280 K, the hydrogen atoms of water are more likely to point towards silica since the hydrophobic

195

property of its surface. Hence, the hydrogen bond network between water molecules is difficult to

196

establish, which can be seen from Figure 6. A hydrogen bond is identified when the distance of two

197

oxygen atoms of water molecules (or the distance of a carbon atom of –CH3 groups on the silica

198

surface and an oxygen atom of water) is within 3.5 Å and the H−O···O angle (or the H−C···O angle)

199

is less than 30°.41 When temperature decreases, more hydrogen bonds formed between water

200

molecules to reduce the energy of the system. Thus, the orientation of water relative to the silica

201

particle becomes disordered. The establishment of hydrogen bond network of water at the interface

202

area is a key factor to capture large amount of CO2: if the CO2 molecules close to the interface are

203

surrounded by local-structured water molecules with hydrogen bonds, they are easier to be captured

204

to form cage-like precursor structures of gas hydrates.42

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

205 206

Figure 5. (a)Orientation distribution of water molecules in the first hydration layer of silica particles at the initial

207

stage of the simulations. (b)Radial distribution function of silica sphere centroid-oxygen of water.

208 209

Figure 6. Hydrogen bonds per water molecule in the first hydration layer of silica particles at the initial stage of

210

the simulations.

211

For the direct gas-liquid interface located at the right side of the simulation box, we calculated

212

the orientation distribution of both water and CO2 molecules near the interface at different

213

temperatures. Here the orientation of water molecule is defined as the angle between the water

214

dipole and the direction of the normal vector of the interface pointing into liquid water, while the

ACS Paragon Plus Environment

Page 10 of 21

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

215

orientation of CO2 is defined as the C-O vector relative to the normal vector of the interface pointed

216

into bulk CO2. The results are presented in Figure 7(a) and (b), which are very similar to the study

217

by Singer et al.43 Different from the dipole distribution induced by local curvature at the silica

218

interface, the dipole of water molecules is almost parallel to the gas-liquid interface. Note that for

219

CO2 in bulk water phase, owing to the electrostatic interactions, a preferred configuration between

220

CO2 and its neighbor water molecule is shown in Figure 7(c). However, at gas-liquid interface, the

221

arrangement of water and CO2 is different from that in bulk system, which may prevent the capture

222

of CO2 molecules.

223 224

Figure 7. Orientation distribution of water and CO2 molecules in the gas-liquid interface at the final stage of the

225

simulations.

226

To investigate the temperature sensitivity on the structure of the silica interface, we performed

227

two additional simulations with temperature variation: (i) an NVT simulation of 200 ns at 275 K is

228

followed by a 100 ns simulation at 270 K, then temperature rises back to 275 K for 400 ns; (ii) same

229

as simulation (i) but the position of CO2 is fixed and the interactions between CO2 and other

230

molecules are set to zero during the first 350 ns (equivalent to no CO2). Then, at 350 ns, to avoid

231

the overlaps of atoms, a relaxing process of 10 ps with soft pairwise interaction was performed by

232

using the command “pair_style soft” in LAMMPS. The amount of CO2 adsorbed in region A and

233

the orientation distribution of water molecules are shown in Figure 8. It can be seen that the dipole

234

distribution of water induced by interface curvature is sensitive to temperature: the hydrophobic

235

effect is the main factor at high temperature while it is controlled by the hydrogen-bond network

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

236

when temperature decreased. If there are no CO2 molecules, part of hydrogen bonds is broken, and

237

water molecules return to the hydrophobic-controlled state when temperature rises. Once the CO2

238

is captured at low temperature, it will be surrounded by water molecules, which will change the

239

hydrogen-bond distribution of water. When temperature increases, the orientation distribution of the

240

interfacial water cannot be reversibly returned back, so that the gas uptake can carry on continuously.

241

242 243

Figure 8. Time evolution of the number of CO2 adsorbed in region A (a) and the orientation distribution of water

244

molecules (b) in temperature-variation simulations. In panel (a), the hollow dots indicate that the temperature is

245

275 K, and the solid dots indicate 270 K.

246 247

4.3. Cooling Stimulation for Gas Uptake The temperature sensitivity of the interfacial properties (the dynamic barrier) in temperature-

ACS Paragon Plus Environment

Page 12 of 21

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

248

variation simulations provides us with a way to enhance gas capture. We consider that accelerating

249

the nucleation through cooling stimulation may be a good method to eliminate the barrier. So, based

250

on the constant temperature uptake experiment (scheme I in Figure 9), we designed the following

251

experimental schemes. Firstly, same as scheme I, sample DSW of 25.0 g was cooled down to 279

252

K and kept at this temperature, and the reactor was then pressurized with pure CO2 to 3.5 MPa and

253

sealed. At this point, gas uptake is not occurring in most cases. Subsequently, the temperature is

254

adjusted to 277 K (cooling stimulation) and remained in a period of time for the starting of CO2

255

adsorption. Finally, the temperature is adjusted back to 279 K for gas uptake. The modified process

256

is marked as scheme II.

257 258

Figure 9. Experimental schemes used in our experiments.

259

The gas uptake kinetics by the cooling stimulation method is shown in Figure 10(a). At the

260

initial stage, almost no gas adsorption was observed (less than 1 V/V), showing the existence of the

261

barrier. When temperature decreased from 279 K to 277 K at 45 min, gas uptake occurs immediately.

262

When the CO2 uptake reaches ~15 V/V, we adjust the temperature back to 279 K. The uptake

263

process of CO2 goes on continuously, and the final gas uptake capacity reached 25 V/V, same as

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

264

the capacity of the experiments carried out at the constant temperature of 279 K (Figure 2), although

265

the two processes are different.

266

267 268

Figure 10. CO2 uptake kinetics in DSW by cooling stimulation. In the upper panel of (a) and (b), the hollow dots

269

indicate that the temperature is 279 K, and the solid dots indicate 277 K, which can also be seen in the lower panel.

270

As the simulation results suggested, the cooling stimulation changed the local structure of

271

DSW, thus the adsorption process is initiated. In scheme II, if the adsorption quantity is high enough

272

(more than 25 V/V) during the stimulation period of 277 K, we found that the gas will be desorbed

273

from the DSW when temperature is back to 279 K. The result is shown in Figure 10(b). As we

274

expected, the desorption process reaches an equilibrium after the volumetric gas capacity decreased

275

to ~25 V/V. Therefore, we considered that 25 V/V is the saturated uptake capacity of DSW at 279

276

K, which is controlled by thermodynamics. If there is no cooling stimulation, however, the uptake

ACS Paragon Plus Environment

Page 14 of 21

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

277

process will be difficult to start because of the barrier controlled by dynamics. Besides, in several

278

cooling stimulation experiments, we found that there is an obvious induction time for the CO2

279

uptake when temperature decreased to 277 K, and the results are shown in Figure 11. This is the

280

same as the viewpoint of classical nucleation theory: macroscopic transition can take place only

281

when the critical nuclei are formed, which requires a certain amount of time to cross the energy

282

barrier.

283 284

Figure 11. CO2 uptake kinetics in DSW with induction time by cooling stimulation. The hollow dots indicate that

285

the temperature is 279 K, and the solid dots indicate 277 K. Note that the three experiments in each panel are under

286

the same conditions.

287

To verify the simulation results of scheme (ii) in Figure 8, i.e. the effect of temperature on the

288

interface of DSW is reversible when there is no gas, we designed a similar experiment according to

289

scheme III shown in Figure 9: after the sample DSW underwent the cooling stimulation, pure CO2

290

is added to the reactor. However, no adsorption occurred, indicating that the structure of DSW is

291

very sensitive to the temperature. Only the special structure of DSW at low temperature can initiate

292

the adsorption process when it is exposed to CO2. In other words, the gas adsorption is strongly

293

related to the morphology of water in the inner core of DSW.

294

4.4. Applicability of Cooling Stimulation Method

295

Based on the cooling stimulation method, we further performed several experiments at other

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

296

temperatures. For example, at 280 K, DSW can uptake ~19 V/V of CO2 by cooling stimulation to

297

277 K. The undercooling of the system to 277 K is sufficient to provide the driving force for CO2

298

hydrate nucleation. After the nuclei formed, the gas uptake will continue even if the temperature

299

rises to 280 K. When temperature is higher than 283 K, however, the cooling stimulation method is

300

invalid: the gas captured by stimulation will be fully desorbed after the temperature rises to 283 K.

301

We note that 283 K is the upper limit of temperature that can be used by the cooling stimulation

302

method. The CO2 uptake capacity obtained by constant temperature method (scheme I) and cooling

303

stimulation method (scheme II) at different temperatures is shown in Figure 12.

304 305

Figure 12. Gas uptake capacity by scheme I and II at different temperatures.

306

Finally, by using cooling stimulation method, one can save a certain amount of energy. We

307

should note that cooling stimulation at higher temperature (e.g. at 280 K) can save more energy, but

308

the final uptake capacity of gas will be reduced. The optimal operating condition should be at a

309

temperature which has the lowest energy consumption per unit mass of CO2 adsorbed. The energy

310

consumption is related to the stimulation method used. In our experiments, we used a simple

311

circulator bath for temperature controlling, and 279 K is the best operating condition for cooling

312

stimulation method.

313

5. Conclusion

314

In this paper, we explored the capture process of small insoluble gas molecules such as CO2

315

by DSW. The DSW can be prepared by simple process with cheap raw materials. We found that

ACS Paragon Plus Environment

Page 16 of 21

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

316

there is a dramatic change in the gas uptake capacity at temperature higher than 279 K: DSW has

317

almost no gas adsorption capacity by constant temperature experiments. This is because the low

318

undercooling increases the induction time of nucleation. MD simulation results show that near the

319

silica surface, the orientation of water and the hydrogen-bond network structure are very sensitive

320

to temperature. With the temperature decreased, the water layer near the interface becomes local-

321

structured by forming more hydrogen bonds, which facilitated the capture of CO2 by forming

322

hydrate precursor easily. The change of the water layer structure is irreversible once the CO2 uptake

323

process occurred at low temperature, hence the gas uptake can carry on continuously even though

324

the temperature rises back. By using the cooling stimulation, gas adsorption can be induced because

325

the energy barrier of nucleation is reduced. However, the saturated uptake capacity of DSW is

326

controlled by thermodynamics, which is only related to the operating temperature, such as ~25 V/V

327

in our experiments at 279 K.

328

Finally, since the formation of gas hydrate precursor at low temperature is very common for

329

conventional insoluble small gas molecules, we note that the cooling stimulation technology can

330

also be used for the storage of a large class of gas. And compared with the direct cooling adsorption

331

process, it can save a certain amount of energy used to keep the low temperature. The DSW prepared

332

by NaCl solution (rather than by pure water) shows a good uptake capacity, which means the direct

333

usage of non-purified water containing electrolytes is possible, and hence the cost on water

334

purification will be also saved. Therefore, it might be a potential industrial method for efficient gas

335

storage. The key to the application of this method is how to intensify hydrate formation at higher

336

temperatures. Adding suitable additives might be a possible way, which will be explored in our

337

subsequent research.

338 339

Conflicts of interest

340

There are no conflicts to declare.

341 342

Acknowledgement

343

This work is supported by National Natural Science Foundation of China (No. 21403009).

344 345

References

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

346

(1) Billemont, P.; Coasne, B.; De Weireld, G. Adsorption of Carbon Dioxide, Methane, and Their Mixtures in Porous

347

Carbons: Effect of Surface Chemistry, Water Content, and Pore Disorder. Langmuir 2013, 29, 3328-3338.

348

(2) Pantatosaki, E.; Pazzona, F. G.; Megariotis, G.; Papadopoulos, G. K. Atomistic Simulation Studies on the

349

Dynamics and Thermodynamics of Nonpolar Molecules within the Zeolite Imidazolate Framework-8. J. Phys. Chem.

350

B 2010, 114, 2493-2503.

351

(3) Xiang, Z.; Leng, S.; Cao, D. Functional Group Modification of Metal–Organic Frameworks for CO2 Capture. J.

352

Phys. Chem. C 2012, 116, 10573-10579.

353

(4) Martin, R. L.; Shahrak, M. N.; Swisher, J. A.; Simon, C. M.; Sculley, J. P.; Zhou, H.-C.; Smit, B.; Haranczyk,

354

M. Modeling Methane Adsorption in Interpenetrating Porous Polymer Networks. J. Phys. Chem. C 2013, 117,

355

20037-20042.

356

(5) Dawson, R.; Stevens, L. A.; Williams, O. S. A.; Wang, W.; Carter, B. O.; Sutton, S.; Drage, T. C.; Blanc, F.;

357

Adams, D. J.; Cooper, A. I. ‘Dry bases’: Carbon Dioxide Capture Using Alkaline Dry Water. Energy Environ. Sci.

358

2014, 7, 1786-1791.

359

(6) Binks, B. P.; Murakami, R. Phase Inversion of Particle-Stabilized Materials from Foams to Dry Water. Nat.

360

Mater. 2006, 5, 865-869.

361

(7) Wang, W.; Bray, C. L.; Adams, D. J.; Cooper, A. I. Methane Storage in Dry Water Gas Hydrates. J. Am. Chem.

362

Soc. 2008, 130, 11608-11609.

363

(8) Ueno, K.; Hamasaki, S.; Wanless, E. J.; Nakamura, Y.; Fujii, S. Microcapsules Fabricated from Liquid Marbles

364

Stabilized with Latex Particles. Langmuir 2014, 30, 3051-3059.

365

(9) Farhang, F.; Nguyen, T. D.; Nguyen, A. V. Non-Destructive High-Resolution X-ray Micro Computed

366

Tomography for Quantifying Dry Water Particles. Adv. Powder Technol. 2014, 25, 1195-1204.

367

(10) Park, J.; Shin, K.; Kim, J.; Lee, H.; Seo, Y.; Maeda, N.; Tian, W.; Wood, C. D. Effect of Hydrate Shell

368

Formation on the Stability of Dry Water. J. Phys. Chem. C 2015, 119, 1690-1699.

369

(11) Carter, B. O.; Wang, W.; Adams, D. J.; Cooper, A. I. Gas Storage in "Dry Water" and "Dry Gel" Clathrates.

370

Langmuir 2010, 26, 3186-3193.

371

(12) Hu, G.; Ye, Y.; Liu, C.; Meng, Q.; Zhang, J.; Diao, S. Direct Measurement of Formation and Dissociation Rate

372

and Storage Capacity of Dry Water Methane Hydrates. Fuel Process. Technol. 2011, 92, 1617-1622.

373

(13) Fan, S.; Yang, L.; Wang, Y.; Lang, X.; Wen, Y.; Lou, X. Rapid and High Capacity Methane Storage in Clathrate

374

Hydrates Using Surfactant Dry Solution. Chem. Eng. Sci. 2014, 106, 53-59.

375

(14) Drage, T. C.; Snape, C. E.; Stevens, L. A.; Wood, J.; Wang, J.; Cooper, A. I.; Dawson, R.; Guo, X.; Satterley,

ACS Paragon Plus Environment

Page 18 of 21

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

376

C.; Irons, R. Materials Challenges for the Development of Solid Sorbents for Post-Combustion Carbon Capture. J.

377

Mater. Chem. 2012, 22, 2815-2823.

378

(15) Ribeiro, C. P.; Lage, P. L. C. Modelling of Hydrate Formation Kinetics: State-of-the-Art and Future Directions.

379

Chem. Eng. Sci. 2008, 63, 2007-2034.

380

(16) Zhang, J.; Lee, J. W. Enhanced Kinetics of CO2 Hydrate Formation under Static Conditions. Ind. Eng. Chem.

381

Res. 2009, 48, 5934-5942.

382

(17) Seo, Y. T.; Lee, H.; Moudrakovski, I.; Ripmeester, J. A. Phase Behavior and Structural Characterization of

383

Coexisting Pure and Mixed Clathrate Hydrates. Chemphyschem 2003, 4, 379.

384

(18) Zhong, Y.; Rogers, R. E. Surfactant Effects on Gas Hydrate Formation. Chem. Eng. Sci. 2000, 55, 4175-4187.

385

(19) Ding, A.; Yang, L.; Fan, S.; Lou, X. Reversible Methane Storage in Porous Hydrogel Supported Clathrates.

386

Chem. Eng. Sci. 2013, 96, 124-130.

387

(20) Rong, X.; Yang, H.; Zhao, N. Rationally Turning the Interface Activity of Mesoporous Silicas for Preparing

388

Pickering Foam and "Dry Water". Langmuir 2017, 33, 9025-9033.

389

(21) Ji, H.; Chen, D.; Wu, G. Molecular Mechanisms for Cyclodextrin-Promoted Methane Hydrate Formation in

390

Water. J. Phys. Chem. C 2017, 121, 20967-20975.

391

(22) Hou, J.; Zhou, W.; Bai, D.; Li, S.; Han, M. Interfacial Effect of Cyclodextrin Inclusion Complex on Gas

392

Adsorption Kinetics of Dry Water Emulsion. Colloids Surf., A 2018, 544, 8-14.

393

(23) Liu, P.; Georgiadis, M. C.; Pistikopoulos, E. N. Advances in Energy Systems Engineering. Ind. Eng. Chem.

394

Res. 2011, 50, 4915-4926.

395

(24) Wang, S.; Baldea, M. Temperature Control and Optimal Energy Management using Latent Energy Storage.

396

Ind. Eng. Chem. Res. 2013, 52, 3247-3257.

397

(25) Chave, K. E. Chemical Reactions and the Composition of Sea Water. J. Chem. Educ. 1971, 48, 148-151.

398

(26) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1-19.

399

(27) Abascal, J. L. F.; Vega, C. A General Purpose Model for the Condensed Phases of Water: TIP4P/2005. J. Chem.

400

Phys. 2005, 123, 234505-234516.

401

(28) Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C. Numerical Integration of the Cartesian Equations of Motion

402

of a System with Constraints: Molecular Dynamics of N-Alkanes. J. Comput. Phys. 1977, 23, 327-341.

403

(29) Harris, J. G.; Yung, K. H. Carbon Dioxide's Liquid-Vapor Coexistence Curve And Critical Properties as

404

Predicted by a Simple Molecular Model. J. Phys. Chem. 1995, 99, 12021-12024.

405

(30) Lopes, P. E. M.; Murashov, V.; Tazi, M.; Demchuk, E.; MacKerell, A. D. Development of an Empirical Force

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

406

Field for Silica. Application to the Quartz−Water Interface. J. Phys. Chem. B 2006, 110, 2782-2792.

407

(31) Jorgensen, W. L.; Madura, J. D.; Swenson, C. J. Optimized Intermolecular Potential Functions for Liquid

408

Hydrocarbons. J. Am. Chem. Soc. 1984, 106, 6638-6646.

409

(32) Smith, D. E.; Dang, L. X. Computer Simulations of NaCl Association in Polarizable Water. J. Chem. Phys.

410

1994, 100, 3757-3766.

411

(33) Hockney, R. W.; Eastwood, J. W. Computer Simulation Using Particles; Taylor & Francis: New York, 1989.

412

(34) Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A: At. Mol. Opt. Phys.

413

1985, 31, 1695-1697.

414

(35) Hoover, W. G. Constant-Pressure Equations of Motion. Phys. Rev. A: At. Mol. Opt. Phys. 1986, 34, 2499-2500.

415

(36) Melchionna, S.; Ciccotti, G.; Lee Holian, B. Hoover NPT Dynamics for Systems Varying in Shape and Size.

416

Mol. Phys. 1993, 78, 533-544.

417

(37) Prasad, P. S. R. Methane Hydrate Formation and Dissociation in the Presence of Hollow Silica. J. Chem. Eng.

418

Data 2015, 60, 304-310.

419

(38) Chari, V. D.; Raju, B.; Prasad, P. S. R.; Rao, D. N. Methane Hydrates in Spherical Silica Matrix: Optimization

420

of Capillary Water. Energy Fuels 2013, 27, 3679-3684.

421

(39) Wang, J.; Wang, R.; Yoon, R.-H.; Seol, Y. Use of Hydrophobic Particles as Kinetic Promoters for Gas Hydrate

422

Formation. J. Chem. Eng. Data 2015, 60, 383-388.

423

(40) Peng, D.-Y.; Robinson, D. B. A New Two-Constant Equation of State. Ind. Eng. Chem. Res. 1976, 15, 59-64.

424

(41) Zhang, J.; Hawtin, R. W.; Yang, Y.; Nakagava, E.; Rivero, M.; Choi, S. K.; Rodger, P. M. Molecular Dynamics

425

Study of Methane Hydrate Formation at a Water/Methane Interface. J. Phys. Chem. B 2008, 112, 10608-10618.

426

(42) Khurana, M.; Yin, Z.; Linga, P. A Review of Clathrate Hydrate Nucleation. ACS Sustainable Chem. Eng. 2017,

427

5, 11176-11203.

428

(43) Zhang, H.; Singer, S. J. Analysis of the Subcritical Carbon Dioxide−Water Interface. J. Phys. Chem. A 2011,

429

115, 6285-6296.

430 431

ACS Paragon Plus Environment

Page 20 of 21

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

432 433

TOC Image

434 435

436

ACS Paragon Plus Environment