Stable Isotope Geochemistry: Theory

Lecture 35

Fractionations

• • • Isotope fractionation can originate from both kinetic effects and equilibrium effects. Equilibrium Fractionations: o Quantum mechanics predicts that the mass of an atom affects its vibrational motion, and therefore the strength of chemical bonds. It also affects rotational and translational motions. From an understanding of these effects of atomic mass, it is possible to predict the small differences in the chemical properties of isotopes quite accurately. Kinetic Fractionations o Lighter isotopes form weaker bonds and therefore react faster. They also diffuse more rapidly. These also lead to isotopic differences.

Equilibrium Fractionations

• • • • •

Equilibrium fractionations arise from translational, rotational and vibrational motions of molecules in gases and liquids and atoms in crystals.

Isotopes will be distributed so as to minimize the vibrational, rotational, and translational energy of a system. All these motions are quantized (but quantum steps in translation are very small).

Of the three types of energies, vibrational energy makes by far the most important contribution to isotopic fractionation. Vibrational motion is the only mode of motion available to atoms in a solid. These effects are small. For example, the equilibrium constant for the reaction: ½C 16 O 2 + H 2 18 O ⇄ ½C 18 O 2 + H 2 16 O is only about 1.04 at 25°C and the ∆G of the reaction, given by -RT ln K, is only -100 J/mol

• • • •

Predicting Fractionations Boltzmann distribution law states that the probability of a

molecule having internal energy E

i

is:

P i

= å

g i e

-

E i

/k

T g n e

-

E i

/k

T

o

n

g is a weighting factor to account for degenerate states.

The denominator is the

partition function

:

Q

= å

n g n e

-

E i

/k

T

From partition functions we can calculate the equilibrium constant: K = Õ

i Q i

n

i

We can divide the partition function into three parts:

Q

total

= Q

trans

Q rot Q vib o The separate partition functions can be calculated separately.

• • • •

Partition Functions

The following applies to diatomic molecules, for which the sums in the partition function fortunately have simple solutions. Principles are the same for multiatomic molecules and crystals, but the equations are more complex because there are many possible vibrations and rotations. Vibrational Partition Function: Note error in book: o

Q vib

= 1

e

-

h e

n /2 k

h

n

T

/k

T

remove 2 in exponential term in denominator where h is planks constant (converting frequency to energy) and ν is the vibrational frequency of the bond.

Rotational Partition Function

Q rot

= 8 p 2

I

k

T

s

h

2 o where I is the moment of inertia I = µr 2 ; i.e., the reduced mass (µ) of atoms times bond length.σ is a symmetry factor; σ= 1 for a non-symmetric molecule ( 18 O 16 O) and 2 for a symmetric one ( 16 O 16 O) Translational partition function for molecule of mass m is (derived from Schrödinger’s equation for particle in a box):

Q trans

= (2 p

m

k

T

) 3/2

V h

3

• • • • •

Partition Function Ratios

Since K = Õ

Q i

n

i

what we really want is the ratios of partition functions for isotope exchange reactions. Most terms cancel (including bond length).

This ratio for two isotopic species (isotopologues) of the same diatomic molecules, e.g., 16 O 18 O and 16 O 16 O, will be:

Q

16

O

18

O Q

16

O

2 =

h

n 16

O

18

O

s 16

O

18

O kT h

n 16

O

2 s 16

O

2

kT e

-

h

n 16

O

18

O

/2

kT

1 -

e

-

h

n 16

O

18

O

/

kT e

-

h

n 16

O

2 /2

kT

1 -

e

-

h

n 16

O

2 /

kT µ

16

O

18

O m

3/2 16

O

18

O µ

16

O

2

m

3/2 16

O

2 = n 16

O

18

O

n 16

O

2

e

-

h

n 16

O

18

O

/2

kT

(1 -

e

-

h

n 16

O

2 /2

kT

(1 -

e

-

h

n

H

2 16

O

/

kT

)

µ

16

O

18

O

s

e

-

h

n 16

O

18

O

/

kT

)

µ

16

O

2 s 16

O

2 3/2

m

16

O

18

O

16

O

18

O

3/2

m

16

O

2 We see the partition function ratio is temperature dependent (which arises only from the vibrational contribution: temperature canceled in other modes).

We also see that we can predict fractionations from measured vibrational frequencies and atomic and molecular masses.

Temperature Dependence

• • • • • • The temperature dependence is: At low-T (~surface T and below), the exponential term is small and the denominator approximates to 1. Hence

Q vib Q vib

»

e

=

h

n

e

-

h

n /2 k

T

1 -

e

-

h

n /k

T

/2k

T

= 1 -

h

n expressed as 2k

T

a =

A

+

B T

At higher temperature, however, this approximation no longer holds and α varies with the inverse square of T: a =

A

+

B T

2 The temperature dependence leads to important applications in geothermometry & paleoclimatology O isotope fractionation between CO 2 and H 2 O

• • • •

Kinetic Fractionations: Reaction Rates

Looking again at the hydrogen molecular bond, we see it takes less energy to break if it is H-H rather than D-H.

This effectively means the activation energy is lower and the rate constant, k, will be higher: DH will react faster than H 2 .

We can calculate a kinetic fractionation factor from the ratio of rate constants: a =

k k D H

=

e

( e 1/2

h

n

D

)/k

T e

( e 1/2

h

n

H

)/k

T

This will make no difference if the reaction goes to completion, but will make a difference for incomplete reactions. (good example is photosynthesis, which does not convert all CO organic carbon).

2 to

•

