Intracavity Conformational Changes in Clathrate Hydrates - The

Jul 25, 2016 - (7-10) There are three widely known crystal structures of clathrate hydrates, .... line-broadening factor, as implemented in EasySpin 5...
0 downloads 0 Views 638KB Size
Subscriber access provided by Northern Illinois University

Article

Intracavity Conformational Changes in Clathrate Hydrates Yun-Ho Ahn, Hyung-Kyu Lim, Hyery Kang, Hyungjun Kim, Minjun Cha, Kyuchul Shin, and Huen Lee J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b04511 • Publication Date (Web): 25 Jul 2016 Downloaded from http://pubs.acs.org on August 2, 2016

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 C 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 27

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

Intracavity Conformational Changes in Clathrate Hydrates Yun-Ho Ahn1, Hyung-Kyu Lim2, Hyery Kang1, Hyungjun Kim2, Minjun Cha3, Kyuchul Shin*4 and Huen Lee*1 1

Department of Chemical and Biomolecular Engineering (BK21+ program), Korea Advanced

Institute of Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 305-701, Republic of Korea. 2

Graduate School of EEWS, Korea Advanced Institute of Science and Technology (KAIST),

291 Daehak-ro, Yuseong-gu, Daejeon 305-701, Republic of Korea. 3

Department of Energy and Resources Engineering, Kangwon National University, 1

Kangwondaehak-gil, Chuncheon-si, Gangwon-do 200-701, Republic of Korea. 4

Major in Applied Chemistry, School of Applied Chemical Engineering, Kyungpook National

University, 80 Daehak-ro, Buk-gu, Daegu 702-701, Republic of Korea.

ACS Paragon Plus Environment

1

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 2 of 27

ABSTRACT

Acyclic hydrocarbon molecules favor the gauche- or cis- conformation, more stable in terms of molecular geometry when they are enclathrated in clathrate hydrates. Once they are captured, they keep their conformation in the hydrate cavities. However, on the basis of Raman spectra and density functional theory (DFT) calculations, we observed the conformational changes of acyclic guest molecule (3-buten-2-one) occurring in the hydrate cavities induced by the intercavity electron transfer after γ-irradiation. This phenomenon is also accompanied by the enhanced thermodynamic stability of clathrate hydrate phases. These results contribute to understand the complex physicochemical properties of clathrate hydrates and to suggest practical application for hydrate-based gas storage and separation.

ACS Paragon Plus Environment

2

Page 3 of 27

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 hydrates, which are composed of host water molecules and small gaseous or organic guest molecules,1 have been extensively investigated because of their potential for practical applications in energy and environmental engineering such as gas separation,2 methane production with carbon dioxide sequestration,3-6 and hydrogen storage.7-10 There are three widely known crystal structures of clathrate hydrates, structure-I (sI, cubic Pm-3n), structure-II (sII, cubic Fd-3m), and structure-H (sH, hexagonal P6/mmm), containing several types of large and small cavities.1 The sI clathrate hydrate is composed of six 51262 and two 512 cavities with 46 water molecules in a unit cell. On the other hand, the sII and sH clathrate hydrates are composed of eight 51264 and sixteen 512 cavities with 136 H2Os and one 51268, two 435663, and three 512 cavities with 34 H2Os, respectively.1 The hydrophobic guest-host interactions occurring in relatively larger cavities of sII (the 51264) and sH (the 51268) hydrates than sI (the 51262) hydrates often lead to higher thermodynamic stability of the clathrate structures.11-13 Many organic molecules having suitable molecular sizes to those cavities thus have been studied to form more stable clathrate hydrates with methane, carbon dioxide, or hydrogen gases for applications to gas storage and transportation.11-13 Several acyclic molecules having two or three rotational conformers are known to prefer the stable gauche conformation when those molecules are enclathrated in hydrate cavities. Isopentane, forming the sH hydrate with methane, was reported to show a gauche-gauche position in the 51268 cavity.14 Lee et al. also investigated several sH hydrate formers and their preferred conformer in the hydrate phases.15 One of other examples is n-butane, forming the sII hydrate with methane, hydrogen sulfide, or xenon.16-18 Although a trans conformer of n-butane is

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

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

Page 4 of 27

