Thermodynamics the branch of science that deals with

Thermodynamics • “the branch of science that deals with energy levels and the transfer of energy between systems and between different states of matter”

What is Energy? ? ? • “It is important to realize that in physics today, we have no knowledge of what energy is. We do not have a picture that energy comes in little blobs of a definite amount. It is not that way. ” –Richard Feynman • HOWEVER Feynman goes on to elaborate that energy has meaning as a way to define, and quantify, changes which bring about reactions – between systems, energy levels, or states of matter. • “How seriously must we take the physical existence of this energy? No more and no less than any other bookkeeping practices. ” –Richard Feynman

Thermodynamics • Thermodynamics answers the following question: • For any reaction - defined by a set of reactants and products set in exactly defined conditions (temperature, pressure, concentration, etc. ) will that reaction go forward spontaneously or not? ? • For ANY geochemical reaction, if thermo says NO, rest assured the reaction will not proceed. • Thermo says NOTHING about the speed a reaction occurs!!

Equilibrium • Anything at equilibrium is theoretically undergoing forward and reverse reactions: • A+B↔C – A + B C AND C A +B • Equilibrium has 2 criteria: – Reaction does not appreciably change in time – Perturbation of that equilibrium will result in a return to the equilibrium

Equations can be ‘added’ together, equilibrium constants also get ‘added’ together! Convenient way to rewrite reactions (to look at more appropriate reactions or to use things you’ve more directly measured…. ) log Keq • Ca. CO 3(calcite) = Ca 2+ + CO 32 - -8. 48 • CO 2(g) + H 2 O = H 2 CO 30 -1. 47 • H 2 CO 30 = H+ + HCO 3 - -6. 35 • H+ + CO 32 - = HCO 3 - +10. 33 Ca. CO 3(calcite) + CO 2(g) + H 2 O = Ca 2+ + 2 HCO 3 - -5. 97

Assessing equilibrium Q reaction quotient, aka Ion Activity Product (IAP) is calculated from knowing activity of all components of a reaction K aka Keq, we get from thermodynamic data – it is one number defined AT EQUILIBRIUM Equilibrium for any reaction is when Q = K

Where do K’s come from? • Measure directly – experimental determination of conditions at equilibrium • Use thermodynamic data – K is directly related to free energy of reaction – DGR

Hydroxylapatite • Ca 5(PO 4)3(OH) = OH- + 3 PO 43 - + 5 Ca 2+ • Log K = -59. 0351 at 25ºC

Aqueous Complexes • • Combinations of ions to form dimers, trimers, etc. , are complexes Why do we care? ? 1. Complexation of an ion also occurring in a mineral increases solubility 2. Some elements occur as complexes more commonly than as free ions 3. Adsorption of elements greatly determined by the complex it resides in 4. Toxicity/ bioavailability of elements depends on the complexation

How do we know about all those species? ? • Based on complexation how any ion interacts with another ion to form a molecule, or complex (many of these are still in solution) • Yet we do not measure how much Ca. NO 3+, Ca. F+, or Ca. PO 4 - there is in a particular water sample • We measure Ca 2+ But is that Ca 2+ really how the Ca exists in a water? ?

Defining Complexes • Use equilibrium expressions: • c. C + l. HL CL + l. H+ • Where B is just like Keq!

Equilibrium • Equilibrium Constant, K (or Keq) describes conditions AT equilibrium (where DGR=0) DGR -DG 0 R = RTln. K DG 0 R = -RTln. K DG 0 R= SDG 0 R products – SDG 0 R reactants DGR= SDGR products – SDGR reactants

DG 0 • Energy at STANDARD STATE • = 25°C, 1 bar Pressure, 1 molal concentration for each

Speciation • Any element exists in a solution, solid, or gas as 1 to n ions, molecules, or solids • Example: Ca 2+ can exist in solution as: Ca++ Ca(H 3 Si. O 4)2 Ca(O-phth) Ca. B(OH)4+ Ca. CH 3 COO+ Ca. CO 30 Ca. Cl+ Ca. F+ Ca. H 2 Si. O 4 Ca. H 3 Si. O 4+ Ca. HCO 3+ Ca. NO 3+ Ca. OH+ Ca. PO 4 Ca. SO 4 Ca. HPO 40 • Plus more species gases and minerals!!

Mass Action & Mass Balance • m. Ca 2+=m. Ca 2++MCa. Cl+ + m. Ca. Cl 20 + Ca. CL 3 - + Ca. HCO 3+ + Ca. CO 30 + Ca. F+ + Ca. SO 40 + Ca. HSO 4+ + Ca. OH+ +… • Final equation to solve the problem sees the mass action for each complex substituted into the mass balance equation

Mineral dissolution/precipitation • To determine whether or not a water is saturated with hydroxyapatite, we could write a dissolution reaction such as: Ca 5(PO 4)3(OH) = OH- + 3 PO 43 - + 5 Ca 2+ • We could then determine the equilibrium constant: • If K = -59. 04, can determine how much Ca 2+ and PO 43 might dissolve at any p. H and T

Activity • Sometimes called ‘effective concentration’, which is misleading and reflects a poor understanding of the property… • Think of more of the effect the rest of a solution has on how easily two ions come together. .

Activity • For solids or liquid solutions: ai=Xigi • For gases: ai=Pigi = fi Xi=mole fraction of component i Pi = partial pressure of component i mi = molal concentration of component i • For aqueous solutions: ai=migi

Activity Coefficients • Where do they come from? ? • The standard state for dissolved ions is actually an infinitely dilute solution… • Activity of phases - gases, minerals, and bulk liquids (H 2 O) are usually pretty close to 1 in waters • Dissolved molecules/ ions have activity coefficients that change with concentration (ions are curved lines relating concentration and activity coefficients, molecules usually more linear relation)

Application to ions in solution • Ions in solutions are obviously nonideal mixtures! a i = g i mi • The activity coefficient, gi, is found via some empirical foundations • Dependent on the other ions in water…

Dissolved species gi • First must define the ionic strength (I) of the solution the ion is in: Where mi is the molar concentration of species i and zi is the charge of species I

Activity Coefficients • Debye-Huckel approximation (valid for I: • Where A and B are constants (depending on T), z is charge, I is ionic strength, and å is a measure of the effective diameter of the ion

Different ways to calculate gi • • Limiting law Debye-Huckel Davies TJ, SIT models • Pitzer, HKW models

Neutral species • Setchnow equation: • Logan=ks. I For activity coefficient (see table 4 -2 for selected coefficients)

Mass Action & Mass Balance • m. Ca 2+=m. Ca 2++MCa. Cl+ + m. Ca. Cl 20 + Ca. CL 3 - + Ca. HCO 3+ + Ca. CO 30 + Ca. F+ + Ca. SO 40 + Ca. HSO 4+ + Ca. OH+ +… • Final equation to solve the problem sees the mass action for each complex substituted into the mass balance equation • Equations for each ion – iterative solution…

Speciation Models • PHREEQC, or Web. PHREEQ, is a USGS program that solves, simultaneously (iteratively really), all of the mass action and mass balance equations for a water’s chemical composition