Hydrogen Bonded Networks in Supercritical Water - The Journal of

Sep 3, 2014 - In this study, the structure of supercritical water (SCW) is investigated by Raman spectroscopy and molecular dynamics simulations. It w...
0 downloads 8 Views 636KB Size
Subscriber access provided by Colorado State University | Libraries

Article

The Hydrogen Bonded Networks in Supercritical Water Qiang Sun, Qianqian Wang, and Dongye Ding J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/jp503474s • Publication Date (Web): 03 Sep 2014 Downloaded from http://pubs.acs.org on September 9, 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 26

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

The Hydrogen Bonded Networks in Supercritical Water

2 3

Qiang Sun* and Qianqian Wang

4

Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, The School of

5

Earth and Planetary Sciences, Peking University, Beijing 100871, China

6

Dongye Ding

7

School of Resources and Environmental Engineering, Shandong University of Technology, Zibo

8

255049, China

9 10

ABSTRACT: In this study, the structure of supercritical water (SCW) is investigated by Raman

11

spectroscopy and molecular dynamics simulations. It was found that the hydrogen bonding in

12

water is closely related to temperature and pressure (or water density). According to the Raman

13

spectroscopic study of SCW, the existence of tetrahedral hydrogen bonds in SCW is also affected

14

by the density of SCW. In addition, for SCW with critical density (0.322 g·cm-3), we suggest that

15

the tetrahedral hydrogen bonding is absent at water critical point (647 K and 22.1 MPa) based on

16

Raman evidence. From the dependence of νmax of the Raman OH stretching bands on temperature

17

and pressure, the structure of SCW can be divided into three-dimensional and chain (or string)

18

hydrogen bonded networks, which correspond to liquid- and gas-like phases, respectively.

19

Keywords: Critical point, Tetrahedral hydrogen bonding, Raman spectroscopy, Molecular

20

simulations

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

23

INTRODUCTION

24

Recently, there has been increasing interest in the structure of water at high temperature and

25

pressure, especially in the supercritical domain. This is because the physical and chemical

26

properties of sub- and supercritical water (SCW) are remarkably different from those under

27

ambient conditions. Additionally, the density of SCW can be controlled between gas- and

28

liquid-like values by varying the pressure and temperature. Then, many fundamental properties,

29

such as the dielectric constant, viscosity, and ion product, can be continuously varied in SCW.1

30

Therefore, SCW can be treated as a highly reactive, polar, and hydrogen bonded solvent that plays

31

a crucial role in many chemical, biological, and industrial processes.2,3

32

Many experimental and theoretical studies have been conducted to investigate the structure of

33

SCW.4–29 According to the neutron scattering data of SCW at 673 K, Postorino et al.4 observed the

34

complete disappearance of the peak located at 1.9 Å, which is a direct evidence of hydrogen

35

bonding in the oxygen–hydrogen radial distribution function (gOH(r)). Therefore, the authors

36

claimed that almost all the hydrogen bonds between water molecules break at this temperature.

37

Later, according to the newly derived site–site pair correlation functions from the old data, the

38

above view was revised by Soper et al.,5 who suggested a highly distorted arrangement of water

39

molecules with a largely reduced degree of hydrogen bonding compared with water at ambient

40

conditions. Because O–H vibrations are sensitive to hydrogen bonding, Raman and infrared (IR)

41

spectroscopy are widely applied to investigate the structure of SCW. These methods indicate that a

42

significant fraction of network water is still present even at high temperatures, pressures, and

43

solute concentrations,13 and that supercritical water can be considered as an ideal mixture of small

44

water clusters ((H2O)n, n = 1–3) at the chemical equilibrium.20 In addition, according to NMR 2

ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26

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

45

studies of SCW, hydrogen bonding persists at supercritical temperatures, and the average number

46

of hydrogen bonds per water molecule is greater than one at the supercritical densities.21

47

According to the X-ray Raman spectrum of supercritical water at 380 °C and density 0.54 g⋅cm−3,

48

Wernet et al.23 suggested that supercritical water is inhomogeneous and can be regarded as a

49

heterogeneous system with small patches of bonded molecules in various tetrahedral

50

configurations and surrounding non-bonded gas-phase-like molecules. Recently, it has been

51

suggested that SCW can be separated into a gas-like phase with short chains (or sheets with local

52

planar geometries) and a liquid-like phase with three-dimensional percolating hydrogen bonded

53

networks according to whether it is below or above the percolation threshold, respectively.24,25

54

Parallel with those experimental efforts, there have been extensive theoretical works to present a

55

