Selectivity of Cobalt Corrole for CO vs. O2 and N2 in Indoor Pollution

Coal combustion causes indoor pollution of CO. In this work, DFT calculations on cobalt corrole (Co(Cor)) with three most common indoor gas molecules (N2, O2 and CO) were performed. The Mulliken spin densities show that the ground states of Co(N2)(Cor), Co(CO)(Cor) and Co(OC)(Cor) have an anti-ferromagnetic coupling fashion of the electrons on the Co 3d z 2 orbital and the π orbital of the corrole ring. However, Co(O2)(Cor) has a triplet ground state. With the spin contamination corrections, the Co(N2)(Cor) binding energy was obtained at −50.6 kcal mol−1 (B3LYP-D3). While CO can interact with Co(Cor) in two different ways, and their binding energies were −22.8 and −10.9 kcal mol−1 (B3LYP-D3) for Co(CO)(Cor) and Co(OC)(Cor), respectively. The natural bond orbital charges on the axial ligands (NO, CO, OC) are increased upon the chemical bond formation. These are the cause of the shorten metal-ligand bond and the increase of the wavenumber of the metal-ligand bond vibrational transitions. While the charges for O2 are decreased, leading to bond elongation as well as the decrease of the wavenumber upon complexation. Overall, O2 was found to be hardly coordinated with Co(Cor). This study provides a detailed molecular understanding of interactions between a gas sensor and gaseous indoor air-pollutants.

some capabilities to coordinate with CO, NO and O 2 , but cannot be used for gas sensing purposes as the former derivatives. Recently several papers have demonstrated that Co(III)(Cor) exhibited an infinite selectivity for CO with respect to N 2 and O 2 12,13 . This specific property enables them to be used as sensing components for gas detectors from the selective coordination of CO to the central Co(III) ion by a chemisorption process where N 2 and O 2 hardly coordinates.
We undertook a quantum chemical study of Co(Cor) interactions with diatomic gas donor molecules (O 2 , N 2 and CO) from a molecular level. Quantum theory of atoms in molecules (AIM) and natural bond orbital (NBO) analysis was performed to understand the nature of weak interaction in the five-coordinate cobalt corrole complexes. This work is imperative to understand the nature of the adsorption bonding interactions.

Results and Discussion
In this study, we reported investigations of the molecular interactions of Co(Cor) with a series of common atmospheric diatomic ligands by different density functional methods. The unsubstituted corrole ligand was used for the theoretical calculations. The optimized structures of the Co(L)(Cor) (L = O 2 , N 2 , CO, OC) obtained at the BP86/def2-TZVP level are provided in Figures S1-S4 (Supporting Information). Only the most stable structures are discussed in the paper. The important parameters affecting indoor air pollution from a molecular level by analyzing the results of AIM, NBO as well as the binding energies were given in details.

