Measurement of a Doubly Substituted Methane Isotopologue

Jun 4, 2014 - 13CH3D, by Tunable Infrared Laser Direct Absorption Spectroscopy ... continuous wave quantum cascade lasers; one tuned at 8.6 μm to...
2 downloads 0 Views 2MB Size
Subscriber access provided by Oregon State University

Article

Measurement of a Doubly-Substituted Methane Isotopologue, CHD, by Tunable Infrared Laser Direct Absorption Spectroscopy

13

3

Shuhei Ono, David T Wang, Danielle S. Gruen, Barbara Sherwood Lollar, Mark Zahniser, Barry J. McManus, and David D. Nelson Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/ac5010579 • Publication Date (Web): 04 Jun 2014 Downloaded from http://pubs.acs.org on June 10, 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.

Analytical Chemistry 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 37

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

Analytical Chemistry

1

Measurement of a Doubly-Substituted Methane

2

Isotopologue, 13CH3D, by Tunable Infrared Laser

3

Direct Absorption Spectroscopy

4

Shuhei Ono,*1, David T. Wang1, Danielle S. Gruen1, Barbara Sherwood Lollar2, Mark S.

5

Zahniser3, Barry J. McManus3, David D. Nelson3 1

6

Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of

7

Technology, Cambridge, MA. [email protected] 2

8 9

3

Department of Earth Sciences, University of Toronto, Toronto, ON, Canada

Center for Atmospheric and Environmental Chemistry, Aerodyne Research, Inc., Billerica,

10

Massachusetts, USA

11 12

Keywords: Tunable Laser Spectroscopy, Quantum Cascade Laser, Methane, Isotopologue,

13

Clumped Isotope

14

* Corresponding author, [email protected]

15

A manuscript prepared for submission to Analytical Chemistry March 18, 2014

16

5520 words, 7 figures (equivalent to ~1200 words), 2 tables (250 to 400 words)

17 -1ACS Paragon Plus Environment

Analytical 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 37

18 19

Abstract

20

Methane is an important energy resource and significant long-lived greenhouse gas. Carbon

21

and hydrogen isotope ratios have been used to better constrain the sources of methane, but

22

interpretations based on these two parameters alone can often be inconclusive. The precise

23

measurement of a doubly-substituted methane isotopologue, 13CH3D, is expected to add a critical

24

new dimension to source signatures by providing the apparent temperature at which methane was

25

formed or thermally equilibrated. We have developed a new method to precisely determine the

26

relative abundance of

27

(TILDAS). The TILDAS instrument houses two continuous wave quantum cascade lasers; one

28

tuned at 8.6 µm to measure

29

13

30

‰ (2ߪ) was achieved for the measurement of

31

mmol) pure methane samples. Smaller quantity samples (ca. 0.5 mL STP) can be measured at

32

lower precision.

33

thermally-equilibrating methane with a range of δD values.

34

corresponds to uncertainties of ±7°C at 25°C and ±20°C at 200°C for estimates of apparent

35

equilibrium temperatures.

36

determine 13CH3D in natural methane samples to distinguish geological and biological sources of

37

methane in the atmosphere, hydrosphere, and lithosphere.

13

CH3D by using tunable infrared laser direct absorption spectroscopy

13

CH3D,

12

CH3D and

12

CH4, and the other at 7.5 µm to measure

CH4. Using an astigmatic Herriott cell with an effective pathlength of 76 m, a precision of 0.2 13