microscopic understanding of SCW,26-29 which are devoted to investigate the hydrogen bonding

56

network and its influence on physicochemical properties of SCW. It was found that hydrogen

57

bonds still exist in SCW, but hydrogen bond number decreases with the increase of temperature or

58

the decrease of the density. From the above studies, the hydrogen bonding in SCW still remains

59

elusive.

60

In this study, in situ Raman spectra of SCW are measured, and the Raman OH stretching bands

61

of water are applied to investigate the hydrogen bonding in water. The hydrogen bonding between

62

water molecules is closely related to the electron delocalization. In comparison with molecular

63

dynamics simulations, quantum chemistry calculations provide a more reasonable approach to

64

investigate the water structure. In this work, due to high-efficiency, molecular dynamics

65

simulations are conducted to study the structure of SCW.

66

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

67

EXPERIMENTAL

68

In this study, deionized water was produced by a Milli-Q water purification system, and its

69

conductivity was 5.5 × 10−6 S/m. The experiments were conducted in an externally heated

70

hydrothermal diamond anvil cell (HDAC)30 equipped with low-fluorescence type Ia diamonds

71

with culet faces 1 mm across. The temperature was measured by a pair of type-K thermocouples

72

attached to both the upper and lower diamond anvils. The pressure at the elevated temperatures

73

was determined from the measured liquid–vapor homogenization temperature using the

74

International Association for the Properties of Water and Steam 1995 equation of state of water.31

75

The errors in the pressure and temperature were estimated to be 5% and ±2 K, respectively.

76

Additionally, to investigate the effects of pressure on the structure of ambient water, high-pressure

77

experiments were conducted up to 250 MPa at 295 K. The experimental pressure was calculated

78

according to the Raman shift of the 464 cm−1 peak of quartz. According to Schmidt and

79

Ziemann,32 this pressure estimation method had an uncertainty of ±50 MPa.

80

The unpolarized Raman spectra were recorded using a confocal micro-Raman system in a

81

backscattering geometry. An Ar+ ion laser was used with an excitation wavelength of 514.5 nm

82

and operated at 25 mW. The spectrometer with an entrance slit of 50 µm was used to collect the

83

signals. The resolution was about 1 cm−1. Each Raman spectrum was recorded for 30 s with a 20×

84

objective.

85

The Raman spectra were analyzed by the Jandel Scientific Peakfit v4.04 program. The Raman

86

spectra were firstly smoothed until the noise diminished, and the baselines were corrected. Then,

87

Gaussian functions were used to fit the Raman OH stretching bands of water.

88

4

ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26

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

89

MOLECULAR DYNAMICS SIMULATIONS

90

In this study, owing to its simplicity and popularity, the well-established extended simple point

91

charge (SPC/E) potential33 was used to model the water–water interaction energy. The SPC/E

92

potential is known to predict well the temperature and pressure dependent number density in

93

liquid to supercritical water while reproducing the structural properties from ambient to

94

supercritical conditions.34 The water–water interaction energy is given by

95

12 6  σ   σ wt − wt   wt − wt VW −W (r ) = ∑ + 4ε wt − wt   −    rOO   rOO   ij 4πε 0 rij

qi q j

96

In this work, the molecular dynamics simulations were performed using the GROMACS

97

package (version 4.07).35,36 The simulations were carried out in the NPT ensemble. The simulation

98

time was 6000 ps with a time step of 2 fs. Additionally, periodic boundary conditions were applied

99

in all three directions. The Nose–Hoover thermostat and Parrinello–Rahman pressure coupling

100

were used in the simulations. The Lennard–Jones interactions were truncated at 1.0 nm. The

101

particle mesh Ewald method was used to calculate the long-range electrostatics forces.

102 103

RESULTS AND DISCUSSION

104

For a single H2O molecule, the vibrational normal modes are 2A1 (including a symmetric

105

stretching vibration ν1 at 3657.05 cm−1 and a bending vibration ν2 near 1595 cm−1) and B1

106

(anti-symmetric stretching vibration ν3 at 3755.97 cm−1),37 which are all Raman active. When a

107

hydrogen bond forms between two water molecules, electron redistribution occurs. For water

108

molecular clusters, hydrogen bond formation increases the O–H bond lengths, and it causes a

109

20-fold greater reduction in the H–O and O–O distances.38 In addition, it causes a red shift of the

110

hydrogen-bonded OH stretching frequency, and the magnitude of the red shift increases with 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

111

cluster size.39 Therefore, the hydrogen bonding strength is closely related to the OH stretching

112

vibrations.

113

According to our recent study on the dependence of the OH vibration frequency on the size of

