Next Article in Journal
Application of Lecitase® Ultra-Catalyzed Hydrolysis to the Kinetic Resolution of (E)-4-phenylbut-3-en-2-yl Esters
Next Article in Special Issue
Theoretical Study on the Quantum Capacitance Origin of Graphene Cathodes in Lithium Ion Capacitors
Previous Article in Journal
Theoretical Studies on the Direct Propylene Epoxidation Using Gold-Based Catalysts: A Mini-Review
Previous Article in Special Issue
Theoretical Study on the Hydrogenation Mechanisms of Model Compounds of Heavy Oil in a Plasma-Driven Catalytic System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Straightforward Design for Phenoxy-Imine Catalytic Activity in Ethylene Polymerization: Theoretical Prediction

1
School of Molecular Science and Engineering, Vidyasirimedhi Institute of Science and Technology, Rayong 21210, Thailand
2
Research Center of Nanoscience and Technology, Shanghai University, Shanghai 200444, China
*
Authors to whom correspondence should be addressed.
Catalysts 2018, 8(10), 422; https://doi.org/10.3390/catal8100422
Submission received: 4 September 2018 / Revised: 25 September 2018 / Accepted: 25 September 2018 / Published: 28 September 2018
(This article belongs to the Special Issue Heterogeneous Catalysis for Energy Conversion)

Abstract

:
The quantitative structure-activity relationship (QSAR) of 18 Ti-phenoxy-imine (FI-Ti)-based catalysts was investigated to clarify the role of the structural properties of the catalysts in polyethylene polymerization activity. The electronic properties of the FI-Ti catalysts were analyzed based on density functional theory with the M06L/6-31G** and LANL2DZ basis functions. The analysis results of the QSAR equation with a genetic algorithm showed that the polyethylene catalytic activity mainly depended on the highest occupied molecular orbital energy level and the total charge of the substituent group on phenylimine ring. The QSAR models showed good predictive ability (R2) and R2 cross validation (R2cv) values of greater than 0.927. The design concept is “head-hat”, where the hats are the phenoxy-imine substituents, and the heads are the transition metals. Thus, for the newly designed series, the phenoxy-imine substituents still remained, while the Ti metal was replaced by Zr or Ni transition metals, entitled FI-Zr and FI-Ni, respectively. Consequently, their polyethylene polymerization activities were predicted based on the obtained QSAR of the FI-Ti models, and it is noteworthy that the FI-Ni metallocene catalysts tend to increase the polyethylene catalytic activity more than that of FI-Zr complexes. Therefore, the new designs of the FI-Ni series are proposed as candidate catalysts for polyethylene polymerization, with their predicted activities in the range of 35,000–48,000 kg(PE)/mol(Cat.)·MPa·h. This combined density functional theory and QSAR analysis is useful and straightforward for molecular design or catalyst screening, especially in industrial research.

Graphical Abstract

1. Introduction

As an inexpensive material, polyethylene (PE) has drawn significant attention due to its excellent physical properties, such as high mechanical strength, high elasticity and high resistivity towards corrosion, light weightness, and reusability [1,2,3,4]. Normally, the polyethylene materials are prepared via the coordination polymerization using different categories of catalysts, among which metallocenes are very promising candidates that provide high activity and selectivity [5,6,7,8,9,10,11]. Many experimental and theoretical methodologies show that the structures of the metallocenes determine their macroscopic properties and catalytic behaviors. Understanding these structures can thus lead to the development of new catalysts with improved properties.
Recently, the development of single site olefin polymerization is one of the hot topics in industrial polyolefins. Generally, to improve the polymerization catalytic performance, there are two desirable methods: modification of ligands or transition metal substitutions [11,12,13,14,15,16,17]. Each fine-tuned ligand unit can offer various structural and electronic properties that influence the catalytic activity [13,18]. In particular, phenoxy-imine can efficiently provide a suitable bulky structure for the ethylene polymerization process. The bidentate ligand was successfully studied by Fujita and co-workers through the generation of a bis(phenoxy-imine) transition metal complex [19,20,21,22,23,24]. The formation of this organic bulky substituent enhances the catalytic activity, thus initiating the polymerization. Therefore, nowadays, many researchers aim to develop efficient catalysts through ligand design and functionality modification [25,26,27,28,29]. In addition, the metal center was also investigated by substitution between early transition metals and late transition metals, such as Ti and Ni, respectively [2,5,30,31].
Since current technologies offer high efficiency and performance in simulations and calculations, computational studies of catalyst chemistry are becoming more and more applied to clarify phenomena in molecules. Significant polyolefin catalyst research has been adopted in order to understand the molecular details. As with polyolefin Ziegler–Natta catalysts, Ratanasak and co-worker declared that the calculation results could reasonably explain the role of both electric and steric parameters on catalytic activity and stereoselectivity via density functional theory (DFT) calculations [32,33]. Moreover, the group of Ojwach reported the DFT study of transition metal complex catalysts, which justified that the molecular structure directly affected the catalytic activity [2,34]. The theoretical study of polyethylene polymerization via a titanium phenoxy-imine catalyst was further carried out by Nikitin and co-workers [35]. Their DFT calculations showed that both the structural and electronic properties were related to the experimental catalytic activities. The key parameter was the energy band gap of the frontier molecular orbitals, which depended on the active site of the catalyst. By tuning this parameter, it was showed that the structures and type of substituent groups on the ligands had a significant impact on the energy band gap and catalytic performance. Nevertheless, there was still a lack of information regarding the metal center effect, although this was directly positioned at the active site.
The quantitative structure–activity relationship (QSAR) method, an approach for examining the activities of these catalysts in relation to their structures, is widely applied in catalyst research [28,32,36,37,38,39]. In particular, this technique can also provide a better explanation of the different types of electronic and steric factors, which not only affect the bond lengths and angles, but also the highest occupied molecular orbital (HOMO) energy, the lowest unoccupied molecular orbital (LUMO) energy, and the atomic charge. As a result, this quantitative analysis can give an effective insight into the efficiency of different catalysts based on their structure and electronic properties. Thus, in this work, the DFT optimization method and QSAR analysis were applied to investigate the electronic descriptors of Ti-phenoxyimine (FI-Ti) metallocene catalysts related to their polyethylene polymerization activity. A genetic algorithm with multiple linear regression (GA-MLR) was selected to derive the QSAR equations [40,41,42,43]. Based on the best obtained QSAR equation, it was then used to predict the polyethylene polymerization activity of the newly designed Zr- and Ni-phenoxyimine catalysts.

2. Results and Discussion

2.1. Electronic Properties

