Wetting Properties of the CO2–Water–Calcite System via Molecular Simulations: Shape and Size Effects

13 September 2019, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO2 and the pore fluids. In this context, the interfacial tension (IFT) between CO2 and the aqueous fluid and the contact angle, theta, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO2–water IFT and the CO2–water–calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO2 model, the SPC/Fw water model indeed underestimates the IFT in the low pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ~5 mN/m. Literature data regarding the water contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO2 atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii, to study the size dependence of the water contact angle on calcite in the presence of CO2. Our results suggest that the contact angle of a cylindrical water droplet on calcite {10.4}, in the presence of CO2, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38 degrees, indicating that calcite is strongly wetted by water.


Interfacial tension
Contact angle
Molecular Dynamics Simulations
CO2 geological storage


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.