114

water clusters,40,41 when three-dimensional hydrogen bonded networks form, the OH stretching

115

vibrations are mainly dependent on the local hydrogen-bonded networks (local hydrogen bonding

116

refers to the interactions between a water molecule with neighboring molecules or with the

117

hydrogen-bonded networks in the first coordination shell of the molecule). This has been

118

confirmed by a Raman spectroscopic study of the effects of increasing pressure on water

119

structures.42 From our recent study,42 with increasing pressure up to 400 MPa at 293 K, the Raman

120

OH stretching bands slightly move to lower wavenumber. However, high pressure does not

121

obviously change the profile of the Raman OH stretching band. Based on our interpretation of the

122

Raman OH stretching bands, this indicates that high pressure has no obvious effects on the first

123

shell of a water molecule. These results are in agreement with other studies.43,44

124

From the above analysis, it is reasonable to assign different OH vibrational frequencies to

125

different local hydrogen bonding motifs. For a water molecule, the local hydrogen-bonded

126

networks can be differentiated by whether the molecule forms hydrogen bonds as a proton donor

127

(D), proton acceptor (A), or a combination of both with neighboring molecules. Under ambient

128

conditions, the main local hydrogen bonding motifs for a water molecule can be classified as

129

double donor–double acceptor (DDAA), double donor–single acceptor (DDA), single

130

donor–double acceptor (DAA), and single donor–single acceptor (DA).40,41 At ambient conditions,

131

the Raman OH stretching band of water can be deconvoluted into five sub-bands located around

132

3041, 3220, 3430, 3572, and 3636 cm−1, which can be assigned to the νDAA-OH, νDDAA-OH, νDA-OH, 6

ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26

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

133

νDDA-OH, and free OH symmetric stretching vibrations, respectively.41

134

To date, many experimental and theoretical studies have been performed to investigate the

135

hydrogen bonding in ambient water. However, there still remains strong debate on the structure of

136

liquid water.45–47 This study focuses on the structure of SCW. Therefore, much attention is paid to

137

the dependence of vmax on temperature and pressure and the Raman spectral changes of tetrahedral

138

hydrogen bonding (3220 cm−1 sub-band for ambient water).

139

For water from ambient to supercritical conditions, the in situ Raman OH stretching bands of

140

water are shown in Figure 1. With increasing temperature, the Raman OH stretching bands move

141

to higher wavenumbers. Simultaneously, the increase of temperature results in the decrease of the

142

intensity of the 3220 cm−1 sub-band (tetrahedral hydrogen bonding), and the Raman OH stretching

143

bands become increasingly symmetric. These results are in agreement with other Raman

144

spectroscopic studies of water at high temperature.13–16 This indicates that increasing temperature

145

results in the decrease of hydrogen bonding strength. Because tetrahedral hydrogen bonding is the

146

main three-dimensional hydrogen bonding, and taking into consideration the dependence of OH

147

vibrations on the size of water clusters,40,41 it is expected that the disappearance of tetrahedral

148

hydrogen bonds may lead to the structural transformation from three-dimensional to chain or

149

string (νDA-OH) hydrogen bonding in SCW.

150

In this study, the νmax is used to investigate the spectral changes of the Raman OH stretching

151

bands of water. The effect of temperature on νmax of the Raman spectra is shown in Figure 2. With

152

increasing temperature along the gas–liquid equilibrium line, νmax moves to high wavenumber. For

153

SCW, the shift of νmax of the Raman spectra is affected by both pressure and temperature.

154

Therefore, the hydrogen bonding in SCW should be closely related to the temperature and 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

155

pressure (or SCW density). Additionally, in reference with νmax in the equilibrium region, the shift

156

of νmax of SCW at 25 MPa is more obvious than that at 100 MPa (Figure 2). This is in agreement

157

with other Raman spectroscopic studies of νmax of SCW.13 From this, the structural changes in

158

water critical region may be different from those of SCW apart from the critical region.

159

For water at ambient temperature, the primary effect of increasing pressure is to modify the

160

second water shell around a central water molecule, pulling it inwards, while changes to the first

161

shell are relatively small.43,44 To investigate the effect of pressure on the structure of SCW, the

162

dependence of νmax on pressure is shown in Figure 3. From the figure, the effect of increasing

163

pressure on νmax is closely related to the density of SCW, and SCW can be divided into a gas-like

164

phase and a liquid-like phase. For high-density (liquid-like) SCW, there is a slight decrease of νmax

165

with increasing pressure. This is similar to the effect of pressure on water at ambient temperature.

166

