Fluid Phase Equilibria, Vol.387, 24-31, 2015
Molecular dynamics simulation investigation of hexanoic acid adsorption onto calcite (10(1)over-bar4) surface
In this paper we report the results of classical molecular dynamics (MD) simulation of hexanoic acid adsorption on calcite (10 (1) over bar4) surface plane using Pavese and AMBER force fields for calcite and hexanoic acid, respectively. Pair correlation function, strictly suggests a well-structured adsorption. Density profile indicates the adsorption occurs through double-bonded O atom of the acid head group by direct interaction with Ca atom at calcite surface. Adsorption orientation of H and double-bonded O atoms was found to be as lock and key with respect to calcite surface Ca and O atoms, facilitating an effective adsorption. Adsorption time evolution indicates that O atom adsorption is accompanied by wobbling between two Ca sites. Hence, in spite of H-bonding by acid group, the surface atoms matrix enforces adsorption with some dynamics. Hexanoic acid molecules adopt a well-defined adsorption layer immediate to the surface patterned by the calcite surface structure. The adsorption energy was estimated to be almost 187 times that of octane, a relevant nonpolar fraction of crude oil. These are useful information on the microscopic behavior and adsorption mechanism of crude oil polar fraction in calcite reservoir. The density profile over simulation time steps is proposed as a practical tool to study the dynamics of adsorbed layer in the same way as the most surface spectroscopic method like surface scattering spectroscopy do but in dynamic mode sub-picosecond scale. The intrinsic orientation of double-bonded O and H atoms of the head-group are presented by bivariate maps which helps establishing the key factors of calcite surface activity. (C) 2014 Elsevier B.V. All rights reserved.
Keywords:Calcite surface;Hexanoic acid adsorption energy;Molecular dynamics simulation;Orientational bivariate map;Self-diffusion