In this work, we studied the relationship between the PE activity and the basic electronic properties of the FI-Ti catalysts. Eighteen FI-Ti catalysts were used as a training set for manipulating the QSAR equation. The structure details of each candidate and their experimental catalyst activity that were collected from Nikitin’s work [35] are shown in Table 1 (for more details regarding experimental data, see Materials and Methods). The range of PE catalytic activities (A) was found to be 4100–54,200 kg(PE)/mol(Cat.)·MPa·h. Therefore, the activities were converted into ln(A) values (listed in Table 1), which is the same unit as the electronic properties. Based on the QSAR equation, the experimental activity would serve as the dependent variable, while the electronic properties would be independent ones. Thus, all catalysts were built and then optimized using DFT with the M06L exchange-correlation functional [44] via the 6-31G** and LANL2DZ basis functions. The M06L, an empirical functional, is accurate for noncovalent interactions, main group thermochemistry, and transition metals [45,46]. In this work, we applied the M06L functional with 6-31G** for C, H, O, N, and Br, and LANL2DZ for Ti, Ni, and Zr. The selection of the M06L functional is based on previous benchmark studies [46,47,48]. Zhao and Truhlar [46] reported that M06L has reasonably good overall performance, and it is the best functional for transition metal energetics, and for the study of organometallic and inorganometallic thermochemistry. Furthermore, M06L was reported as a good choice for calculations on some small non-covalent dimers in both geometric and electronic studies [47]. In addition, the M06L functional was successfully carried out on dihydropyrazine series annulated linear polyacene systems to study the aromaticity and HOMO–LUMO energy [48].
Based on the optimized FI-Ti structures, with consideration of a symmetric skeleton structure of metallocene catalysts (see inserted figure of Table 2), there is only one Ti atom, thus, the partial charge of the Ti atom was collected. Alternatively, a pair of O, N, and Cl atoms coordinated to Ti metal; therefore, the charge of O, N, and Cl was calculated from the average charge of (eA + eB)/2, where the eA and eB are the partial charge of OA, NA, or ClA, and OB, NB, or ClB, respectively. For the phenoxy-imine ligands, there are three substituents attached on phenyl rings; R1 and R2 were o- and p-substitutions on phenoxy ring I, respectively, whereas R3 was a m-, or p-substituent on phenylimine ring II. Therefore, to elucidate the electronic effects of ligand substitution, the partial charge of all atoms in each substituent group was combined, which recalled the total charge of the R1, R2 and R3 groups. However, based on the bis-phenoxyimine structure, there are pairs of R1, R2, and R3; therefore, the electronic charge of ligand substitutions was calculated from (tA + tB)/2, where tA is the total charge of the R1A, R2A, or R3A group, and tB is the total charge of the R1B, R2B, or R3B group. In addition, the HOMO energy, LUMO energy, and the HOMO–LUMO energy gap (E-gap) were also included, as listed in Table 2.

2.2. QSAR Model Relationship

The main steps in the QSAR study were the calculations and selections of structural descriptors as numerical variables representing each chemical structure. Firstly, the ln(A) activity and all electronic descriptors were calculated, and the pairwise correlation-located correlation coefficient was a number that was a measure of the strength and the direction of the correlation between two variables. It was be expressed using the variable r, where r was between 1 and −1. Thus, the data points with a perfect strength line are r = −1 or r = 1. When r < 0 or r > 0, the data would have negative or positive associations, respectively. The obtained correlation coefficients of greater than ±0.8 were classified as being intercorrelated, and only one of them with a high correlation with the activity data was considered in developing the model. Based on our correlation matrix in Table 3, the high correlation coefficient was only found on HOMO and LUMO pairwise, with a 0.890 correlation coefficient value. However, both of them resulted in similar correlation coefficients to ln(A); therefore, these two descriptors would be further selected by the occurrences of the population in the genetic algorithm (GA) process. Then, all 10 electronic descriptors were used as input descriptors for manipulating the QSAR equation of the FI-Ti catalysts.
GAs are a well-known and widely employed variable selection method that is used to solve the optimization problems defined by fitness criteria, applying the evolution hypothesis of Darwin and various genetic functions, i.e., cross-over and mutation [49,50,51,52]. GAs have been applied in various fields for many applications involving regression problems, and they offer a new approach to the problem of QSAR models. By using the major characteristics of Darwinian evolution, GAs can construct complicated models that not only use linear polynomials, but also higher-order polynomials, splines and Gaussians, allowing for the production of a more widespread model [43,53,54,55,56,57]. Generally, GAs work with a set of strings that are called a population. This population is evolved in such a manner that leads it toward the objective of the search. There are three operations of selection, crossover and mutation, which are iteratively performed, and then newly added members are scored according to a fitness criterion that is related to the quality of the regression fit to the data. Thus, the selection probabilities must be re-evaluated each time a new member is added to the population, and then the convergence is triggered by a lack of progress in the highest and average scores of the population [57]. Therefore, in this work, to select the most relevant descriptors, the evolution of the population was simulated and the occurrences of populations of 10 descriptors are depicted in Figure 1. As the GA population of each parameter was considered to be the basic value of fitness for the QSAR model, it was found that the HOMO energy level exhibits the highest value, followed by the R3 total charge, and then the Ti charge with occurrences in population of 2097, 1067, and 351, respectively, whereas the rest of the descriptors only resulted in a range of 10–67 populations, as shown in Figure 1. This indicates that the FI-Ti catalytic activity is closely related to the HOMO energy level, the R3 total charge, and the Ti charge. Looking into the details of both correlation coefficient and occurrences of populations, we observed that the complex with the lower HOMO energy level combined with the higher negative charge of the R3 substituent group showed the better catalytic activity (see in Table 2). For example, the complex 3Ti with a low HOMO energy of −0.194 eV and high R3 negative charge of −1.230 achieved the highest PE activity value of 54,200 kg(PE)/mol(Cat.)·MPa·h, while complex 1Ti with R3 positive charge gave the lowest catalytic activity of 4100 kg(PE)/mol(Cat.)·MPa·h. Consequently, when the substituents differ, the Ti center atom charge was also affected, and complexes with a Ti charge of ~1.24–1.25 showed higher catalytic activity than the others. Thus, these three descriptors would further be set as independent variables in the QSAR equation investigation.
From the GA population results, we found the promising parameter candidates for QSAR model construction, then, those three electronic parameters of the HOMO energy level, R3 total charge and Ti atomic charge were used to generate Equation (1) see Table 4. With consideration of QSAR Equation (1), the linear function was found to describe the relationship between the PE catalytic activity in terms of ln(A) and three electronic properties, the R3 charge, HOMO energy, and Ti atomic charge. It consisted of ramp functions that are continuous functions, segmented in three parts: xfit(t) = x1 for tt1, x2 for tt2, and linearly connected between t1 and t2. This type of function is very useful and has many applications in engineering and mathematics [58,59]. As shown in Table 4, Equation (1) exhibited an excellent correlation with an R-squared (R2) value of 0.992. The obtained model also showed a high prediction ability with a cross-validated technique R-squared (R2cv) value of 0.936. In addition, the data distribution was evaluated by the F value, which resulted in 168.295. Although it seems that this model is a good model, the population data of the GA imposed a question on us that Ti charge population was actually very small compared with the other two descriptors. Furthermore, based on the “head-hat” design concept, the Ti center would be replaced with other metals; thus, to avoid the direct effect of the metal’s partial charge, the elimination of the Ti charge parameter was suitable to accomplish a better QSAR model.