Kinetic Fractionation: Diffusion Lighter isotopic species will diffuse more rapidly.

o Energy is equally partitioned in a gas (or liquid). The translational kinetic energy is simply E = ½mv 2 .

o o o Consider two molecules of carbon dioxide, 12 C 16 O 2 and 13 C 16 O 2 , in a gas. If their energies are equal, the ratio of their velocities is (45/44) 1/2 , or 1.011. Thus 12 C 16 O 2 can diffuse 1.1% further in a given amount of time at a given temperature than 13 C 16 O 2 . But this applies to ideal gases (i.e., low pressures where collisions between molecules are infrequent. For the case of air, where molecular collisions are important, the ratio of the diffusion coefficients of the two CO 2 species is the ratio of the square roots of the reduced masses of CO 2 and air (mean molecular weight 28.8):

D

12

CO D

13

CO

2 2 =

µ

12

CO

2 -

air µ

13

CO

2 -

air

= 4.1906

4.1721

= 1.0044

o leading to a 4.4‰ fractionation (actually observed).

• • • • •

Rayleigh Distillation/Condensation

Different

isotopologues

of water evaporate at different rates and have different condensation temperatures.

We can imagine two ways in which condensation occurs: o droplets remain in isotopic

equilibrium

with vapor o droplets do not remain in equilibrium:

fractional

condensation If the fractionation between vapor and liquid is α, for fractional condensation, the fractionation, ∆, varies with fraction of vapor remaining, ƒ , as: ∆ = 1000(ƒ α-1 -1) For equilibrium condensation it is: ∆ = æ ææ 1 (1 1

f

) / a + ƒ æ ææ æ 1000 Fractional condensation can lead to quite extreme compositions of remaining vapor.

Isotopic composition of vapor when the fraction of original vapor, ƒ, remains.

• • • •

Isotope Fractionations

As a rule, heavy isotopes partition preferentially into phases in which they are most strongly bound (because this results in the greatest reduction in system energy).

o Covalent bonds, and bonds to heavier atoms, are generally strongest and hence will most often incorporate the heavy isotope.

Largest fractionations will occur where the atomic environment or bond energy differences are greatest o So, for example, fractionation of O between silicates are not large, because the O is mainly bound to Si.

Fractionations tend to be large between different oxidation states of an element (e.g., for C, N, S).

Lighter isotopes are likely to be enriched in the products of incomplete reactions and also reactions where diffusion is important.

• • • • •

Mass Dependent Fractionation

if a 4‰ fractionation of δ 18 O is observed in a particular sample, what value of δ 17 O do we predict? We might guess it would ½ as much. Mass occurs in a variety of ways in the partition function, as m ratio of fractionation of 17 O/ 16 O to that of 18 O/ 16 3/2 , as reduced mass, and in the exponential term. Consequently, the O in most cases is about 0.52. Nevertheless, the fractionation between isotopes predicted by this equation is proportional to the difference in mass – this is referred to as mass-dependent fractionation. There are some exceptions where the ratio of fractionation of 17 O/ 16 O to that of 18 O/ 16 O is ≈1. Since the extent of fractionation in these cases seems independent of the mass difference, this is called mass-independent fractionation. Examples o Oxygen in meteorites o o Sulfur in Archean sulfides Oxygen in stratospheric gases

• • • • • •

Mass Independent Fractionation

Most examples of ‘MIF’ seem to be related to photochemical reactions.

Formation of ozone in the stratosphere involves the energetic collision of monatomic and molecular oxygen: O + O 2 → O 3 The ozone molecule is in a vibrationally excited state and subject to dissociation if it cannot lose this excess energy. The excess vibrational energy can be lost either by collisions with other molecules, or by partitioning to rotational energy. In the stratosphere, collisions are infrequent, hence repartitioning of vibrational energy represents an important pathway to stability.

Because there are more possible energy transitions for asymmetric species such as 16 O 16 O 18 O and 16 O 16 O 17 O than symmetric ones such as 16 O 16 a stable molecule.

O 16 O, the former can more readily repartition its excess energy and form In the troposphere, collisions more frequent, reducing this effect.

Isotope Geothermometry

• • One of the principal uses of stable isotopes is geothermometry. Stable isotope geothermometers are based on the temperature dependence of the fractionation factor or equilibrium constant, which can generally be expressed as: a =

A

+

B T

2 o (at low temperatures, the form of changes to α ∝

1/T).

Temperature dependence can be theoretically calculated or experimentally measured. Measuring the isotopic composition of two phases allows us to calculate the temperature at which they equilibrated (assuming, of course, that they did equilibrate).

• • • •

Isotope “Clumping”

Consider the distribution of 9.1). CO and O molecules or “isotopologues”, such as 16 O 17 O, etc. 2 18 O between CO and O 12 C 16 O, 12 C 17 2 (Example O, 13 C 18 O, The distribution of isotopes between these species will not be random but rather some of these isotopologues will be thermodynamically favored. o Essentially, grouping the heavy isotopes in one molecule, e.g., 13 C 18 O, reduces bond energy by a bit more than twice the reduction of putting one heavy isotope in the molecule. Thus “clumping” of heavy isotopes reduces system energy.

This ‘clumping’ depends on temperature (greater at low T). By analyzing the various isotopologues of the species, one can calculate equilibrium temperatures.

o The advantage is that we need analyze just one phase involved in the reaction, for example, carbonate precipitated from water.

o In addition to calculating temperature, one can also calculate the isotopic composition of the water.

This is a very new field, but holds great promise in isotope geothermometry.