However, for low-density (gas-like) SCW, νmax significantly decreases with increasing pressure.

167

This means that the hydrogen bonded networks in SCW are closely related to the SCW density,

168

and there is an obvious structural difference between low- and high-density SCW.

169

From the above dependence of νmax on pressure, different hydrogen bonded characteristics can

170

be expected in SCW with different density. For SCW with high density, similar to ambient water,

171

its structure can be regarded as three-dimensional hydrogen bonded networks. However, from the

172

OH vibrations of molecular clusters, low-density SCW may be composed of dimer (H2O)2 and

173

trimer (H2O)3 clusters, as described by Tassaing et al.17,20 Therefore, for gas-like SCW,

174

three-dimensional hydrogen bonding is lost, and hydrogen bonding is present in the small

175

molecular clusters as chains (or strings) of hydrogen bonds. From the above explanation of the

176

Raman OH stretching bands, the tetrahedral hydrogen bonded network is the predominant 8

ACS Paragon Plus Environment

Page 8 of 26

Page 9 of 26

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

177

three-dimensional hydrogen bonding motif in water. Therefore, it is reasonable to regard

178

tetrahedral hydrogen bonding as the structural characteristic to differentiate the different structural

179

types of SCW. Additionally, the critical region of SCW located at the transition between the liquid-

180

and gas-like phases may correspond to the structural transformation from three-dimensional to

181

chain (or string) hydrogen bonding.

182

According to the dependence of OH vibrations on water molecular clusters (H2O)n,40,41 OH

183

vibrations are closely related to the structural characteristics of hydrogen bonding. When the size

184

of the clusters is less than the hexamer (n < 6), only chain (or string) hydrogen bonding is

185

expected in the small clusters, and the increase of cluster size leads to an obvious decrease of

186

DA-OH stretching vibrations. However, when three-dimensional hydrogen bonding occurs (n ≥ 6),

187

different local hydrogen bonding motifs correspond to different OH vibrational frequencies, and

188

the increase of cluster size does not affect the OH vibration frequency. Thus, the disappearance of

189

three-dimensional hydrogen bonded networks undoubtedly plays an important role in the νmax of

190

the Raman OH vibrations. Therefore, the obvious shift of νmax in the water critical region, as

191

shown in Figure 2, can be attributed to the structural transformation from three-dimensional to

192

chain (or string) hydrogen bonding. This is also in agreement with the above discussion on the

193

water critical region, and corresponds to structural transformation from three-dimensional to chain

194

(or string) hydrogen bonding in SCW.

195

From the Raman OH stretching bands of water, tetrahedral hydrogen bonding (νDDAA-OH

196

sub-band) is the main three-dimensional hydrogen bonded network. Therefore, it is important to

197

investigate the structural changes of tetrahedral hydrogen bonds in SCW. The Raman OH

198

stretching bands of SCW with densities from 0.11 to 0.62 g⋅cm−3 at 673 K are shown in Figure 4. 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

199

No tetrahedral hydrogen bonding sub-band is observed in SCW with low density (0.11 and 0.19

200

g⋅cm−3). However, the tetrahedral hydrogen bonding sub-band is observed in SCW when the

201

density is higher than the critical density. Additionally, the increase of the SCW density increases

202

the intensity of the tetrahedral hydrogen bonding sub-band. Therefore, as with hydrogen bonding

203

in SCW, the existence of tetrahedral hydrogen bonded networks in SCW is also closely related to

204

the water density.

205

From the above, for SCW in critical region, the obvious shift of νmax should be closely related to

206

the disappearance of tetrahedral hydrogen bondings, which leads to the structural transformation

207

from three-dimensional to chain (or string) hydrogen bonds. Therefore, it is necessary to

208

investigate the structural behaviors of tetrahedral sub-band at water critical point. Figure 4 shows

209

the Raman OH stretching bands of SCW at 646 K and 0.39 g⋅cm-3, and 647 K and 0.11 g⋅cm-3.

210

From the figure, no tetrahedral hydrogen bonds can be found in the gas-like SCW, and only weak

211

tetrahedral sub-band can be observed in the SCW with density slightly higher than critical density.

212

After carefully reexamining other Raman spectra of SCW at critical point,13 and considering the

213

dependence of hydrogen bonding on SCW density, we suggest that the critical point may

214

correspond to the disappearance of tetrahedral hydrogen bondings in SCW. Therefore, the critical

215

point should connect to a known continuous phase transition in water just described in Alphonse et

216

al. study.48

217

In this work, to investigate the structure of SCW, molecular dynamics simulations were also

218