2.3. QSAR Model Improvement via the Elimination of Insignificant Descriptors and Outliers

As discussed above, the Ti charge exhibited a significantly lower occurrence population when compared with the two properties in Equation (1). Thus, an elimination of the weak parameter strategy was applied. The two-descriptor models of the HOMO energy levels, and the R3 total charge was generated. As shown in Table 4, the QSAR Equation (2) still obtained a good predictive ability with an R2 value of 0.992 and rising on R2cv of 0.984. Furthermore, the data distribution F value of Equation (2) was also increased to 230.917, which was higher than that in Equation (1). It is noteworthy that with only two electronic property variables, the obtained QSAR model resulted in higher reliability and better predictive ability. The analysis results confirmed that the partial Ti atom charge plays a less important role in PE catalytic activity prediction.
Generally, with the exception of the predictive ability, the model with the higher F value was considered to be the better model. Although Equation (2) showed very good correlation and predictable ability, if there were some choices for more accurate predicted results, it would be completed to improve the model. Thus, based on Equation (2), the accuracy was considered via the residual value that was calculated from the experimental and predicted ln(A) of PE polymerization, as illustrated in Figure 2. The graph shows that there were two outliers exhibiting more than a |2| scaled residual value, which are the complexes 6Ti (2.128) and 12Ti (−2.168). These outliers show high residual results for the predicted activity compared to the others, and this reason led to the error of the whole model. Thus, to further improve the QSAR model, the data of 6Ti and 12Ti complexes were eliminated, and then QSAR Equation (3) was generated. The results indicated that the new model of Equation (3) achieved a superior higher F value of 564.305, which was twice that of the previous one, while still exhibiting high prediction ability with an R2 of 0.997 and an R2cv of 0.927. Therefore, both Equations (2) and (3) were selected to predict the PE activity of FI-Ti catalysts using the GA-MLR of the ramp function, and their graphical plots between the experimental and predicted PE activities are depicted in Figure 3a,b, respectively.

2.4. Prediction of PE Polymerization Activity for Newly Designed Zr and Ni Phenoxyimine Catalysts

We propose the design concept of “head-hat”, where the hats are phenoxy-imine substituents and the heads are transition metals. In this work, based on the training set of the FI-Ti catalysts, the Ti metal center was then replaced by Zr or Ni transition metals, while keeping the same phenoxy-imine substituents, namely, the FI-Zr and FI-Ni series, respectively, as shown in Table 5. The new FI-Zr and FI-Ni complex structures were systemically built and optimized using the same criteria as the FI-Ti catalysts. Their HOMO energy levels and R3 total charges were then extracted. The prediction of the PE activities of the new FI-Zr and FI-Ni catalysts was carried out by applying QSAR model Equations (2) and (3), as listed in Table 5. Furthermore, to investigate the PE polymerization tendency performance, the results of the overall PE catalytic activities versus the Ti, Zr, and Ni phenoxy-imine complexes were plotted using Equations (2) and (3), as depicted in Figure 4a,b, respectively. Based on the graphical plots, the trend of higher PE activity was found to follow the order of FI-Ni > FI-Ti > FI-Zr, which agreed well with both QSAR Equations (2) and (3).
In Table 5, assuming the purposed metallocenes for PE catalysts could obtain the theoretical PE activities prediction of higher than 35,000 kg(PE)/mol(Cat.)·MPa·h. Thus, an insight into the details of predicted PE activities for the FI-Zr and FI-Ni series based on Equation (2) consisted of 7Zr, 2Ni, 3Ni, 4Ni, 6Ni, 7Ni, 10Ni, 12Ni, 13Ni, and 15Ni. Whereas the predictions from Equation (3) were found at 9Zr, 2Ni, 3Ni, 10Ni and 11Ni. Therefore, by the intersection results of Equations (2) and (3), the new proposed catalysts for PE polymerization are 2Ni, 3Ni, and 10Ni and they are expected to be potential catalysts with predicted PE catalytic activities (Equation (3)) of 47324, 37871, and 41383 kg(PE)/mol(Cat.)·MPa·h, respectively. When looking into the details of these proposed catalysts, our QSAR finding suggested that the higher HOMO energy level and the negative total charge of R3 substituents tend to increase the PE activity.

3. Materials and Methods

3.1. Source of Experimental Datasets

For the experimental data of PE activity over the FI-Ti catalysts, the time-averaged activity of the FI-Ti catalysts in ethylene polymerization is defined as:
activity = M C 2 H 4 0 τ W p ( t ) d t C T i p C 2 H 4 τ
where M(C2H4) and p(C2H4) denote the molecular weight and the pressure of ethylene, respectively; CTi is the initial concentration of a catalyst, τ is polymerization time, and Wp(t) is the propagation rate given by:
W p = k e f f [ T i ] a c t [ C 2 H 4 ]
where keff is the effective rate constant (activity of a single active site) and [Ti]act and [C2H4] denote the concentrations of active sites and ethylene, respectively. For the QSAR study, the assumption of the same conditions for obtaining PE polymerization activity is required; thus, the catalytic activity should be determined by a composite function of keff [Ti]act, as reported by Nikitin et al. [35]. The FI-Ti catalytic activities for PE polymerization were given in kg(PE)/mol(Cat.)·MPa·h units, and the temperature conditions were chosen at 40 °C except 1Ti, 4Ti, 7Ti, 10Ti, 13Ti and 16Ti, which were estimated by their 40 °C performance condition by (activity(30 °C) + activity(50 °C))/2. Therefore, the FI-Ti catalytic activities for PE polymerization in this work were collected from Nikitin’s work and they are listed in Table 1.

3.2. Structural Optimization

All the calculations were performed with the M06L density functional method using the Gaussian 09 software package (Gaussian, Inc., Wallingford, CT, USA, 2009) [60]. The 6-31G** basis set was used for C, O, N, H, and Cl, while LANL2DZ was used for Ti, Zr, and Ni. The structures of the metal complexes were optimized and their energies were calculated to extract electronic and steric properties. Ten electronic properties were collected: the charges of metal centers with Ti/Zr/Ni atoms, the average charge of O, N, and Cl, charges of substitution groups on phenoxy-imine ligands of R1, R2, and R3, HOMO energy, LUMO energy, and bond distances of metal centers to surrounding O, N, and Cl atoms.

3.3. QSAR Model