CH3D abundance in ca. 10 mL STP (i.e., 0.42

The accuracy of the ∆13CH3D measurement is 0.7‰ (2ߪ), evaluated by The precision of ±0.2‰

The TILDAS instrument offers a simple and precise method to

-2ACS Paragon Plus Environment

Page 3 of 37

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

Analytical Chemistry

Introduction

38 39

Methane is the second most important long-lived greenhouse gas, and an increasingly

40

important energy resource.1 Abrupt changes in reservoirs and fluxes of methane in the Earth’s

41

surface may have been responsible for triggering abrupt climate change in the past2. The

42

atmospheric methane level increased through much of the twentieth century, but the growth rate

43

slowed and the concentration appeared to have plateaued between 1999 and 2007, at which point

44

growth resumed.3 The main drivers for this secular trend have been debated, with explanations

45

ranging from changes in source (e.g., wetlands vs. fossil fuel emissions) and/or sink (oxidation

46

by OH radical) strengths.3-5

47

The majority of environmental methane is of biogenic (originating from microbial

48

methanogenesis) or thermogenic (produced by thermal cracking of higher molecular weight

49

hydrocarbons, kerogen and coal) origin, with contributions from biomass burning and, in certain

50

geologic environments, putative abiotic sources (e.g., serpentinite and deep crustal fluids).6-8

51

Carbon (13C/12C) and hydrogen (D/H) isotope ratios of methane and associated short-

52

chain hydrocarbons have been used to elucidate their sources9,10 but such data alone can often be

53

inconclusive because of significant overlap in isotopic signatures associated with microbial,

54

thermogenic and abiogenic gases. This is because the

55

carbon source (e.g., CO2, acetate) and its isotopic composition,11 the isotope effects associated

56

with microbial methanogenesis or thermal cracking, and the effects of secondary processes, such

57

as oxidation and mixing.12,13 The D/H ratio of methane is also a complex function of reaction

58

pathways, the D/H ratio of environmental water, and the D/H ratio of precursor compounds.14,15

59

Therefore, the development of a new constraint, such as an isotopic thermometer, could unlock

60

critical information to constrain the source of methane.

13

C/12C ratio of methane depends on the

-3ACS Paragon Plus Environment

Analytical 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 37

61

In this study, we report a methane isotope thermometry method based on the

62

measurements of the doubly-substituted (“clumped”) methane isotopologue, 13CH3D. Following

63

earlier studies of doubly-substituted carbon dioxide (13C16O18O),19 the precise measurements of

64

four or more isotopologues of methane (12CH4,

65

estimation of the temperature at which a sample of methane was formed or thermally

66

equilibrated.20,21

13

CH4,

12

CH3D, and

13

CH3D) may allow

67

Carbon (13C/12C) and hydrogen (D/H) isotope ratios of methane have been routinely measured

68

by gas-source isotope-ratio mass spectrometry (IRMS). Conventional IRMS techniques for

69

measuring carbon- and hydrogen-isotope ratios involve combustion or pyrolysis of methane,

70

respectively, and measurement of the isotopologue ratios of the product CO2 or H2.9,22 Direct

71

measurement of

72

13

73

3 to 8 ppm) and interferences from adduct (13CH5+) and

74