thermodynamically preferred in the gas phase, only the gauche form is observed in its clathrate hydrate.19-22 At this stage, a question arises: Is it possible that a guest molecule changes its conformation in a hydrate cavity? Isopentane and n-butane are converted to suitable conformations for hydrate cavities when they are enclathrated.14, 16-18 After enclathration, guest molecules keep the conformations and no cases of intracavity changes of guest conformation in clathrate hydrates have been reported to date. Recently, Ohgaki et al. reported that the interconversions between npropyl and i-propyl radicals or between i-butyl and tert-butyl radicals induced by hydrogen transfer to an adjacent cavity were observed in the annealing process of γ-irradiated propane and isobutane hydrates,23-24 but there is no conformational change of the propane and isobutane molecules. In this study, we report for the first time intracavity conformational changes of guest molecules in clathrate hydrates. For the study, 3-buten-2-one (also called methyl vinyl ketone), having a suitable molecular size to fit into the 51264 cavity, was chosen. 3-buten-2-one is a α,β-unsaturated carbonyl compound containing both a carbon-carbon double bond and a carbonyl group. The two conjugated double bonds can rotate 180° on the axis of a carbon-carbon single bond resulting in two conformers, s-cis and s-trans, in nature (Figure S1).25-26 It was reported that the existence ratio of s-cis and s-trans conformers is temperature dependent and that s-trans is the favored conformer at a low temperature condition on the basis of IR and Raman spectroscopic studies.2527

In this study, we attempted to form the clathrate hydrate of 3-buten-2-one with methane as a

help guest and investigate the conformational changes of 3-buten-2-one in the hydrate cavity by Raman spectroscopy.

ACS Paragon Plus Environment

4

Page 5 of 27

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

EXPERIMENTAL SECTION Synthesis of binary (3-buten-2-one + CH4) hydrate samples Deionized water of ultrahigh purity was supplied from Merck (Germany) and 3-buten-2-one was supplied by Sigma-Aldrich Inc. CH4 gas was purchased from Special Gas (Korea) with stated minimum purities of 99.99 mol %. An amount of 10 g of 3-buten-2-one solution having 5 mol % of concentration was initially loaded into the stirring cell. The stirring cell has a vertical magnetic drive agitator so that continuous stirring is possible when making binary hydrates. A four-wire type Pt-100Ω probe was used for the temperature sensing, and a pressure transducer (Druck, PMP5073) was also used for the pressure sensing inside the cell. After loading the solution, air existed inside the cell was flushed out by continuous injecting of CH4. After few seconds, CH4 was pressurized up to around 120 bar by using microflow syringe pump (Teledyne, ISCO 260D). Then, temperature of stirring cell was lowered with a circulating bath (Jeio Tech., RW-2025G) until binary hydrates were formed. After binary hydrate formation detected with increased temperature and decreased pressure, the cell was kept in a 263 K bath for sufficient time (3 days) to produce binary hydrate without any icy impurities. To recover the binary hydrate sample, the stirring cell was quenched in liquid nitrogen and after a few seconds, the pressure was released to atmospheric pressure. The binary hydrate particles were ground (~200 µm) in a mortar and pestle in liquid nitrogen. These binary hydrate samples were used for HRPD and Raman analysis without γ-ray irradiation. A part of sample was irradiated at 90 kGy dose (15 kGy per 1 h) by a 60Co γ-ray source at KAERI in Jeongeup, Korea. The sample was being immersed in liquid nitrogen during the irradiation. These γ-ray irradiated binary hydrate samples were used for Raman and ESR analysis.

ACS Paragon Plus Environment

5

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 6 of 27

Structure Analysis Using Synchrotron High Resolution Powder Diffraction (HRPD) The HRPD patterns were measured using the High-Resolution Powder Diffraction beamline (9B) at Pohang Accelerator Laboratory (PAL) in Korea. During the measurements, the θ/2θ scan mode with a fixed time of 2 s, and a step size of 0.005° for 2θ = 0-120° and the beamline with a wavelength of 1.5472 Å were used for each samples. The binary hydrate powder stored in liquid nitrogen was quickly transferred to the sample stage cooled down to 80 K in air, and experiment was carried out at around 77 K to minimize the possible sample damage. Conformational Characterization Using Raman Spectroscopy Raman spectra were measured using the Horiba Jobin Yvon LabRAM HR UV/Vis/NIR high resolution dispersive Raman microscope equipped with CCD detector and cooled by liquid nitrogen. The resolution of recorded spectra was 0.75 cm-1 at 1800 grating mode. The excitation source was an Ar-ion laser emitting a 514.53 nm line. The intensity of laser was typically 30 mW. For variable temperature experiment, LINKAM unit was used to control the sample temperature. At each temperature, sufficient time was provided to reach the equilibrium state. Radical Characterization Using Electron Spin Resonance (ESR) The electron spin resonance (ESR) spectrum were measured using the JEOL PX2300 equipment. The samples were kept under liquid nitrogen to maintain the radical species during the experiment. Multi-component fitting and simulation of the ESR spectrum were performed using the Garlic function with the Gaussian line broadening factor, implemented in EasySpin 5.0.19 software.28