In QSAR studies, the choices of descriptors are the important steps to achieve the model with high correlation and prediction ability [61]. Hence, the GA method as an effective technique to select the appropriate variables was used in this work. Catalytic activity was defined as the dependent variable, and there were a total of 10 molecular electronic properties (descriptors), including atomistic charges and quantum-chemical descriptors, which were fitted via the leave-one-out cross-validation coefficient value as a fitness function. The genetic function approximation (GFA) algorithm implemented in the Materials Studio 7.0 package (Accelrys Software Inc., San Diego, CA, USA, 2013) was used to construct the QSAR models [62].

3.4. Statistical Terms

R2 is the fraction of the total variance of the Y variable (ln(A)) that is explained by the genetic function approximation equation. R2 is calculated as:
SSR SST
where SSR is the sum of squares of regression and SST is the total sum of squares. Generally, the closer the value is to 1.0, the better the genetic function approximation equation explains the Y variable.
R2cv is derived from cross validation process, in this work, we used the leave-one-out validation. The R2cv is calculated as follows:
1 PRESS SST
where PRESS is the predictive sum of squares of a model. The R2cv represents the predictive power of a model for further prediction on some new catalysts. For a good model, R2cv should be close to R2.
The following is used for the F test:
SSR / ( p 1 ) SSE / ( n p )
where SSE is the sum of squares of errors. This is the sum of squared differences between measurements, and model predictions over the entire data set. SSR is the sum of squares due to regression. n and p represent the number of data points and the number of parameters in a regression model, respectively. Thus, the ratio F = {SSR/(p − 1)}/{SSE/(np)} is then calculated and compared with the tabulated values of the F distribution for different values of n and p. The F value is derived from the F test, indicating the probability of a true relationship, or the significance level of the QSAR model. The F value is the ratio between explained and unexplained variance for a given number of degrees of freedom. Therefore, the larger the F value, the greater the probability that the QSAR equation is significant [62,63,64].

4. Conclusions

The phenoxy-imine catalysts were tailored by varying the substituent groups on the ligand units and also the metal center. Their activities were then calculated via DFT. After analyzing the raw data of the DFT simulations and experimental values of FI-Ti with the QSAR method, the calculated properties from DFT illustrated the significant relation between realistic catalytic activity and the calculated electronic descriptors. For the best obtained QSAR models, consisting of two models of Equations (2) and (3), their R2 and R2cv values are greater than 0.927. From insight into the QSAR models, the HOMO energy level and the total charge of R3 substituent electronic properties could be used to predict the PE activity of the FI-Ti catalysts. The design concept of “head-hat” was applied by replacing the head of Ti metal by Zr or Ni metals, where the hats were still kept as phenoxy-imine substituents. Based on the obtained QSAR models derived from the FI-Ti training set, the PE polymerization activities of the new metallocene catalysts of the FI-Zr and FI-Ni series were then predicted using QSAR Equations (2) and (3). To purpose the candidate catalysts, higher PE polymerization activities than 35,000 kg(PE)/mol(Cat.)·MPa·h on those QSAR models were considered, and it was found that based on all metallocene complexes that consist of Ni as a metal center atom showed better PE performance than that of the Ti and Zr metal center pheoxy-imines. Thus, the new purposed catalysts for PE polymerization are 2Ni, 3Ni and 10Ni metallocene complexes. The present work shows that the combined DFT and QSAR studies are straightforward methods to derive the relationship between PE catalytic activity and electronic descriptors, and the obtained QSAR model could be used for evaluating and designing new catalyst structures in search of promising polyethylene catalysts. This tailor-made study is beneficial to industrial research, and it is a simple method.

Author Contributions

Visualization and Simulations, P.C.; Conceptualization, P.M. and V.P.; Formal analysis, P.M., V.P., D.Z., and L.S.; Writing—original draft, P.C., H.W. and P.M.

Funding

Thailand Research Fund (RTA6080005).

Acknowledgments

