¡Hola!, this site will be dedicated to the last two days of the course and a place to download files.
You've learned quite a lot of things this week so I propose the following game.
After a long field campaign in the Antarctic, a bunch of petrologists from our PGG group brought a collection of 10 unusual samples (samples PGG23-01 to PGG23-10 ) from two sites but with no contextual information (i.e., "xenoliths"). Back in Granada (Armilla), they discovered that most of the samples contained significant amounts of an unknown element that they temporarily called X.
Inspired by the discovery, members of the team wrote a proposal and after a year they got funded by a generous grant to constrain one of the two thermal gradients where the strategic element can be found for future explorations and (of course) revolutionize our current knowledge of the thermal state of the lithosphere.
You will start by looking at the petrological data (bulk rock chemistry, mineralogy and mineral chemistry of the 10 samples).
Site 1 Samples from PGG23-01 to PGG23-05.
Site 2 Samples from PGG23-06 to PGG23-10.
XRF data and EPMA data. (XRF and EPMA with Numbers format for Macs.
No more information is given (but as you know the PGG is excellently equipped with all experimental techniques necessary to accomplish the project). Your challenge is to discuss what kind of experimental data you would need to solve the puzzle.
Good luck!
Here's the data that we managed to collect yesterday using just density measurements and X-ray diffraction data (cell paramters and Z number)
vic EoS = 2 |
X2O3(1)
G0 = -1560000 S0 = 49.9 V0 = 2.658
c1 = 160.0 c2 = .600E-2 c3 = -2560600 c5 = -599.2
b1 = .429E-4 b5 = -.429E-3 b6 = 2620000. b7 = -388 b8 = 4
end
menz EoS = 2
MGO(3)X2O3(1)SIO2(3)
G0 = -5900105 S0 = 280.3 V0 = 8.318
c1 = 600.0 c3 = -5100000 c5 = -4000.2
b1 = .2500E-4 b5 = -.2500E-3 b6 = 2007000. b7 = -300.55 b8 = 4
end
rom EoS = 2
MGO(1)X2O3(1)
G0 = -2106030 S0 = 81.5 V0 = 3.978
c1 = 342.7 c2 = -.5037E-2 c3 = -2205100 c5 = -1678.1
b1 = .431E-4 b5 = -.431E-3 b6 = 1945000. b7 = -291.75 b8 = 4
end
This is the calorimetric data at high temperatures.
This is the X-ray volume data at high temperatures and high pressures.
Low temperature calorimetric data for vicentite.
We didn't really have time to go through the actual computation of Gibbs but if you want to see how to compute it you can have a look at this Jupyter notebook, at the end you will see the algebraic expression for computing the apparent Gibbs energy for any pressure and temperature
Here's a quickly gathered Jupyter notebook to go through the individual solutions of the exercise (but not the final answer), for that, you need more experimental data to constrain the Margules parameter of solid solution for vicentite (you can see those experiments in the GitHub folder "Solid solution".
Here is a summary of the key expression needed to understand what is behind the scenes of a phase diagram (click here).
A few notes about the simplest solid solutions (ideal and regular). (Clic here)
An example of olivine as a simple non-ideal binary solid solution implemented in Python (Clic here)
Oxygen fugacity based on Fe-Mg exchange in olivine-orthopyroxene in magnetite buffered assemblages (Clic here).