TY - JOUR
T1 - Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation
AU - Leyssale, Jean Marc
AU - Galvez, Matthieu E.
AU - Valdenaire, Pierre Louis
AU - Pellenq, Roland
AU - van Duin, Adri C.T.
N1 - Funding Information:
Access to the computer facilities of the Centre Informatique National Supérieur de l'Enseignement Supérieur (grant A0040910065) and of the Mésocentre de Calcul Intensif en Acquitaine are gratefully acknowledged. MG thanks the Branco Weiss/Society in Science fellowship, and the MultiScale Materials Science for Energy and Environment laboratory at MIT for supporting 2 visiting stays during which this project emerged. ACTvD acknowledges funding from the Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials (MUSE) Center, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Science. We thank 3 anonymous reviewers for thorough and constructive reviews. Insightful comments by Xin Zhong also improved the manuscript.
Funding Information:
Access to the computer facilities of the Centre Informatique National Supérieur de l’Enseignement Supérieur (grant A0040910065) and of the Mésocentre de Calcul Intensif en Acquitaine are gratefully acknowledged. MG thanks the Branco Weiss/Society in Science fellowship, and the MultiScale Materials Science for Energy and Environment laboratory at MIT for supporting 2 visiting stays during which this project emerged. ACTvD acknowledges funding from the Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials (MUSE) Center, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Science. We thank 3 anonymous reviewers for thorough and constructive reviews. Insightful comments by Xin Zhong also improved the manuscript.
Publisher Copyright:
© 2022
PY - 2022/7/15
Y1 - 2022/7/15
N2 - We present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth's lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO2 and CH4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth's mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.
AB - We present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth's lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO2 and CH4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth's mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.
UR - http://www.scopus.com/inward/record.url?scp=85131916416&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85131916416&partnerID=8YFLogxK
U2 - 10.1016/j.gca.2022.04.024
DO - 10.1016/j.gca.2022.04.024
M3 - Article
AN - SCOPUS:85131916416
SN - 0016-7037
VL - 329
SP - 106
EP - 118
JO - Geochmica et Cosmochimica Acta
JF - Geochmica et Cosmochimica Acta
ER -