This work was supported by Vidyasirimedhi Institute of Science and Technology (VISTEC) for both funds and equipment. Thanks goes to support from the Thailand Research Fund (RTA6080005) and Research Center of Nano Science and Technology, Shanghai University. We would like to thank NANOTEC, Thailand Sci Park for cluster computer and software supports. PM and LS would like to thank the Shanghai Municipal Science and Technology Commission of Professional and Technical Service Platform for Designing and Manufacturing of Advanced Composite Materials (16DZ2292100).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xie, H.; Hua, X.; Liu, B.; Wu, C.; Cui, D. Phosphinimino-amino supported complex: Synthesis, polymerization of ethylene and dearomatisation of pyridine. J. Organomet. Chem. 2015, 798, 335–340. [Google Scholar] [CrossRef]
  2. Nyamato, G.S.; Ojwach, S.O.; Akerman, M.P. Ethylene oligomerization studies by nickel(II) complexes chelated by (amino)pyridine ligands: Experimental and density functional theory studies. Dalton Trans. 2016, 45, 3407–3416. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, Z.; Do, L.H. Customizing polyolefin morphology by selective pairing of alkali ions with nickel phenoxyimine-polyethylene glycol catalysts. Organometallics 2017, 36, 4691–4698. [Google Scholar] [CrossRef]
  4. Yu, F.; Yang, Y.; He, D.; Gong, D.; Chen, Z.-R. Pressure-sensitive supported FI catalyst for the precise synthesis of uni- and bimodal polyethylene. Ind. Eng. Chem. Res. 2017, 56, 4684–4689. [Google Scholar] [CrossRef]
  5. Goerl, C.; Betthausen, E.; Alt, H.G. Di- and trinuclear iron/titanium and iron/zirconium complexes with heterocyclic ligands as catalysts for ethylene polymerization. Polyhedron 2016, 118, 37–51. [Google Scholar] [CrossRef]
  6. Woo, J.O.; Kang, S.K.; Park, J.-E.; Son, K.-S. Synthesis, characterization, and ethylene polymerization behavior of cr(III) catalysts based on bis(pyrazolylmethyl)pyridine and its derivatives. J. Mol. Catal. A Chem. 2015, 404–405, 204–210. [Google Scholar] [CrossRef]
  7. Steelman, D.K.; Xiong, S.; Pletcher, P.D.; Smith, E.; Switzer, J.M.; Medvedev, G.A.; Delgass, W.N.; Caruthers, J.M.; Abu-Omar, M.M. Effects of pendant ligand binding affinity on chain transfer for 1-hexene polymerization catalyzed by single-site zirconium amine bis-phenolate complexes. J. Am. Chem. Soc. 2013, 135, 6280–6288. [Google Scholar] [CrossRef] [PubMed]
  8. Martínez, S.; Cruz, V.L.; Ramos, J.; Martínez-Salazar, J. Polymerization activity prediction of zirconocene single-site catalysts using 3d quantitative structure–activity relationship modeling. Organometallics 2012, 31, 1673–1679. [Google Scholar] [CrossRef]
  9. Kim, S.; Park, J.; Song, B.; Yoon, S.-W.; Go, M.; Lee, J.; Lee, B. Preparation of thiophene-fused and tetrahydroquinoline-linked cyclopentadienyl titanium complexes for ethylene/α-olefin copolymerization. Catalysts 2013, 3, 104. [Google Scholar] [CrossRef]
  10. Shamiri, A.; Chakrabarti, M.; Jahan, S.; Hussain, M.; Kaminsky, W.; Aravind, P.; Yehye, W. The influence of ziegler-natta and metallocene catalysts on polyolefin structure, properties, and processing ability. Materials 2014, 7, 5069. [Google Scholar] [CrossRef] [PubMed]
  11. Alt, H.G.; Köppl, A. Effect of the nature of metallocene complexes of group iv metals on their performance in catalytic ethylene and propylene polymerization. Chem. Rev. 2000, 100, 1205–1222. [Google Scholar] [CrossRef] [PubMed]
  12. Luo, W.; Li, A.; Liu, S.; Ye, H.; Li, Z. 2-benzimidazol-6-pyrazol-pyridine chromium(III) trichlorides: Synthesis, characterization, and application for ethylene oligomerization and polymerization. Organometallics 2016, 35, 3045–3050. [Google Scholar] [CrossRef]
  13. Bariashir, C.; Wang, Z.; Du, S.; Solan, G.A.; Huang, C.; Liang, T.; Sun, W.-H. Cycloheptyl-fused nno-ligands as electronically modifiable supports for M(II) (M = Co, Fe) chloride precatalysts; probing performance in ethylene oligo-/polymerization. J. Polym. Sci. Part A Polym. Chem. 2017, 55, 3980–3989. [Google Scholar] [CrossRef]
  14. 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]
  15. Park, J.-E.; Kang, S.K.; Woo, J.O.; Son, K.-S. Highly active chromium(III) complexes based on tridentate pyrazolyl pyridyl ligands for ethylene polymerization and oligomerization. Dalton Trans. 2015, 44, 9964–9969. [Google Scholar] [CrossRef] [PubMed]
  16. Nyamato, G.S.; Alam, M.G.; Ojwach, S.O.; Akerman, M.P. (pyrazolyl)-(phosphinoyl)pyridine iron(II), cobalt(II) and nickel(II) complexes: Synthesis, characterization and ethylene oligomerization studies. J. Organomet. Chem. 2015, 783, 64–72. [Google Scholar] [CrossRef]
  17. Falivene, L.; Cavallo, L.; Talarico, G. Buried volume analysis for propene polymerization catalysis promoted by group 4 metals: A tool for molecular mass prediction. ACS Catal. 2015, 5, 6815–6822. [Google Scholar] [CrossRef]
  18. Nomura, K.; Mitsudome, T.; Igarashi, A.; Nagai, G.; Tsutsumi, K.; Ina, T.; Omiya, T.; Takaya, H.; Yamazoe, S. Synthesis of (adamantylimido)vanadium(V) dimethyl complex containing (2-anilidomethyl)pyridine ligand and selected reactions: Exploring the oxidation state of the catalytically active species in ethylene dimerization. Organometallics 2017, 36, 530–542. [Google Scholar] [CrossRef]
  19. Makio, H.; Terao, H.; Iwashita, A.; Fujita, T. Fi catalysts for olefin polymerization—A comprehensive treatment. Chem. Rev. 2011, 111, 2363–2449. [Google Scholar] [CrossRef] [PubMed]
  20. Mitani, M.; Furuyama, R.; Mohri, J.-I.; Saito, J.; Ishii, S.; Terao, H.; Nakano, T.; Tanaka, H.; Fujita, T. Syndiospecific living propylene polymerization catalyzed by titanium complexes having fluorine-containing phenoxy−imine chelate ligands. J. Am. Chem. Soc. 2003, 125, 4293–4305. [Google Scholar] [CrossRef] [PubMed]
  21. Mitani, M.; Mohri, J.-I.; Yoshida, Y.; Saito, J.; Ishii, S.; Tsuru, K.; Matsui, S.; Furuyama, R.; Nakano, T.; Tanaka, H.; et al. Living polymerization of ethylene catalyzed by titanium complexes having fluorine-containing phenoxy-imine chelate ligands. J. Am. Chem. Soc. 2002, 124, 3327–3336. [Google Scholar] [CrossRef] [PubMed]
  22. Saito, J.; Mitani, M.; Matsui, S.; Tohi, Y.; Makio, H.; Nakano, T.; Tanaka, H.; Kashiwa, N.; Fujita, T. A new titanium complex having two phenoxy-imine chelate ligands for ethylene polymerization. Macromol. Chem. Phys. 2002, 203, 59–65. [Google Scholar] [CrossRef]
  23. Terao, H.; Ishii, S.; Mitani, M.; Tanaka, H.; Fujita, T. Ethylene/polar monomer copolymerization behavior of bis(phenoxy-imine)ti complexes: Formation of polar monomer copolymers. J. Am. Chem. Soc. 2008, 130, 17636–17637. [Google Scholar] [CrossRef] [PubMed]
  24. Makio, H.; Kashiwa, N.; Fujita, T. Fi catalysts: A new family of high performance catalysts for olefin polymerization. Adv. Synth. Catal. 2002, 344, 477–493. [Google Scholar] [CrossRef]
  25. Wang, Z.; Solan, G.A.; Mahmood, Q.; Liu, Q.; Ma, Y.; Hao, X.; Sun, W.-H. Bis(imino)pyridines incorporating doubly fused eight-membered rings as conformationally flexible supports for cobalt ethylene polymerization catalysts. Organometallics 2018, 37, 380–389. [Google Scholar] [CrossRef]
  26. Ishii, S.; Nakano, T.; Kawamura, K.; Kinoshita, S.; Ichikawa, S.; Fujita, T. Development of new selective ethylene trimerization catalysts based on highly active ethylene polymerization catalysts. Catal. Today 2018, 303, 263–270. [Google Scholar] [CrossRef]
  27. Wang, Z.; Liu, Q.; Solan, G.A.; Sun, W.-H. Recent advances in ni-mediated ethylene chain growth: Nimine-donor ligand effects on catalytic activity, thermal stability and oligo-/polymer structure. Coord. Chem. Rev. 2017, 350, 68–83. [Google Scholar] [CrossRef]
  28. Preston, A.Z.; Kim, J.; Medvedev, G.A.; Delgass, W.N.; Caruthers, J.M.; Abu-Omar, M.M. Steric and solvation effects on polymerization kinetics, dormancy, and tacticity of Zr-salan catalysts. Organometallics 2017, 36, 2237–2244. [Google Scholar] [CrossRef]
  29. Mitchell, N.E.; Anderson, W.C., Jr.; Long, B.K. Mitigating chain-transfer and enhancing the thermal stability of co-based olefin polymerization catalysts through sterically demanding ligands. J. Polym. Sci. Part A. Polym. Chem. 2017, 55, 3990–3995. [Google Scholar] [CrossRef]
  30. D’Auria, I.; Milione, S.; Caruso, T.; Balducci, G.; Pellecchia, C. Synthesis of hyperbranched low molecular weight polyethylene oils by an iminopyridine nickel(II) catalyst. Polym. Chem. 2017, 8, 6443–6454. [Google Scholar] [CrossRef]
  31. Tuskaev, V.A.; Gagieva, S.C.; Kurmaev, D.A.; Kolosov, N.A.; Fedyanin, I.V.; Bulychev, B.M. Titanium (IV) complexes stabilized by 2,6-bis-(α,α′-diphenyl(hydroxy)methyl)-pyridine: Catalytic activity in olefin polymerization and impact of lithium and magnesium chlorides. Inorg. Chim. Acta 2015, 425, 275–281. [Google Scholar] [CrossRef]
  32. Ratanasak, M.; Rungrotmongkol, T.; Saengsawang, O.; Hannongbua, S.; Parasuk, V. Towards the design of new electron donors for ziegler–natta catalyzed propylene polymerization using qspr modeling. Polymer 2015, 56, 340–345. [Google Scholar] [CrossRef]
  33. Ratanasak, M.; Parasuk, V. Roles of malonate donor on activity and stereoselectivity of ziegler-natta catalyzed propylene polymerization. J. Organomet. Chem. 2015, 775, 6–11. [Google Scholar] [CrossRef]
  34. Nyamato, G.S.; Ojwach, S.O.; Akerman, M.P. Potential hemilabile (imino)pyridine palladium(II) complexes as selective ethylene dimerization catalysts: An experimental and theoretical approach. Organometallics 2015, 34, 5647–5657. [Google Scholar] [CrossRef]
  35. Nikitin, S.V.; Nikitin, V.V.; Oleynik, I.I.; Oleynik, I.V.; Bagryanskaya, E.G. Activity of phenoxy-imine titanium catalysts in ethylene polymerization: A quantum chemical approach. J. Mol. Catal. A Chem. 2016, 423, 285–292. [Google Scholar] [CrossRef]
  36. Yang, W.; Yi, J.; Ma, Z.; Sun, W.-H. 2d-qsar modeling on the catalytic activities of 2-azacyclyl-6-aryliminopyridylmetal precatalysts in ethylene oligomerization. Catal. Commun. 2017, 101, 40–43. [Google Scholar] [CrossRef]
  37. Manz, T.A. Deactivation of Ti and Zr half-metallocene complexes activated with B(C6F5)3: A case study in constructing DFT-based qsars to predict unimolecular rate constants. RSC Adv. 2015, 5, 48246–48254. [Google Scholar] [CrossRef]
  38. Cruz, V.L.; Martinez, S.; Ramos, J.; Martinez-Salazar, J. 3D-QSAR as a tool for understanding and improving single-site polymerization catalysts. Organometallics 2014, 33, 2944–2959. [Google Scholar] [CrossRef]
  39. Cruz, V.L.; Martinez, S.; Martinez-Salazar, J.; Polo-Cerón, D.; Gómez-Ruiz, S.; Fajardo, M.; Prashar, S. 3D-QSAR study of ANSA-metallocene catalytic behavior in ethylene polymerization. Polymer 2007, 48, 4663–4674. [Google Scholar] [CrossRef]
  40. Melanie, M. Genetic algorithms: An overview. Complexity 1995, 1, 31–39. [Google Scholar] [Green Version]
  41. McCall, J. Genetic algorithms for modelling and optimisation. J. Comput. Appl. Math. 2005, 184, 205–222. [Google Scholar] [CrossRef]
  42. Fisz, J.J. Combined genetic algorithm and multiple linear regression (GA-MLR) optimizer: Application to multi-exponential fluorescence decay surface. J. Phys. Chem. A 2006, 110, 12977–12985. [Google Scholar] [CrossRef] [PubMed]
  43. Rogers, D.; Hopfinger, A.J. Application of genetic function approximation to quantitative structure-activity relationships and quantitative structure-property relationships. J. Chem. Inf. Comp. Sci. 1994, 34, 854–866. [Google Scholar] [CrossRef]
  44. Zhao, Y.; Truhlar, D.G. A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. J. Chem. Phys. 2006, 125, 194101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Della Sala, F.; Fabiano, E.; Constantin, L.A. Kinetic-energy-density dependent semilocal exchange-correlation functionals. Int. J. Quantum Chem. 2016, 116, 1641–1694. [Google Scholar] [CrossRef] [Green Version]
  46. Zhao, Y.; Truhlar, D.G. The m06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four m06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar]
  47. Remya, K.; Suresh, C.H. Which density functional is close to ccsd accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using gaussian09. J. Comput. Chem. 2013, 34, 1341–1353. [Google Scholar] [CrossRef] [PubMed]
  48. Rakhi, R.; Suresh, C.H. A dft study on dihydropyrazine annulated linear polyacenes: Aromaticity, stability and homo-lumo energy modulation. PCCP 2016, 18, 24631–24641. [Google Scholar] [CrossRef] [PubMed]
  49. Ahmad, S.; Gromiha, M.M. Design and training of a neural network for predicting the solvent accessibility of proteins. J. Comput. Chem. 2003, 24, 1313–1320. [Google Scholar] [CrossRef] [PubMed]
  50. Hunger, J.; Huttner, G. Optimization and analysis of force field parameters by combination of genetic algorithms and neural networks. J. Comput. Chem. 1999, 20, 455–471. [Google Scholar] [CrossRef]
  51. Waller, C.L.; Bradley, M.P. Development and validation of a novel variable selection technique with application to multidimensional quantitative structure-activity relationship studies. J. Chem. Inf. Comp. Sci. 1999, 39, 345–355. [Google Scholar] [CrossRef]
  52. Aires-de-Sousa, J.; Hemmer, M.C.; Gasteiger, J. Prediction of 1H NMR chemical shifts using neural networks. Anal. Chem. 2002, 74, 80–90. [Google Scholar] [CrossRef] [PubMed]
  53. Thakur, P.; Srivastava, D.C.; Gupta, P.K. The genetic algorithm: A robust method for stress inversion. J. Struct. Geol. 2017, 94, 227–239. [Google Scholar] [CrossRef]
  54. Douguet, D.; Thoreau, E.; Grassy, G. A genetic algorithm for the automated generation of small organic molecules: Drug design using an evolutionary algorithm. J. Comput. Aided Mol. Des. 2000, 14, 449–466. [Google Scholar] [CrossRef] [PubMed]
  55. Kubiny, H. Variable selection in QSAR studies. I. An evolutionary algorithm. Quant. Struct. Act. Relat. 1994, 13, 285–294. [Google Scholar] [CrossRef]
  56. Leardi, R. 3-genetic algorithms in feature selection. In Genetic Algorithms in Molecular Modeling; Devillers, J., Ed.; Academic Press: Cambridge, MA, USA, 1996; pp. 67–86. [Google Scholar]
  57. Wold, S.; Eriksson, L.; Clementi, S. Statistical validation of QSAR results. In Chemometric Methods in Molecular Design; VCH, P.O. Box: Weinheim, Germany, 2008; pp. 309–338. [Google Scholar]
  58. Balakrishnan, N.; Colton, T.; Everitt, B.; Piegorsch, W.; Ruggeri, F.; Teugels, J.L. Ramp function. Wiley StatsRef Stat. Ref. Online 2014. [Google Scholar] [CrossRef]
  59. Mudelsee, M. Ramp function regression: A tool for quantifying climate transitions. Comput. Geosci. 2000, 26, 293–307. [Google Scholar] [CrossRef]
  60. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, revision B.01; Gaussian Inc.: Wallingford UK, 2009. [Google Scholar]
  61. Pourbasheer, E.; Aalizadeh, R.; Ganjali, M.R.; Norouzi, P.; Shadmanesh, J. QSAR study of ACK1 inhibitors by genetic algorithm-multiple linear regression (GA-MLR). J. Saudi Chem. Soc. 2014, 18, 681–688. [Google Scholar] [CrossRef]
  62. Materials Studio Modeling; release 7.0; Accelrys Software Inc.: San Diego, CA, USA, 2013.
  63. Smith, H.J.; Williams, H. Introduction to the Principles of Drug Design, 2nd ed.; The stonebrige press: Butterworth, Malaysian, 1988; pp. 240–264. [Google Scholar]
  64. Leardi, R. Genetic algorithms in chemometrics and chemistry: A review. J. Chemom. 2001, 15, 559–569. [Google Scholar] [CrossRef]
