Next Article in Journal
Cobalt Based Catalysts on Alkali-Activated Zeolite Foams for N2O Decomposition
Next Article in Special Issue
Synthesis of Ethylene/1-Octene Copolymers with Ultrahigh Molecular Weights by Zr and Hf Complexes Bearing Bidentate NN Ligands with the Camphyl Linker
Previous Article in Journal
Biodiesel: Modified Mixed Oxides as Catalyst for Transesterification of Rapeseed Oil
Previous Article in Special Issue
Adjusting Ortho-Cycloalkyl Ring Size in a Cycloheptyl-Fused N,N,N-Iron Catalyst as Means to Control Catalytic Activity and Polyethylene Properties
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Catalysts
Volume 10
Issue 12
10.3390/catal10121396
catalysts-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Reactivity and Structures for the Neutral and Cationic Bis(imino)pyridyl Iron and Cobalt Species by DFT Calculations

by
Zilong Li
1,2,
Yanping Ma
1,* and
Wen-Hua Sun
1,2,3,*
1
Key Laboratory of Engineering Plastics and Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
2
CAS Research/Education Center for Excellence in Molecular Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
3
State Key Laboratory for Oxo Synthesis and Selective Oxidation, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou 730000, China
*
Authors to whom correspondence should be addressed.
Catalysts 2020, 10(12), 1396; https://0-doi-org.brum.beds.ac.uk/10.3390/catal10121396
Submission received: 9 November 2020 / Revised: 27 November 2020 / Accepted: 29 November 2020 / Published: 30 November 2020
(This article belongs to the Special Issue Catalysis in Plastics for the 21st Century)
Browse Figures
Versions Notes

Abstract

:
Density Functional Theory (DFT) method was adopted to investigate and compare the reaction mechanisms of ethylene polymerization catalyzed by neutral, cationic bis(imino)pyridyl (PDI) iron and cobalt derivatives. The electronic structure and the oxidation states of the metal center and the PDI ligand were analyzed by taking spin states, natural bond orbital (NBO) charge distribution, etc. into consideration, revealing that the reactivity is closely related to the valence electron numbers instead of the charge numbers. The neutral Co(0) had the lowest reactivity as it possessed the most electrons. During the formation of the cationic Co(+)/Fe(+), one electron was mainly lost from PDI ligand rather than the metal center while the metal center maintained +II valence state through the process. Moreover, a special unsymmetrically bidentate N^N coordination manner was found to provide the deficient metal surroundings with 14e, which may initiate the reactivity of some unsymmetrical species with rich electrons. Finally, an anion [AlMe4]− participating process was proposed to explain the presence of the experimentally observed LCo(+)B(C2H4). A special intermediate, Co(+)B(C2H4) [AlMe4]− with Co in +I and absence of Co–C σ bond, was obtained. These calculation results may provide fundamental information for further understanding and designing the ethylene polymerization catalysts.

Graphical Abstract

1. Introduction

Much attention was focused on late transition metal complexes by researchers and engineers since Brookhart’s and Gibson’s research groups independently found that bis(imino)pyridine cobalt and iron complexes had high activities toward ethylene polymerization in 1998 [1,2]. Against this backdrop, the late transition metal complexes have made tremendous progress on ethylene polymerization by modifying the catalyst structures [3,4,5,6,7,8,9,10]. The introduction of a carbocyclic ring to the pyridine backbone [11,12,13,14,15] and the variation of substitution pattern in the N-aryl groups [16,17,18,19,20,21,22,23,24,25] were found as effective ways to improve the catalytic performance. Numerous studies focusing on the relationship between catalyst structure and their catalytic activity have gradually been reported [26,27,28,29,30,31,32,33,34,35,36,37].
Spectroscopic techniques such as Electron Paramagnetic Resonance (EPR) and Nuclear Magnetic Resonance (NMR) were developed to help understand the nature of the active species in details [26,27,28,29,30,31,32,33,34]. Gibson et al. revealed that the iron centers in bis(imino)pyridyl iron(II) pre-catalysts(LFeCl2) were oxidized upon treatment with methylaluminoxane (MAO) to provide a Fe(III) active species(LFe) [26]. In contrast, for the cobalt catalytic systems [27,28,29], the (iPrPDI)Co(II)Cl2 (iPrPDI = (1E,1′E)-1,1′-(pyridine-2,6-diyl)bis(N-(2,6-diisopropylphenyl)ethan-1 -imine)) was firstly reduced to (iPrPDI)Co(I)Cl by an alkylating agent (MAO), then Co(I) alkyl species (LCoCH3) was obtained. Talsi and coworkers [30,31,32] remarked the catalytic behavior of the bis(imino)pyridine iron and cobalt dihalide precatalyst upon addition of MAO or AlMe3, demonstrating the different types of active sites precursors present in different catalytic systems. Erker and colleagues [33] found that in the (ArPDI)CoCH3/Li[B(C6F5)4]/ethylene system (ArPDI = [bis(N-aryliminoethyl-κN,N′)pyridine-κN]), Li+ ions may interact with the “arms” of chelated donor ligands forming a reactive cationic reactant on the open coordination site. According to Talsi’s results [34], the interaction of L3MeCo(II)Cl2 and AlMe3 at −10 °C reduced cobalt (II) (L3MeCo(II)Cl2) into a cobalt (I) complex (L3MeCo(II)(μ-Me)(μ-Cl)AlMe2). Furthermore, in the presence of ethylene, newly formed complexes, such as (iPrPDI)Co(I)C2H4+, were also observed. However, up to now, the nature of the active species and the reaction mechanism of this special PDI-based iron or cobalt species were still in dispute.
In addition to experimental researches on the catalytic mechanism of the cobalt and iron catalysts based on a PDI backbone, many studies on theoretical calculations that overlay the mechanism were also reported [35,36,37]. In 1999, Gibson et al. [38] performed full ab initial computations on the bis(imino)pyridyl iron catalytic system for ethylene polymerization, and identified the key structure of (ArPDI)Fe(II)CH3C2H4 as the product of the first monomer insertion into the cationic alkyl using the BP86 functionals. While just the singlet multiplicities was considered during calculation, this was proved to not be accurately described [39]. Later, Morokuma etc. theoretically investigated the mechanisms of chain propagation and β-hydride transfer (BHT) by considering three potential energy surfaces of the catalyst model, in which two cationic bis(imino)pyridyl-Fe(II) systems were explored with different steric effect [40]. They found that the BHT pathway could be suppressed by increasing the steric bulk and the chain propagation took place mainly on quintet or triplet surfaces. In order to explore the effect of cocatalyst in the catalytic process, Burin and coworkers [41] chose different additives including (AlMeO)6(AlMe3), (PhOAlMe2)2 and TMA and calculated the activation and termination process of ethylene polymerization catalyzed by bis(imino)pyridyl iron catalysts using the unrestricted M06 functional. Altogether, the “alkyl-free” theory of polymerization catalyzed by bis(imino)pyridyl cobalt and polymerization mechanism have not been finalized and further researches are still necessary.
The reactivity differed due to the different charge or spin states of the active center. Herein, we used Density Functional Theory (DFT) to investigate and compare the reactivities of the cationic (ArPDI)MCH3+ and neutral (ArPDI)MCH3 (Ar = 2,6-Me2Ph, M = Fe, Co) of the iron and cobalt active species for the ethylene polymerization. Due to the limitation of the DFT calculations, just the following reaction processes were investigated: the first and second ethylene coordination and insertion, β-hydrogen transfer (BHT) to monomer and β-hydrogen elimination (BHE) reaction. Further, two kinds of spin states were considered. The differences of the electronic properties such as the NBO charge distribution, valence electron number and oxidation states and so on between the neutral and cationic iron or cobalt species were studied in detail. Moreover, in the case of the cobalt complex, a new reaction mechanism involving the cocatalyst ions (AlMe4−) was also simply analyzed.

2. Results and Discussion

To gain insight into the nature of the active species formed upon treatment of the corresponding metal dihalide complexes with MAO, various alkyl species based on bis(imino)pyridyl iron and cobalt complexes have been isolated in the experiments [2,27,42]. However, the activities and structures of the isolated species are related to the methodologies and additives which the researchers used during their alkylation synthesis. That means the actual electronic structures of the active centers varies in the charge distribution, spin states and the oxidation levels. Since the LMCH3 (M = Fe and Co, L = ArPDI) species are always considered as the initiated species towards ethylene polymerization, we select it as the model to explore the comparison between the iron and cobalt species by considering their neutral and cationic structures with different spin states. As shown in Scheme 1, the following reaction process are studied: the first coordination of the ethylene (Coord I), first insertion (Ins I), the second coordination (Coord II), second insertion (Ins II), β-hydrogen transfer (BHT) and β-hydrogen elimination (BHE). Comparison of energy curves (such as ethylene insertion and β-hydrogen transfer) demonstrated that there are some differences between the behavior of iron and cobalt complexes in polymerization, and also the spin states. The energy interval between the energy level of π-complex and sum of ethylene and catalyst energies could give the monomer coordination energy. Thus, the bond order of the special LCoCH3 and LCoCH3-C2H4(2) and orbitals of the reaction intermediates were analyzed. Additionally, the charge distribution of the system and the oxidation state for the PDI ligand were discussed. Finally, the effect of cocatalyst ions on the ethylene reaction process by cobalt system was also considered and compared.

