Structure, Dynamics, and Hydration Free Energy of Carbon Dioxide in Aqueous Solution: A Quantum Mechanical/Molecular...

The solvation of carbon dioxide in solution represents a key step for the capture and fixation CO₂ in nature, which may be further influenced by the formation of (bi)carbonate species and/or the formation of CO₂ clusters in solution. The latter processes are strongly dependent on the exact environment of the liquid state (e.g., pH value, solvated ions, etc.) and may interfere with the experimental determination of structural, dynamical, and thermodynamic properties. In this work a hybrid quantum mechanical/molecular mechanical (QM/MM) simulation approach at correlated ab initio level of theory resolution-of-identity second-order Møller-Plesset Perturbation Theory (RI-MP2) has been applied in the framework of thermodynamic integration (TI) to study structure, dynamics, and the hydration free energy of a single carbon dioxide molecule in aqueous solution. A detailed analysis of the individual QM/MM potential energy contributions demonstrate that the overall potential remains highly consistent over the entire sampling phase and that no artificial contributions are influencing the determination of the hydration free energy. The latter value of 0.01 ± 0.92 kcal/mol was found in very good agreement with the values of 0.06 and 0.24 kcal/mol obtained via quasi-chemical theory and experimental measurements, respectively. In order to obtain detailed information about the C- and O _C-water interaction, conically restricted regions with respect to the main axis of the CO₂ molecule have been employed in structural analysis. The presented data not only provide detailed information about the hydration properties of CO₂ but act as a critical validation of the simulation technique, which will be beneficial in the study of nonaqueous solvents such as pure and aqueous NH₃ solutions, which have been suggested as potential candidates to capture CO₂ from anthropogenic sources.