Figure 1. Occurrences of populations of 10 electronic descriptors by the genetic algorithm (GA) method.
Figure 1. Occurrences of populations of 10 electronic descriptors by the genetic algorithm (GA) method.
Catalysts 08 00422 g001
Figure 2. A scaled residual plot of training set obtained from Equation (2) for FI-Ti catalysts.
Figure 2. A scaled residual plot of training set obtained from Equation (2) for FI-Ti catalysts.
Catalysts 08 00422 g002
Figure 3. The plot of the experimental PE activities versus the predicted ones obtaining from Equation (2) (a), and Equation (3) (b) for FI-Ti catalysts.
Figure 3. The plot of the experimental PE activities versus the predicted ones obtaining from Equation (2) (a), and Equation (3) (b) for FI-Ti catalysts.
Catalysts 08 00422 g003
Figure 4. Trend of the predicted PE activity of all FI catalysts via Equation (2) (a) and Equation (3) (b).
Figure 4. Trend of the predicted PE activity of all FI catalysts via Equation (2) (a) and Equation (3) (b).
Catalysts 08 00422 g004aCatalysts 08 00422 g004b
Table 1. Chemical structure of FI-Ti catalysts with different substituents R1, R2, and R3 (data cited from reference [35]).
Table 1. Chemical structure of FI-Ti catalysts with different substituents R1, R2, and R3 (data cited from reference [35]).
Catalysts 08 00422 i001
ComplexR1R2R3Activity (A) (kg(PE)/mol(Cat.)·MPa·h)ln(A)
1Tit-BuHH41003.613
2Tit-BuHm-OAll17,8004.250
3Tit-BuHp-OAll54,2004.734
4Tit-BuMeH17,1004.233
5Tit-BuMem-OAll13,6004.134
6Tit-BuMep-OAll39,0004.591
7Tit-But-BuH19,4004.288
8Tit-But-Bum-OAll12,4004.093
9Tit-But-Bup-OAll33,0004.519
10TiCumylHH27,4004.438
11TiCumylHm-OAll32,6004.513
12TiCumylHp-OAll34,5004.538
13TiCumylMeH34,0504.532
14TiCumylMem-OAll19,8004.297
15TiCumylMep-OAll29,4004.468
16TiCumylCumylH18,5504.268
17TiCumylCumylm-OAll26,0004.415
18TiCumylCumylp-OAll25,7004.410
Table 2. FI-Ti electronic descriptor data.
Table 2. FI-Ti electronic descriptor data.
Catalysts 08 00422 i002
ComplexTi aO bN bCl bR1 cR2 cR3 cHOMO dLUMO dE-Gap d
1Ti1.235−0.704−0.586−0.374−1.0620.1360.136−0.199−0.1311.840
2Ti1.232−0.703−0.585−0.374−1.0620.345−1.117−0.197−0.1291.858
3Ti1.247−0.706−0.588−0.376−1.0570.130−1.230−0.194−0.1271.843
4Ti1.234−0.707−0.589−0.377−1.058−0.4500.134−0.193−0.1291.736
5Ti1.227−0.705−0.589−0.378−1.062−0.446−1.109−0.191−0.1261.770
6Ti1.250−0.709−0.592−0.380−1.055−0.442−1.134−0.191−0.1241.812
7Ti1.252−0.713−0.610−0.372−1.056−1.0830.134−0.197−0.1291.863
8Ti1.235−0.707−0.588−0.380−1.033−1.054−1.097−0.191−0.1241.817
9Ti1.236−0.707−0.586−0.381−1.030−1.054−1.116−0.187−0.1241.720
10Ti1.243−0.708−0.607−0.362−1.1470.1390.136−0.198−0.1331.765
11Ti1.273−0.706−0.603−0.372−1.3030.137−1.129−0.199−0.1301.867
12Ti1.232−0.694−0.578−0.374−1.2540.122−1.131−0.189−0.1241.765
13Ti1.242−0.696−0.584−0.376−1.248−0.0260.133−0.189−0.1251.739
14Ti1.269−0.707−0.606−0.376−1.302−0.445−1.125−0.191−0.1271.764
15Ti1.224−0.700−0.579−0.376−1.257−0.432−1.121−0.185−0.1231.688
16Ti1.286−0.703−0.538−0.385−1.222−1.3530.120−0.189−0.1211.856
17Ti1.236−0.705−0.602−0.370−1.216−1.350−1.101−0.189−0.1231.794
18Ti1.270−0.709−0.598−0.384−1.262−1.305−1.102−0.190−0.1251.750
a = partial charge (e), b = average charge (e), c = average of total charge of substituent group (e), d = in Ha unit.
Table 3. Correlation matrix of polyethylene (PE) activity in term of ln(A) and 10 electronic properties. High correlations (≥±0.8) are labelled in gray cells.
Table 3. Correlation matrix of polyethylene (PE) activity in term of ln(A) and 10 electronic properties. High correlations (≥±0.8) are labelled in gray cells.
Des.ln(A)Ti aO bN bCl bR1 cR2 cR3 cHOMOLUMOE-Gap
ln(A)1.0000.1550.119−0.0840.033−0.2980.041−0.3800.3550.327−0.225
Tia0.1551.000−0.3120.119−0.294−0.468−0.2900.114−0.1560.0400.391
Ob0.119−0.3121.0000.175−0.003−0.2680.027−0.1690.4490.296−0.472
Nb-0.0840.1190.1751.000−0.568−0.032−0.1280.1590.3970.5500.043
Clb0.033−0.294−0.003−0.5681.000−0.0650.5010.219−0.524−0.6360.087
R1c-0.298−0.468−0.268−0.032−0.0651.0000.0020.112−0.280−0.2250.230
R2c0.041−0.2900.027−0.1280.5010.0021.0000.017−0.480−0.5670.103
R3c-0.3800.114−0.1690.1590.2190.1120.0171.000−0.347−0.3960.096
HOMO0.355−0.1560.4490.397−0.524−0.280−0.480−0.3471.0000.890−0.684
LUMO0.3270.0400.2960.550−0.636−0.225−0.567−0.3960.8901.000−0.275
E-Gap-0.2250.391−0.4720.0430.0870.2300.1030.096−0.684−0.2751.000
a = partial charge (e), b = average charge (e), c = average of total charge of substituent group (e).
Table 4. Statistical results of quantitative structure-activity relationship (QSAR) equations on FI-Ti series.
Table 4. Statistical results of quantitative structure-activity relationship (QSAR) equations on FI-Ti series.
NoEquation cR2R2cvF Value
Equation (1)ln(Activity) = 26.09 × (R3 Charge + 0.22) − 7.47 × (R3 Charge + 1.18) + 497.87 × (HOMO + 0.19) − 494.72 × (HOMO + 0.19) − 7.81 × ( 1.24 − Ti-Charge) + 1943.36 × (−0.20 - HOMO) − 3345.07 × (−0.20 − HOMO) + 4.730.9920.936168.295
Equation (2) aln(Activity) = −4.91 × (R3 Charge + 1.22) + 23.54 × (R3 Charge + 0.13) + 33949.84 × (−0.19 − HOMO) − 2361.91 × (−0.20 − HOMO) + 743.80 × (−0.20 − HOMO) − 33941.61 × (−0.19 − HOMO) + 5.000.9920.984230.917
Equation (3) bln(Activity) = −9.15 × (R3 Charge + 1.17) + 30.70 × (R3 Charge + 0.24) − 257.00 × (HOMO + 0.189) + 225.65 × (HOMO + 0.19) −1938.84 × (−0.20 − HOMO) + 533.89 * (−0.20 − HOMO) + 4.740.9970.927564.305
a reduced Ti-charge descriptor, b deleted outliers, c linearity ramp function equation.
Table 5. Electronic descriptors and predicted activities from Equations (2) and (3) of newly designed FI-Zr and FI-Ni catalysts. The higher predicted activities of PE polymerization more than 35,000 kg(PE)/mol(Cat.) MPa·h are highlighted in red numbers.
Table 5. Electronic descriptors and predicted activities from Equations (2) and (3) of newly designed FI-Zr and FI-Ni catalysts. The higher predicted activities of PE polymerization more than 35,000 kg(PE)/mol(Cat.) MPa·h are highlighted in red numbers.
ComplexSubstituentsDescriptorsPredicted PE Activity a
R1R2R3R3-ChargeHOMO
Equation (2)Equation (3)
2-Zrt-BuHm-OAll−1.115−0.19912,26011,474
3-Zrt-BuHp-OAll−1.124−0.19418,13419,368
4-Zrt-BuMeH0.135−0.19419,14117,009
5-Zrt-BuMem-OAll−1.105−0.19715,33512,979
6-Zrt-BuMep-OAll−1.119−0.19031,76424,327
7-Zrt-But-BuH0.135−0.19843,69732,147
8-Zrt-But-Bum-OAll−1.096−0.19212,70910,737
9-Zrt-But-Bup-OAll−1.117−0.18930,87936,923
10-ZrCumylHH0.136−0.19925,3102709
11-ZrCumylHm-OAll−1.128−0.19720,01521,072
12-ZrCumylHp-OAll−1.124−0.19117,23919,574
13-ZrCumylMeH0.134−0.19133,90121,331
14-ZrCumylMem-OAll−1.123−0.19317,62418,965
15-ZrCumylMep-OAll−1.116−0.18730,70531,103
16-ZrCumylCumylH0.119−0.19017,81712,411
17-ZrCumylCumylm-OAll−1.096−0.18924,49024,515
18-ZrCumylCumylp-OAll−1.101−0.18925,76827,286
2-Nit-BuHm-OAll−1.132−0.18836,58747,324
3-Nit-BuHp-OAll−1.136−0.18338,49737,871
4-Nit-BuMeH0.138−0.18340,24528,649
5-Nit-BuMem-OAll−1.104−0.18226,80917,645
6-Nit-BuMep-OAll−1.135−0.17737,84922,786
7-Nit-But-BuH0.138−0.18340,24527,974
10-Nit-But-Bum-OAll−0.140−0.18743,84941,383
11-Nit-But-Bup-OAll−1.126-0.18634,18735,222
12-NiCumylHH−1.145−0.18042,38034,431
13-NiCumylHm-OAll0.137−0.17938,55521,053
14-NiCumylHp-OAll−1.124−0.18133,42224,122
15-NiCumylMeH−1.132−0.17636,79419,909
16-NiCumylMem-OAll0.130−0.17828,55713,369
17-NiCumylMep-OAll−1.101−0.17625,91510,772
18-NiCumylCumylH−1.098−0.17825,05011,616
a in units of kg(PE)/mol(Cat.)·MPa·h.