2.1. Energy Profiles and Structures

For clear description, the symbol sM(q)-i was used to represent the intermediates in Scheme 1, where the ‘M’ could be the iron or cobalt metals, the ‘s’ is the spin quantum number, ‘q’ denotes charges and the ‘i’ corresponds to the intermediates in Scheme 1. For example, the initiated LCoCH3 will be 1Co(0)-1 when it is neutral in the triplet.

2.1.1. Energy Profiles and Structures of the Neutral LCoCH3 System

Figure 1 shows the relative zero-point energy profile of the corresponding reaction in Scheme 1 for the neutral cobalt species, sCo(0). The optimized structures of all the points in Figure 1 are listed in Figures S1 and S2. The energies of the reactants: 1Co(0)-1 + C2H4 were taking as the starting point, as the triplet 1Co(0)-1 has lower energy than the singlet 0Co(0)-1. Although both the singlet and triplet states were considered, the energies of most of the intermediates and transition state (TS) structures in triplet state 1Co(0) are lower than the corresponding singlet 0Co(0). The diffidence is 13.5~15.9 kcal/mol for the intermediates and 3.2~8.9 kcal/mol for the TS states. Although the energy differences between 0Co(0)-1 and triplet 1Co(0)-1 (1′ in Figure 1) is not much, it seems that the behavior of neutral cobalt species tend to be in high spin states instead of the low spin state. However, the BHE process is a bit different. The transition state Co(0)-TS3/9 for the hydrogen transferred to the cobalt metal center, is lower in energy for its singlet state than triplet state. Thus, spin crossover behavior exists in the BHE process of the neutral cobalt alkyl complexes. This phenomenon has been reported by Holland in his special concentration on the BHE process [43].
The binding energy for the initiated 1Co(0)-1 with the ethylene is −2.86 kcal/mol, resulting 1Co(0)-2. Then, by overcoming a considerable overall energy barrier of 21.25 kcal/mol, the coordinated ethylene inserts into the Co–C bond in 1Co(0)-2, producing the cobalt alkyl complex 1Co(0)-3 (LCoC3H7) with releasing energy of −23.62 kcal/mol. The second ethylene insertion (Ins II) starts from 3 reacting with the second ethylene, which was separated from the Ins I in Figure 1 for clarity. The energy barrier of Ins II is much higher than the Ins I. However, it could be compensated partly by the releasing energy from the formation of 1Co(0)-3, then the overall barrier for Ins II will be 9.50 kcal/mol. Additionally, in comparison with the first one (−2.86 kcal/mol), the coordination of the second ethylene to 1Co(0)-3 is unfavorable (binding energy of 7.47 kcal/mol). These results indicate that, in the ethylene polymerization catalyzed by the active species with similar electronic properties to the neutral Co(0)-1, the first ethylene insertion process could be the rate-determining step with exothermic energies ready for further reactions.
As shown in Scheme 1, after the formation of 1Co(0)-3, three reaction pathways could proceed: second insertion (Ins II) for chain growth, β-hydrogen transfer (BHT) and β-hydrogen elimination (BHE) for chain terminations. BHE just requires transferring the β-hydrogen to the metal center while Ins II and BHT need to coordinate another ethylene firstly. Obviously, the BHE process is more advantageous than BHT and Ins II. The overall barrier of BHE in the triplet state is 4.37 kcal/mol, while it is only 1.99 kcal/mol in the singlet TS via a spin crossover. Then the singlet BHE complex 0Co(0)-9 is provided. Since the overall barrier of the BHT is much higher than that of BHT and Ins II reaction paths, it will not be considered as a termination process for the neutral cobalt alkyls species. So BHE process is kinetically favored over the thermodynamically favored Ins II to form the propagated cobalt alkyl complexes by 7.5 kcal/mol. For the bis(imino)pyridyl metal species with low steric substituents, the BHE is competitive chain termination mechanism for producing short oligomer chains with unsaturated end group [44].
Figure 2 shows the binding information of the ethylene coordinated 0Co(0)-2, 1Co(0)-2, 1Co(0)-4 and 1Co(0)-6. The binding energy ΔE is defined as ΔE = E2—(E1 + EC2H4) where E2 refers to the zero point energies of the coordinated species 1Co(0)-2 or 1Co(0)-4 etc., and E1 is for the precoordinated 1Co(0)-1 or 1Co(0)-3. It is found that the distance between the cobalt atom and the coordinated ethylene in 0Co(0)-2 is smaller than 1Co(0)-2 complex (2.08, 2.08 Å vs 2.43, 2.33 Å). The double bond of the coordinated ethylene is longer in the 0Co(0)-2 (1.39 Å) than 1Co(0)-2 (1.35 Å). These exhibited stronger interactions in 0Co(0)-2 than in 1Co(0)-2. However, the binding energy of 1Co(0)-2 has negative values (−2.86 kcal/mol) while it is positive in 0Co(0)-2 (3.63 kcal/mol). Asymmetrical coordination was found between the bis(imino)pyridyl ligand and the cobalt metal center in 1Co(0)-2, as may contribute to the more binding energies than 0Co(0)-2 The asymmetric coordination could be easily observed from the distances between cobalt metal center and the nitrogen atoms in the bis(imino)pyridyl ligand. They are 2.24, 1.95 and 2.24 Å in 1Co(0)-1, while they turn to 2.71, 2.03 and 2.08 Å in 1Co(0)-2. The same phenomena was also observed in the ethylene coordinated 1Co(0)-4 and 1Co(0)-6, which is the intermediates for the further Ins II or BHT reactions. It is special to find this kind of asymmetrical coordination between the PDI ligand and the metal atom in the neutral triplet cobalt species. Therefore, the PDI ligand is bidentate rather than tridentate when interacts with the cobalt (Figure 2). As reported for the structure of the iron species, it is generally considered that the active center is 14-electron alkyl divalent cation structure [42,45]. It can be speculated that the cobalt species may also become a 14-electron structure by forming the bidentate coordination structure during the process of catalyzing ethylene polymerization. Chirik et al. [10] speculate that when bis(imino)pyridyl cobalt catalyzes the [2 + 2] cycloadditions, the cobalt will de-coordinate with one of the nitrogen on both sides, which is consistent with the optimized intermediate structure 1Co(0)-2 and 1Co(0)-4. This special intermediate may provide new electronic properties for the interactions between the cobalt complexes and the coordinated ethylene.
In the transition state (TS) structures of ethylene insertion 0Co(0)-2/3 and 1Co(0)-2/3 (see Supporting Information Figures S1 and S2), the distance between the cobalt atom and γ-hydrogen on the connected polymer chain indicates an agnostic interaction leading to stabilization of structure. Further, the bidentate interactions recover to the tridentate manner during the transition state process.
To further clarify the special unsymmetrically bidentate NN coordination character in the π-complexes 1Co(0)-2, Mayer bond order analysis was performed by Multiwfn [46]. The bond order of the bonds in the active center for 1Co(0)-2 and 0Co(0)-2 are listed in Table 1 where the concerned atom numbers are illustrated in Scheme 2. The greater the value of bond order, the stronger the bond between atoms. Usually, the bond order for one π-bond equals two, and for one σ-bond equals one.
As shown in Table 1, the bond order value of most carbon–carbons bond is greater than 1.4, but differs between the 0Co(0)-2 and 1Co(0)-2. The bond orders are 1.41/1.42 in the aromatic pyridine ring (C(1)~C(5)) for 0Co(0)-2, indicating well delocalized π-bonds. In contrast, less conjugation degree exist in the pyridine ring of 1Co(0)-2, which could be derived from the broad range of the bond order value (1.33~1.49) in the ring.
As shown in Scheme 2, the C(11)–N(12) and C(10)–N(13) bonds could be considered as two ‘arms’ in bis(imino)pyridine skeleton. Both ‘arms’ and pyridine ring in 0Co(0) have symmetrical bonding properties and exhibit good delocalized π-bond properties. However, unsymmetrical and localized bonding character were found in 1Co(0)-2. The bond order of C(1)–C(10) (1.23) and C(10)–N(13) (1.44) on one ‘arm’ are different from those of C(5)–C(11) (1.08) and C(11)–N(12) (1.82) on the other ‘arm’ in 1Co(0)-2, while they are almost the same in the 0Co(0)-2.
Further, the interactions between the cobalt metal and the three nitrogen atoms are intramolecular coordination bonds in which the bond order is mostly less than 0.5. Among the three Co---N interactions, stronger values were found between the central nitrogen on pyridine and the cobalt metal center. Additionally, they are weaker in 1Co(0)-2 than in 0Co(0)-2. The binding levels of N(6)---Co(66), N(12)---Co(66) and N(13)---Co(66) are 0.42, 0.11 and 0.39 in 1Co(0)-2, but they are 0.61, 0.47 and 0.48 in 0Co(0)-2. The weak bond of the N(12)---Co(66) (0.11) on one arm of the 1Co(0)-2 suggests that there is almost no coordination interaction between the N(12) and the cobalt metal. As a result, the cobalt atom has more interaction with one of the ‘arms’ (C(5)–C(11)–N(12)) than the other, and the delocalization degree of this arm is smaller. Thus, the PDI ligand interacts with the cobalt in 1Co(0)-2 by a bidentate coordination instead of the normal tridentate. By this way, the deficient 14-electron metal center may be created to inspire the activity of neutral cobalt LCo-R towards ethylene insertion.
In order to clearly demonstrate and compare the electronic properties of the 1Co(0)-2 and 0Co(0)-2, their HOMO orbitals were obtained by Multiwfn (Figure 3). The distribution of the HOMO orbitals of singlet 0Co(0)-2 are well symmetrically delocalized in the structures. In contrast, the distribution in triplet 1Co(0)-2 is asymmetric about the axis of the cobalt atom and the central nitrogen atom on the pyridine. This may be due to the fact the two unpaired electrons in 1Co(0)-2 arise the unsymmetrically bidentate NN coordination character.

