First-principles Modeling of Minerals, Melts and Fluids at High Pressures and High Temperatures
Institute of Geology and Mineralogy, University of Cologne (Germany)
Local Project ID:
HPC Platform used:
JUWELS of JSC
The solid Earth is mainly composed of rocks consisting of both minerals, i.e. crystalline phases of natural occurrence, and amorphous or glassy material. In geologically active regions, e.g., continental margins, subduction zones or volcanic environments, high temperatures and pressures and/or volatile chemical components lead to the formation of molten rock, mineral alteration processes, or the release of magmatic-hydrothermal fluids. Such processes have shaped our present-day Earth and triggered chemical differentiation that is reflected in various rock types, ranging from common basalts or granites to technologically important ore deposits. The solid, liquid and gaseous phases found in the Earth’s interior are generally chemically complex, and their properties are often not directly accessible at the relevant thermodynamic conditions. Therefore, interpretation of field observations and understanding of the underlying mechanisms of geological processes require insights from laboratory experiments and numerical simulations.
In the context of this long-term supercomputing project, our research focused on the following two goals:
To achieve these goals, proven computational chemistry methods were used, namely molecular dynamics simulations in combination with electronic structure or classical interaction potential methods. This numerical modeling approach makes it possible to simulate the motion of atoms and molecules at defined volume (or pressure) and temperature in real time, i.e. on the time scales of femto- to nanoseconds. Depending on the system, the particle interactions are either described by rather accurate but computationally expensive first-principles methods based on density-functional theory, or by more efficient force fields, which usually have a limited scope. In both cases, the system sizes and simulation times are limited by the available computational resources. The allocated computational time allowed us to perform extensive simulations with dedicated MPI-OpenMP-parallelized codes (e.g. CPMD, CP2K, ABINIT) on the JURECA and JUWELS CPU supercomputers at JSC.
One of our research interests is the study of speciation and thermodynamic behavior of metal ions in supercritical aqueous fluids. Such fluids are geologically important in various environments such as magmatic-hydrothermal systems, hot geothermal reservoirs, or subduction zones. We considered pure H2O and chloride and fluoride solutions as simplified solvents for the metal cations (e.g. Sahle et al. 2013; Jahn et al. 2015; Sahle et al. 2017; Stefanski et al. 2020). Significant structural changes in solvents are observed when pressure and temperature vary from near-surface conditions to conditions in the lower crustal or upper mantle (e.g. Sahle et al. 2018; Stefanski et al. 2018; Elbers et al. 2021). Solvent properties determine speciation and solubility of metal ions. For example, the low dielectric constant of water destabilizes highly charged species at high temperatures and pressures, leading to the formation of (near) neutral contact ion pairs.
Fig. 1 shows, as an example, structural information of Y speciation in mixed fluoride and chloride supercritical solutions from ab initio molecular dynamics simulations using the CPMD code. In order to obtain the stability constants of different complexes, the dissociation constants were derived using the potential of mean force approach of thermodynamic integration. This involves determining the change in free energy during stepwise removal of individual ligands (e.g. Cl or F). In practice, this is done by performing constraint molecular dynamics simulations in which the distance between the cation and the dissociating anion is held constant and the average force between these two ions is determined.
Repeating this procedure for different distances between the equilibrium distance of the associated complex and a distance at which the complex is considered dissociated, the potential of mean force is obtained. The free energy of dissociation is given by the integral of this potential. Stability constants describe the reverse process of stepwise association, but they are easily derived from the (standardized) dissociation energies. Finally, species distribution diagrams are constructed based on the computed species stability constants (Fig. 2). In the example presented here, the fluoride species become significant and increasingly dominant at fluoride concentrations above a few millimole per kg H2O. Further details of this study are published in Stefanski and Jahn (2020).
Molecular simulations are an indispensable complement to in-situ experimental studies of the vibrational and electronic properties of aqueous fluids at high pressure and temperature. Once reliable structural models are established, the thermodynamic properties derived from the simulations provide missing parameters for thermodynamic and reactive flow simulations in a wide range of geological settings. Furthermore, the simulations can help to develop new types of thermodynamic models, e.g., for the development of superhot geothermal resources (e.g. www.geoproproject.eu).
The results obtained in this project (chpo15) have been published in about 30 peer-reviewed articles, many of them in collaboration with experimentalists. Another highlight from this project, a pioneering study of the structural evolution of SiO2 glass up to extreme pressures of 172 GPa, was presented a previous report. We gratefully acknowledge project funding from the DFG (JA1469/10-1).
Elbers, M., Schmidt, C., Sternemann, C., Sahle, C., Jahn, S., Albers, C., Sakrowski, R., Gretarsson, H., Sundermann, M., Tolan, M., Wilke, M., 2020. Ion association in hydrothermal aqueous NaCl solutions: Implications for the microscopic structure of supercritical water, Phys. Chem. Chem. Phys. 23, 14845-14856. https://doi.org/10.1039/D1CP01490K
Jahn, S., Dubrail, J., Wilke, M., 2015. Complexation of Zr and Hf monomers in supercritical aqueous solutions: Insights from ab initio molecular dynamics simulations, Chem. Geol. 418, 30-39. https://doi.org/10.1016/j.chemgeo.2014.10.012
Sahle, C.J., Sternemann, C., Schmidt, C., Lehtola, S., Jahn, S., Simonelli, L., Huotari, S., Hakala, M., Pylkkänen, T., Nyrow, A., Mende, K., Tolan, M., Hämäläinen, K., Wilke, M., 2013. Microscopic structure of water at elevated pressures and temperatures, PNAS 110, 6301-6306. https://doi.org/10.1073/pnas.1220301110
Sahle, C.J., Niskanen, J., Schmidt, C., Stefanski, J., Gilmore, K., Forov, Y., Jahn, S., Wilke, M., Sternemann, C., 2017. Cation hydration in supercritical NaOH and HCl aqueous solutions, J. Phys. Chem. B 121, 11383-11389. https://doi.org/10.1021/acs.jpcb.7b09688
Sahle, C.J., Niskanen, J., Gilmore, K., Jahn, S., 2018. Exchange-correlation functional dependence of the O 1s excitation spectrum of water, J. Electron. Spectrosc. Relat. Phenom. 222, 57-62. https://doi.org/10.1016/j.elspec.2017.09.003
Stefanski, J., Schmidt, C., Jahn, S., 2018. Aqueous sodium hydroxide (NaOH) solutions at high pressure and temperature: insights from in situ Raman spectroscopy and ab initio molecular dynamics simulations, Phys. Chem. Chem. Phys. 20, 21629-21639. https://doi.org/10.1039/c8cp00376a
Stefanski, J., Jahn, S., 2020. Yttrium speciation in subduction-zone fluids from ab initio molecular dynamics simulations, Solid Earth 11, 767-789. https://doi.org/10.5194/se-11-767-2020
Prof. Dr. Sandro Jahn, Dr. Johannes Stefanski, Maximilian Schulze (all: Institute of Geology and Mineralogy, University of Cologne)
Prof. Dr. Sandro Jahn
Institute of Geology and Mineralogy
University of Cologne
Zülpicher Str. 49b, D-50674 Köln (Germany)
e-mail: s.jahn [at] uni-koeln.de
Local project ID: chpo15