Empirical and Ab Initio Thermodynamic Models of Minerals and Melts
This is the home page for the five day MELTS/FPMD training school to be held on the island of Milos; for the two to three day workshops held at Caltech, and elsewhere, please see the alphaMELTS workshops page; for MELTS short courses held at Goldschmidt and other conferences see the Goldschmidt 2020 page.
Early career scientists are invited to attend a training school on thermodynamic modelling of natural silicate liquids and their equilibria with minerals and fluids, from atomistic to macroscopic scales. The methods to be taught are the foundations of modern attempts to systemize igneous/metamorphic petrology and the Earth’s deep volatile (water and carbon) cycles and to make complex natural phenomena computable. Few young scientists today are trained in both approaches, but the intent of this course is to bring them together in order to highlight the strengths and weaknesses of each, the rapid advances being made in both areas, and the ways that convergence of the two approaches can lead to robust predictive understanding of natural magmatic phenomena. Students and young scientists from across the EU, and beyond, will learn the theory underlying both the macroscopic and atomistic approaches and be trained in practical use of both sets of modelling tools.
The macroscopic or classical portion will introduce the fundamental thermodynamics of multicomponent phase equilibria and the ideas behind the design of solution models before outlining the specific details of the MELTS family of models and the associated software packages (Ghiorso and Sack, 1995; Asimow and Ghiorso, 1998; Smith and Asimow, 2005; Antoshechkina and Ghiorso, 2018). The current incarnations of this approach — pMELTS (Ghiorso et al., 2002), rhyolite-MELTS (Gualda et al., 2012), and the H2O-CO2-bearing extension of rhyolite-MELTS (Ghiorso and Gualda, 2015) — enable calculation from subsolidus to superliquidus conditions from ambient pressure up to about 3 GPa in systems from ultramafic to highly felsic with mixed volatiles and a wide range of oxidation states. The lecture portion will be followed by practical exercises demonstrating and teaching the use of the software, including new capabilities to fully and seamlessly integrate MELTS calculations into MATLAB, Python, and other programming environments, enabling for the first time a range of large-scale calculations and code integration.
The atomistic portion will centre around the application of first-principles molecular dynamics (FPMD) to understanding static, dynamic, thermodynamic, and transport properties of multicomponent silicate liquids. In MD the atoms move according to Newtonian dynamics under the action of interatomic forces. In FPMD the ionic forces are computed directly from the electronic structure of the system, which is obtained by solving an approximate form of the Schrödinger equation. We will present different ensembles, statistical analysis of the MD runs, and applications in geophysics. Again, the lecture series will be followed by practical tutorials. We will employ the UMD package, developed at Lyon (Caracas, R., Kobsch, A., Solomatova, N. V., Li, Z., Soubiran F. & Hernandez, J.-A., UMD – an open-source python-based package to analyze ab initio molecular dynamics simulations, submitted.) to perform post-processing tasks. We will exemplify with obtaining a series of structural, transport and thermodynamic properties of various silicate melts.
We will complement the theoretical part with lectures on experimental apparatus and in situ measurements. Finally, we will bring the two approaches together with a final practical exercise in which the heat capacity of mixing is evaluated both by application of the MELTS model and from a set of FPMD calculations.
The novelty of the proposed event, which builds upon a very successful recent series of NSF-funded MELTS workshops in the USA and previous workshops sponsored by CECAM and EGU in Europe, is to combine and showcase complementary progress in both empirical and ab initio thermodynamic modelling. The two approaches have only recently begun to capture the complex behavior of multicomponent fluid solutions relevant to magmatic phenomena at realistic and experimentally challenging physico-chemical conditions of the deep Earth. Consequently, in the last decade a large part of the geoscientific community focused on their critical role in the evolution of terrestrial planets, from the early solar system, dominated by Giant Impacts and large-scale magma oceans, all the way to the current evolved state of solidified planets.
Workshop programme
- Using MELTS code the participants will be able to do the following:
- To model magmatic evolution scenarios as a series of steps in temperature and pressure (Gibbs energy minimization), temperature and volume (Helmholtz energy minimization), enthalpy and pressure (entropy maximization) or entropy and pressure (enthalpy minimization).
- To apply these scenarios to exploring open- and closed-system magmatic processes such as energy constrained assimilation, adiabatic decompression melting, or post-entrapment crystallization in phenocryst-hosted melt inclusions.
- To compute equilibrium states in systems constrained to follow oxygen fugacity buffers.
- To simulate forward, down-temperature, fractional crystallization and to learn what is possible in terms of inverse (up-temperature) fractionation modeling.
- To compute complete models of the melting regime underlying a mid-ocean ridge.
- To access all these calculations from within MATLAB or Python in order to enable seamless coupling to geodynamic codes or large-scale modeling efforts.
All workshop participants will leave with the necessary software installed and configured on their own computers and with membership in the users forum for ongoing communication among users and developers of the software.
- Using the UMD package the participants will be able to do the following:
- To extract all relevant results from simulations of ab intio molecular dynamics, and construct UMD (=”Universal Molecular Dynamics”) files.
- To apply the individual components of the UMD package to analyze these results.
- To calculate the pair distribution functions, determine the bond length, the size of the coordination sphere, the average coordination numbers.
- To build the connectivity matrix and from there to obtain the chemical speciation, including the population analysis, the polymerization of the melt, and the lifetimes of the different chemical species.
- To determine the mean-square displacements, and from there to extract the diffusion coefficients.
- To compute the self-correlation of the atomic velocities from which to obtain the vibrational spectrum of the fluid and the diffusion coefficients.
- To calculate the self-correlation of the stress tensor from which to estimate the viscosity of the fluid.
- To monitor the formation of gas bubbles as a model for magma degasing and devolatilization.
In addition, invited experimental lectures will illustrate modern methods for determining liquid properties at elevated temperature and pressure in the lab, including:
- In situ measurements in static presses: piston-cylinder apparatus, multi-anvil presses, diamond anvil cells; applications of synchrotron radiation.
- Shock wave methods for equations of state and phase transition determination.
- Criteria for attainment of equilibrium in phase relations and uncertainty assessment in transport property measurements.
Schedule:
The approximate schedule will be:
- 13th September 2020: Travel to Milos
- 14th September 2020: Introduction to Melt Thermodynamics and MELTS + Practice MELTS
- 15th September 2020: Introduction to Ab Initio Molecular Dyanmics + Practice UMD
- 16th September 2020: Practice MELTS
- 17th September 2020: Practice UMD
- 18th September 2020: Fieldtrip
Workshop Organizers
Paul Asimow is Professor of Geology and Geochemistry at Caltech. He developed pHMELTS and many of the tools since incorporated into alphaMELTS. He has supervised a number of students and postdocs with MELTS-related projects, both to extend and improve the underlying thermodynamic models and to apply the software to real world magmatic processes. He is also involved in shock-wave experiments and theoretical calculations on silicate liquid and minerals.
Razvan Caracas is a Senior Researcher at the École Normal Supérieure de Lyon, France. He is a computational mineral physicist interested in the functioning of planetary interiors. He use ab initio methods based on the density-functional theory and the density-functional perturbation theory to understand the behavior of minerals in a variety of planetary environments.
Ioannis Baziotis is Assistant Professor in Mineralogy and Petrology at the Agricultural University of Athens, Greece. He studies the mineralogy, petrology and geochemistry of terrestrial (mantle xenoliths, high and ultra-high pressure metamorphic rocks) and extra-terrestrial materials (Martian meteorites and chondrites), including metamorphism and impact-related processes inferred from Martian meteorites.
Workshop Sponsor
The European Geosciences Union is Europe’s premier geosciences union, dedicated to the pursuit of excellence in the Earth, planetary, and space sciences for the benefit of humanity, worldwide. It was established in 2002 as a merger of the European Geophysical Society (EGS) and the European Union of Geosciences (EUG), and has headquarters in Munich, Germany.