2.1.2. Energy Profiles and Structures of the Cationic LCo+CH3 System

As shown in Figure 4, zero point energy profiles of the cationic quartet (3/2Co(+)) and the doublet (1/2Co(+)) cobalt species are expressed by setting the energies of the reactants: 1/2Co(+)-1 and ethylene as the zero point. The optimized structures of all the points in Figure 4 are listed in Figures S3 and S4. For the LCo(+)CH3 (Co(+)-1) species, the doublet 1/2Co(+)-1 is more stable in energy than the quartet 3/2Co(+)-1 by 1.16 kcal/mol. Differently, the ethylene coordinated Co(+)-2 is more stable in the quartet state than the doublet. It is interesting to find that for the cationic cobalt alkyl species, all the coordinated intermediates such as Co(+)-2, Co(+)-4, Co(+)-6, Co(+)-7 have lower energy in quartet than in doublet state. However, the transition process of Ins I and II prefer to proceed in doublet state instead of quartet. Overall, except for the two transition states Co(+)TS4/5 and Co(+)TS6/7, the energy differences between the quartet and the doublet states are not much within 3.5 kcal/mol. This means that for the cationic cobalt alkyl species, the spin crossover may occur throughout the reaction pathway considered in Scheme 1.
It is obvious that the overall energy barrier of the Ins I (19.31 kcal/mol) is higher than Ins II (4.24 kcal/mol), BHT (7.59 kcal/mol) and BHE (2.39 kcal/mol), predicting the Ins I to be the rate determining step.
In comparison with the neutral cobalt species (Figure 1), the overall barrier for the process considered in Scheme 1 has lower energies for the corresponding cationic species (Figure 4). However, the energy released from the formation of the cobalt alkyl species is greater in the neutral system. For example, the energies of the overall transition barriers for Ins I and II in the neutral system are 21.25 and 9.5 kcal/mol, respectively, while they are 19.31 and 4.24 kcal/mol in the cationic. The formation of the neutral 1Co(0)-3 and 1Co(0)-5 are −23.62 and −42.74 kcal/mol, respectively, while they are −15.10 and −34.53 kcal/mol for 3/2Co(+)-3 and 3/2Co(+)-5. So, the cationic system may have better reactivities than the neutral one. The energy difference between the favorable BHE and the second favorable process Ins II is just 1.85 kcal/mol in Co(+) species while it is 7.51 kcal/mol in Co(0), indicating that the cationic species are more favorable for the chain propagation to generate polymers. The energetically unfavorable BHE in cationic could also facilitate the chain growth to polymers. For either the neutral or the cationic system, the Ins I is the rate-determining step.
Similar to the neutral π-coordinated species, the 3/2Co(+)-4 (quartet 4′ in Figure 4) and 3/2Co(+)-6 (Figure S4) exhibits slightly asymmetry between metal center and nitrogen atoms in the bis(imino)pyridyl ligand.

2.1.3. Energy Profiles and Structures of the Cationic LFe+CH3 System

Previous theoretical works on the PDI-derived catalysts are mostly for cationic iron species. Ziegler examined the reactivity of the bis-ortho aryl substituted [(PDI)Fe(C3H7)]+ complex towards ethylene polymerization by BP86 [35]. Along with transition states for Ins, BHT and BHE, metal−alkyl and ethylene π coordination structures were reported. However, only singlet spin state energy profile was reported [40]. Later, Musaev reported the same structures and catalytic behaviors by B3LYP, in which high spin-state structures are more favorable. They also found that the BHT is favorable than the Ins process when the ligand steric is low. By B3LYP methods, the reactivity of the [(PDI)Fe(CH3)]+ was studied by de Bruin [47]. It was suggested that the termination via the BHT transition state has around 10 kcal/mol lower activation enthalpy than ethylene insertion. Further, the [(PDI)Fe(CH3)]2+ with iron in +III oxidation state has a much lower insertion barrier than BHT.
In order to compare the different charge effect between iron and cobalt complexes the computational study was carried out on both of the neutral and cationic bis(imino)pyridine (PDI) iron complexes. The energy profiles for the cationic PDI iron species were depicted in Figure 5. Further, the optimized structures were listed in Figures S5 and S6. For most of the points shown in Figure 5, the quintet state is more stable in energy than the triplet state. The transition state of Ins II located on the triplet spin state is consistent with the results reported by Musaev who has examined the ethylene insertion to the [(PDI)-FeC3H7]+ structure [47].
Similar to the cobalt species, Ins I process is also the rate determining step. Close overall transition barriers were found for the three reaction pathways of Ins II, BHT, and BHE in the cationic iron species, in which the energy difference is just 3.14 kcal/mol. It seems that the BHE has lower transition energy. However, BHE is an energetically unfavorable process producing 2Fe(+) and C3H6, with endothermal energy of 2.75 kcal/mol.
Despite the transition barrier is not high, the BHE process of both cationic cobalt and iron species occurs with a bit endothermal energy. The overall energy of the rate-determining transition barrier is 19.35 kcal/mol in Fe(+) while it is 19.31 kcal/mol in Co(+). Identical reactivity may exist in the cationic iron and cobalt systems. However, it is difficult for the active species to possess strictly positive charge in the actual catalytic process, since the anion ions interact with them more or less in the real experiments.

2.1.4. Energy Profiles and Structures of the Neutral LFeCH3 System

As shown in Figure 6, the energy profiles for the neutral iron species Fe(0) were illustrated and the optimized structures were listed in Figures S7 and S8. Almost all the quartet state species is more stable in energy than doublet state species, indicating that for neutral bis(imino)pyridyl (PDI) iron complexes, doublet state species have an overwhelming superiority.
The barrier for Ins I is also the rate determining step for Fe(0). The value for it is 16.13 kcal/mol, as it is lower than either neutral Co(0) or cationic Co(+) and Fe(+). Different from Fe(+) and Co(0), the BHT process for Fe(0) is more favorable than the BHE. Overall, in comparison with Ins II and BHE, the BHT is the most preferred reaction paths. Similar with neutral Co(0), the neutral iron alkyl intermediates Fe(0)-3 and Fe(0)-5 are very stable with considerable exothermic energies indicating that the neutral systems tend to form stable metal alkyls. Moreover, the binding energies of 3/2Fe(0)-1 with ethylene to form 3/2Fe(0)-2 is −5.74 kcal/mol, which is more stable than that of 1Co(0)-2. However, it may be due to the one less electron in the 3/2Fe(0)-2 than the 1Co(0)-2, the unsymmetrical coordination has not been observed in the neutral iron system.

2.2. General Electronic Properties

2.2.1. Oxidation Level of the PDI Ligand