carried out to obtain the radial distribution functions (RDFs) gHH, gOH, and gOO, which are shown

219

in Figure S1 in the Supporting Information. For ambient water, the simulated radial distribution

220

functions are in good agreement with the experimental results of Soper.49 With increasing 10

ACS Paragon Plus Environment

Page 10 of 26

Page 11 of 26

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

221

temperature from 295 to 647 K, the nearest neighbor hydrogen-bonded OH and HH peaks in the

222

gOH and gHH functions become gradually smeared, clearly indicating that the hydrogen bonding

223

interactions between nearest neighbor water molecules considerably weaken. In addition, the

224

increase of temperature also leads to broadening of the first OO peak, and lowers the intensity of

225

the second peak at 4.5 Å, which is the characteristic peak of the tetrahedral arrangement of the

226

water molecules at ambient conditions. At the critical point, a well-resolved peak at 1.9 Å can still

227

be detected in gOH, which is the signature of the presence of hydrogen bonding, and the tetrahedral

228

hydrogen bonding in gOO is hardly distinguishable. For SCW at 673 K, inspection of the site–site

229

radial distribution functions confirms that SCW at 0.11 g·cm−3 is gas-like, with all three RDFs

230

slowly converging to their asymptotic values without oscillations. This is different from the RDFs

231

of liquid-like SCW.

232

To measure the local structure in SCW, the hydrogen bonding number (average number of

233

hydrogen bonds per water molecule, nHB) was calculated from the molecular dynamics simulations.

234

Despite the fundamental role of hydrogen bonding in molecular processes, there is still not a

235

robust, self-consistent, and universally accepted definition of a hydrogen bond.50,51 In this work,

236

the geometrical definition of hydrogen bonding is used to determine the hydrogen bonds in water.

237

This definition makes use of both the oxygen–oxygen distance (rOO) and the ∠OOH angle

238

between two water molecules. For two water molecules to be considered hydrogen bonded, we

239

used 3.5 Å and 30° for the upper bounds of rOO and ∠OOH, respectively, based on a previous

240

study.52 From Figure 5, with increasing temperature, the decrease of nHB is closely related to the

241

water density. At the critical point, the hydrogen bonding number was determined to be 1.536.

242

From the above discussion on the structure of SCW, this value may be used to discriminate the 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

243

type of SCW (gas-like or liquid-like).

244

From the above, the structure of SCW can be divided into three-dimensional and chain (or

245

string) hydrogen bonded networks, which correspond to liquid- and gas-like phases, respectively.

246

The critical region of SCW located at the transition between the liquid- and gas-like phases may

247

correspond to the structural transformation from three-dimensional to chain (or string) hydrogen

248

bonding. Additionally, the tetrahedral hydrogen bonding is the primary three-dimensional

249

hydrogen bonding in water, and can be regarded as the structural characteristic to differentiate the

250

different structural types of SCW.

251 252

CONCLUSIONS

253

In this work, Raman spectroscopy and molecular dynamics simulations were used to investigate

254

the structure of supercritical water. We found that the hydrogen bonding in water is closely related

255

to temperature and pressure (or water density). According to the Raman spectroscopic study of

256

SCW, we suggest that the water critical point may correspond to the disappearance of tetrahedral

257

hydrogen bonding in SCW. From the dependence of νmax of the Raman OH stretching bands on

258

temperature and pressure, the structure of SCW can be divided into three-dimensional and chain

259

(or string) hydrogen bonded networks, which correspond to liquid- and gas-like SCW, respectively.

260

Therefore, the critical region of SCW is located at the structural transition between the

261

three-dimensional and chain (or string) hydrogen bonded networks.

262 263

ASSOCIATED CONTENT

264

* Supporting Information 12

ACS Paragon Plus Environment

Page 12 of 26

Page 13 of 26

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

265

The fitted equations of the νmax of supercritical water and the site–site radial distribution functions

266

of simulated water. This information is available free of charge via the Internet at

267

http://pubs.acs.org.

268 269

AUTHOR INFORMATION

270

Corresponding Author

271

*E-mail: [email protected]

272

Tel: +86 136 9358 6882

273

Notes

274

The authors declare no competing financial interest.

275 276

ACKNOWLEDGEMENTS

277

The reviewers are greatly appreciated for providing good suggestions to revise the paper. This

278

work is supported by the National Natural Science Foundation of China (Grant Nos. 41073048

279

and 41373057).

280 281 282 283 284 285 286 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

287

REFERENCES

288 289 290 291 292 293 294 295 296

(1) Shaw, R.W.; Brill, B.; Clifford, A.; Eckert, A.; Franck, E.U. Supercritical Water-A Medium for Chemistry. Chem. Eng. News 1991, 69, 26–39 (2) Akiya, N.; Savage, P. Roles of Water for Chemical Reactions in High-temperature Water. Chem. Rev. 2002, 102, 2725–2750 (3) Bröll, D.; Kaul, C.; Krämer, A.; Krammer, P.; Richter, T.; Jung, M.; Vogel, H.; Zehner, P. Chemistry in Supercritical Water. Angew. Chem. Int. Edit. 1999, 38, 2998–3014 (4) Postorino, P.; Tromp, R.H.; Ricci, M.A.; Soper, A.K.; Neilson, G.W. The Interatomic Structure of Water at Supercritical Temperatures. Nature 1993, 366, 668–670

297

(5) Soper, A.K.; Bruni, F.; Ricci, M.A. Site-site Pair Correlation Functions of Water from 25 to

298

400 °C: Revised Analysis of New and Old Diffraction Data. J. Chem. Phys. 1997, 106,

299

247–254

300

(6) Bellissent-Funel, M.-C.; Tassaing, T.; Zhao, H.; Beysens, D.; Guillot, B.; Guissani, Y. The

301

Structure of Supercritical Heavy Water as Studied by Neutron Diffraction. J. Chem. Phys.

302

1997, 107, 2942–2949

303 304 305 306

(7) Tassaing, T.; Bellissent-Funel, M.-C.; Guillot, B.; Guissani, Y. The Partial Pair Correlation Functions of Dense Supercritical Water. Europhys. Lett. 1998, 42, 265–270 (8) Soper, A.K. The Radial Distribution Functions of Water and Ice from 220 to 673 K and at Pressures up to 400 MPa. Chem. Phys. 2000, 258, 121–137

307

(9) Otomo, T.; Iwase, H.; Kameda, Y.; Matubayasi, N.; Itoh, K.; Ikeda, S.; Nakahara, M. Partial

308

Pair Correlation Functions of Low-density Supercritical Water Determined by Neutron 14

ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26

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

309

Diffraction with the H/D Isotopic Substitution Method. J. Phys. Chem. B 2008, 112,

310

4687–4693

311

(10) Chialvo, A.A..; Cummings, P.T. Microstructure of Ambient and Supercritical Water: Direct

312

Comparison between Simulation and Neutron Scattering Experiments. J. Phys. Chem. 1996,

313

100, 1309–1316

314 315

(11) Ohtaki, H.; Radnai, T.; Yamaguchi, T. Structure of Water under Subcritical and Supercritical Conditions Studied by Solution X-ray Diffraction. Chem. Soc. Rev. 1997, 26, 41–51

316

(12) Sahle, C.J.; Sternemann, C.; Schmidt, C.; Lehtola, S.; Jahn, S.; Simonelli, L.; Huotari, S.;

317

Hakala, M.; Pylkkänen, T.; Nyrow, A.; et al. Microscopic Structure of Water at Elevated

318

Pressures and Temperatures. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, 6301–6306

319

(13) Ikushima, Y.; Hatakeda, K.; Saito, N.; Arai, M. An in situ Raman Spectroscopy Study of

320

Subcritical and Supercritical Water: The Peculiarity of Hydrogen Bonding near the Critical

321

Point. J. Chem. Phys. 1998, 108, 5855–5860

322

(14) Yui, K.; Uchida, H.; Itatani, K.; Koda, S. Raman OH Stretching Frequency Shifts in

323

Supercritical Water and in O2 and Acetone-aqueous Solutions near the Water Critical Point.

324

Chem. Phys. Lett. 2009, 477, 85–89

325

(15) Hu, Q.; Lu, X.; Lu, W.; Chen, Y.; Liu, H. An Extensive Study on Raman Spectra of Water

326

from 253 to 753 K at 30 MPa: A New Insight into Structure of Water. J. Mol. Spectrosc.

327

2013, 292, 23–27

328

(16) Sugimoto, K.; Koshi, M.; Koda, S. Micro Solution-structure in O2-H2O Mixtures near the

329

Critical Point Studied by Raman Spectroscopic Method. 14th International Conference on the

330

Properties of Water and Steam in Kyoto 2004, 182–187 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

331 332

Page 16 of 26

(17) Tassaing, T.; Danten, Y.; Besnard, M. Supercritical Water: Local Order and Molecular Dynamics. Pure Appl. Chem. 2004, 76, 133–139

333

(18) Ricci, M.A.; Nardone, M.; Fontana, A.; Andreani, C.; Hahn, W. Light and Neutron Scattering

334

Studies of the OH Stretching Band in Liquid and Supercritical Water. J. Chem. Phys. 1998,

335

108, 450–454

336

(19) Kandratsenka, A.; Schwarzer, D.; Vöhringer, P. Relating Linear Vibrational Spectroscopy to

337

Condensed-phase Hydrogen-bonded Structures: Liquid-to-supercritical Water. J. Chem. Phys.

338

2008, 128, 244510

339

(20) Tassaing, T.; Garrain, P.A.; Bégué, D.; Baraille, I. On the Cluster Composition of

340

Supercritical Water Combining Molecular Modeling and Vibrational Spectroscopic Data. J.

341

Chem. Phys. 2010, 133, 034103

342 343 344 345 346

(21) Matubayasi, N.; Wakai, C.; Nakahara, M. Structural Study of Supercritical Water. I. Nuclear Magnetic Resonance Spectroscopy. J. Chem. Phys. 1997, 107, 9133–9140 (22) Hoffmann, M.M.; Conradi, M.S. Are There Hydrogen Bonds in Supercritical Water? J. Am. Chem. Soc. 1997, 119, 3811–3817 (23) Wernet, Ph.; Testemale, D.; Hazemann, J.L.; Argoud, R.; Glatzel, P.; Pettersson, L.G.M.;

347

Nilsson,

348

Hydrogen-bonding Disparities in Supercritical Water. J. Chem. Phys. 2005, 123, 154503

349

(24) Bernabei, M.; Botti, A.; Bruni, F.; Ricci, M.A; Soper, A.K. Percolation and

350 351 352

A.;

Bergmann,

U.

Spectroscopic

Characterization

of

Microscopic

Three-dimensional Structure of Supercritical Water. Phys. Rev. E 2008, 78, 021505 (25) Bernabei, M.; Ricci, M.A. Percolation and Clustering in Supercritical Aqueous Fluids. J. Phys.: Condens. Matter 2008, 20, 494208 16

ACS Paragon Plus Environment

Page 17 of 26

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

353 354

(26) Boero, M.; Terakura, K.; Ikeshoji, T.; Liew, C.C.; Parrinello, M. Water at Supercritical Conditions: A First Principles Study. J. Chem. Phys. 1999, 115, 2219–2227

355

(27) Skarmoutsos, I.; Dellis, D.; Samios, J. The Effect of Intermolecular Interactions on Local

356

Density Inhomogeneities and Related Dynamics in Pure Supercritical Fluids: A Comparative

357

Molecular Dynamics Simulation Study. J. Phys. Chem. B 2009, 113, 2783–2793

358

(28) Ma, H.; Ma, J. Density Dependence of Hydrogen Bonding and the Translational-orientational

359

Structural Order in Supercritical Water: A Molecular Dynamics Study. J. Chem. Phys. 2011,

360

135, 054504

361 362

(29) Kang, D.; Dai, J.; Yuan, J. Changes of Structure and Dipole Moment of Water with Temperature and Pressure: A First Principles Study. J. Chem. Phys. 2011, 135, 024505

363

(30) Bassett, W.A.; Shen, A.H.; Bucknum, M.; Chou, I.M. A New Diamond Anvil Cell for

364

Hydrothermal Studies to 2.5 GPa and from −190 to 1200 °C. Rev. Sci. Instrum. 1993, 64,

365

2340–2345

366

(31) Wagner, W.; Prub, A. The IAPWS Formulation 1995 for the Thermodynamic Properties of

367

Ordinary Water Substance for General and Scientific Use. J. Phys. Chem. Ref. Data 2002, 31,

368

387–535

369

(32) Schmidt, C.; Ziemann, M. In-situ Raman Spectroscopy of Quartz: A Pressure Sensor for

370

Hydrothermal Diamond-anvil Cell Experiments at Elevated Temperatures. Am. Mineral.

371

2000, 85, 1725–1734

372 373 374

(33) Berendsen, H.J.C.; Grigera, J.R.; Straatsma, T.P. The Missing Term in Effective Pair Potentials. J. Phys. Chem. 1987, 91, 6269–6271 (34) Hura, G.; Russo, D.; Glaeser, R.M.; Head-Gordon, T.; Krack, M.; Parrinello, M. Water 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

375

Structure as a Function of Temperature from X-ray Scattering Experiments and Ab Initio

376

Molecular Dynamics. Phys. Chem. Chem. Phys. 2003, 5, 1981–1991

377

(35) Berendsen, H.J.C.; van der Spoel, D.; van Drunen, R. GROMACS: A Message-passing

378

Parallel Molecular Dynamics Implementation. Comput. Phys. Commun. 1995, 91, 43–56

379

(36) Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4: Algorithms for Highly

380

Efficient, Loaded-balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput.

381

2008, 4, 435–447

382 383 384 385 386 387

(37) Fraley, P.E.; Rao, K.N. High Resolution Infrared Spectra of Water Vapor ν1 and ν3 Bands of H2O. J. Mol. Spectrosc. 1969, 29, 348–364 (38) Ludwig, R. The Effect of Hydrogen Bonding on the Thermodynamic and Spectroscopic Properties of Molecular Clusters and Liquids. Phys. Chem. Chem. Phys. 2002, 4, 5481–5487 (39) Xantheas, S.S. Cooperativity and Hydrogen Bonding Network in Water Clusters. Chem. Phys. 2000, 258, 225–231

388

(40) Sun, Q. The Raman OH Stretching Bands of Liquid Water. Vib. Spectrosc. 2009, 51, 213–217

389

(41) Sun, Q. Local Statistical Interpretation for Water Structure. Chem. Phys. Lett. 2013, 568,

390 391 392 393 394 395 396

90–94 (42) Sun, Q. Raman Spectroscopic Study of the Effects of Dissolved NaCl on Water Structure. Vib. Spectrosc. 2012, 62, 110–114 (43) Saitta, A.M.; Datchi, F. Structure and Phase Diagram of High-density Water: The Role of Interstitial Molecules. Phys. Rev. E 2003, 67, 20201 (44) Soper, A.K.; Ricci, M.A. Structures of High-density and Low-density Water. Phys. Rev. Lett. 2000, 84, 2881 18

ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26

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

397 398

(45) Nilsson, A.; Pettersson, L.G.M. Perspective on the Structure of Liquid Water. Chem. Phys. 2011, 389, 1–34

399

(46) Soper, A.K. Recent Water Myths. Pure Appl. Chem. 2010, 82, 1855–1867

400

(47) Clark, G.N.I.; Cappa, C.D.; Smith, J.D.; Saykally, R.J.; Head-Gordon, T. The Structure of

401

Ambient Water. Mol. Phys. 2010, 108, 1415–1433

402

(48) Alphonse, N.K.; Dillon, S.R.; Dougherty, R.C.; Galligan, D.K.; Howard, L.N. Direct Raman

403

Evidence for a Weak Continuous Phase Transition in Liquid Water. J. Phys. Chem. A 2006,

404

110, 7577–7580

405

(49) Soper, A.K. The Radial Distribution Functions of Water as Derived from Radiation Total

406

Scattering Experiments: Is There Anything We Can Say for Sure? ISRN Phys. Chem. 2013,

407

2013, 1−67

408 409 410 411 412 413

(50) Prada-Gracia, D.; Shevchuk, R.; Rao, F. The Quest for Self-consistency in Hydrogen Bond Definitions. J. Chem. Phys. 2013, 139, 084501 (51) Kumar, R.; Schmidt, J.R.; Skinner, J.L. Hydrogen Bonding Definitions and Dynamics in Liquid Water. J. Chem. Phys. 2007, 126, 204107 (52) Luzar, A.; Chandler, D. Effect of Environment on Hydrogen Bond Dynamics in Liquid Water. Phys. Rev. Lett. 1996, 76, 928

414 415 416 417 418

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

419

Figure 1. Raman OH stretching bands of water up to 1073 K and 383 MPa. The homogeneous

420

density is determined from the measurement of the homogeneous temperature, and the pressure is

421

determined by the equation of state of water.31

422 423 424

Figure 2. Dependence of νmax of the Raman OH stretching bands of water on temperature. The

425

fitted lines are shown to guide the eye.

426 427 428

Figure 3. The effect of pressure on νmax of SCW at 673, 773, and 873 K. The νmax of ambient

429

water up to 250 MPa is also shown. The solid lines are fitted lines, and the fitted equations of

430

SCW are listed in the Supporting Information.

431 432 433

Figure 4. Raman OH stretching bands of SCW at 673 K. The inset shows the Raman spectra of

434

water at the critical temperature.

435 436 437

Figure 5. Hydrogen bonding number per water molecule (nHB) as a function of density derived

438

from the simulations. The fitted lines are shown to guide the eye.

439 440 20

ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26

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

80x80mm (300 x 300 DPI)

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

80x80mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 26

Page 23 of 26

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

80x80mm (300 x 300 DPI)

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

80x80mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26

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

80x80mm (300 x 300 DPI)

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

50x50mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 26