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,