Prasetyo Niko, Hofer Thomas S
Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry , University of Innsbruck , Innrain 80-82A , A-6020 Innsbruck , Austria.
Austria-Indonesia Centre (AIC) for Computational Chemistry , Universitas Gadjah Mada , Sekip Utara, Yogyakarta 55281 , Indonesia.
J Chem Theory Comput. 2018 Dec 11;14(12):6472-6483. doi: 10.1021/acs.jctc.8b00557. Epub 2018 Dec 3.
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 -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.
二氧化碳在溶液中的溶剂化是自然界中捕获和固定二氧化碳的关键步骤,这可能会受到(双)碳酸盐物种的形成和/或溶液中二氧化碳簇的形成的进一步影响。后一过程强烈依赖于液态的确切环境(例如,pH值、溶剂化离子等),并且可能会干扰结构、动力学和热力学性质的实验测定。在这项工作中,一种在相关的从头算理论分辨率恒等式二阶莫勒-普莱塞特微扰理论(RI-MP2)水平上的混合量子力学/分子力学(QM/MM)模拟方法已被应用于热力学积分(TI)框架中,以研究水溶液中单个二氧化碳分子的结构、动力学和水合自由能。对各个QM/MM势能贡献的详细分析表明,在整个采样阶段,总势能保持高度一致,并且没有人为贡献影响水合自由能的测定。发现后者的值为0.01±0.92千卡/摩尔,与通过准化学理论和实验测量分别获得的0.06和0.24千卡/摩尔的值非常吻合。为了获得有关C-水和O-水相互作用的详细信息,在结构分析中采用了相对于二氧化碳分子主轴的锥形受限区域。所呈现的数据不仅提供了有关二氧化碳水合性质的详细信息,而且作为模拟技术的关键验证,这将有利于研究非水溶剂,如纯氨溶液和氨水溶液,它们已被建议作为从人为来源捕获二氧化碳的潜在候选物。