ACS Paragon Plus Environment

6

Page 7 of 27

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

Calculation Methods Density Functional Theory (DFT) calculations for assigning Raman peaks have been carried out by using Jaguar 8.4.29 We used the Becke three-parameter functional (B3) combined with the correlation functional of Lee, Yang, and Parr (LYP) 30, and the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional 31 along with 6-31+G** basis sets. After fully optimizing the molecular structures of 3-buten-2-one in gas-phase state, we calculated the Hessians for vibrational frequencies by evaluating analytic gradients. From these gas-phase vibrational frequency data, we derived the scaling factors for each C=C and C=O peaks to match with the experimental gas-phase Raman peaks (Table S1). In order to estimate the peak shifts of vibrational modes when the guest molecules are in aqueous and enclathrated phase, we employed Poisson-Boltzmann (PB) implicit solvation method for aqueous phase and explicit 51264 water cage which is sampled from sII clathrate structure with experimentally refined lattice parameter in this work. Before inserting guest molecules into the 51264 cage, the positions of hydrogen atoms in cage were fully optimized while the oxygen atoms were fixed. We note that the frequency calculations of enclathrated guest molecule, we diagonalized the partial Hessian matrix of guest molecule with fixing the position of host cage structure. For direct and quantitative comparison with inducing an appropriate error cancellation, we scaled the calculated Raman frequencies a priori. The scaling factors were slightly differently chosen depending on peaks to align gas-phase experimental frequencies and DFT calculated frequencies, ranging from 0.96 to 1.01 as listed in Table S1. The final scaled and unscaled DFT calculated Raman frequencies are listed in Table S2.

ACS Paragon Plus Environment

7

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 27

RESULTS AND DISCUSSION From a synchrotron high resolution powder diffraction (HRPD) analysis, the binary 3-buten-2one + CH4 clathrate hydrate sample was proven to be sII clathrate hydrate of a cubic Fd-3m space group with a small amount of ice impurity (Figure 1). The lattice parameter obtained from pattern refinement by the Le Bail fitting method32 was 17.1959(3) Å at 80 K. Because we used a 5 mol % 3-buten-2-one solution, which is slightly less concentrated than the stoichiometric 5.6 mol % concentration for the complete occupation of 51264 cavities in the sII hydrate phase, so as to minimize unreacted 3-buten-2-one, methane molecules might occupy a fraction of the 51264 cavity. Although the molecular size of 3-buten-2-one is slightly larger than the mean-free diameter of 51264 cavities, it possibly occupies large cages of sII hydrates as other large molecules do due to the flexibility of water frameworks.33-34 As the next step, we checked the conformations of 3-buten-2-one in the aqueous solution and in the clathrate hydrate by Raman spectroscopy. The Raman spectra of the binary 3-buten-2-one + CH4 clathrate hydrate were measured with various initial temperatures from 93 to 213 K shown in Figure 2. The C‒H stretching modes of methane occupying the 512 and 51264 cavities of sII hydrate were observed at 2914 cm-1 and 2904 cm-1, respectively, and the sII 3-buten-2-one + CH4 hydrate kept the clathrate structure stable until 133 K. However, as the temperature increased, the hydrate started to dissociate and a large amount of methane already disappeared at 153 K. By the way, considering both the C=C band (1623 cm-1, s-trans; 1627 cm-1, s-cis) and the C=O band (1705 cm-1, s-trans; 1715 cm-1, s-cis) of neat 3-buten-2-one,26, 35 two bands of the binary 3-buten-2-one + CH4 clathrate hydrate observed at 1626 cm-1 and 1692 cm-1 imply that they are from the C=C bond and C=O bond of one 3-buten-2-one conformer, which is dominantly favored in the cavity (shown in the spectra recorded at 93 K, 113 K, and 133 K of

ACS Paragon Plus Environment

8

Page 9 of 27

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

The Journal of Physical Chemistry

Figure 2). The two tiny shoulders appearing at 1618 and 1672 cm-1 are considered to be from the unreacted 3-buten-2-one molecules that exist outside the hydrate media, or possibly from another geometrically less favored conformer in the cavity. As the hydrate started to dissociate from 153 K, two tiny peaks at 1618 cm-1 (C=C stretching mode) and 1672 cm-1 (C=O stretching mode) rapidly grew, implying that they are from the 3-buten-2-one released from the clathrate cavity. In Figure S2, the Raman spectra of the frozen 5 mol % 3-buten-2-one solution in the region from 1600 to 1725 cm-1 show two same peaks, at 1618 cm-1 for C=C and at 1672 cm-1 for C=O stretching modes, indicating that only one conformer exists from 93 K to 253 K. At this stage, we checked the conformation of 3-buten-2-one in the aqueous solution and in the clathrate hydrate by assigning the Raman band with the help of density functional theory (DFT) calculation. We adopted two different B3LYP and PBE exchange-correlation functionals for cross-checking the results, and the detailed calculation methods and assigned Raman frequencies are in the calculation methods of experimental section and the supporting information, respectively. First, both the B3LYP and PBE methods predicted that the Raman peaks will show a slight blue-shift (2~10 cm-1) from the aqueous solution to the clathrate phase for both s-cis and s-trans forms. However, both the C=C and C=O stretching modes of the Raman spectra (Figure 2) show larger peak shifts than the predicted shifts for the enclathration while retaining conformations of guests (8 cm-1 (1618→1626 cm-1) and 20 cm-1 (1672→1692 cm-1), respectively). The calculation results for the peak shifts from the aqueous s-trans form to the enclathrated s-cis form predicted a similar degree of peak shifts to the experimental observations. Therefore, it is reasonable that assigning s-cis form to observed peaks in clathrate phase. Considering the calculated Raman peaks in aqueous phase, it is also confirmed that solvated 3buten-2-one by water molecules prefers the s-trans form and the band observed in the aqueous

