Molecular Dynamics Simulations of CO2/Water/Quartz Interfacial Properties: Impact of CO2 Dissolution in Water

The safe trapping of carbon dioxide (CO2) in deep saline aquifers is one of the major concerns of CO2 sequestration. The amount of capillary trapping is dominated by the capillary pressure of water and CO2 inside the reservoir, which in turn is controlled by the interfacial tension (IFT) and the contact angle (CA) of CO2/water/rock systems. The measurement of IFT and CA could be very challenging at reservoir conditions, especially in the presence of toxic cocontaminants. Thus, the ability to accurately predict these interfacial properties at reservoir conditions is very advantageous. Although the majority of existing molecular dynamics (MD) studies of CO2/water/mineral systems were able to capture the trends in IFT and CA variations with pressure and temperature, their predictions often deviated from experimental data, possibly due to erroneous models and/or overlooked chemical reactions. The objective of this study was to improve the MD predictions of IFT and CA of CO2/water/quartz systems at various pressure and temperature conditions by (i) considering the chemical reactions between CO2 and water and (ii) using a new molecular model for α-quartz surface. The results showed that the presence of carbonic acid at the CO2/water interface improved the predictions of IFT, especially at low temperature and high pressure where more CO2 dissolution occurs. On the other hand, the effect on CA was minor. The slight decrease in CA observed across the pressure range investigated could be attributed to an increase in the total number of H-bonds between fluid molecules and quartz surface.