Thermally accessible triplet state of π-nucleophiles does exist. Evidence from first principles study of ethylene interaction with copper species†
Abstract
Three different models of ethylene interaction with copper species, namely, the Cu(100) surface, odd-numbered copper clusters C2H4/Cun (where n = 3, 7, 11, 15, 17, 19, 21, 25 and 27) and atomic copper C2H4/Cu were studied theoretically. It was found that the ethylene molecule possesses three different types of bonding depending on the presence of the unpaired spin on the reacting copper atom. These bonding structures demonstrate different types of band gap (bulk) or SOMO–LUMO gaps (cluster/atom), where SOMO stands for the singly occupied and LUMO means the lowest unoccupied molecular orbitals of the copper species. The obtained results are in good agreement with the previous experimental and computational results on the structural, spectral and energetic properties of the studied species. The bulk copper and sub-nanosized clusters (n > 7) build up the mono-π-bonded ground state complexes with ethylene where the latter species possesses the C2v symmetry. The single-atom complex C2H4/Cu forms the CS-symmetrical ground state 2A′ and the excited 2B2 and 4B state complexes of the C2v and C2 symmetry, respectively. The 2A′ state complex is mono-σ-bonded and involves the singlet ethylene moiety. The more tightly bound excited 2B2 complex has the di-σ-bonded structure and corresponds to the triplet ethylene. The adiabatic energy difference between the 2B2 and 2A′ states is equal to 10.8 kcal mol−1 and can be ascribed to the singlet–triplet splitting of the ethylene moiety interacting with copper. The QTAIM analysis supports the coordination type of the Cu–C bonds in all the studied complexes. Formation of the C2H4/Cu(100), C2H4/Cun and C2H4/Cu species is in accord with the well-known Dewar–Chatt–Duncanson model, in such a way that the opposing σ-donation step yields the ground state complex (2A′), while the subsequent more expensive supporting π*-back donation step provides the excited 2B2 state complex. In the present paper we have developed a computational procedure to optimize the latter complex.