A molecular dynamics study of carbon dioxide in water: diffusion,structure and thermodynamics

Atomistic simulations are reported of a model of CO₂ in water. CO₂ is modelled by partial charges and Lennard-Jones interaction sites on each atom; the SPC/E model for water is used. Good agreement with experiment is found for the translational diffusion constants. The variation of the dynamics with the potential parameter was investigated. As expected, the orientational correlation times increase as the magnitude of the quadrupole moment is increased, but the translational diffusion constants are found to be surprisingly insensitive to the magnitude of the CO₂ quadrupole moment. The translational friction coefficient was resolved into electrostatic, Lennard-Jones and cross-terms; the Lennard-Jones contribution is found to be the largest. Varying the Lennard-Jones size parameter affects both translational and reorientational motion. In order to try to understand these results further, the variation of solvation free energy was investigated and the solvent structure around carbon dioxide was examined as the electrostatic and Lennard-Jones parameters were changed. The temperature dependence of the self-diffusion constant of pure SPC/E water was determined.