ACS Paragon Plus Environment

9

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 10 of 27

solution showed redshift due to the solvation effect compared to the bands of gas-phase 3-buten2-one.26, 35 Therefore, from the Raman spectra and DFT calculation results shown in Figure 2, the conformation change of 3-buten-2-one from the s-cis to s-trans is due to the hydrate dissociation and there is no intracavity conformational transformation. By the way, recent studies reported hydrogen atom transfer (referred to as the ‘hydrogenpicking’ phenomenon) through a hexagonal face of a 51264 cavity between hydrocarbon guests of γ-irradiated hydrates.23-24 In order to provide driving force for conformational transformation in a hydrate cavity, the 3-buten-2-one + CH4 clathrate hydrate was γ-irradiated to 90 kGy at 77 K and analysed again by Raman spectroscopy. At 93 K, the obtained Raman spectrum of the γ-irradiated sample is similar to that of the nonirradiated sample (Figure 3). As the temperature increases, however, two peaks at 1620 cm-1 and at 1674 cm-1 newly appear and increase until 133 K even though the hydrate phase does not start to dissociate. As shown in Figure 3, most of the methane stays in the 512 and the 51264 cavities of sII hydrate until the temperature increases to 173 K. Thus, newly appearing two peaks at 1620 cm-1 and at 1674 cm-1 strongly imply the existence of the s-trans form of 3-buten-2-one in the 51264 cavity of sII hydrate, corresponding to the slight calculated peak shifts of s-trans form from aqueous to hydrate phase. Meanwhile, calculation results predict that enclathrated s-cis and strans forms are slightly distorted around 10 and 5 degree, respectively in order to have geometrically optimized configuration in the cavity (Figure 4). However, this slight distortion does not induce high energy increment comparable with the third conformer of 3-buten-2-one, having a nonplanar or gauche structure.26

ACS Paragon Plus Environment

10

Page 11 of 27

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

One interesting aspect in Figure 3 is that the Raman signals of the s-trans form disappear at 153 K and reappear as the temperature rises to 173 K. Here, we note that the spectroscopic results shown in Figure 3 were reproducible for the other samples obtained from two additionally repeated experiments (Figure S3). Reappearing s-trans signals over 173 K are due to 3-buten-2one molecules released from the hydrate cavities because the intensities of C‒H stretching modes at 2904 and 2914 cm-1, originating from enclathrated methane molecules, start to decrease from 173 K. Therefore, we can conclude that there would be a defined temperature window (~ 133 K) where the s-trans conformer could exist stably in the hydrate phase, and the intracavity conformational transformation of 3-buten-2-one occurs in this defined temperature window. In order to identify the mechanism underlying this phenomenon, we obtained the electron spin resonance (ESR) spectrum of γ-irradiated 3-buten-2-one + CH4 hydrate at 77 K (Figure 5). The spectrum shows strong 1:3:3:1 hyperfine signals, indicating methyl radicals generated from the dissociation of CH4 molecules (coupling constant: 2.35 mT).36 The presence of small peaks merged in the signal of the methyl radical in the spectrum imply that there exist at least two different radicals in the irradiated sample, but the relatively low intensities of these peaks make identification of them difficult. One possible radical for the irradiated (3-buten-2-one + CH4) hydrate system is a ketyl radical anion (Figure 5), which is generated when the 3-buten-2one molecule captures a solvated electron in liquid ammonia.37 In order to check whether this ketyl radical was formed by the γ-irradiation, we simulated the ESR spectrum containing both a methyl radical and a ketyl radical anion. The EasySpin program was used for this simulation with the coupling constants of radicals obtained from the literature (2.3 mT for methyl radical36 and 0.92 mT for ‒CH3, 1.17 mT for =CH2, and 0.05 mT for =CH‒ groups of ketyl radical37) for the initial input. As shown in Figure 5, the simulated spectrum is well-matched with the

ACS Paragon Plus Environment

11

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 12 of 27

experimental spectrum and the refined values of the coupling constants after the simulation (2.35 mT for methyl radical and 1.11 mT / 1.42 mT / 0.07 mT for ‒CH3 / =CH2 / =CH‒ groups of ketyl radical) were in good agreement with the values in the literature, although there are small uncertain shoulders. Slight differences in the coupling constants between the literature and the simulation are ascribed to the enclathration effect of the ketyl radical. Therefore, we conclude that γ-irradiation generated methyl radicals from CH4 and ketyl radical anions from 3-buten-2one guest molecules, respectively. From the ESR spectroscopic results, a possible mechanism for the conformational change is as follows (Scheme 1): (1) A few 3-buten-2-one molecules are changed to ketyl radicals by γirradiation. (2) The carbon atom of the carbonyl group has a pyramidal geometry rather than triangular and freely rotates. (3) The electron of the ketyl radical transfers to neighboring 3buten-2-one as the temperature increases and forms a new ketyl radical. (4) After the electron transfer, the old ketyl radical is returned to 3-buten-2-one with either the s-cis or s-trans form. Because the mechanism proposed above is a kind of chain reaction, a few ketyl radical anions only can initiate a considerable number of conformational changes as observed in the Raman spectroscopy (Figure 3). Another remarkable feature of this phenomenon is that the dissociation temperature of 3-buten2-one + CH4 hydrate increases by over 20 degrees after γ-irradiation. As shown in Figure 2, most methane molecules of 3-buten-2-one + CH4 hydrate were already released at 153 K. On the other hand, after γ-irradiation, a large amount of methane was still inside the cages even at 173 K (Figure 3). Although the exact mechanism has not been revealed yet, this enhanced thermodynamic stability of the 3-buten-2-one hydrate by γ-irradiation can provide an important contribution to the modification of clathrate materials.

ACS Paragon Plus Environment

12

Page 13 of 27

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

CONCLUSION The guest conformation in the hydrate cavity often provides abundant information for the guesthost interactions. This work demonstrates the first example of the intracavity conformational change of an acyclic guest molecule enclathrated in a hydrate cage. The 3-buten-2-one molecule favors enclathration in the s-cis form in the 51264 cavity, similarly to the gauche form of n-butane in the sII hydrate,22 while the neat or solvated 3-buten-2-one prefers the s-trans over the s-cis conformation. The γ-irradiation can induce the intracavity conformational change of 3-buten-2one, i.e., from the s-cis to the s-trans form in the hydrate cavity. The γ-irradiation causes a few 3buten-2-one molecules to form ketyl radicals and the s-trans conformation appears when the ketyl radical returns to the 3-buten-2-one by electron transfer to the neighboring 3-buten-2-one guest in the 51264 cavity. This work also shows the possibility that the conformational transformation enhances the thermodynamic stability of clathrate hydrate phases. The present findings provide an important key to understand the radiolysis chemistry occurring in the clathrate phase and to modify the clathrate hydrate materials.