CH3+, CH2+) also complicate mass-16 (12CH4+) and mass-17 (12CH3D+ and

75

Resolving these isobars requires the use of double-focusing high resolution gas-source

76

IRMS.21,23 Stolper et al. (2014)21 report the first precise measurements of methane clumped

77

isotopologue abundance (combined

78

(M/∆M) from 16,000 to 25,000; contributions of adduct ions (13CH5+ on 13CH3D+, and 12CH5+ on

79

12

80

13

CH3D+ ion is required, however, for mass-spectrometric determination of

CH3D. This is technically challenging because of the low fractional abundance of 13CH3D (ca.

13

12

CH2D2+ ions. Fragment ions (e.g., 13

CH4+) signals.

CH3D+12CH2D2) by using medium mass-resolving power

CH3D+) were corrected as a function of mass-16 ion current. In this study, we report the application of tunable infrared laser direct absorption spectroscopy 13

81

(TILDAS) for the fully-resolved measurement of a clumped isotopologue of methane,

CH3D.

82

TILDAS is a virtually non-destructive technique, and offers a promising approach for precise

83

and direct measurements of methane isotopologues.

Bergamaschi et al. (1994)24 reported

-4ACS Paragon Plus Environment

Page 5 of 37

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

Analytical Chemistry

13

measurements of

85

(ߥ3 band). Precisions of 0.44 and 5.1 ‰ were achieved for δ13C and δD measurements,

86

respectively (δ values are defined as the deviation of the isotope ratios

87

with respect to reference materials). Recently, a tunable laser spectroscopy instrument

88

with an interband cascade laser tuned at 3.27μm was used for the detection of methane on

89

the Mars Curiosity rover mission.25 Tsuji et al., (2012)26 reported the first optical measurement

90

of

91

analytical precision of 20‰ (1σ), however, was too large to resolve the thermodynamically-

92

predicted natural range of variation of 13CH3D abundances (ca. 6‰).

13

CH4 and

12

84

CH3D using a lead salt diode laser tuned in the ~3.3 µm regions

13C/12C

and D/H

CH3D using a difference-frequency-generation laser in the 3.4µm region. Their reported

93

The TILDAS instrument used in this study takes advantage of quantum cascade lasers

94

(QCLs) that allow for continuous wave (cw) operation near room temperature, with high power

95

(tens of mW) and narrow linewidths (1,000K), and referred to

110

as a stochastic distribution.

111

predicts that the equilibrium for (1) lies slightly toward the right (K≃1.006 at room temperature,

112

see Supporting Information, SI).

113

At lower temperatures, however, statistical mechanics theory

Defining ∆13CH3D as the natural logarithm for equation (2) gives: ଵଷ

Δ CHଷ D = ln(‫ = )ܭ‬ln

ൣ 13CH3 D൧ ൣ 12CH3 D൧

− ln

ൣ 13CH4 ൧ ൣ 12CH4 ൧

(3)

114

The definition of ∆13CH3D in equation (3) yields values that are practically identical to those

115

calculated using the definition used by Stolper et al. (2014)21 (i.e., the deviation from stochastic

116

distribution. see SI).

117 118

In this study, we use conventional delta notation to express ratios of isotopologue abundance in a sample with respect to those of a reference. For example, ߜ 13CH3 D =

119 120

(ሾ 13CH3 Dሿ/ሾ 12CH4 ሿ)sample −1 (ሾ 13CH3 Dሿ/ሾ 12CH4 ሿ)reference

(4)

For practical purposes (under expected ranges of multiply-substituted isotopologue abundance for natural samples), δ13CH4 is identical to δ13C and δ12CH3D is identical to δD. -6ACS Paragon Plus Environment

Page 7 of 37

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

Analytical Chemistry

Experimental Section

121

122

Synthesis of Methane Isotopologues

123

Methane isotopologues,

13

13

CH3D, CH3D (natural abundance

C/12C), and CH4 (D-depleted)

124

were synthesized for spectral characterization and to produce a series of working standard gases.

125

The clumped methane isotopologue,

126

13

127

iodomethane and D2O. D-depleted methane was synthesized from aluminum carbide (Al4C3,

128

natural C isotope abundance) and D-depleted water (D content 2-3 ppm) (See SI for detailed

129

description of synthesis procedures).

13

CH3D, was synthesized by the Grignard reaction from

C-iodomethane and D2O. Similarly, CH3D (with natural-abundance 13C) was synthesized from

130 131 132

Selection of Absorption Lines Among several possible methane infrared band systems, a spectral region around 8.6 µm was Tsuji et al. (2012)26 reported measurements of

13

133

selected.

134

However, they used dry ice to cool the absorption cell in order to reduce hot bands and high-J

135

fundamental bands from major isotopologues. The absorption band system at 8.6 µm is highly

136

advantageous because of less interference from hot bands in this region, and the availability of

137

line parameters in the HITRAN database.29 A single deuterium substitution in CH4 reduces its

138

symmetry from Td to C3v. As a result, the triply degenerate C-H bending mode of CH4 centered

139

at 1300 cm-1 splits into C-D bending of 13CH3D and 12CH3D. The resultant C-D bending mode is

140

about 150 cm-1 lower in frequency than the C-H bending vibration.

141

frequency shift by deuterium substitution allows us to measure 13CH3D and 12CH3D lines that are

142

free from strong lines originating from the much more abundant 12CH4 (Figure 1). -7ACS Paragon Plus Environment

CH3D in the 3.4 µm region.

This relatively large

Analytical 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 37

The spectral window between 1168.885 and 1168.995 cm-1 was chosen to access three

143

13

CH3D,

12

CH3D and

12

144

isotopologue lines of

CH4 (Figure 1 and 2). These three lines are

145

virtually-free from spectral interference, and have relatively low lower-state energies (59, 470,

146

and 861 cm-1 for 13CH3D, 12CH3D, and 12CH4, respectively), such that temperature-dependence is

147

expected to be minor (