Wieghardt and colleagues found a structural parameter (Δ) to determine electron exchange between the metal and the bis(imino)pyridyl ligand along with the oxidation level of metal according to the following Equation (1) [29];
Δ = [(d2 + d2’) /2 − (d1 + d1’ + d3 + d3’) /4]
Each factor was defined as depicted in Scheme 3. Based on the equation, we found the changes linearly with the reduction of (PDI0) ligand. This finding could be helpful for the electronic exchange definitions between the metal and bis(imino)pyridine skeleton structure.
To figure out how the electronic exchange occurred between the metal and the bis(imino)pyridine skeleton structure, we calculated the Δ value of the structures for the Ins I and II process and listed in Figure 7 (for more details see Supporting Information in Tables S1–S8).
It is necessary to have an overall cognition about the structure parameter Δ of iron and cobalt species before analyzing the data in Figure 7. It is characteristic of neutral (PDI)0 ligands when the values are in the range of 0.15–0.19 Å (av 0.17 Å). The Δexp values in the range 0.09~0.12 Å (average, 0.102 Å) are typical for a π-radical anion (PDI)1−. Structures with an average Δexp value of 0.062 ± 0.02 Å have been found for complexes containing a dianion (PDI)2−. Having this in mind, the calculated Δ in this work will reflect the oxidation level of the PDI ligand.
As illustrated in Figure 7, it is clear that for most cationic systems, the calculated values Δ are located in the range of 0.15~0.19 Å, which indicates that the PDI ligands are neutral. Then, the oxidation state of the metal center will be +II, which is normally accepted in the related reports. There are some differences between the high and low spin states for the cationic species, in which the latter has lower values than the former. This means that the interactions between the PDI ligand and metal center are stronger in the low spin state than the high spin state. Further, it could also be supported from the shorter bond in optimized structures in Figures S1–S8. In comparison with the iron system, the values of the cationic cobalt species are always greater. Especially for the low spin cationic species, the values of iron system are in the range of 0.144~0.159 while they are 0.165~0.174 for cobalt. The values for the high spin states are almost the same for the iron and cobalt species.
For most of the neutral iron and cobalt structures, the Δ values are within 0.09~0.13 Å. Theoretically, that the PDI ligand is an anion revealing the +II oxidation level of the metal center in the neutral LMCH3 structure, which is exhibits as +I in appearance. According to the results mentioned above, the ground state of these neutral species is high spin state. It is similar to the cationic species, the Δ values for high spin states are almost the same for iron and cobalt system. However, lower values are found in 0Co(0) structures. Particularly, the Δ values of the transition barrier for the Ins I and Ins II in 0Co(0) are rather low. They are 0.061 and 0.058, respectively, suggesting that PDI is dianion during the insertion step.
Overall, in comparison with the cationic species, the one more extra electron in the neutral species may mainly transfer to the PDI ligand instead of the metal center. By this way, the metal center does always keep the +II valence in either cationic or neutral species.

2.2.2. Charge Distribution

In order to figure out the difference of the charge distribution among sM(q)-1 and the ethylene coordinated sM(q)-2, we do the electronic population analysis by natural bond orbital (NBO) method [48]. Firstly, the model molecule in Scheme 2 is divided into four moieties: PDI ligand, M center, CH3 and two Ars (aromatic rings connected to the two nitrogen atoms in PDI ligand). The NBO charge distribution of sM(q)-1 was illustrated in Figure 8a and Table S9. The charge changes aroused by the coordination of the ethylene in sM(q)-2 were listed in Figure 8b.
As shown in Figure 8a, for both of the cationic and neutral species, the metal center and the two Ars possess positive charges while the CH3 moiety that binds to the metal has negative values. The positive values for the two Ars decrease slightly from the cationic species (0.492~0.515) to the neutral ones (0.327~0.385). Similarly, the fluctuation of the negative charges in the CH3 moiety is also not great. Additionally, it is obvious that the cationic complex has more positive charges than the neutral. The NBO charge distributions here could strongly support the conclusion that the PDI ligand is almost neutral in the cationic species while it is anion in the neutral species. The charge values of PDI in both of the neutral iron and cobalt species are close to zero while they are the most negative moiety in the neutral species. Therefore, the loss of one electron from the neutral system is mainly from the PDI ligand instead of the metal center.
After coordination of ethylene molecules, the charge distribution changes slightly in the cationic species while it changes more in the neutral species. In sM(q)-2, the C2H4 moiety has positive NBO charges in both cationic and neutral systems, indicating that some of the electrons from C2H4 are transferred to other moieties, mainly to the metal center. Additionally, the electrons transferred from C2H4 moiety are less in the neutral sM(0)-2 than the cationic sM(+)-2. Besides the metal center, there are significant changes of charges on the CH3 and PDI moieties in the neutral due to the coordination. For 3/2Fe(0)-2, the iron atom become more positively charged while CH3 moiety is more negatively. In contrast, for the 1/2Fe(0)-2, 1Co(0)-2 and 0Co(0)-2, the metal atoms receive less positive charges while CH3 less negative.

2.2.3. General Comparison on Valence Electrons

For all the four species, Co(+), Co(0), Fe(+), Fe(0), the first ethylene insertion step is the rate determining one. Additionally the energies of the transition barriers increase in the sequence of Fe(0) (16.13 kcal/mol) < Co(+)(19.31 kcal/mol) < Fe(+)(19.35 kcal/mol) < Co(0)(21.25 kcal/mol), though the difference is not much, within 5.12 kcal/mol. The neutral cobalt system Co(0) has the lowest reactivity because it possesses the most more electrons; however, the Fe(+), which lacks the electrons most, is not the most active one here. As for the spin states of the Co(+), Fe(0), and Fe(+), Ins I and Ins II process have lower energy in the high spin states, but spin crossover may exist for the BHT and BHE process in most cases. Additionally, the spin crossover is almost throughout the Co(+) species. BHT or BHE is the kinetically favorable process in comparison with the Ins II as they have a low steric substituent, but they are always energetically unfavorably. The Co(0) and Fe(+) tend to undergo hydrogen transfer through BHE for the possible chain termination, while Fe(0) and Co(+) through BHT.
As demonstrated in Scheme 4, the valence electron numbers were considered generally. The unpaired electron numbers are the same for Co(+) and Fe(0), which is one electron less, and better reactivity, than the Co(0). Additionally, the unpaired electrons are located mainly at the cobalt metal center for Co(+) whereas those in Fe(0) are distributed in both of the iron metal and the PDI ligand. This may explain the occurrence of more spin crossover in Co(+) system. Additionally, the loss of one electron from the Co(0) to form the Co(+) is mainly from the PDI ligand instead of the cobalt metal, as is the same for Fe(0) to Fe(+). The better reactivity of neutral Fe(0) may be attributed to the almost 14e deficient electrons around the iron atom. For the cobalt species, it could also reach the 14e center by losing the CH3 anion group from Co(0) to afford the Co(+)B (the symbol B to mark the bare nature or no alkyl coordination of the cobalt center) which could be supported by the observation of the alkyl free intermediate in the experiment [27,28]. In the next part, the coordination of ethylene or the cocatalyst anion groups to the Co(+)B will be analyzed for the possible reaction mechanism.
It is worth mentioning that the numbers of the valence electrons in Co(0) is 16e, but when one of the nitrogen atoms in the PDI lost interactions with the metal center (as shown in Scheme 4), the 14e deficient center may be formed by an unsymmetrical manner. The good activity in the single cyclopenta ring fused PDI cobalt system towards the ethylene oligomerization in the experiments [49] could be attributed to the NNN unsymmetrically loose coordination manner in the metal center.

2.3. Possible Catalytic Process for the Cobalt Species

The electron deficient metal alkyl species [LM-R]+ was well accepted mode for both early and late transition metal complexes when LnMX2 activated by MAO. However, Gibson [27,28] reported that the active species for cobalt species is anticipated to be the [LCo]+ instead of the cobalt alkyl cation [LCo-R]+. Strong coordinating interactions were found for [LM]+ cation with the anions [Cl-MAO]− or [Me-MAO]− derived from cocatalysts. In the NMR investigation by Talsi et. al [34], neutral complexes of LCo(μ-Me)2AlMe2 were observed in the catalyst systems LCoCl2/AlMe3. However, in the presence of ethylene, the proposed structure of [LCo(η-C2H4)]+[AlMe3Cl]− is the major cobalt species. The polymerization manner catalyzed by the cobalt complex initiated with cocatalysts is still unclear. Even the exact nature of MAO and MMAO remains somewhat unclear, previous computational work showed that explicit inclusion of a MAO model did not significantly alter energy and structures for modeling ethylene oligomerization [50]. Thus, the [AlMe4]− was considered as the anion ions that function as the cocatalysts.
To further explore the intrinsic mechanism of ethylene reactions catalyzed by the cobalt species, as depicted in Figure 9 and Scheme 5, we contrast the energy profiles of the reactions of LsCo(+)B with C2H4 and anion ions [CH3]− and [AlMe4]−, respectively. Only the first ethylene insertion process was considered, as it appears to be the rate-determining step in the above calculated results. Both of the single and triplet states were considered with the optimized structures listed in Figures S9 and S10. The triplet is lower in energy than the singlet spin state. It was found that the reactions of both of the anion ions [CH3]− and [AlMe4]− with the LsCo(+)B release plenty of energy, generating the neutral species 1A and 1C. The coordination of the [AlMe4]− to L1Co(+)B is more stable in energy than the [CH3]− by 15.48 kcal/mol. This indicates that the similar anion ions [Cl-MAO]− or [Me-MAO]− derived from the cocatalysts may prefer to interact with the cobalt center closely to form the neutral species, instead of resulting the ion pairs with weak interactions.
In comparison with anion ions [CH3]− and [AlMe4]−, the coordination of C2H4 to the bare cationic L1Co(+)B produces less energy, but the formation of the L1Co(+)C2H4 (1B) is still considerably stable in energy by 16.27 kcal/mol. The electron deficient 14e L1Co(+)B which is electron deficient has better binding abilities than either the cationic L3/2Co(+)CH3 (3.62 kcal/mol) or the neutral L1Co(0)CH3 (−2.86 kcal/mol) to coordinate the ethylene monomers. This may explain the experimental observation of the LCo(η-C2H4)]+. Additionally, for the process of the ethylene insertion, the energy barrier for the [AlMe4]− participating one is lower than the [CH3]− by 5.46 kcal/mol.
As mentioned above, the loss of electrons between the neutral and cationic species mainly caused the reduction of the oxidation state of the PDI ligand instead of the metal center. This means the metal center keeps the +II valence state in both the neutral and cationic systems. The oxidation state of the PDI ligand in the processes in Figure 9 was monitored and compared by calculating the structure parameter Δ. As listed in Table 2, similar to the cationic sCo(+)-1, the Δ value of L1Co(+)B is 0.181 Å, indicating the neutral character of the PDI ligand. Consequently, the valence state of the cobalt metal center in L1Co(+)B is +I instead of +II. The neutral PDI ligand and the +I valence cobalt metal center were retained in the ethylene coordinated cationic L1Co(+)C2H4 (1B in Figure 9). However, the Δ values of 0.123 Å, 0.119 Å and 0.114 Å for 1A, 1C and 3A/3B indicated the presence of a monoanionic PDI ligand (more details see Tables S10 and S11). 2B could be generated from the 1B (neutral PDI and +I cobalt center) and 1C (monoanionic PDI and +II Cobalt center). Particularly, the calculated Δ value of the structure 2B appears to be 0.168 Å which is in the range of the neutral PDI. Additionally, there is no Co–C σ bond formed in 2B, as is different from 2A. In 2B, the carbon atom in methyl group forms σ bond neither with Cobalt nor with Aluminum. Thus, the bonding in 2B appears to be a 3-center 2-electron character. Additionally, the valence state of cobalt should be +I in 2B. Moreover, the interaction of the cobalt atom with the [AlMe4]− (Co---C([AlMe4]−), 2.39 Å) is less strong than with the coordinated ethylene (Co---C(C2H4), 2.02 Å). The chain growth in the cobalt species may proceed by attacking the carbon atom of the binding C2H4 from the methyl group in [AlMe4]−.

