A swift introduction to thermodynamic modelling (Part II)

¡Hola!, this site will be dedicated to the last two days of the course and a place to download files.

Course logo

The last challenge

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!

NEWS!!! (1)

Here's the data that we managed to collect yesterday using just density measurements and X-ray diffraction data (cell paramters and Z number)

NEWS!! (2)

This is the calorimetric data at high temperatures.

NEWS!! (3)

This is the X-ray volume data at high temperatures and high pressures.

NEWS!! (4)

Low temperature calorimetric data for vicentite.

NEWS!! (5)

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 g(P,T,mineral) that is:

(1)g(P,T)=g(Pr,Tr)+s(Pr,Tr)(T+Tr)+c1(TTr)+c22(T2Tr2)c3(1Tr1T)+2c5(TTr)T(c1ln(TTr)+c2(TTr)+c32(1Tr21T2)+2c5(1Tr1T))13(v0((b6+b7(TTr)(4P4P+b7(TTr)+1)3/41)(b1(TTr)+2b5(TTr)+1)

 

SOLUTIONS

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".

Resources

Other unrelated resources