Electronic configurations and spin densities. Previous study on Co(Cor) indicates the ground state is a
Co(III) triplet ground state with an electronic configuration of (d ) x y 2 2 14 . The ground states of both the CO and the N 2 free gases are singlet states with closed-shell electronic configurations, but for the O 2 free gas the ground state is a triplet state with two unpaired electrons occupying two degenerate anti-bonding molecular orbitals. Upon complexation with the ligands, the spin state of the central cobalt atom is changed from an intermediate-spin (S = 1) state to a low-spin (S = 0) state. The Co(Cor) system has C 2v symmetry, and the five-coordinate Co(L)(Cor) systems possess either C s or C 1 symmetry. In this study, the Co-L bond distances were predicted much large by B3LYP-D3/def2-TZVP. In some cases, the bond was unexpected broken, for instance, the Co-O 2 bond distance in the triplet state was predicted at 3.250 Å (B3LYP-D3/def2-TZVP), and the O-O distance (1.205 Å, B3LYP-D3/def2-TZVP) is also close to the free O 2 value (1.204 Å, B3LYP-D3/def2-TZVP). This implies that the B3LYP-D3 functional is not suitable in predicting the structures for this type of weak interaction. On the other hand, the Co-L bonds predicted by BP86/def2-TZVP looks much reasonable, i.e., 1.727-1.902 Å for the Co-CO distances; 1.964 Å for the Co-OC distance; 1.788-2.106 Å for the Co-N 2 distances; 1.802-1.854 Å for the Co-O 2 distances. Furthermore, the single point calculations with the B3LYP/def2-TZVP, OLYP/def2-TZVP, B3LYP-D3/def2-TZVP levels of theory were carried out, making use of the BP86/def2-TZVP optimized structures. Since, CO has two binding sites: the oxygen atom and the carbon atom. Both of the two possible binding sites were considered in this study.
The corrole ligand belongs to the point group C 2v . The two HOMO corrolate π orbitals are near-degenerate, which are represented by either a 2 or b 2 symmetry. Similarly, the two LUMO corrolate π orbitals are represented by either a 2 * or b 2 * symmetry. The a 2 and b 2 HOMOs of corrole are close analogues of the a 1u and a 2u porphyrin HOMOs, respectively. They are known as the Gouterman's four-orbital to explain the electronic absorption spectra of porphyrins 15 . Co(Cor) tends to bind various axial ligands, thus the electrons in the HOMO might be readily excited into the cobalt 3d orbitals to form a Co(II) π-cation radical state. The Mulliken spin densities of the Co(L)(Cor) singlet and triplet states are presented in Tables 1 and 2. A quantity that is useful in this respect is the so-called effective number of unpaired electrons (NUE) and denoted as 16,17 : where n i (0 ≤ n i ≤ 2) is the occupation number of the natural orbital; < > S 2 is the expectation value of the total spin-squared operator; S is the total spin quantum number. NUE may be used to diagnose the character of the DFT solution obtained for any open-shell state.
As shown from Table 1, the two pure functionals, OLYP and BP86 obtained a restricted solution (<S 2 > = 0). It is due to the degree of separation of the positive and negative spin densities, which depends on the xc functional. The separation by pure functionals is relatively small 18 . By mixing HF with DFT, the hybrid functionals, B3LYP (20% HF) and B3LYP-D3 (20% HF) obtain a Broken Symmetry (BS) solution. The two pure functionals predict the closed-shell singlet state for the Co(L)(Cor) (L = O 2 , N 2 , CO, OC), which are non-broken symmetry solutions and can be described as low-spin Co(III)-like system. Contrary to this expectation, the hybrid functionals show that the most stable state is a BS solution. For the Co(L)(Cor) (L = N 2 , CO, OC) BS states, the electron in the π orbital ("a 2u " HOMO) in the corrole HOMO is anti-ferromagnetically coupled with the unpaired electron in the d 3 z 2 orbital. The Co ion in the complex is better described as a partial Co(II) character and the corrole ring is a radical with a charge of "−2" 19 . In the BS singlet Co(II)(S = 1/2)(L)(Cor •2− ), the negative Mulliken spin populations at the meso carbons ( Fig. 1) further proves the biradical character of the singlet system. As seen from the spin density plot (Fig. 1), the unpaired electron in the Co d 3 z 2 orbital is anti-ferromagnetically coupled with another unpaired electron in the π orbital on the corrole ring to form a singlet state.
For the Co(O 2 )(Cor) BS state (B3LYP-D3/def2-TZVP), the spin population of O 2 is 1.3975. In contrast, the free O 2 gas molecule has two unpaired electrons. There is no significant bonding between O 2 and Co(Cor). The electronic configuration for the Co(Cor) fragment is described as one electron in the metal anti-ferromagentically coupled one electron in the corrole. The expectation value, 1.3247, also shows that there is a large spin contamination. Because it is a Cor 2− π radical state with electron transfer from corrole π orbital to one of the metal orbitals.
The ground state was found to be the triplet state, which is the same as BP86. All of the above models are different from the models for iron porphyrins: Weiss model where Fe(III) is anti-ferromagnetically coupled with O 2 − to form a singlet ground state and Pauling model where the system was described as Fe(II)(S = 0)O 2 (S = 0) 20 .   The BS singlet state contains some amount of triplet, quintet, …, states. The contamination was approximately corrected by the standard spin projection technique (Eqs 2-4) 21,22 . The BP86 and OLYP functionals could not find the BS solution. The correction to the energy was done for B3LYP and B3LYP-D3, and shown in Table 3. From the expectation values, we can notice that the BS singlet states are strongly contaminated by the triplet state, for Co(O 2 )(Cor) it contains 34% singlet and 66% triplet character for Co(O 2 )(Cor). For the other three system, the triplet character is much smaller, while 65% singlet and 35% triplet character for Co(N 2 )(Cor), 89% singlet and 11% triplet character for Co(CO)(Cor), and 55% singlet and 45% triplet character for Co(OC)(Cor). Meanwhile, the correction for Co(O 2 )(Cor) is positive, for other three systems are negative. This is due to the different ground spin state as shown in Table 4, where it shows the relative energies of the singlet and triplet electronic states of Co(L)(Cor) (L = O 2 , N 2 , CO, OC) with respect to the ground state. The results obtained with OLYP and B3LYP are listed in the supporting information (Table S1). As seen from Table 4, the pure functional delivers singlet ground states for Co(N 2 )(Cor), Co(CO)(Cor) and Co(OC)(Cor), but a triplet ground state for Co(O 2 )(Cor). While, B3LYP-D3 gives a different ground state for Co(N 2 )(Cor). However, the singlet state of Co(N 2 )(Cor) becomes the ground state after the spin contamination corrections. The singlet-triplet splitting is relatively different to the method and also the systems. The BP86 methods give a 1.1, −6.3, and −17.5 kcal mol −1 for Co(O 2 )(Cor), Co(N 2 )(Cor) and Co(CO)(Cor), respectively. In contrast, the B3LYP-D3 method obtained    relatively different singlet-triplet splitting energies, where the corresponding data for the three systems are 8.2, −6.5, and −9.8 kcal mol −1 . However, the trends observed between them are the pure functional overstabilized low-spin relative to high-spin. In contrast, the high-spin state becomes more favored over the low-spin state with hybrid functionals.
Geometric parameters. In Table 5 the mean bond lengths, cobalt out-of-plane displacement and ligand stretching frequency for the selected systems are given, together with the ligand bond length and the ligand stretching frequencies of the free gases. The mean bond lengths of Co-N in the four systems are very close to each other. The axil ligand bond lengths and stretching frequencies were slightly larger than the ones in free gas phase. Large deviation in bond lengths was found in the Co-L bond length, while 1.  19 , slightly longer than the values in our study. Pure functionals always obtain a substantially shorter metal-ligand distance, which are in line with the similar studies 23 . The closest analogue is the binding with cobalt porphyrins. The CO interacts with Co(II) porphyrins show to have much weaker Co-CO bond distances at 2.01 Å which is also larger than the one studied here 24 . In Co(Cor) (without axil ligand), the cobalt out-of-plane displacement is calculated at 0.065 Å (BP86/def2-TZVP). After binding with the gas molecules, Co moves out towards the axial ligand. The Co out-of-plane displacements differs for the four systems: the Co(CO)(Cor) singlet state gives the largest displacement about 0.436 Å; in the other systems, the displacements are less than 0.32 Å. N 2 , CO and OC bind with Co(Cor) in a linear mode with σ-donation and π-back donation in the singlet states, but they are bonded in bent-modes with angles of 164.4°-168.7° in the triplet states. Meanwhile, O 2 binds in a bent-mode with an angle of about 120.4°-121.4° in both the singlet and triplet states. The CO and O 2 binding angles are in good agreement with the ones presented in the iron porphyrin histidine complex 25 . In the linear Co-L (L = N 2 , CO, OC) bond, the Co d 3 z 2 orbital is destabilized by the σ-donation from the axial ligand into the Co d 3 z 2 orbital. Similarly, the Co 3 d xy orbital is destabilized by the σ-antibonding interactions with the four corrole N atoms. While, the d x y 2 2 − orbital is an essentially non-bonding orbital. In the triplet state, one of the electron from the bonding orbitals (Co d 3 z 2 or 3 d xy ) excites into its corresponding σ-antibonding orbital. In the case of occupying the d 3 z 2 antibonding orbital, the overlap of the Co d 3 z 2 and the ligand σ orbital is much reduced. Thus, a bending of the Co-L bond is obtained 26 Tables 1-2. ZPVE obtained with BP86/def2-TZVP varies from about 1.0 to 3.7 kcal mol −1 depending on system, and they were included into all binding energies. In the previous study, the OLYP functional has been found to give a better match to the experimentally available binding energies for porphyrin systems than other  Table 5. Selected bond lengths (Å), cobalt out-of-plane displacements (Å), angles (degrees) and ligand stretching frequencies (cm −1 ) for the lowest Co(L)(Cor) states and the free gas molecules at the BP86/def2-TZVP level a . a r (Co-N) is the mean bond length between Co and the four nitrogen atoms. r (Co-L) is the distance between Co and the ligand. θ (Co-L) is the angle between Co and the ligand. d (Co-N4) is the cobalt out-of-plane (plane is defined by the four nitrogen atoms) displacement. r (L) is the bond length of the ligand. υ (L) is the ligand stretching vibrational frequency.
functionals (PBE0, B3LYP, BP86) 20 . It clearly shows the DFT binding energies are strongly functionals dependent (Tables 1 and 2), which conforms with previous studies 14,20,26 . The strongest BEs are predicted by BP86, then the weakest BEs are obtained by B3LYP. This further proves the hybrid functionals are not suitable to calculate this type of weak interactions. On the other hand, as compared with the results obtained by OLYP, B3LYP-D3 gives a better match than BP86 and B3LYP. The energy differences between B3LYP and B3LYP-D3 are less than 5 kcal mol −1 . This implies that the dispersion corrections are about 5 kcal mol −1 . Thus, they are quite important in calculating binding energies. To the best of our knowledge, there is no experimental measurement on the binding energies for these type of corrole systems. The close study is the binding energy with CO and O 2 to heme (porphyrin) system. In the computational study of the binding energies of CO and O 2 to the four-coordinate iron porphyrin, the binding energies were calculated as −52.7 kcal mol −1 for CO and −43.6 kcal mol −1 for O 2 (Fe: OLYP/QZVPP; ligand: OLYP/TZVPP) 20 . This is the same binding trend as our DFT results: CO forms a stronger bond than O 2 does. Since N 2 is isoelectronic with CO, the bonding in dinitrogen complexes is closely allied to that in carbonyl compounds, although N 2 is a weaker σ-donor and π-acceptor than CO. However, we found that N 2 forms a stronger bond than CO does. Moreover, the O 2 hardly coordinates to Co(Cor). This is a different conclusion as the previous studies, where CO is a significantly selected by cobalt corrole with respect to N 2 and O 2 12,13 . The different conclusion is probably due to the selected system, where Barbe et al. synthesized six different corroles: 5,15-dimesityl-10-(4-aminophenyl) corrole, 5,15-dimesityl-10-(4-chloroacetamidophenyl) corrole, 5,15-dimesityl-10-(2,4,6-trimethoxyphenyl) corrole, 5,15-dimesityl-10-(4-acetamidophenyl) corrole, 5,15-dimesityl-10-(2,4,6-trimethoxyphenyl) corrole, and 5,10,15-tris(2,6-dichlorophenyl) corrole. In this study, a very simplified corrole was used. The groups on the meso carbons may modify the chemical properties of the corrole system. AIM and NBO analysis. Topological analysis by using quantum theory of atoms in molecules (QTAIM) is commonly used to analyze structure, bonding and chemical reactivity. In this study, AIM analysis was used to evaluate the bond strength. AIM analysis was performed using the wavefunctions calculated at the B3LYP-D3/ def2-TZVP level. The AIM plots of the complexes with bond critical points (BCPs), ring critical points (RCPs) and electron density paths are shown in Fig. 2. The AIM molecular graph shows the BCPs along the lines joining the Co atom and the ligand, which clearly prove the presence of chemical bonds. The topological parameters, including electron densities ρ(r), Laplacian ∇ 2 ρ(r) at the BCPs with the B3LYP-D3/def2-TZVP method are listed in Table 6. It shows that the values of ρ(r) at BCP are in the range of 0.0514-0.1246 a.u. for the Co(L)(Cor) complexes. The positive values of ∇ 2 ρ(BCP) imply a closed-shell interaction between two bonded atoms 27 . In addition, the values of ∇ 2 ρ(BCP) are all positive as well, ranging from 0.3216-1.1472 a.u. It is well known that the higher values of ρ(BCP) indicate stronger molecular interactions. Therefore, as the results in Table 6  Electron density redistribution always accompanies the formation of a chemical bond. NBO analysis provides an efficient method for studying inter-molecular interactions and a convenient basis for investigating charge transfer (CT) in a molecular system. The charge donating centers and the charge accepting centers are often the lone pairs of the acceptor and the anti-bonding orbitals of the donor, respectively. Values listed in Table 6