3. Computational Methods

Density functional theory (DFT) calculation was performed using the Gaussian 09 software package [51]. The geometric structures of the relevant species were optimized in the gas phase with the hybrid B3LYP which was used to produce reasonable results by other researchers [43,52]. In the investigation performed by Busico et al. [53], the B3LYP demonstrated moderate accuracy among the 22 functionals on transition states energies, agostic interactions, reaction energies and π-coordination energies for transition metal catalysts. The standard 6-311G(d) basis set was used for nonmetal atoms, while the SDD was chosen for the metal atom [54]. Vibration frequencies were calculated to check that the reaction intermediates were all positive and the TSs had only one imaginary frequency. The zero-point energy changes (ΔE) were used to discuss the catalytic reaction process. Additionally, possible spin states were considered for all the investigated systems. Bruin et al [41] have confirmed that the inclusion of solvent or expansion of the wave function in a larger basis set did not affect the calculations performance in the gas phase. Then, in this theoretical calculation solvent effect was not considered. Intrinsic reaction coordinate (IRC) calculations [55,56] were also performed to check a TS connects two appropriate local minima in some of the reaction paths. In addition, natural bond orbital (NBO) analysis [46,49,57,58] was conducted on the optimized structures to acquire a further insight into the charge distributions of the system [59]. Some of the bond order was also analyzed by the Mayer bond order analysis by Multiwfn [46].

4. Conclusions

The reactivities and electronic structures of the neutral, cationic bis(imino)pyridyl iron and cobalt species were explored and compared in this work by DFT calculations. Different spin states were considered for the four species Co(+), Co(0), Fe(+), Fe(0). The results revealed that for all the species, the rate-determining step is the first ethylene insertion. The neutral cobalt system Co(0) had the lowest reactivity as it possessed the most electrons. As for the spin states in Co(+), Fe(0) and Fe(+), the Ins I and Ins II processes prefer to take place in high spin states, but the BHT and BHE tend to occur along with spin crossover. Additionally, the spin crossover is almost throughout the Co(+) species. Overall, due to the low steric substituent in the investigated species here, the BHT or BHE is the kinetically favorable process in comparison with the Ins II, but they are always energetically unfavorable. The Co(0) and Fe(+) tend to perform chain termination through BHE, while Fe(0) and Co(+) through BHT.
Additionally, the electron numbers are the same for Co(+) and Fe(0), which is one electron less and better reactivity than the neutral Co(0). The unpaired electrons are located mainly on the cobalt metal center in Co(+) while they are on both of the iron metal center and the PDI ligand in Fe(0). The more electrons in the cobalt metal center may explain the occurrence of more spin crossover in Co(+) system. Additionally the loss of one electron from the Co(0) to form the Co(+) is mainly from the PDI ligand instead of the cobalt metal, as is the same for Fe(0) to Fe(+).
Moreover, a special unsymmetrically NN coordination manner was found between the PDI ligand and the cobalt metal in the neutral cobalt species Co(0). This may result in 14e deficient center. Through this center, the good activity in some asymmetric structures of electron-rich cobalt in the experiments may be explained. Furthermore, a process involving [AlMe4]− was proposed to explain the experimentally observed LCo(+)B(C2H4)(1B) which could be generated from the LCo(+)B with 14e deficient center by ethylene coordination. A special intermediate, Co(+)B(C2H4)[AlMe4]− with Co in +I and absence of Co–C σ bond, was obtained. In this intermediate, the interaction of the cobalt metal atom with the [AlMe4]− was less strong than with the coordinated ethylene. The chain growth in the cobalt species may proceed by attacking the carbon atom of the binding C2H4 from the methyl group in [AlMe4]−. These calculations results may provide fundamental information for understanding and designing the catalysts of ethylene polymerization.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4344/10/12/1396/s1, Figure S1 The DFT optimized structures of the neutral cobalt system 0Co(0) in Figure 1, Figure S2 The DFT optimized structures of the neutral cobalt system 1Co(0) in Figure 1, Figure S3 The DFT optimized structures of the cationic cobalt system 1/2Co(+) in Figure 3, Figure S4 The DFT optimized structures of the cationic cobalt system 3/2Co(+) in Figure 3, Figure S5 The DFT optimized structures of the cationic iron system 1Fe(+) in Figure 5, Figure S6 The DFT optimized structures of the cationic iron system 2Fe(+) in Figure 5, Figure S7 The DFT optimized structures of the neutral iron system 1/2Fe(0) in Figure 6, Figure S8 The DFT optimized structures of the neutral iron system 3/2Fe(0) in Figure 6, Figure S9 The DFT optimized structures of the cobalt system 0LCo(+)B in Figure 9, Figure S10 The DFT optimized structures of the cobalt system 1LCo(+)B in Figure 9, Table S1 The factors for calculating the structure parameter Δ in 0Co(0), Table S2 The factors for calculating the structure parameter Δ in 1Co(0), Table S3 The factors for calculating the structure parameter Δ in 1/2Co(+), Table S4 The factors for calculating the structure parameter Δ in 3/2Co(+), Table S5 The factors for calculating the structure parameter Δ in 1Fe(+), Table S6 The factors for calculating the structure parameter Δ in 2Fe(+), Table S7 The factors for calculating the structure parameter Δ in 1/2Fe(0), Table S8 The factors for calculating the structure parameter Δ in 3/2Fe(0), Table S9 The NBO charge distribution for sM(q)-1 and sM(q)-2, M includes iron and cobalt, s means the spin states and the q refers to the charge number, Table S10 The factors for calculating the structure parameter Δ in 0Co(+)B, Table S11 The factors for calculating the structure parameter Δ in 1Co(+)B.

Author Contributions