ACS Paragon Plus Environment

13

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 27

Figure 1. Synchrotron HRPD patterns of the 3-buten-2-one + CH4 hydrate recorded at 80 K. (Red circles, observed HRPD pattern; Black solid line, calculated HRPD pattern; tick marks, cubic Fd-3m (top) and hexagonal P63/mmc phases (bottom)).

ACS Paragon Plus Environment

14

Page 15 of 27

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

The Journal of Physical Chemistry

Figure 2. Raman spectra of 3-buten-2-one + CH4 hydrate measured at various temperatures. Enclathrated s-cis 3-buten-2-one (1626, 1692 cm-1) and CH4 molecule (2904, 2914 cm-1) were observed from 93 K to 153 K. Green and blue tick marks indicate the calculated Raman shifts for s-trans and s-cis conformers in the aqueous solution (top) and in the clathrate hydrate (bottom), respectively.

ACS Paragon Plus Environment

15

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 16 of 27

Figure 3. Raman spectra of the γ-irradiated 3-buten-2-one + CH4 hydrate measured at various temperatures. Enclathrated s-trans 3-buten-2-one was detected (1620, 1674 cm-1) from 93 K to 133 K. Green and blue tick marks indicate the calculated Raman shifts for s-trans and s-cis conformer in the aqueous solution (top) and in the clathrate hydrate (bottom), respectively.

ACS Paragon Plus Environment

16

Page 17 of 27

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. Geometrically optimized configurations of s-cis and s-trans forms in the clathrate hydrate cavities. (distortion angles: 9.5º and 5.1 º for s-cis and s-trans, respectively, using PBE; 10.6 º and 5.4 º for s-cis and s-trans, respectively, using B3LYP)

ACS Paragon Plus Environment

17

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 18 of 27

Figure 5. Experimental and simulated ESR spectra of the γ-irradiated 3-buten-2-one + CH4 hydrate measured at 77 K. The inset figure shows simulated ESR patterns from each guest molecule in the hydrate system.

ACS Paragon Plus Environment

18

Page 19 of 27

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

Scheme 1. A schematic illustration of the conformational change mechanism in large cavities of sII hydrates.

ACS Paragon Plus Environment

19

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 27

ASSOCIATED CONTENT Supporting Information. Structures of 3-buten-2-one conformers and their size, Raman spectra of 3-buten-2-one solution, Scaling factors used in DFT calculation, and Detailed results of Raman shifts from both experiment and calculations. This material is available free of charge via the Internet at http://pubs.acs.org.

AUTHOR INFORMATION Corresponding Author * (H.L.) E-mail: [email protected], Phone: +82-42-350-3917 * (K.S.) E-mail: [email protected], Phone: +82-53-950-5587 Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes Any additional relevant notes should be placed here.

ACKNOWLEDGMENT This research was funded by the the Nuclear Energy Research Infrastructure Program (NRF2015M2B2A4031401) through a National Research Foundation Korea (NRF) grant funded by

ACS Paragon Plus Environment

20

Page 21 of 27

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 Ministry of Science, ICT and Future Planning (MSIP). HRPD experiments at PLS (Beamline 9B) were supported by POSTECH.

ABBREVIATIONS sI, Structure-I; sII, Structure-II; sH, Structure-H; HRPD, High Resolution Powder Diffraction; ESR, Electron Spin Resonance; DFT, Density Functional Theory; B3LYP, Becke 3-Parameter (Exchange), Lee, Yang and Parr; PBE, Perdew-Burke-Ernzerhof.

ACS Paragon Plus Environment

21

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 22 of 27

REFERENCES (1)

Sloan, E. D., Fundamental Principles and Applications of Natural Gas Hydrates. Nature

2003, 426, 353-359. (2)

Kang, S. P.; Lee, H., Recovery of CO2 from Flue Gas Using Gas Hydrate:

Thermodynamic Verification through Phase Equilibrium Measurements. Environ. Sci. Technol. 2000, 34, 4397-4400. (3)

Lee, H.; Seo, Y.; Seo, Y. T.; Moudrakovski, I. L.; Ripmeester, J. A., Recovering Methane

from Solid Methane Hydrate with Carbon Dioxide. Angew. Chem. Int. Edit. 2003, 42, 50485051. (4)

Park, Y.; Kim, D. Y.; Lee, J. W.; Huh, D. G.; Park, K. P.; Lee, J.; Lee, H., Sequestering

Carbon Dioxide into Complex Structures of Naturally Occurring Gas Hydrates. P. Natl. Acad. Sci. USA 2006, 103, 12690-12694. (5)

Shin, K.; Park, Y.; Cha, M. J.; Park, K. P.; Huh, D. G.; Lee, J.; Kim, S. J.; Lee, H.,

Swapping Phenomena Occurring in Deep-Sea Gas Hydrates. Energ. Fuel 2008, 22, 3160-3163. (6)

Koh, D. Y.; Kang, H.; Kim, D. O.; Park, J.; Cha, M.; Lee, H., Recovery of Methane from

Gas Hydrates Intercalated within Natural Sediments Using CO2 and a CO2/N2 Gas Mixture. Chemsuschem 2012, 5, 1443-1448. (7)

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, 434, 743-746.

ACS Paragon Plus Environment

22

Page 23 of 27

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

(8)

Shin, K.; Kim, Y.; Strobel, T. A.; Prasad, P. S. R.; Sugahara, T.; Lee, H.; Sloan, E. D.;

Sum, A. K.; Koh, C. A., Tetra-N-Butylammonium Borohydride Semiclathrate: A Hybrid Material for Hydrogen Storage. J. Phys. Chem. A 2009, 113, 6415-6418. (9)

Koh, D. Y.; Kang, H.; Lee, H., Multiple Guest Occupancy in Clathrate Hydrates and Its

Significance in Hydrogen Storage. Chem. Commun. 2013, 49, 6782-6784. (10) Koh, D. Y.; Kang, H.; Jeon, J.; Ahn, Y. H.; Park, Y.; Kim, H.; Lee, H., Tuning Cage Dimension in Clathrate Hydrates for Hydrogen Multiple Occupancy. J. Phys. Chem. C 2014, 118, 3324-3330. (11) Mehta, A. P.; Sloan, E. D., Structure-H Hydrate Phase-Equilibria of Paraffins, Naphthenes, and Olefins with Methane. J. Chem. Eng. Data 1994, 39, 887-890. (12) Tohidi, B.; Danesh, A.; Todd, A. C.; Burgass, R. W.; Ostergaard, K. K., Equilibrium Data and Thermodynamic Modelling of Cyclopentane and Neopentane Hydrates. Fluid Phase Equilibr. 1997, 138, 241-250. (13) de Deugd, R. M.; Jager, M. D.; Arons, J. D., Mixed Hydrates of Methane and WaterSoluble Hydrocarbons Modeling of Empirical Results. Aiche J. 2001, 47, 693-704. (14) Shin, K.; Park, Y.; Hong, J. H.; Lee, H., Gauche Conformation of Acyclic Guest Molecules Appearing in the Large Cages of Structure-H Clathrate Hydrates. Korean J. Chem. Eng. 2007, 24, 843-846.

ACS Paragon Plus Environment

23

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 24 of 27

(15) Lee, J. W.; Lu, H. L.; Moudrakovski, I. L.; Ratcliffe, C. I.; Ohmura, R.; Alavi, S.; Ripmeester, J. A., 13C NMR Studies of Hydrocarbon Guests in Synthetic Structure H Gas Hydrates: Experiment and Computation. J. Phys. Chem. A 2011, 115, 1650-1657. (16) Wu, B. J.; Robinson, D. B.; Ng, H. J., 3-Phase and 4-Phase Hydrate Forming Conditions in Methane + Isobutane + Water. J. Chem. Thermodyn. 1976, 8, 461-469. (17) Davidson, D. W.; Garg, S. K.; Gough, S. R.; Hawkins, R. E.; Ripmeester, J. A., Characterization of Natural-Gas Hydrates by Nuclear Magnetic-Resonance and DielectricRelaxation. Can. J. Chem. 1977, 55, 3641-3650. (18) Ripmeester, J. A.; Ratcliffe, C. I., 129Xe NMR Studies of Clathrate Hydrates: New Guests for Structure II and Structure H. J. Phys. Chem. 1990, 94, 8773-8776. (19) Davidson, D. W.; Gough, S. R.; Lee, F.; Ripmeester, J. A., Molecular-Reorientation in 1,3-Dioxane Hydrate and Compositions of Structure II Clathrates. Rev. Chim. Miner. 1977, 14, 447-463. (20) Subramanian, S.; Sloan, E. D., Trends in Vibrational Frequencies of Guests Trapped in Clathrate Hydrate Cages. J. Phys. Chem. B 2002, 106, 4348-4355. (21) Luzi, M.; Schicks, J. M.; Naumann, R.; Erzinger, J.; Udachin, K.; Moudrakowski, I.; Ripmeester, J. A.; Ludwig, R., Investigations on the Influence of Guest Molecule Characteristics and the Presence of Multicomponent Gas Mixtures on Gas Hydrate Properties. Proceedings of the 6th International Conference on Gas Hydrates, Vancouver, Canada, July 6-10, 2008.