Share and Cite

MDPI and ACS Style

Chasing, P.; Maitarad, P.; Wu, H.; Zhang, D.; Shi, L.; Promarak, V. Straightforward Design for Phenoxy-Imine Catalytic Activity in Ethylene Polymerization: Theoretical Prediction. Catalysts 2018, 8, 422. https://doi.org/10.3390/catal8100422

AMA Style

Chasing P, Maitarad P, Wu H, Zhang D, Shi L, Promarak V. Straightforward Design for Phenoxy-Imine Catalytic Activity in Ethylene Polymerization: Theoretical Prediction. Catalysts. 2018; 8(10):422. https://doi.org/10.3390/catal8100422

Chicago/Turabian Style

Chasing, Pongsakorn, Phornphimon Maitarad, Hongmin Wu, Dengsong Zhang, Liyi Shi, and Vinich Promarak. 2018. "Straightforward Design for Phenoxy-Imine Catalytic Activity in Ethylene Polymerization: Theoretical Prediction" Catalysts 8, no. 10: 422. https://doi.org/10.3390/catal8100422

APA Style

Chasing, P., Maitarad, P., Wu, H., Zhang, D., Shi, L., & Promarak, V. (2018). Straightforward Design for Phenoxy-Imine Catalytic Activity in Ethylene Polymerization: Theoretical Prediction. Catalysts, 8(10), 422. https://doi.org/10.3390/catal8100422

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