Data curation, Z.L.; Formal analysis, Z.L. and Y.M.; Funding acquisition, W.-H.S.; Investigation, Z.L.; Project administration, W.-H.S.; Supervision, Y.M.; Writing—original draft, Z.L.; Writing—review & editing, Z.L., Y.M. and W.-H.S. All authors have read and agreed to publish this version of manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 21871275).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Small, B.L.; Brookhart, M.; Bennett, A.M.A. Highly active iron and cobalt catalysts for the polymerization of ethylene. J. Am. Chem. Soc. 1998, 120, 4049–4050. [Google Scholar] [CrossRef]
  2. Britovsek, G.J.P.; Gibson, V.C.; Kimberley, B.S.; Maddox, P.J.; McTavish, S.J.; Solan, G.A.; White, A.J.P.; Williams, D.J. Novel olefin polymerization catalysts based on iron and cobalt. Chem. Commun. 1998, 849–850. [Google Scholar] [CrossRef]
  3. Zhang, W.; Wang, S.; Du, S.; Guo, C.-Y.; Hao, X.; Sun, W.-H. 2-(1-(2,4-Bis((di(4-fluorophenyl)methyl)-6-methylphenylimino)ethyl)-6-(1-(arylimino)ethyl)pyridylmetal (iron or cobalt) Complexes: Synthesis, Characterization, and Ethylene Polymerization Behavior. Macromol. Chem. Phys. 2014, 215, 1797–1809. [Google Scholar] [CrossRef]
  4. Ahmed, S.; Yang, W.; Ma, Z.; Sun, W.-H. Catalytic Activities of Bis(pentamethylene)pyridyl(Fe/Co) Complex Analogues in Ethylene Polymerization by Modeling Method. J. Phys. Chem. A 2018, 122, 9637–9644. [Google Scholar] [CrossRef] [PubMed]
  5. Huang, C.; Zakharov, V.A.; Semikolenova, N.V.; Matsko, M.A.; Solan, G.A.; Sun, W.-H. A comparative kinetic study of ethylene polymerization mediated by iron, cobalt and chromium catalysts bearing the same N,N,N-bis(imino)trihydroquinoline. J. Catal. 2019, 369, 1–9. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, Z.; Solan, G.A.; Zhang, W.; Sun, W.-H. Carbocyclic-fused N,N,N-pincer ligands as ring-strain adjustable supports for iron and cobalt catalysts in ethylene oligo-/polymerization. Coord. Chem. Rev. 2018, 363, 92–108. [Google Scholar] [CrossRef] [Green Version]
  7. Phillips, A.M.F.; Suo, H.; Guedes da Silva, M.d.F.C.; Pombeiro, A.J.L.; Sun, W.-H. Recent developments in vanadium-catalyzed olefin coordination polymerization. Coord. Chem. Rev. 2020, 416, 213332. [Google Scholar] [CrossRef]
  8. Bariashir, C.; Huang, C.; Solan, G.A.; Sun, W.-H. Recent advances in homogeneous chromium catalyst design for ethylene tri-, tetra-, oligo- and polymerization. Coord. Chem. Rev. 2019, 385, 208–229. [Google Scholar] [CrossRef] [Green Version]
  9. Suo, H.; Solan, G.A.; Ma, Y.; Sun, W.-H. Developments in compartmentalized bimetallic transition metal ethylene polymerization catalysts. Coord. Chem. Rev. 2018, 372, 101–116. [Google Scholar] [CrossRef] [Green Version]
  10. Chirik, P.J. Carbon-Carbon Bond Formation in a Weak Ligand Field: Leveraging Open-Shell First-Row Transition-Metal Catalysts. Angew. Chem. Int. Ed. 2017, 56, 5170–5181. [Google Scholar] [CrossRef]
  11. Zhang, W.; Chai, W.; Sun, W.-H.; Hu, X.; Redshaw, C.; Hao, X. 2-(1-(Arylimino)ethyl)-8-arylimino-5,6,7-trihydroquinoline Iron(II) Chloride Complexes: Synthesis, Characterization, and Ethylene Polymerization Behavior. Organometallics 2012, 31, 5039–5048. [Google Scholar] [CrossRef]
  12. Huang, F.; Xing, Q.; Liang, T.; Flisak, Z.; Ye, B.; Hu, X.; Yang, W.; Sun, W.-H. 2-(1-Aryliminoethyl)-9-arylimino-5,6,7,8-tetrahydrocycloheptapyridyl iron(II) dichloride: Synthesis, characterization, and the highly active and tunable active species in ethylene polymerization. Dalton Trans. 2014, 43, 16818–16829. [Google Scholar] [CrossRef] [PubMed]
  13. Sun, W.-H.; Kong, S.; Chai, W.; Shiono, T.; Redshaw, C.; Hu, X.; Guo, C.; Hao, X. 2-(1-(Arylimino)ethyl)-8-arylimino-5,6,7-trihydroquinolylcobalt dichloride: Synthesis and polyethylene wax formation. Appl. Catal. A 2012, 447–448, 67–73. [Google Scholar] [CrossRef]
  14. Huang, F.; Zhang, W.; Sun, Y.; Hu, X.; Solan, G.A.; Sun, W.-H. Thermally stable and highly active cobalt precatalysts for vinyl-polyethylenes with narrow polydispersities: Integrating fused-ring and imino-carbon protection into ligand design. New J. Chem. 2016, 40, 8012–8023. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, Y.; Suo, H.; Huang, F.; Liang, T.; Hu, X.; Sun, W.-H. Thermo-stable 2-(arylimino)benzylidene-9-arylimino-5,6,7,8-tetrahydro cyclohepta[b]pyridyliron(II) precatalysts toward ethylene polymerization and highly linear polyethylenes. J. Polym. Sci. Part A Polym. Chem. 2017, 55, 830–842. [Google Scholar] [CrossRef]
  16. Appukuttan, V.K.; Liu, Y.; Son, B.C.; Ha, C.-S.; Suh, H.; Kim, I. Iron and Cobalt Complexes of 2,3,7,8-Tetrahydroacridine-4,5(1H,6H)-diimine Sterically Modulated by Substituted Aryl Rings for the Selective Oligomerization to Polymerization of Ethylene. Organometallics 2011, 30, 2285–2294. [Google Scholar] [CrossRef]
  17. Du, S.; Wang, X.; Zhang, W.; Flisak, Z.; Sun, Y.; Sun, W.-H. A practical ethylene polymerization for vinyl-polyethylenes: Synthesis, characterization and catalytic behavior of α,α′-bisimino-2,3:5,6-bis(pentamethylene)pyridyliron chlorides. Polym. Chem. 2016, 7, 4188–4197. [Google Scholar] [CrossRef]
  18. Zada, M.; Guo, L.; Ma, Y.; Zhang, W.; Flisak, Z.; Sun, Y.; Sun, W.-H. Activity and Thermal Stability of Cobalt(II)-Based Olefin Polymerization Catalysts Adorned with Sterically Hindered Dibenzocycloheptyl Groups. Molecules 2019, 24, 2007. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Z.; Ma, Y.; Guo, J.; Liu, Q.; Solan, G.A.; Liang, T.; Sun, W.-H. Bis(imino)pyridines fused with 6- and 7-membered carbocylic rings as N,N,N-scaffolds for cobalt ethylene polymerization catalysts. Dalton Trans. 2019, 48, 2582–2591. [Google Scholar] [CrossRef] [Green Version]
  20. Guo, L.; Zada, M.; Zhang, W.; Vignesh, A.; Zhu, D.; Ma, Y.; Liang, T.; Sun, W.-H. Highly linear polyethylenes tailored with 2,6-bis[1-(p-dibenzo-cycloheptylarylimino)ethyl]pyridylcobalt dichlorides. Dalton Trans. 2019, 48, 5604–5613. [Google Scholar] [CrossRef]
  21. Chen, Q.; Zhang, W.; Solan, G.A.; Zhang, R.; Guo, L.; Hao, X.; Sun, W.-H. CH(phenol)-Bridged Bis(imino)pyridines as Compartmental Supports for Diiron Precatalysts for Ethylene Polymerization: Exploring Cooperative Effects on Performance. Organometallics 2018, 37, 4002–4014. [Google Scholar] [CrossRef] [Green Version]
  22. Chen, Q.; Zhang, W.; Solan, G.A.; Liang, T.; Sun, W.-H. Methylene-bridged bimetallic bis(imino)pyridine-cobaltous chlorides as precatalysts for vinyl-terminated polyethylene waxes. Dalton Trans. 2018, 47, 6124–6133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Khoshsefat, M.; Dechal, A.; Ahmadjo, S.; Mortazavi, S.M.M.; Zohuri, G.; Soares, J.B.P. Amorphous to high crystalline PE made by mono and dinuclear Fe-based catalysts. Eur. Polym. J. 2019, 119, 229–238. [Google Scholar] [CrossRef]
  24. Huang, F.; Zhang, W.; Yue, E.; Liang, T.; Hu, X.; Sun, W.-H. Controlling the molecular weights of polyethylene waxes using the highly active precatalysts of 2-(1-aryliminoethyl)-9-arylimino-5,6,7,8-tetrahydrocycloheptapyridylcobalt chlorides: Synthesis, characterization, and catalytic behavior. Dalton Trans. 2016, 45, 657–666. [Google Scholar] [CrossRef] [PubMed]
  25. Du, S.; Zhang, W.; Yue, E.; Huang, F.; Liang, T.; Sun, W.-H. α,α′-Bis(arylimino)-2,3:5,6-bis(pentamethylene)pyridylcobalt Chlorides: Synthesis, Characterization, and Ethylene Polymerization Behavior. Eur. J. Inorg. Chem. 2016, 2016, 1748–1755. [Google Scholar] [CrossRef]
  26. Britovsek, G.J.P.; Clentsmith, G.K.B.; Gibson, V.C.; Goodgame, D.M.L.; McTavish, S.J.; Pankhurst, Q.A. The nature of the active site in bis(imino)pyridine iron ethylene polymerisation catalysts. Catal. Commun. 2002, 3, 207–211. [Google Scholar] [CrossRef]
  27. Humphries, M.J.; Tellmann, K.P.; Gibson, V.C.; White, A.J.P.; Williams, D.J. Investigations into the mechanism of activation and initiation of ethylene polymerization by bis(imino)pyridine cobalt catalysts: Synthesis, structures, and deuterium labeling studies. Organometallics 2005, 24, 2039–2050. [Google Scholar] [CrossRef]
  28. Gibson, V.C.; Humphries, M.J.; Tellmann, K.P.; Wass, D.F.; White, A.J.; Williams, D.J. The nature of the active species in bis(imino)pyridyl cobalt ethylene polymerisation catalysts. Chem. Commun. 2001, 2252–2253. [Google Scholar] [CrossRef]
  29. Roemelt, C.; Weyhermueller, T.; Wieghardt, K. Structural characteristics of redox-active pyridine-1,6-diimine complexes: Electronic structures and ligand oxidation levels. Coord. Chem. Rev. 2019, 380, 287–317. [Google Scholar] [CrossRef]
  30. Bryliakov, K.P.; Semikolenova, N.V.; Zakharov, V.A.; Talsi, E.P. Active intermediates of ethylene polymerization over 2,6-bis(imino)pyridyl iron complex activated with aluminum trialkyls and methylaluminoxane. Organometallics 2004, 23, 5375–5378. [Google Scholar] [CrossRef]
  31. Bryliakov, K.P.; Talsi, E.P.; Semikolenova, N.V.; Zakharov, V.A. Formation and Nature of the Active Sites in Bis(imino)pyridine Iron-Based Polymerization Catalysts. Organometallics 2009, 28, 3225–3232. [Google Scholar] [CrossRef]
  32. Bryliakov, K.P.; Talsi, E.P. Frontiers of mechanistic studies of coordination polymerization and oligomerization of α-olefins. Coord. Chem. Rev. 2012, 256, 2994–3007. [Google Scholar] [CrossRef]
  33. Kleigrewe, N.; Steffen, W.; Blomker, T.; Kehr, G.; Frohlich, R.; Wibbeling, B.; Erker, G.; Wasilke, J.C.; Wu, G.; Bazan, G.C. Chelate bis(imino)pyridine cobalt complexes: Synthesis, reduction, and evidence for the generation of ethene polymerization catalysts by Li+ cation activation. J. Am. Chem. Soc. 2005, 127, 13955–13968. [Google Scholar] [CrossRef] [PubMed]
  34. Soshnikov, I.E.; Semikolenova, N.V.; Bushmelev, A.N.; Bryliakov, K.P.; Lyakin, O.V.; Redshaw, C.; Zakharov, V.A.; Talsi, E.P. Investigating the Nature of the Active Species in Bis(imino)pyridine Cobalt Ethylene Polymerization Catalysts. Organometallics 2009, 28, 6003–6013. [Google Scholar] [CrossRef]
  35. Deng, L.Q.; Margl, P.; Ziegler, T. Mechanistic aspects of ethylene polymerization by iron(II)-bisimine pyridine catalysts: A combined density functional theory and molecular mechanics study. J. Am. Chem. Soc. 1999, 121, 6479–6487. [Google Scholar] [CrossRef]
  36. Kwon, D.H.; Small, B.L.; Sydora, O.L.; Bischof, S.M.; Ess, D.H. Challenge of Using Practical DFT to Model Fe Pendant Donor Diimine Catalyzed Ethylene Oligomerization. J. Phys. Chem. C 2019, 123, 3727–3739. [Google Scholar] [CrossRef]
  37. Cruz, V.L.; Ramos, J.; Martinez-Salazar, J.; Gutierrez-Oliva, S.; Toro-Labbe, A. Theoretical Study on a Multicenter Model Based on Different Metal Oxidation States for the Bis(imino)pyridine Iron Catalysts in Ethylene Polymerization. Organometallics 2009, 28, 5889–5895. [Google Scholar] [CrossRef]
  38. Griffiths, E.A.H.; Britovsek, G.J.P.; Gibson, V.C.; Gould, I.R. Highly active ethylene polymerisation catalysts based on iron: An ab initio study. Chem. Commun. 1999, 1333–1334. [Google Scholar] [CrossRef]
  39. Tondreau, A.M.; Milsmann, C.; Patrick, A.D.; Hoyt, H.M.; Lobkovsky, E.; Wieghardt, K.; Chirik, P.J. Synthesis and Electronic Structure of Cationic, Neutral, and Anionic Bis(imino)pyridine Iron Alkyl Complexes: Evaluation of Redox Activity in Single-Component Ethylene Polymerization Catalysts. J. Am. Chem. Soc. 2010, 132, 15046–15059. [Google Scholar] [CrossRef]
  40. Khoroshun, D.V.; Musaev, D.G.; Vreven, T.; Morokuma, K. Theoretical study on bis(imino)pyridyl-Fe(II) olefin poly- and oligomerization catalysts. Dominance of different spin states in propagation and beta-hydride transfer pathways. Organometallics 2001, 20, 2007–2026. [Google Scholar] [CrossRef]
  41. Boudene, Z.; Boudier, A.; Breuil, P.-A.R.; Olivier-Bourbigou, H.; Raybaud, P.; Toulhoat, H.; de Bruin, T. Understanding the role of aluminum-based activators in single site iron catalysts for ethylene oligomerization. J. Catal. 2014, 317, 153–157. [Google Scholar] [CrossRef] [Green Version]
  42. Bouwkamp, M.W.; Bart, S.C.; Hawrelak, E.J.; Trovitch, R.J.; Lobkovsky, E.; Chirik, P.J. Square planar bis(imino)pyridine iron halide and alkyl complexes. Chem. Commun. 2005, 3406–3408. [Google Scholar] [CrossRef] [PubMed]
  43. Bellows, S.M.; Cundari, T.R.; Holland, P.L. Spin Crossover during β-Hydride Elimination in High-Spin Iron(II)– and Cobalt(II)–Alkyl Complexes. Organometallics 2013, 32, 4741–4751. [Google Scholar] [CrossRef]
  44. Huang, C.; Huang, Y.; Ma, Y.; Solan, G.A.; Sun, Y.; Hu, X.; Sun, W.-H. Cycloheptyl-fused N,N,N′-chromium catalysts with selectivity for vinyl-terminated polyethylene waxes: Thermal optimization and polymer functionalization. Dalton Trans. 2018, 47, 13487–13497. [Google Scholar] [CrossRef] [Green Version]
  45. Bouwkamp, M.W.; Lobkovsky, E.; Chirik, P.J. Bis(imino)pyridine iron(II) alkyl cations for olefin polymerization. J. Am. Chem. Soc. 2005, 127, 9660–9661. [Google Scholar] [CrossRef] [PubMed]
  46. Lu, T.; Chen, F.W. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  47. Raucoules, R.; de Bruin, T.; Raybaud, P.; Adamo, C. Theoretical Unraveling of Selective 1-Butene Oligomerization Catalyzed by Iron−Bis(arylimino)pyridine. Organometallics 2009, 28, 5358–5367. [Google Scholar] [CrossRef]
  48. Reed, A.E.; Weinstock, R.B.; Weinhold, F. Natural population analysis. J. Chem. Phys. 1985, 83, 735–746. [Google Scholar] [CrossRef]
  49. Ba, J.; Du, S.; Yue, E.; Hu, X.; Flisak, Z.; Sun, W.-H. Constrained formation of 2-(1-(arylimino)ethyl)-7-arylimino-6,6-dimethylcyclopentapyridines and their cobalt(ii) chloride complexes: Synthesis, characterization and ethylene polymerization. RSC Adv. 2015, 5, 32720–32729. [Google Scholar] [CrossRef]
  50. Pasha, F.A.; Basset, J.-M.; Toulhoat, H.; de Bruin, T. DFT Study on the Impact of the Methylaluminoxane Cocatalyst in Ethylene Oligomerization Using a Titanium-Based Catalyst. Organometallics 2015, 34, 426–431. [Google Scholar] [CrossRef]
  51. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision E.01; Gaussian, Inc.: Wallingford, CT, USA, 2013. [Google Scholar]
  52. Minaev, B.; Baryshnikova, A.; Sun, W.-H. Spin-dependent effects in ethylene polymerization with bis(imino)pyridine iron(II) complexes. J. Organomet. Chem. 2016, 811, 48–65. [Google Scholar] [CrossRef]
  53. Ehm, C.; Budzelaar, P.H.M.; Busico, V. Calculating accurate barriers for olefin insertion and related reactions. J. Organomet. Chem. 2015, 775, 39–49. [Google Scholar] [CrossRef]
  54. Dolg, M.; Wedig, U.; Stoll, H.; Preuss, H. Energy-adjusted abinitio pseudopotentials for the 1st-row transition-elements. J. Chem. Phys. 1987, 86, 866–872. [Google Scholar] [CrossRef]
  55. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Gonzalez, C.; Schlegel, H.B. Reaction-path following in mass-weighted internal coordinates. J. Phys. Chem. 1990, 94, 5523–5527. [Google Scholar] [CrossRef]
  57. Reed, A.E.; Weinhold, F. Natural localized molecular-orbitals. J. Chem. Phys. 1985, 83, 1736–1740. [Google Scholar] [CrossRef]
  58. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  59. Reed, A.E.; Schleyer, P.V. chemical bonding in hypervalent molecules—The dominance of ionic bonding and negative hyperconjugation over d-orbital participation. J. Am. Chem. Soc. 1990, 112, 1434–1445. [Google Scholar] [CrossRef]
Catalysts 10 01396 sch001 550
Scheme 1. The reaction pathway of ethylene polymerization catalyzed by LMCH3 (M = Fe and Co) species.
Scheme 1. The reaction pathway of ethylene polymerization catalyzed by LMCH3 (M = Fe and Co) species.
Catalysts 10 01396 sch001
Catalysts 10 01396 g001 550
Figure 1. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the neutral cobalt complex in singlet (dash line) and triplet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LCoCH3 (1) (triplet state) and C2H4. The number n corresponds to Co(0)-n, the TS m/n refers to the transition state between the two intermediates Co(0)-m and Co(0)-n. The n’ was used to distinguish structures between the two spin states, as is the same in the following Figures.
Figure 1. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the neutral cobalt complex in singlet (dash line) and triplet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LCoCH3 (1) (triplet state) and C2H4. The number n corresponds to Co(0)-n, the TS m/n refers to the transition state between the two intermediates Co(0)-m and Co(0)-n. The n’ was used to distinguish structures between the two spin states, as is the same in the following Figures.
Catalysts 10 01396 g001
Catalysts 10 01396 g002 550
Figure 2. Optimized geometries of the ethylene coordinated cobalt alkyl complexes 1Co(0)-2, 1Co(0)-4 and 1Co(0)-6. (Distances are given in angstrom (Å), bonding energy is given in the bracket in kcal/mol).
Figure 2. Optimized geometries of the ethylene coordinated cobalt alkyl complexes 1Co(0)-2, 1Co(0)-4 and 1Co(0)-6. (Distances are given in angstrom (Å), bonding energy is given in the bracket in kcal/mol).
Catalysts 10 01396 g002
Catalysts 10 01396 sch002 550
Scheme 2. Atom labels around the active centers and fragmental units for the ethylene coordinated sM(q)-2 π-complexes. The hydrogen atoms were omitted for charity.
Scheme 2. Atom labels around the active centers and fragmental units for the ethylene coordinated sM(q)-2 π-complexes. The hydrogen atoms were omitted for charity.
Catalysts 10 01396 sch002
Catalysts 10 01396 g003 550
Figure 3. HUMO orbitals of 1Co(0)-2 and 0Co(0)-2 (green; positive orbital phase and blue; negative orbital phase, isovalue of orbital is 0.05).
Figure 3. HUMO orbitals of 1Co(0)-2 and 0Co(0)-2 (green; positive orbital phase and blue; negative orbital phase, isovalue of orbital is 0.05).
Catalysts 10 01396 g003
Catalysts 10 01396 g004 550
Figure 4. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the cationic cobalt complex in quartet (dash line) and doublet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LCo+CH3 (1) (doublet state) and C2H4. The number n corresponds to Co(+)-n, the TS m/n refers to the transition state between the two intermediates Co(+)-m and Co(+)-n.
Figure 4. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the cationic cobalt complex in quartet (dash line) and doublet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LCo+CH3 (1) (doublet state) and C2H4. The number n corresponds to Co(+)-n, the TS m/n refers to the transition state between the two intermediates Co(+)-m and Co(+)-n.
Catalysts 10 01396 g004
Catalysts 10 01396 g005 550
Figure 5. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the cationic iron complex in triplet (dash line) and quintet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LFeCH3 (+) (1)(quintet state) and C2H4. The number n corresponds to Fe(+)-n, the TS m/n refers to the transition state between the two intermediates Fe(+)-m and Fe(+)-n.
Figure 5. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the cationic iron complex in triplet (dash line) and quintet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LFeCH3 (+) (1)(quintet state) and C2H4. The number n corresponds to Fe(+)-n, the TS m/n refers to the transition state between the two intermediates Fe(+)-m and Fe(+)-n.
Catalysts 10 01396 g005
Catalysts 10 01396 g006 550
Figure 6. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the neutral iron complex in doublet (dash line) and quartet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LFeCH3 (1, quartet state) and C2H4. The number n corresponds to Fe(0)-n, the TS m/n refers to the transition state between the two intermediates Fe(0)-m and Fe(0)-n.
Figure 6. Energy profiles (kcal/mol) for the first and second ethylene coordination and insertion (black), β-H transfer (BHT, in blue) and β-H elimination (BHE, in red) catalyzed by the neutral iron complex in doublet (dash line) and quartet (solid line) spin states. The energies are relative zero-point energies with respect to the reactants: LFeCH3 (1, quartet state) and C2H4. The number n corresponds to Fe(0)-n, the TS m/n refers to the transition state between the two intermediates Fe(0)-m and Fe(0)-n.
Catalysts 10 01396 g006
Catalysts 10 01396 sch003 550
Scheme 3. Definition of the factors to calculating the structure parameter Δ in the bis(imino)pyridine (PDI) skeleton structure.
Scheme 3. Definition of the factors to calculating the structure parameter Δ in the bis(imino)pyridine (PDI) skeleton structure.
Catalysts 10 01396 sch003
Catalysts 10 01396 g007 550
Figure 7. The structure parameter Δ calculated for both iron and cobalt species.
Figure 7. The structure parameter Δ calculated for both iron and cobalt species.
Catalysts 10 01396 g007
Catalysts 10 01396 g008 550
Figure 8. (a) The natural bond orbital (NBO) charge distribution for sM(q)-1 and (b) the charge changes aroused by ethylene coordination in sM(q)-2, where M is iron or cobalt, q is zero or positive, with high and low spin states.
Figure 8. (a) The natural bond orbital (NBO) charge distribution for sM(q)-1 and (b) the charge changes aroused by ethylene coordination in sM(q)-2, where M is iron or cobalt, q is zero or positive, with high and low spin states.
Catalysts 10 01396 g008
Catalysts 10 01396 sch004 550
Scheme 4. The oxidation states of the metal center and the PDI ligand, the numbers of the valence electrons around the metals.
Scheme 4. The oxidation states of the metal center and the PDI ligand, the numbers of the valence electrons around the metals.
Catalysts 10 01396 sch004
Catalysts 10 01396 sch005 550
Scheme 5. The reaction pathway of ethylene polymerization catalyzed by LsCo(+)B on addition of [CH3]− or [AlMe4]− or C2H4.
Scheme 5. The reaction pathway of ethylene polymerization catalyzed by LsCo(+)B on addition of [CH3]− or [AlMe4]− or C2H4.
Catalysts 10 01396 sch005
Catalysts 10 01396 g009 550
Figure 9. Comparison of the energy profiles for LsCo(+)B (L = ArPDI, s means spin states, here is singlet and triplet, the B means the cobalt metal having no other coordination besides the L) on addition of [CH3]− or [AlMe4]− or C2H4.
Figure 9. Comparison of the energy profiles for LsCo(+)B (L = ArPDI, s means spin states, here is singlet and triplet, the B means the cobalt metal having no other coordination besides the L) on addition of [CH3]− or [AlMe4]− or C2H4.
Catalysts 10 01396 g009
Table
Table 1. Bond order of the bonds in the active center for 1Co(0)-2 and 0Co(0)-2 by Mayer bond order analysis (“---” represents coordination interaction and “–” represents a chemical bond).
Table 1. Bond order of the bonds in the active center for 1Co(0)-2 and 0Co(0)-2 by Mayer bond order analysis (“---” represents coordination interaction and “–” represents a chemical bond).
BondsBond OrderBondsBond Order
1Co(0)-20Co(0)-2 1Co(0)-20Co(0)-2
C(1)-–C(2)1.361.42C(10)–N(13)1.441.47
C(2)–C(3)1.491.41C(11)–N(12)1.821.48
C(3)–C(4)1.331.41N(12)---Co(66)0.110.47
C(4)–C(5)1.471.42N(13)---Co(66)0.390.48
C(1)–N(6)1.131.15C(56)–Co(66)0.730.79
C(5)–N(6)1.251.15C(60)–C(63)1.631.29
C(1)–C(10)1.231.24C(60)---Co(66)0.240.54
C(5)–C(11)1.081.24C(63)---Co(66)0.180.49
N(6)---Co(66)0.420.61
Table
Table 2. The structure parameter Δ calculated for the species in Figure 9.
Table 2. The structure parameter Δ calculated for the species in Figure 9.
0Co(+)B1A1B1C2A2BTS2A/3ATS2B/3B3A/3B
Singlet0.150.10.1580.1140.0850.0840.0610.0770.095
Triplet0.1810.1230.180.1190.1230.1680.1120.1160.114
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Ma, Y.; Sun, W.-H. Comparison of the Reactivity and Structures for the Neutral and Cationic Bis(imino)pyridyl Iron and Cobalt Species by DFT Calculations. Catalysts 2020, 10, 1396. https://0-doi-org.brum.beds.ac.uk/10.3390/catal10121396

AMA Style

Li Z, Ma Y, Sun W-H. Comparison of the Reactivity and Structures for the Neutral and Cationic Bis(imino)pyridyl Iron and Cobalt Species by DFT Calculations. Catalysts. 2020; 10(12):1396. https://0-doi-org.brum.beds.ac.uk/10.3390/catal10121396

Chicago/Turabian Style

Li, Zilong, Yanping Ma, and Wen-Hua Sun. 2020. "Comparison of the Reactivity and Structures for the Neutral and Cationic Bis(imino)pyridyl Iron and Cobalt Species by DFT Calculations" Catalysts 10, no. 12: 1396. https://0-doi-org.brum.beds.ac.uk/10.3390/catal10121396

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop