Ab initio Study of the Interactions between CO2 and N-Containing Organic Heterocycles

• Published in: Chemphyschem 2009-02, Vol.10 (2), p.374-383
• Main Authors: Vogiatzis, Konstantinos D., Mavrandonakis, Andreas, Klopper, Wim, Froudakis, George E.
• Format: Article
• Language: English
• Subjects: Ab initio calculations; Atomic and molecular collision processes and interactions; Atomic and molecular physics; Calculations and mathematical techniques in atomic and molecular physics (excluding electron correlation calculations); carbon dioxide; Coupled cluster theory; Density-functional theory; donor-acceptor systems; Electronic structure of atoms, molecules and their ions: theory; Exact sciences and technology; hydrogen bonds; Interatomic and intermolecular potentials and forces, potential energy surfaces for collisions; Intermolecular and atom-molecule potentials and forces; Physics; weak interactions
• ISSN: 1439-4235; 1439-7641
• DOI: 10.1002/cphc.200800583

In the garden of dispersion: High‐accuracy ab initio calculations are performed to determine the nature of the interactions and the most favorable geometries between CO2 and heteroaromatic molecules containing nitrogen (see figure). Dispersion forces play a key role in the stabilization of the dimer, because correlation effects represent about 50 % of the total interaction energy. The interactions between carbon dioxide and organic heterocyclic molecules containing nitrogen are studied by using high‐accuracy ab initio methods. Various adsorption positions are examined for pyridine. The preferred configuration is an in‐plane configuration. An electron donor–electron acceptor (EDA) mechanism between the carbon of CO2 and the nitrogen of the heterocycle and weak hydrogen bonds stabilize the complex, with important contributions from dispersion and induction forces. Quantitative results of the binding energy of CO2 to pyridine (C5H5N), pyrimidine, pyridazine, and pyrazine (C4H4N2), triazine (C3H3N3), imidazole (C3H4N2), tetrazole (CH2N4), purine (C5H4N4), imidazopyridine (C6H5N3), adenine (C5H5N5), and imidazopyridamine (C6H6N4) for the in‐plane configuration are presented. For purine, three different binding sites are examined. An approximate coupled‐cluster model including single and double excitations with a perturbative estimation of triple excitations (CCSD(T)) is used for benchmark calculations. The CCSD(T) basis‐set limit is approximated from explicitly correlated second‐order Møller–Plesset (MP2‐F12) calculations in the aug‐cc‐pVTZ basis in conjunction with contributions from single, double, and triple excitations calculated at the CCSD(T)/6‐311++G** level of theory. Extrapolations to the MP2 basis‐set limit coincide with the MP2‐F12 calculations. The results are interpreted in terms of electrostatic potential maps and electron density redistribution plots. The effectiveness of density functional theory with the empirical dispersion correction of Grimme (DFT‐D) is also examined. In the garden of dispersion: High‐accuracy ab initio calculations are performed to determine the nature of the interactions and the most favorable geometries between CO2 and heteroaromatic molecules containing nitrogen (see figure). Dispersion forces play a key role in the stabilization of the dimer, because correlation effects represent about 50 % of the total interaction energy.