ACS Paragon Plus Environment

24

Page 25 of 27

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

(22) Takeya, S.; Fujihisa, H.; Hachikubo, A.; Sakagami, H.; Gotoh, Y., Distribution of Butane in the Host Water Cage of Structure II Clathrate Hydrates. Chem. Eur. J. 2014, 20, 17207-17213. (23) Ohgaki, K.; Nakatsuji, K.; Takeya, K.; Tani, A.; Sugahara, T., Hydrogen Transfer from Guest Molecule to Radical in Adjacent Hydrate Cages. Phys. Chem. Chem. Phys. 2008, 10, 8082. (24) Kobayashi, N.; Minami, T.; Tani, A.; Nakagoshi, M.; Sugahara, T.; Takeya, K.; Ohgaki, K., Intermolecular Hydrogen Transfer in Isobutane Hydrate. Energies 2012, 5, 1705-1712. (25) Noack, K.; Jones, R. N., Conformational Equilibria in Open-Chain α, β-Unsaturated Ketones. Can. J. Chem. 1961, 39, 2225-2235. (26) Durig, J. R.; Little, T. S., Conformational Barriers to Internal-Rotation and Vibrational Assignment of Methyl Vinyl Ketone. J. Chem. Phys. 1981, 75, 3660-3668. (27) Bowles, A. J.; George, W. O.; Cunliffe.Db, Conformations of some αβ-unsaturated carbonyl compounds. Part II. Infrared and Raman spectra of methyl and ethyl acrylates and trans-crotonates. J. Chem. Soc. B 1970, 1070-1075. (28) Stoll, S.; Schweiger, A., Easyspin, a Comprehensive Software Package for Spectral Simulation and Analysis in EPR. J. Magn. Reson. 2006, 178, 42-55. (29) Bochevarov, A. D.; Harder, E.; Hughes, T. F.; Greenwood, J. R.; Braden, D. A.; Philipp, D. M.; Rinaldo, D.; Halls, M. D.; Zhang, J.; Friesner, R. A., Jaguar: A High-Performance Quantum Chemistry Software Program with Strengths in Life and Materials Sciences. Int. J. Quantum Chem. 2013, 113, 2110-2142.

ACS Paragon Plus Environment

25

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 26 of 27

(30) Lee, C. T.; Yang, W. T.; Parr, R. G., Development of the Colle-Salvetti CorrelationEnergy Formula into a Functional of the Electron-Density. Phys. Rev. B 1988, 37, 785-789. (31) Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (32) Rodriguezcarvajal, J., Recent Advances in Magnetic-Structure Determination by Neutron Powder Diffraction. Physica B 1993, 192, 55-69. (33) Cha, M.; Shin, K.; Lee, H., Spectroscopic Identification of Amyl Alcohol Hydrates through Free OH Observation. J. Phys. Chem. B 2009, 113, 10562-10565. (34) Cha, M.; Kwon, M.; Youn, Y.; Shin, K.; Lee, H., Phase Equilibria and Spectroscopic Identification of (2-Methylpropane-2-Peroxol + Gaseous Guests) Hydrates. J. Chem. Eng. Data 2012, 57, 1128-1133. (35) Nyquist, R. A., Interpreting Infrared, Raman, and Nuclear Magnetic Resonance Spectra. Academic Press: San Diego 2001, 1, 292-293. (36) Bunce, N. J., Introduction to the Interpretation of Electron-Spin Resonance-Spectra of Organic Radicals. J. Chem. Educ. 1987, 64, 907-914. (37) Elson, I. H.; Kemp, T. J.; Greatore.D; Jenkins, H. D. B., Electron spin resonance studies of reduction by solvated electrons in liquid ammonia. Part 6.-α,β-Unsaturated ketones. J. Chem. Soc., Faraday Trans. 2 1973, 69, 665-683.

ACS Paragon Plus Environment

26

Page 27 of 27

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

Table of Contents Graphic

ACS Paragon Plus Environment

27