Conclusions
The main aim of this study was to investigate the parameters affecting indoor air pollution from a molecular level. Theoretical study of the molecular interactions of cobalt corrole with a series of common atmospheric diatomic ligands (O 2 , N 2 , and CO) was performed with several DFT methods (BP86, OLYP, B3LYP, B3LYP-D3). The accuracy of the theoretical methods was corrected with spin contaminations. For N 2 , CO and OC, the ground state is a singlet state, in which the unpaired electron in the Co center is anti-ferromagnetically coupled with the unpaired electron in the corrole π cation radical. The complexes are better described as Co(II)(S = 1/2)(L)(Cor •2− ). However, for Co(O 2 )(Cor), the ground state is a triplet state. AIM analysis revealed that molecular interactions were present in the complexes, forming a chemical bond. NBO analysis clearly indicated the existence of charge transfer. The calculated binding energies show that N 2 forms the strongest complex with Co(Cor), while O 2 hardly coordinates with Co(Cor). The reason may be due to the simplified corrole in this study. In our study, the results provide a better understanding of the molecular interaction between a gas sensor and indoor pollutants.

Methods
DFT calculations on Co(Cor) systems ( Fig. 3) with several axial ligands, denoted as Co(L)(Cor) (L = O 2 , N 2 , CO, OC), were carried out with Gaussian 09 (Revision E.01) program package 28 . Geometries were optimized for all the possible low-lying spin states (Fig. 3, ligands at positions A, B and C) with BP86 plus the def2-TZVP basis set on all atoms. Co(Cor) system has a C 2v symmetry. All Co(L)(Cor) model systems possess either a C s or C 1 symmetry. In order to find the lowest-lying states, the C 1 symmetry was used in all calculations. Frequency calculations to approximate zero-point vibrational energy (ZPVE) were performed at the BP86/def2-TZVP level. A "verytight" optimization convergence criteria and an "ultrafine" numerical integration grid were used for the DFT calculations to obtain accurate results 29,30 . Two more functionals were used as well. One hybrid functional was employed: B3LYP has 20% Hatree-Fock (HF) exchange and 80% Slater (S) exchange. Besides BP86, another pure Generalized Gradient Approximation (GGA) functional was used: OLYP. Moreover, B3LYP-D3 is a hybrid functional, containing Grimme's D3 empirical dispersion correction 31 . The empirical dispersion correction has been shown to give an important contribution for general main group thermochemistry, kinetics, and noncovalent interaction 32 . The spin-unrestricted DFT formalism were used in all cases. For the (S = 0) state, both B3LYP and B3LYP-D3 give an unrestricted solution. For the two pure functionals, the unrestricted calculation converges to a restricted solution. The binding energies (BEs) are defined as the difference between the energies of the complex and the sum of the monomers, and then they are corrected with ZPVEs.
For the singlet states, a singlet biradical model (Broken Symmetry (BS) approach) was obtained with hydride functional, where there are two singly occupied molecular orbitals (SOMOs) of equal energy, and the spins of the two electrons are weakly anti-ferromagnetically coupled. However, the wavefunction is no longer an eigenfunction of the total spin, <S 2 >, thus introducing some error. This error is called spin contamination. The contamination could be approximately corrected by a standard spin projection technique 21,22 . The correction to the energy of the BS state is are the expectation values of the total spin-squared operator for the two calculations, S is the total spin quantum number, E BS is the energy of the BS state, E HS is the energy of the high-spin (HS) state which is obtained from a separate energy calculation at the equilibrium geometry of the BS state.