Researcher profile

B. Guillot

B. Guillot contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - Baseline
5works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2012arXiv

Viscosity and viscosity anomalies of model silicates and magmas: a numerical investigation

We present results for transport properties (diffusion and viscosity) using computer simulations. Focus is made on a densified binary sodium disilicate 2SiO$_2$-Na$_2$O (NS2) liquid and on multicomponent magmatic liquids (MORB, basalt). In the NS2 liquid, results show that a certain number of anomalies appear when the system is densified: the usual diffusivity maxima/minima is found for the network-forming ions (Si,O) whereas the sodium atom displays three distinct régimes for diffusion. Some of these features can be correlated with the obtained viscosity anomaly under pressure, the latter being be fairly well reproduced from the simulated diffusion constant. In model magmas (MORB liquid), we find a plateau followed by a continuous increase of the viscosity with pressure. Finally, having computed both diffusion and viscosity independently, we can discuss the validity of the Eyring equation for viscosity which relates diffusion and viscosity. It is shown that it can be considered as valid in melts with a high viscosity. On the overall, these results highlight the difficulty of establishing a firm relationship between dynamics, structure and thermodynamics in complex liquids.

preprint2011arXiv

Carbon dioxide in silicate melts: A molecular dynamics simulation study

The distribution, recycling and storage of carbon in the Earth are of fundamental importance to understand the global carbon cycle between the deep Earth and near surface reservoirs. Degassing of CO2 at mid-ocean ridges may give information on the source region but the very low solubility of CO2 in tholeitic basalts has for consequence that near all Mid-Ocean Ridge Basalts glasses exsolve their CO2 rich vapor at shallow depth as they approach the ocean floor. Hence their CO2 contents mostly represent the pressure at eruption and not the source region. Recent petrological investigations have shown that the presence of carbonates at depth in the upper mantle has a large effect on the solidus of carbonated silicates by inducing incipient melting at much lower temperature. So the role of carbon-rich melts at great depth is now becoming a credible scenario to explain the extraction of CO2 from the source region to the surface. During the last three decades many studies have been devoted to measure the solubility of CO2 in silicate melts of various composition. But due to experimental difficulties these studies were generally restricted to low and moderate pressures (below ~20 kbar). By performing a series of molecular dynamics simulation where a supercritical CO2 phase is in contact with a silicate melt of various composition (from felsic to ultrabasic) at different temperatures (1473-2273K) and pressures (20-150kbar), we have been able to evaluate the solubility of CO2, the population of molecular and carbonate species, their diffusivity through the melt and the local structure. We show that this kind of molecular simulation is a useful theoretical guide to better understand the behavior of CO2 in magmas at depth.

preprint2011arXiv

Electronic redistribution around oxygen atoms in silicate melts by ab initio molecular dynamics simulation

The structure around oxygen atoms of four silicate liquids (silica, rhyolite, a model basalt and enstatite) is evaluated by ab initio molecular dynamics simulation. Thanks to the use of maximally localized Wannier orbitals to represent the electronic ground state of the simulated system, one is able to quantify the redistribution of electronic density around oxygen atoms as a function of the cationic environment and melt composition. It is shown that the structure of the melt in the immediate vicinity of the oxygen atoms modulates the distribution of the Wannier orbitals associated with oxygen atoms. In particular the evaluation of the distances between the oxygen-core and the orbital Wannier centers and their evolution with the nature of the cation indicates that the Al-O bond in silicate melts is certainly less covalent than the Si-O bond while for the series Mg-O, Ca-O, Na-O and K-O the covalent character of the M-O bond diminishes rapidly to the benefit of the ionic character. Furthermore it is found that the distribution of the oxygen dipole moment coming from the electronic polarization is only weakly dependent on the melt composition, a finding which could explain why some empirical force fields can exhibit a high degree of transferability with melt composition.

preprint2011arXiv

Noble gases in high-pressure silicate liquids: A computer simulation study

The test particle method has been used in conjunction with molecular dynamics simulations to evaluate the solubility of noble gases in silicate melts of various compositions. At low pressure the calculated solubility constants (the inverse of the Henry's constant) are in excellent agreement with data of the literature. In particular it is found that the solubility constant (i) decreases when the size of the noble gas increases, (ii) decreases from silica-rich to silica-poor composition of the melt, and (iii) is positively correlated with the temperature. Moreover it is shown that the solubility is governed primarily by the entropic cost of cavity formation for inserting the noble gas into the melt and secondarily by its solvation energy. Interestingly, the behaviour of these two contributions differ from each other as the entropic cost of cavity formation increases strongly with the size of the solute atom to insert whereas large atoms are better solvated than small ones. With all melt composition investigated here (silica, rhyolite, MORB and olivine), the calculated solubility curves exhibit the same qualitative behaviour with pressure; a steep rise culminating in a broad maximum followed by a gradual decrease of the solubility at higher pressure. At variance with LHDAC experiments (Chamorro et al. (1996, 1998) and Bouhifd et al. (2006, 2008)) where a Ar solubility drop is observed at about 50 kbar in silica and molten olivine and in the pressure range ~100-160 kbar with other melt composition, we do not find such a sudden change of the solubility.

preprint2009arXiv

Computer modeling of natural silicate melts: what can we learn from ab initio simulations

The structural and dynamical properties of four silicate liquids (silica, rhyolite, a model basalt and enstatite) are evaluated by ab initio molecular dynamics simulation using the density functional theory and are compared with classical simulations using a simple empirical force field. For a given composition, the structural parameters of the simulated melt vary little between the two calculations (ab initio versus empirical) and are in satisfactory agreement with structure data available in the literature. In contrast, ionic diffusivities and atomic vibration motions are found to be more sensitive to the details of the interactions. Furthermore, it is pointed out that the electronic polarization, as evaluated by the ab initio calculation, contributes significantly to the intensity of the infrared absorption spectra of molten silicates, a spectral feature which cannot be reproduced using nonpolarizable force field. However the vibration modes of TO4 species and some structural details are not accurately reproduced by our ab initio calculation, shortcomings which need to be improved in the future.