Hydrogen Isotope Fractionation in the Epidote–Hydrogen and Epidote

With the objective of development of a theoretical methodology for calculation of ... hydrogen isotope fractionation in the epidote–water and epidot...
0 downloads 0 Views 494KB Size
Subscriber access provided by University of South Dakota

Letter

Hydrogen Isotope Fractionation in the Epidote–Hydrogen and Epidote–Water Systems: Theoretical Study and Implications Abu Md. Asaduzzaman, and Jibamitra GANGULY ACS Earth Space Chem., Just Accepted Manuscript • DOI: 10.1021/ acsearthspacechem.8b00076 • Publication Date (Web): 06 Sep 2018 Downloaded from http://pubs.acs.org on September 12, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

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

ACS Earth and Space Chemistry

1

1

Hydrogen Isotope Fractionation in the Epidote–Hydrogen and Epidote–Water Systems:

2

Theoretical Study and Implications

3

Abu Asaduzzaman1,2 and Jibamitra Ganguly*,3

4 5

1

6

Middletown, PA 17057, USA

7

2

Materials Science and Engineering, University of Arizona, Tucson, AZ 85721, USA

8

3

Department of Geosciences, University of Arizona, Tucson, AZ 85721, USA

9

*Corresponding author

School of Science, Engineering and Technology, Pennsylvania State University - Harrisburg,

10 11 12

Key words: Hydrogen Isotope, Density Functional Theory, Partition Function, Epidote, Equilibrium Constant

13

ACS Paragon Plus Environment

ACS Earth and Space 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 17

2

14

ABSTRACT Experimental data for hydrogen isotope and in general for stable isotope

15

fractionation between mineral and water often tend to be widely divergent primarily because of

16

solution and re-precipitation of finely ground mineral grains in the experimental charges. With

17

the objective of development of a theoretical methodology for calculation of stable isotope

18

fractionation that may be suitably integrated with experimental data to produce better constraints

19

on the fractionation behavior, we have calculated hydrogen isotope fractionation in the epidote-

20

water and epidote-hydrogen systems for which reliable experimental data are available. The

21

calculations were carried out using a combination of classical and statistical thermodynamics and

22

quantum chemical methods. The calculated fractionation factors vs. temperature in both systems

23

have been found to be in very good agreement with the experimental data and provide a

24

theoretical framework for extrapolation of experimental data as a function of temperature.

25

1. INTRODUCTION

26

Quantitative data about deuterium and hydrogen isotope fractionation between hydrous minerals

27

and their fluid environments as function of temperature and pressure permit indirect

28

determination of the D/H ratio of the fluid phase from the measured isotopic ratios of the

29

hydrous minerals that have important implications for a variety of geological processes. Thus,

30

there has been a sustained effort in the geological community to experimentally determine the

31

hydrogen isotopic composition of hydrous minerals in equilibrium with water. Additionally,

32

knowledge of the hydrogen isotopic fractionation between the solar nebula and nominally

33

anhydrous minerals, such as olivine, is of critical importance to understand if the nebular water

34

was responsible for a significant component of the Earth’s water via adsorption mediated

35

accretion on mineral surfaces during nebular condensation process1,2. The experimental studies

36

in the mineral-water systems have, however, suffered from several problems such as failure to

37

attain equilibrium fractionation, especially at relatively low temperatures, and solution-

38

precipitation that often resulted in wide divergence of the published data. This has led to revival

39

of efforts for theoretical calculations using statistical thermodynamics, following the groundwork

40

laid by Urey3 and subsequent developments4,5,6,7,8,9.

41

The primary purpose of this paper is to calculate hydrogen isotope fractionation between

42

epidote and water/molecular hydrogen using a combination of thermodynamics and quantum

43

chemical methodologies, and to compare the results with reliable experimental data that are

ACS Paragon Plus Environment

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

ACS Earth and Space Chemistry

3

44

available for these systems in order to assess the validity of the theoretical method. Additionally,

45

we suggest a method of fitting of the experimental lnα vs. 1/T2 relation for mineral-water

46

systems (where α is the bulk isotopic fractionation factor) that is better suited to extrapolation

47

than the conventional practice of linear regression and extrapolation. The demonstrated

48

agreement between the calculated and experimental data would also lay the groundwork for a

49

follow-up study on the hydrogen isotope fractionation between solar nebula and chemisorbed

50

water on olivine in the nebular environment.

51

2. THERMODYNAMICS OF HYDROGEN ISOTOPE FRACTIONATION IN

52

MINERAL-WATER/H2 SYSTEMS

53

Let us first consider the general case of hydrogen isotope fractionation between a crystal (x’l)

54

and water. The problem may be treated in terms of an exchange reaction as follows.

55

x’l(OH) + ½ D2O = x’l(OD) + ½ H2O

56

for which the equilibrium constant is

57

K =





(a)

(1)

/  

    

58

where the parenthetical terms indicate activities of the specified components. From the available

59

thermodynamic data on the excess Gibbs free energy of mixing in the CH4-CD4 system10, which

60

is found to be symmetrical to composition, the activity coefficient of either species at 50%

61

composition is found to be 1.0004 at 120 K. Thus, the deviation from ideality of mixing of the

62

isotopic species in a given phase should be expected to be negligible at temperatures of

63

geological and planetary interests, and hence the activity of a species may be equated to its mole

64

fraction (using a standard state of pure species at P-T condition of interest).

65

66

The experimentally measured mineral/water bulk fractionation factor is given by α=











(2)

67

This is equivalent to Ka if D2O and H2O are the only species in water since in that case, a(D2O) =

68

(XD)2 and a(H2O)= (XH)2 for the case of ideal mixing, where a and X stand, respectively, for the

ACS Paragon Plus Environment

ACS Earth and Space 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 17

4

69

activity and mole fraction of the specified species. However, in general, this is not the case

70

because of the presence of other H-D-O species in the system.

71

In this work, we are interested in comparing the results of theoretical calculation with the

72

experimentally determined bulk isotopic fractionations. The experiments were carried out at ~2-

73

4 kb pressure in sealed gold capsules that are effectively impervious to hydrogen diffusion at

74

such pressures. If the molar H2 to O2 ratio within the gold capsules remains close to 2 (as would

75

be the case for a closed system in which H2 and O2 are derived only from the stoichiometric

76

dissociation of water), then H2 content of the vapor phase becomes vanishingly small11 and

77

hence may be easily neglected in the mass balance calculations. Thus, we need to consider only

78

three species in the H-D-O system, viz. H2O, D2O and HDO. These species are related by the

79

equilibrium ½ H2O + ½ D2O = HDO

80 81

for which the equilibrium constant is

82

K =

83



  

⁄

(b)

(3)

For the molar abundance of D and H in water we can write

84

D ≈ 2D O + HDO = 2D O + K  $H OD O%&/

(4a)

85

H ≈ 2H O + HDO = 2H O + K  $H OD O%&/

(4b)

86

(D, H, D2O and HDO in these equations represent the molar abundances of the respective

87

species) so that 

 

88

89

90

91

'()*

  &⁄

=   

φ

(5)

where φ

 / , -.  / =  / , -.  /

Thus, α, and Ka are related according to

ACS Paragon Plus Environment

(6)

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

ACS Earth and Space Chemistry

5  8

αx 0 l − water ≡ 

 

92

=

-9 :

(7)

93

It is easy to see from the last two equations that α equals Ka when Kb = 2. From available

94

thermodynamic data in the JANAF Tables12, Kb is found to increase from 1.9239 at 300 K to

95

1.9892 at 1200 K. Thus, since Kb/2 ~1 and (D2O) 0 for product,