Geometry and Stability of Titanium Chloride Species Adsorbed on the

Geometry and Stability of Titanium Chloride Species Adsorbed on the. (100) and (110) Cuts of the MgCl2 Support of the Heterogeneous. ZieglerrNatta Cat...
12 downloads 4 Views 198KB Size
Macromolecules 2000, 33, 8953-8962

8953

Geometry and Stability of Titanium Chloride Species Adsorbed on the (100) and (110) Cuts of the MgCl2 Support of the Heterogeneous Ziegler-Natta Catalysts Guglielmo Monaco,† Massimiliano Toto,† Gaetano Guerra,‡ Paolo Corradini,† and Luigi Cavallo*,† Dipartimento di Chimica, Universita` di Napoli, Complesso Monte S. Angelo, Via Cintia, I-80126 Naples, Italy, and Dipartimento di Chimica, Universita` di Salerno, Via S. Allende, Baronissi, Salerno, Italy Received June 6, 2000; Revised Manuscript Received September 20, 2000 ABSTRACT: Possible structures of TiCl4 molecules and TiCl3 fragments adsorbed on (110) and (100) faces of MgCl2, simulated by clusters of different size and shape, have been studied in the framework of density functional theory. For both monomeric TiCl4 and TiCl3, coordination on the (110) face is favored relative to coordination on the (100) face. TiCl3 fragments can bind together on the (100) face either forming or not Ti-Ti bonds, resulting in the formation of polynuclear TinCl3n species. The steric environment of the ending Ti atoms of such polynuclear species, with n > 2, is extremely similar to that of the C2 symmetric sites proposed several years ago for TiCl3-based catalytic systems and presents a strict analogy with the well-established models for isospecific polymerization with catalytic systems based on C2 symmetric metallocenes. This analogy holds also in the case of a TiCl3 fragment adsorbed on the (110) cut when both of its vicinal positions are occupied. When just one of the vicinal positions is occupied, sites of C1 symmetry can be formed, which have two minimum-energy structures with distinct positions (inward and outward) for the dangling chlorine. These inward and outward geometries can be expected to interconvert easily, the inward arrangement being favored.

1. Introduction The world market for polypropylene is currently over 20 × 106 tons, and more than 80% of the global polymer1 is obtained with heterogeneous Ziegler-Natta catalysts.2-4 Over the years, these catalysts have evolved from simple TiCl3 crystals into the current systems based on MgCl2 as a support for TiCl4. Different routes have been developed for the preparation of the supported catalysts. The mechanical route is based on comilling MgCl2, TiCl4, and a Lewis base (the so-called internal donor) for several hours. Alternatively, chemical routes can be used both to generate active MgCl2 and to incorporate the Ti compound and the Lewis base. Usually, a MgCl2 precursor is treated with the Lewis base and excess TiCl4 above 80 °C. Among internal donors, aromatic esters (ethyl benzoate in particular), alkoxysilanes, and diethers have been shown to be particularly suitable. Finally, the catalyst is activated by addition of alkylating reducing speciessAl(C2H5)3 being the most used onespossibly mixed with a second electron donor (the so-called external donor). The resulting active system is of extreme chemical complexity, and the intimate nature of the active sites is still a matter of debate. The characterization of heterogeneous Ziegler-Natta catalysts has been attempted by means of several experimental approaches,5-26 but definitive answers have not been achieved yet. As for the support structure, the primary particles of activated MgCl2 are composed of a few irregularly stacked Cl-Mg-Cl sandwiches of the kind shown in Figure 1.6 Observations on MgCl2 microcrystals by optical microscopy5 and by HRTEM12 †

Universita` di Napoli. Universita` di Salerno. * Corresponding author. E-mail: [email protected].



Figure 1. Schematic representation of MgCl2. The full line Cl atoms represent the Cl atoms above the plane of the Mg atoms, whereas the dashed Cl atoms represent the Cl atoms below the plane of the Mg atoms. Possible (100) and (110) lateral cuts with 5- and 4-coordinated Mg atoms are also indicated.

confirmed the proposed27 copresence of (100) and (110) lateral cuts. For electroneutrality reasons, these two lateral cuts contain coordinatively unsaturated Mg2+ ions with coordination number 4 and 5 on the (110) and (100) cuts, respectively, as shown in Figure 1. As for treatment with TiCl4, Corradini et al. have proposed that preferential coordination would occur on lateral cuts, leading to adsorbed TiCl4 species on (100) and (110) cuts as well as to adsorbed Ti2Cl8 dimers on the (100) cut.27 In fact, MgCl2 has crystal structures very similar to those of violet titanium trichloride. This dictates the possibility of an epitactical coordination of TiCl4 units (or TiCl3 units, after reduction) on the lateral coordinatively unsatured faces of MgCl2 crystals, giving rise to reliefs crystallographically coherent with the matrix. Moreover, EXAFS studies7,8 are consistent with this model pointing to the coordination of TiCl4 as a monomer7 or a dimer on the (100) face.8 The knowledge of catalyst structure after reduction with alkylaluminum compounds is also incomplete. As

10.1021/ma000988h CCC: $19.00 © 2000 American Chemical Society Published on Web 11/28/2000

8954

Monaco et al.

for the Ti oxidation state, literature reports are often contradictory, owing to the different catalysts and analytical methods used.7,9,13-20 The only reasonable conclusion is that under polymerization conditions a considerable reduction of Ti(IV) takes place, not only to Ti(III) but to Ti(II) as well. However, Ti(II) is usually considered not to be active for propene polymerization.19 The structure of TiCl3 fragments adsorbed on the support is not known. Similarly to isolated TiCl4 molecules, isolated TiCl3 fragments could be epitactically adsorbed on the (110) and (100) cuts, while dimeric Ti2Cl6 species could be grown on the (100) cut. These species could be obtained by direct reduction of the corresponding TiCl4 and Ti2Cl8 adsorbed species.27,28 However, rearrangements of TiCl3 species on the surface are also possible, and ESR experiments suggest the formation of polynuclear species of the type TimCln.17,18 It is also worth noting that these experiments indicate that most of the Ti(III) is ESR silent.17,18,20 Useful information on the nature of active sites can also be obtained by 13C NMR analyses of the polymers produced.28,29 These analyses indicate that some active sites can interconvert in a time shorter than the average time of chain growth21-26 and that at least one of the active sites has an environment of C1 symmetry.24,26 Although the whole framework is extremely complicated, some theoretical efforts have been attempted in this field. Most of the work concentrated on the elementary steps of the polymerization reaction,30-34 and the ab initio investigations of Pakkanen and co-workers on the interactions occurring between TiCl4 and/or Lewis bases and MgCl2 clusters represent an exception.35-39 However, for TiCl4 they investigated small complexes, and thus, although valuable, their study cannot be conclusive. A more detailed study of the interactions between TiCl4 and the (110) MgCl2 face is reported in a paper by Parrinello and co-workers,40 which used first principles molecular dynamics to investigate possible structures of TiCl4 species adsorbed on the (110) face.40 They found that TiCl4 can adopt various geometries, the one corresponding to the TiCl4 epitactically grown on the MgCl2 (110) face being the most stable. However, according to the authors, a higher energy 5-coordinated geometry (similar to that already considered by Pakkanen and co-workers)35 would be kinetically favored. Considering the remarkable contribution given by first principles computational approaches to the understanding of the behavior of the homogeneous ZieglerNatta catalysts, we have been spurred to use these tools to study the industrially relevant heterogeneous catalysts. For this reason, in this paper we attempted a systematic density functional theory (DFT) investigation of possible structures formed by adsorption of TiCl4 on both faces of MgCl2, and we also investigated the possible formation of dimeric Ti2Cl8 species on the (100) face. We also attempted a characterization of adsorbed Ti species after reduction of Ti(IV) to Ti(III). For this reason, we investigated possible structures of TiCl3 fragments adsorbed on the (110) face and of monomeric, dimeric, and tetrameric TiCl3 species on the (100) face.41 Since very little is known about these systems at molecular level, only a systematic comparison between the various species involved can be of some use. The main aim of this study is a simple structural and energetical characterization of possible adsorbed Ti species from a static point of view and a consideration of their possible stereospecific behavior. When analyzing

Macromolecules, Vol. 33, No. 24, 2000

the results, it should be remembered that the evaluated energy differences only represent a contribution to the total free energy differences, since solvent and entropic effects are not considered. Moreover, the role of the alkylaluminum and that of the Lewis bases have not been considered. Of course, the presence of these species is needed to achieve an active and highly stereospecific catalyst.42 However, activation and deactivation (mainly due to reactions of adsorbed Ti species with alkylaluminum compounds and Lewis bases)42 are out of the scope of the present work, which only focuses on the structure of possible precursors of active sites. For the same reason, adsorption of Ti(II) species has not been considered, since Ti(II) is generally thought to be inactive in propene polymerization.19 As for the role played by Lewis bases in modifying the catalysts behavior, it has been the subject of a complementary study which has been reported elsewhere.43 2. Models and Methods Except when specifically indicated, the geometry of the MgCl2 clusters has been kept fixed. The Mg-Cl distances and all the Cl-Mg-Cl angles have been set equal to 2.49 Å and 90°.44 This corresponds to the simplifying assumption that the atoms on the surface present a structure close to that present in the bulk of the crystal. A different approach could consider a total relaxing of all the degrees of freedom. However, due to the dimension of the considered clusters, we think that this last approach would imply deformations that would not be representative of the rearrangements that occur on the MgCl2 surface. To check for the reasonability of these assumptions, test calculations with a partial relaxation of the MgCl2 cluster will also be discussed. With regards to the size of the chosen clusters, we decided to consider only situations in which the Mg atoms that are involved in the coordination of the Ti species should have (previous to Ti adsorption) coordination numbers 4 and 5 on the (110) and (100) lateral cuts, respectively, as requested by electroneutrality reasons (see Figure 1). Moreover, all the Cl atoms of the cluster involved in the coordination of the Ti species should bridge two Mg atoms (see Figure 1). These minimum requirements lead to the Mg9Cl18 and Mg8Cl16 clusters as models for the (100) and (110) lateral cuts, respectively, shown in Figure 2. These clusters represent a reasonable compromise between reliability and computational efficacy. For the sake of consistency, these clusters were used throughout the paper. As we considered adsorption of tetrameric Ti species, however, Mg atoms that do not respect the above requirements are involved. Only for this case, we used the larger Mg13Cl26 cluster. Computational Details. The density functional calculations were carried out by using the Amsterdam Density Functional package ADF, developed by Baerends et al.45-47 The geometry optimization procedure has been developed by Versluis and Ziegler.48 The electronic configurations of the systems were described by an uncontracted triple-ζ STO basis set on Ti(3s, 3p, 3d, 4s, 4p) and Mg(3s, 3p, 3d).49 A double-ζ STO basis sets augmented with a single 3d polarization function was used for Cl(3s, 3p + 3d).49 The 1s22s22p6 configuration on Ti, Mg, and Cl was treated by the frozen-core approximation. Energies and geometries have been obtained by using the local potential by Vosko et al.,50 augmented in a self-consistent manner by the Becke’s exchange gradient correction51 and by the Perdew’s correlation gradient correction.52,53 Whenever systems with unpaired electrons were considered, an unrestricted formalism has been used. Although transition metal-ligand dissociation energetics obtained with the DFT approach outlined above have been shown to be correct within 5 kcal/mol of the experimental result (usually overestimated in terms of absolute size),54,55 for bridged-bonded aluminum compounds DFT methods have been shown to underestimate the heats of formation.56 However, a recent comparative computational study has shown

Macromolecules, Vol. 33, No. 24, 2000

Titanium Chloride Species 8955 EB ) (E{A} + ES - E{A}/S)/m

(1)

where E{A}, ES, and E{A}/S are the electronic energies of the isolated species which compose {A}, of the isolated support (either relaxed or unrelaxed) and of the system formed upon adsorption of the Ti species, respectively, while m is the number of Ti atoms present in {A}. Relative binding energies ∆EB have been computed assuming two reference structures: TiCl4 and TiCl3 adsorbed on a cluster that simulates the (100) face. The sign convention is such that positive ∆EB correspond to energy gain. Often the adsorption leads to a significant change of the geometry. To investigate the influence that such change in geometry has on the values of the adsorption energies, the binding energy has been decomposed as follows:59

EB ) EIntra + EInter

(2)

where EInter ) (E*{A} + E*S - E{A}/S)/m is the energy per Ti atom needed to dissociate the adsorbed species {A} avoiding any relaxation from the geometry found upon adsorption, and EIntra ) EIntra{A} + EIntra(S) is the sum of the energy contributions (E{A} - E*{A})/m and (ES - E*S)/m required to restore the equilibrium geometry of the adsorbed species and that of the support. Excluding two test calculations, the geometry of the support has been kept fixed and thus ES ) E*S.

3. Results and Discussion

Figure 2. Optimized geometries of structures 1a, 1b, 1c, and 1d, corresponding to a TiCl4 molecule adsorbed on the (100) face of a Mg9Cl18 cluster, and to 6-fold, 4-fold and 5-fold TiCl4 molecules adsorbed on the (110) face of a Mg8Cl16 cluster, respectively. Only the positions of the TiCl4 atoms have been optimized. 1d has been obtained without inclusion of exchangecorrelation gradient corrections (see methods). Hollow bonds are used to indicate bonds between the TiCl4 and the Mg8Cl16 cluster. that the DFT functional we have chosen is in excellent agreement with the high-level CSSD(T) coupled cluster calculations performed by Jensen and Børve.57 Energy Decomposition. In this paper we will be concerned with the adsorption of groups of TiCl4 molecules and/ or TiCl3 fragments on supports, S, formed by MgCl2 clusters. We will indicate as binding energy per Ti atom, EB, of a given group of species {A} only the electronic contribution to the binding free energy GB. The free and adsorbed species can be expected to have markedly different nonelectronic contribution to their free energy. Therefore, the values computed for EB cannot be used to determine the percentage of adsorbed species. However the relative binding energies can be expected to be substantially more reliable58 and can be reasonably used to get an idea of the preferred arrangements of the adsorbed species. For a given group of species {A} and a given support S, the binding energy EB per Ti atom is computed as

The starting points of our calculations have been the models of epitactic adsorption of Ti chlorides species proposed long ago by Corradini et al.27 As will be seen, the geometries found in the calculations are generally very similar to epitactic models, although there are deviations that are of interest. These deviations can be attributed to the fact that the length of a Ti-Cl bond (∼2.2 Å) is considerably shorter than the length of a Mg-Cl bond (∼2.5 Å), and thus epitactic placements require a deformation of the adsorbates (negative EIntra). Moreover, unbridged Cl atoms have a stronger trans influence than bridged ones,60 and this is a further source of deformation. To get an easy visualization of the deformations, it is possible to focus only on the shortest Ti-Cl bonds (say shorter than ∼3 Å). Finally, the relative binding energies (∆EB) calculated for Ti(IV) (upper part) and Ti(III) (lower part) species on (100) and (110) lateral cuts are collected in Table 1, together with their intramolecular (∆EIntra) and intermolecular (∆EInter) contributions. Adsorption of Ti(IV) Species. The optimized geometry of a TiCl4 molecule epitactically adsorbed on the (100) face of the MgCl2 cluster 1a is reported in Figure 2. The geometry around the Ti atom is a strongly distorted trigonal bipyramid, with the axis almost perpendicular to the (100) face. The long bond distances Ti-Cl(surface) and Mg-Cl(TiCl4) suggest a weak interaction between the TiCl4 molecule and the MgCl2 cluster. Since there are only four “short” Ti-Cl bonds, an alternative description of the geometry is that of a distorted tetrahedron. In the following, 1a will be the reference state at 0 kcal/mol for the adsorption of TiCl4 molecules. In agreement with previous hypotheses27 and with the findings of Parrinello,40 the most stable geometry of a TiCl4 molecule epitactically adsorbed on the (110) face of the MgCl2 cluster 1b presents the Ti atom octahedrally coordinated. The TiCl4 molecule is considerably closer to the MgCl2 surface relative to 1a, and the better TiCl4/support interaction is confirmed by the high ∆EB (8.3 kcal/mol). We also tried to optimize a TiCl4 molecule

8956

Monaco et al.

Macromolecules, Vol. 33, No. 24, 2000 Table 1. Energy Decomposition of Total Binding Energies, ∆EB

label 1a 1b 1c 2a 2b 3a 3b 4in 4out 5 7-S0 7-S2 9-S0 9-S2 8- S0 8-S4

adsorbed species

unpaired spins

clustera, j

∆Eintra,b kcal/mol

∆Einter,c kcal/mol

∆EBd (EBe), kcal/mol

TiCl4 6-fold TiCl4 6-fold TiCl4 4-fold TiCl4 Ti2Cl8 two TiCl4 TiCl3 TiCl3 TiCl4-TiCl3 TiCl4-TiCl3 Ti2Cl7 Ti2Cl6 Ti2Cl6 Ti2Cl6 Ti2Cl6 Ti4Cl12 Ti4Cl12 Ti4Cl12 Ti4Cl12

0 0 0 0 0 0 1 1 1 1 1 0 2 0 2 0 4 0 4

(100)-1a (110)-1b (110)-2b (110)-1b (100)-1a (110)-2b (100)-1a (110)-1b (110)-4in (110)-4in (100)-1a (100)-1a (100)-1a (100)-1a (100)-1a (100)-1a (100)-1a (100)j (100)j

0 -20.5 -20.4 4.9 -8.4 -21.2 0 -8.7 -19.6,f -9.7g -19.4,f -13.1g -1.7 3.2 6.7 2.3,h -18.6i 3.3,h -21.3i 12.6 12.2 11.5 11.9

0 28.8 30.7 -4.5 9.8 25.6 0 13.6 24.5 23.9 -7.6 -4.6 -7.7 0.7 1.3 -11.5 -10.8 -8.9 -8.9

0 (5.0) 8.3 (13.3) 10.3 (15.3) 0.4 (5.4) 1.4 (6.4) 9.4 (14.4) 0 (24.9) 4.9 (29.8) 5.7 (24.8) 3.5 (22.6) 1.0 (15.9) -1.4 (23.5) -1.0 (23.9) -6.1 (18.2) -5.5 (18.8) 1.1 (26.0) 1.4 (26.3) 2.6 (27.5) 3.0 (27.9)

a The different clusters used to simulate the lateral cuts are labeled with the Miller indices of the considered cut and with the name of first species in which the clusters are present. b ∆EIntra is the energy difference between EIntra per Ti atom of the TimCln species and EIntra of monomeric TiCl4 and/or TiCl3 adsorbed on the (100) cut, which amount to -5.8 and -6.7 kcal/mol, respectively. c ∆EInter is the energy difference between EInter per Ti atom of the TimCln species and EInter of monomeric TiCl4 and/or TiCl3 adsorbed on the (100) cut, which amount to 10.8 and 31.6 kcal/mol, respectively. d ∆EB is the energy difference between EB per Ti atom of the TimCln species and EB of monomeric TiCl4 and/or TiCl3 adsorbed on the (100) cut. e Absolute binding energies per Ti atom with respect to isolated MgCl2 clusters and TiCl4 molecules and/or TiCl3 fragments. f ∆EIntra of the TiCl4 molecule. g ∆EIntra of the TiCl3 fragment. h ∆EIntra per Ti atom of the Ti2Cl6 species. i Unnormalized ∆EIntra of the MgCl2 cluster. j A Mg13C26 cluster, obtained by symmetrical lengthening of the (100)-1a cluster, has been used.

adsorbed on the (110) face with a 5-fold-coordinated Ti atom, similar to those proposed by Pakkanen35 and Parrinello.40 All the performed attempts led to 1c, which essentially corresponds to an undeformed TiCl4 molecule slightly interacting with the MgCl2 cluster (∆EB 0.4 kcal/mol). Since the molecular dynamics simulations of Parrinello were performed without the inclusion of gradient corrections and were checked with gradient corrected static calculations by using the Becke-LeeYang-Parr level of theory,61 we performed a further geometry optimization without including the BeckePerdew gradient corrections.51-53 In this case, the geometry optimization converged to the 5-fold geometry 1d, very similar to that proposed by Parrinello.62 The TiCl4 dimer, Ti2Cl8, coordinated to the (100) face of the MgCl2 cluster, 2a, represents a stable structure of Cs symmetry with two Cl atoms forming bridges between the two octahedrally coordinated Ti atoms, in agreement with models proposed several years ago.27 The Ti2Cl8 fragment is closer to the MgCl2 surface relative to the TiCl4 monomer, and the calculated TiTi and Ti-Cl(surface) distances (3.89 and 2.75 Å, respectively) are in reasonable agreement with the corresponding distances of a model derived from an EXAFS spectrum analysis (3.7 and 3.0 Å, respectively).8 ∆EB amounts to 0.7 kcal/mol, indicating a good interaction between the Ti2Cl8 fragment and the MgCl2 cluster and that the formation of dimeric Ti2Cl8 species is favored relative to the adsorption of isolated TiCl4 molecules. Two such molecules adsorbed on vicinal coordination sites on the (110) face, 2b, also represent a stable situation. In fact, ∆EB per TiCl4 molecule is only 0.9 kcal/mol smaller than ∆EB of an isolated TiCl4 molecule on the same MgCl2 cluster. Adsorption of Ti(III) Species. The lowest energy geometry of the system obtainable after reduction of adsorbed TiCl4 and consisting of a TiCl3 fragment adsorbed on the (100) face of the MgCl2 cluster, 3a, shows the Ti atom tetrahedrically coordinated and short Ti-Cl(surface) and Mg-Cl(TiCl3) distances.41 In the

Figure 3. Optimized geometries of structures 2a and 2b, corresponding to a Ti2Cl8 fragment adsorbed on the (100) face of a Mg9Cl18 cluster and to two TiCl4 molecules adsorbed on vicinal coordination positions on the (110) face of a Mg8Cl16 cluster, respectively. Only the positions of the Ti2Cl8 atoms and of the two TiCl4 molecules have been optimized. Hollow bonds are used to indicate bonds between the Ti2Cl8 fragment or the two TiCl4 molecules and the Mg9Cl18 cluster.

following, 3a will be the reference state at 0 kcal/mol for structures with adsorbed TiCl3 fragments. The structure composed by a TiCl3 fragment adsorbed on the (110) face of the MgCl2 cluster, 3b, has a ∆EB value, 4.9 kcal/mol, which suggests a better interaction between the TiCl3 fragment and the (110) MgCl2 cut.

Macromolecules, Vol. 33, No. 24, 2000

Titanium Chloride Species 8957

Figure 5. Optimized geometries of structures 4in and 4out, corresponding to the adsorption of a TiCl3 fragment and a TiCl4 molecule on the (110) face of a Mg11Cl22 cluster. The main difference between the two structures is the orientation of the dangling chlorine of the TiCl3 fragment. Only the positions of the Ti chloride species have been optimized. Figure 4. Optimized geometries of structure 3a corresponding to a TiCl3 fragment adsorbed on the (100) face of a Mg9Cl18 cluster and of structures 3b and 3c corresponding to C1 and a C2 symmetric TiCl3 fragments adsorbed on the (110) face of a Mg8Cl16 cluster, respectively. Structure 3c corresponds to a transition state. Only the positions of the TiCl3 atoms have been optimized. Hollow bonds are used to indicate bonds between the TiCl3 fragment and the Mg9Cl18 cluster.

Structure 3b is highly asymmetric, and the geometry around the Ti atom is a distorted octahedron with a vacant coordination position. The removal of a chlorine atom from the parent structure 1b results in the loss of a trans influence and thus in the remarkable shortening of a Ti-Cl(surface) distance. Moreover, it can be noted that one of the Cl atoms forming a Ti-Cl-Mg bridge has a position that differs significantly from the epitactical one. This leads to different geometries for the two vicinal magnesium atoms: a square pyramid and a trigonal bipyramid. The asymmetry of 3b implies the existence of another minimum (equivalent by symmetry), in which the only dangling chlorine occupies the vacant coordination position on the Ti atom of 3b. The transition state, 3c, for the switch between these two situations has a C2 symmetry and has a higher energy relative to the minima of only 3.3 kcal/mol, thus suggesting a fast interconversion. When a TiCl4 molecule is coordinated adjacent to 1b, we obtain two minima that are no longer equivalent. In fact, the TiCl4 molecule can be adsorbed on the side of the squarepyramidal magnesium atom or on the side of the trigonal-bipyramidal magnesium atom. In analogy with previous works,28 the two structures will be named inward and outward, 4in and 4out, respectively. Despite the presence of a repulsive interaction between the

TiCl4 molecule and the dangling Cl atom of the TiCl3 fragment, 4in is favored by 4.4 kcal/mol relative to 4out. This preference can be ascribed to the deformation of the TiCl3 fragment, since the deformation energy of the TiCl3 fragment in 4in is 3.4 kcal/mol higher than in 4out, which compares well with the difference in the binding energy reported above. Beside the monomer TiCl3, it is possible to imagine different structures arising from the reduction of adsorbed Ti(IV) species. Starting form the dimer Ti2Cl8, progressive reduction to Ti(III) would lead to Ti2Cl7 and finally to Ti2Cl6. The ∆EB corresponding to the optimized geometry of a Ti2Cl7 species adsorbed on the (100) face, 5, indicates a preference over an isolated TiCl4 molecule and an isolated TiCl3 fragment adsorbed on the same face. The geometry of the pentacoordinated Ti(III) is similar to a trigonal bipyramid; i.e., the dangling chlorine lies in the middle of the two octahedral positions, as expected.30,32-34,63 Upon further reduction, a Ti2Cl6 species can be formed. Although in this paper we are concerned with geometry and stability of adsorbed species, a few remarks on isolated Ti2Cl6 species are useful. Therefore, we considered the structures of free Ti2Cl6 species with zero and two unpaired electron spins, 6-S0 and 6-S2, respectively. The formation energy of 6-S0 and 6-S2 relative to free TiCl3 fragments is -37.0 and -29.3 kcal/ mol, respectively, with 6-S0 favored by 7.7 kcal/mol with respect to 6-S2. The main feature of 6-S0 is the presence of a Ti-Ti bond. In structure 6-S2, instead, the Ti-Ti distance is longer due to the absence of a Ti-Ti bond and to repulsive interactions between the two singly occupied orbitals which are oriented along the Ti-Ti axis.64 The structures composed by Ti2Cl6 coordinated

8958

Monaco et al.

Macromolecules, Vol. 33, No. 24, 2000

Figure 6. Optimized geometry of structure 5, corresponding to the Ti2Cl7 fragment adsorbed on the (100) face of a Mg9Cl18 cluster. Only the positions of the Ti2Cl7 atoms have been optimized. Hollow bonds are used to indicate bonds between the Ti2Cl6 fragment and the Mg9Cl18 cluster.

Figure 8. Optimized geometries of structures 7-S0 and 7-S2, corresponding to the Ti2Cl6 fragment adsorbed on the (100) face of a Mg9Cl18 cluster. Only the positions of the Ti2Cl6 atoms have been optimized. Hollow bonds are used to indicate bonds between the Ti2Cl6 fragment and the Mg9Cl18 cluster.

Figure 7. Optimized geometries of structures 6-S0 and 6-S2, corresponding to the free Ti2Cl6 fragment.

to the (100) face of the MgCl2 cluster, with zero and two unpaired electron spins, 7-S0 and 7-S2, respectively, represent stable structures of Cs symmetry. In structure 7-S0 the Ti-Ti bond distance is stretched by only 0.12 Å with respect to the same distance in the free 6-S0 Ti2Cl6 model species. In structure 7-S2, instead, the TiTi distance is 0.55 Å longer with respect to the same distance in the free 6-S2 Ti2Cl6 species. In both structures the geometry of the Ti is a distorted trigonal bipyramid, as expected.30,32-34,63 It is interesting to note that the deviations of the dangling Cl atoms from the middle positions are opposite in the two structures, being more perpendicular or parallel to the surface. The ∆EB per Ti atom of 7-S0 and 7-S2 amount to -1.4 and -1.0 kcal/mol, respectively, with 7-S2 slightly favored. The preference for the S2 geometry upon coordination is probably due to the absence of a Ti-Ti bond, which allows for a greater flexibility of the Ti2Cl6 fragment, which can easily adopt a geometry apt to be bonded on the surface. For Ti(IV), the dimeric species 2a is the larger cluster that can be grown on the MgCl2 (100) face, whereas after Ti(IV) reduction to Ti(III) also larger clusters could be formed in principle. The optimized geometries of Ti4Cl12 adsorbed on the (100) face of the MgCl2 cluster, with zero and four unpaired electron spins, are denoted 8-S0 and 8-S4, respectively.67 Structure 8-S0 presents an alternance of short-long-short Ti-Ti distances, which are very close to the Ti-Ti distances calculated for the coordinated Ti2Cl6 fragments 7-S0 and 7-S2, respectively. This suggests that 8-S0 can be considered as composed by two Ti2Cl6 fragments presenting a Ti-Ti bond and no unpaired electron spins, connected by two bridging Cl atoms.

Structure 8-S4, instead, presents Ti-Ti distances much more similar to each other, with the central Ti-Ti distance only longer by 0.17 Å. This suggests that 8-S4 can be considered as composed by four TiCl3 fragments connected by bridging Cl atoms. The ∆EB per Ti atom amounts to 1.1 and 1.4 kcal/mol for 8-S0 and 8-S4, respectively. As for the case of the dimeric structures 7-S0 and 7-S2, situations with different number of unpaired spins are of similar energy. Since Ti4Cl12 adsorption in structures 7-S0 and 7-S2 involves Mg atoms that (previous Ti4Cl12 coordination) are not characterized by a coordination number 5, we decided to make test calculations by considering Ti4Cl12 coordination on the bigger Mg13Cl26 cluster and reported in the last two coloumns of Table 1. These calculations indicated that while EB are somewhat dependent from the cluster size, ∆EB are quite more stable. In fact, the adsorbed Ti4Cl12 species with 0 unpaired spins is slightly more stable than the same species with 4 unpaired spins, also with the Mg13Cl26 cluster. The relative binding energies (∆EB) calculated for Ti(IV) (upper part) and Ti(III) (lower part) species on (100) and (110) lateral cuts are collected in Table 1, together with their intramolecular (∆EIntra) and intermolecular (∆EInter) contributions, which generally are considerably larger than the corresponding ∆EB. To check for the effects due to the relaxation of the atoms at the surface of the (100) MgCl2 cut, we partially optimized the geometry of the Mg9Cl18 cluster without, 9, and with coordinated Ti2Cl6 fragments with zero and two unpaired electron spins, 9-S0 and 9-S2, respectively. The optimized atoms are those enclosed in the dashed boxes. In agreement with the results of Parrinello and co-workers,40 we noticed a tendency of the Cl atoms to stretch out of the MgCl2 surface, resulting in a partial screening of the unsaturated Mg atoms. Structures 9-S0 and 9-S2 are quite similar to 7-S0 and

Macromolecules, Vol. 33, No. 24, 2000

Titanium Chloride Species 8959

Figure 9. Optimized geometries of structures 8-S0 and 8-S4, corresponding to the Ti4Cl12 fragment adsorbed on the (100) face of a Mg9Cl18 cluster. Only the positions of the Ti4Cl12 atoms have been optimized. Hollow bonds are used to indicate bonds between the Ti4Cl12 fragment and the Mg9Cl18 cluster.

7-S2, respectively, supporting the reliability of our models. The binding of the Ti2Cl6 fragment to the MgCl2 cluster in 9-S0 and 9-S2 is weaker by roughly 10 kcal/ mol relative to those in 7-S0 and 7-S2. However, also for these more relaxed models the S2 structure is slightly favored, by 0.6 kcal/mol, relative to the S0 structure. This implies that although the absolute EB obtained with fixed MgCl2 clusters are overestimated, the energy differences between different structures seem to be substantially independent from a relaxation of the MgCl2 clusters. Discussion. Although we have calculated only electronic contributions to the binding free energies, it can be reasonably expected that the relative binding energy data reported in Table 1 can be used to gain insights into the species which can be more probably formed during catalyst preparation. First of all, coordination on the (110) face is favored relative to coordination on the (100) face for both TiCl4 molecules and TiCl3 fragments. This is in good qualitative agreement with the higher acidity of the magnesium atoms of the (110) face, as estimated from simple electrostatic calculations.27 In effect, the preferential binding to the (110) face is essentially related to the high attractive intermolecular contributions which more than counterbalance the destabilizing intramolecular contributions (cf. EInter and EIntra values of 1a and 3a with those of 1b and 3b, respectively). As for the (110) lateral cut, Ti(IV) and Ti(III) would be bound preferentially as isolated TiCl4 molecules and TiCl3 fragments, i.e., 1b and 3b, respectively. However, a small energy is required to obtain adjacent monomeric species (like 2b or 4). As for the (100) lateral cut, besides the isolated TiCl4 and TiCl3 species 1a and 3b, dimeric Ti2Cl8 and Ti2Cl6 species (2a, 7-S0, and 7-S2),27 and even longer TinCl3n

Figure 10. Optimized geometries of structures 9, 9-S0, and 9-S2, corresponding to a Mg9Cl18 cluster and to the Ti2Cl6 fragment adsorbed on the (100) face of a Mg9Cl18 cluster, respectively. All the atoms enclosed in the dashed boxes have been optimized. Hollow bonds are used to indicate bonds between the Ti2Cl6 fragment and the Mg9Cl18 cluster.

(e.g., 8-S0 or 8-S2) species17,18 are feasible. In fact, dimeric and tetrameric Ti(III) species are characterized by considerable bonding between the individual TiCl3 fragments (positive EIntra) which more than counterbalance, for this MgCl2 lateral cut, the reduction of attractive interaction with the surface per Ti atom, with respect to the monomeric TiCl3 species (negative ∆EInter in Table 1). As discussed in ref 27, the relative amounts of monomeric and dimeric Ti species on the (100) cut of MgCl2 are not only determined by their relative binding energy values. In fact, each dimeric species (Ti2Cl8 or Ti2Cl6) occupies only three coordination sites, corresponding to three Mg atoms, whereas two monomeric species (2TiCl4 or 2TiCl3) occupy four sites, corresponding to four Mg atoms. Therefore, dimers are favored at high coverage of the surface (which can be achieved by adsorption of Ti species and possibly Lewis bases). This is also true for polynuclear TinCl3n species, which can be formed if atomic rearrangements of Ti(III) species could occur. The prevalence of dimeric and polymeric Ti(III) species would be in line with the ESR silency of a large fraction of Ti(III).13,17,18,20 Stereoselectivity Implications. The Ti(III) monomers and the terminal titanium atoms of the TimCl3m

8960

Monaco et al.

Macromolecules, Vol. 33, No. 24, 2000 Scheme 1

clusters are potential polymerization sites after alkylation of the dangling Ti-Cl bond. Instead, stable Ti(IV) species (excluding the monomer adsorbed on the (100) face) require the abstraction of a chloride anion in order to make place to the incoming monomer. For a given active site, the local symmetry that characterizes the two coordination positions available to the growing chain and to the monomer is extremely relevant to determine the stereospecificity of the insertion reaction. Although alkylation and monomer coordination could lead to geometry changes, some insights into the possible stereoselectivity of the polymerization sites obtained by activation of the various Ti species can be safely obtained by inspection of the models. First of all, we recall that molecular mechanics studies have shown that the stereoselectivity of these catalysts is not due to direct interactions of the chiral catalyst with the monomer, but to interactions between the chiral catalyst (one or two Cl atoms in this case) with the growing chain. These interactions force the chain to assume a chiral orientation, which, in turn, leads to a discrimination between the two prochiral faces of the propene monomer.68,69 In Scheme 1 three different situations (A, B, C) are shown, which differ in the number (0 for A, 1 for B, 2 for C) of conditioning Cl atoms (marked by stars). Ti adsorbed species with structures substantially similar to structure A of Scheme 1 will be substantially aspecific, while those similar to structure B of Scheme 1 will be substantially nonste-

reoselective when the growing chain is outward coordinated, whereas will be substantially stereoselective when the growing chain is inward coordinated. Finally, Ti species similar to structure C of Scheme 1 should be substantially stereospecific, since both coordination positions available to the growing chain will be substantially stereoselective. It is worth noting that structures B and C of Scheme 1 correspond to the Cossee’s C1 model70,71 and to the Allegra’s C2 model,72 first developed for TiCl3-based catalysts. In this framework, a highly isotactic polymer cannot be produced by isolated monomeric species as, for instance, 1a, 1b, 1c, 3a, and 3b as shown by molecular mechanics studies, since they lack conditioning chlorines, as structure A of Scheme 1. Structures 5, 7-S2, and 4, instead, are more similar to the Cosse’s C1 model and resemble structure B of Scheme 1. The stereospecific behavior of structure 7-S2 has been investigated through molecular mechanics studies in the framework of the regular chain migratory mechanism for olefins polymerization.68 These and analogous calculations on the C1 symmetric homogeneous catalytic systems based on group 4 metallocenes73,74 indicated that a regular back-skip of the growing chain to the coordination position previously occupied by the monomer at each step is required to ensure stereospecificity. In fact, by exchanging the relative positions of the growing chain and of the monomer, we switch from a highly stereoselective site (when the growing chain is coordinated close to the MgCl2 surface, inward position) to a substantially nonstereoselective site (when the growing chain is coordinated far from the MgCl2 surface, outward position). The driving force for a possible back-skip of the chain could be an energetic preference of the growing chain to occupy one of the two available coordination positions.73,74 In this respect, it is worth noting that the presence of C1 symmetric sites analogous to the homogeneous ones has been reported for heterogeneous Ziegler-Natta catalysts also.24,26 Among the models discussed in this paper, it can be seen that the terminal Ti of the linear clusters TinCl3n and the monomer adsorbed on the (110) face with both vicinal position saturated are very similar to Allegra’s C2 model, strongly reminds the classical C2 symmetric homogeneous system based on metallocenes of group 4,75 and substantially resembles structure C of Scheme 1. The enantioselective behavior of these systems was rationalized in the framework of the growing chain orientation mechanism.69 The dimeric structure 7-S0 has a C2-like local symmetry around the Ti atoms. This structure has not been investigated through nonbonded interactions, but it is reasonable to think that exchanging the relative positions of the growing chain and of the monomer, we switch from a highly stereoselective site (when the growing chain is inward coordinated) to a poorly stereoselective site (when the growing chain is outward coordinated). In fact, the dangling Cl atom in the almost axial position could induce a chiral orientation of the growing chain, but oscillations of this Cl atom in the relatively open environment could reduce the effectiveness of the steric control. Finally, we want to remark that the stereoselective behavior can be easily modified by the presence of other species (internal and external donors in particular), which can either directly change the steric environment of the active species or modify the surface distri-

Macromolecules, Vol. 33, No. 24, 2000

bution of adsorbed species with different stereospecificity. The latter mechanism favors both dimers27 and linear clusters on the (100) face and lowers the number of monomers with vacant vicinal positions on the (110) face. Moreover, donors could saturate one or both of the vacant positions of a TiCl3 fragment adsorbed on the (110) face and thus can influence directly the stereospecificity. 4. Conclusions In this paper we attempted a systematic DFT investigation of the interactions possibly occurring between TiCl4 molecules and TiCl3 fragments with both (100) and (110) lateral cuts of MgCl2 clusters. Models with a fixed geometry of the MgCl2 support have been considered. However, we think that the energy differences found for different structures would be substantially confirmed by calculations with a more relaxed support. In effect, in the case of the Ti2Cl6 dimer, calculations on structures 9, 9-S0, and 9-S2, in which a partial relaxation of the support was allowed, substantially confirmed the results obtained with the more rigid 7-S0 and 7-S2 structures. Thus, the following conclusions can be drawn: (i) TiCl4 molecules and TiCl3 fragments are preferentially adsorbed on the (110) rather than on the (100) lateral cut of MgCl2. (ii) With respect to monomeric species, dimeric Ti2Cl8 on the (100) lateral cut is energetically favored, while Ti2Cl6 species are unfavored. These species are very similar to the models proposed by Cossee,70,71 suggested by EXAFS studies,8 and calculated to be possible isospecific site by molecular mechanics calculations.27 (iii) Besides monomeric and dimeric species, polynuclear Ti(III) species on the (100) lateral cut of MgCl2 are possible and appear to be energetically favored, independent of the number of unpaired spins. This result supports the formation of polynuclear Ti chloride species first hypothesized by Brant and Speca17,18 and could contribute to explain the ESR silency of a large fraction of Ti(III).13,17,18,20 The characterization of the various adsorbed Ti species is extremely relevant in terms of their potential stereospecific behavior in propene polymerization. It has to be noted that the several structures of similar energy we have found are consistent with a family of polymerization sites with different local symmetries (from C1 to C2-like to almost C2), steric environment, and hence different stereospecific properties (from substantially aspecific sites to isospecific ones). As far as the (100) face is concerned, dimeric Ti2Cl6 sites on the (100) cut are very similar to the C1 symmetric sites proposed by Cossee,70,71 while polynuclear TinCl3n sites with n > 2 are extremely similar to the C2 symmetric sites proposed several years ago by Allegra for TiCl3-based catalytic systems,72 and present a strict analogy with the well-established models for isospecific polymerization with catalytic systems based on C2 symmetric metallocenes. On the (110) face there are many possibilities. In fact, although an isolated TiCl3 molecule adsorbed on the (110) face can be thought to produce a poorly stereoregular polymer with a prevalence of syndiotactic diads,76 C1 and C2 symmetric species, which can produce isotactoid and isotactic polymers, are obtainable by occupation of one and two vicinal positions. An easy switch between these three polymerization modes has

Titanium Chloride Species 8961

been recently proposed on the basis of a high-resolution 13C NMR microstructural analysis.26 This study is far from being comprehensive. In fact, the Ti species we have calculated can only represent a basis for the description of the real effective catalyst, since the relative distribution, interconversion, and stereospecificity of the various Ti species are sensitively modified by the presence of alkylaluminum compounds and, particularly, by the added Lewis bases. A study relative to the interaction between Lewis bases (1,3diethers in particular) and the MgCl2 support is reported elsewhere,43 and further work is in progress. In any case, we believe that the present study can represent a considerable improvement of our understanding of the interaction between the various Ti species and the MgCl2 support, the two fundamental components of supported heterogeneous Ziegler-Natta catalysts. Acknowledgment. We thank Prof. V. Busico of the University of Naples and scientists of Montell Polyolefins for useful discussions. This work has been supported by the Ministero della Ricerca Scientifica e Tecnologica of Italy within PRIN 1998 and 2000 fundings, and by Montell Polyolefins. References and Notes (1) Moore, E. P. J. Polypropylene Handbook: Polymerization, Characterization, Properties, Applications; Hanser Publishers: Munich, 1996. (2) Natta, G. Nobel Lectures in Chemistry, 1963-1970; Elsevier: Amsterdam, 1972; p 27. (3) Ziegler, K. Nobel Lectures in Chemistry, 1963-1970; Elsevier: Amsterdam, 1972; p 6. (4) Kissin, Y. V. Isospecific Polymerization of Olefins; SpringerVerlag: New York, 1985. (5) Giannini, U.; Giunchi, G.; Albizzati, E.; Barbe`, P. C. Frontiers in Polymerization Catalysis 1987, NATO Advanced Research Workshop, Feb 1-6. (6) Auriemma, F.; Talarico, G.; Corradini, P. In Progress and Development of Catalytic Olefin Polymerization; Sano, T., Uozumi, T., Nakatani, H., Terano, M., Eds.; Technology and Education Publisher: Tokyo, Japan, 2000; p 7. (7) Jones, P. J. V.; Oldman, R. J. In Etero EXAFS MgCl2; Kaminsky, W., Sinn, H., Eds.; Springer-Verlag: Berlin, 1988; p 337. (8) Popatov, A. G.; Kriventsov, V. V.; Kochubey, D. I.; Bukatov, G. D.; Zakharov, V. A. Macromol. Chem. Phys. 1997, 198, 3477. (9) Magni, E.; Somorjai, G. A. J. Phys. Chem. B 1998, 102, 8788. (10) Hasebe, K.; Mori, H.; Terano, M. J. Mol. Catal. 1997, 124, L1. (11) Mori, H.; Hasebe, K.; Terano, M. Macromol. Chem. Phys. 1998, 199, 2709. (12) Mori, H.; Sawada, M.; Higuchi, T.; Hasebe, K.; Otsuka, N.; Terano, M. Macromol. Rapid Commun. 1999, 20, 245. (13) Chien, J. C. W.; Wu, J. C. J. Polym. Sci., Part A: Polym. Chem. 1982, 20, 2461. (14) Chien, J. C. W.; Hu, Y. J. Polym. Sci., Part A: Polym. Chem. 1989, 27, 897. (15) Sergeev, S. A.; Polubayarov, V. A.; Zakharov, V. A.; Anufrienko, U. F.; Bukatov, G. D. Makromol. Chem. 1985, 186, 243. (16) Weber, S.; Chien, J. C. W.; Hu, Y. J. Polym. Sci., Part A: Polym. Chem. 1989, 27, 1499. (17) Brant, P.; Speca, A. N. Macromolecules 1987, 20, 2740. (18) Brant, P.; Speca, A. N.; Johnston, D. C. J. Catal. 1988, 113, 250. (19) Kashiwa, N.; Yoshitake, J. Makromol. Chem. 1984, 185, 1133. (20) Fuhrmann, H.; Herrmann, W. Macromol. Chem. Phys. 1994, 195, 3509. (21) Busico, V.; Corradini, P.; De Martino, L. Makromol. Chem., Rapid Commun. 1990, 11, 49. (22) Paukkeri, R. V. T.; Lehtinen, A. Polymer 1993, 34, 2488. (23) Busico, V.; Cipullo, R.; Corradini, P.; Landriani, L.; Vacatello, M.; Segre, A. L. Macromolecules 1995, 28, 1887.

8962

Monaco et al.

(24) Busico, V.; Cipullo, R.; Talarico, G.; Segre, A. L.; Chadwick, J. C. Macromolecules 1997, 30, 4787. (25) Randall, J. C. Macromolecules 1997, 30, 803. (26) Busico, V.; Cipullo, R.; Monaco, G.; Talarico, G.; Vacatello, M.; Chadwick, J. C.; Segre, A. L.; Sudmeijer, O. Macromolecules 1999, 32, 4173. (27) Corradini, P.; Barone, V.; Fusco, R.; Guerra, G. Gazz. Chim. Ital. 1983, 113, 601. (28) Corradini, P.; Busico, V.; Guerra, G. Comprehensive Polymer Science; Allen, G., Bevington, J. C., Eds.; Pergamon Press: Oxford, 1989; Vol. 4, pp 29-50. (29) Zambelli, A. L. P.; Bajo, G.; Bovey, F. A. Macromolecules 1975, 8, 687. (30) Cavallo, L.; Guerra, G.; Corradini, P. J. Am. Chem. Soc. 1998, 120, 2428. (31) Fujimoto, H.; Yamasaki, T.; Mizutani, H.; Koga, N. J. Am. Chem. Soc. 1985, 107, 6157. (32) Jensen, V. R.; Børve, K. J.; Ystenes, M. J. Am. Chem. Soc. 1995, 117, 4109. (33) Sakai, S. J. Phys. Chem. 1994, 98, 12053. (34) Sakai, S. Int. J. Quantum Chem. 1997, 65, 739. (35) Puhakka, E.; Pakkanen, T. T.; Pakkanen, T. A. Surf. Sci. 1995, 334, 289. (36) Puhakka, E.; Pakkanen, T. T.; Pakkanen, T. A.; Iiskola, E. J. Organomet. Chem. 1996, 511, 19. (37) Puhakka, E.; Pakkanen, T. T.; Pakkanen, T. A. J. Mol. Catal. 1996, 120, 143. (38) Puhakka, E.; Pakkanen, T. T.; Pakkanen, T. A. J. Chem. Phys. 1997, 101, 6063. (39) Puhakka, E.; Pakkanen, T. T.; Pakkanen, T. A. J. Mol. Catal. 1997, 123, 171. (40) Boero, M.; Parrinello, M.; Terakura, K. J. Am. Chem. Soc. 1998, 120, 2746. (41) TiCl3 does not exist as an isolated molecule in the conditions used in Ziegler-Natta catalysis. However, Ti(III) species adsorbed on the MgCl2 support are clearly feasible after reduction of adsorbed TiCl4 molecules. (42) Barbe´, C.; Cecchin, G.; Noristi, L. Adv. Polym. Sci. 1987, 81, 1. (43) Toto, M.; Morini, G.; Guerra, G.; Corradini, P.; Cavallo, L. Macromolecules 2000, 33, 1134. (44) Partin, D. E.; O’Keeffe, M. J. Solid State Chem. 1991, 95, 176. (45) ADF 2.3.0, User Guide Vrije Universiteit Amsterdam, Amsterdam, The Netherlands, 1996. (46) Baerends, E. J.; Ellis, D. E.; Ros, P. Chem. Phys. 1973, 2, 41. (47) te Velde, B.; Baerends, E. J. J. Comput. Phys. 1992, 99, 84. (48) Versluis, L.; Ziegler, T. J. Chem. Phys. 1998, 88, 322. (49) Snijders, J. G.; Vernooijs, P.; Baerends, E. J. At. Nucl. Data Tables 1981, 26, 483. (50) Vosko, S. H.; Wilk, L.; Nusair, M. Can. J. Phys. 1980, 58, 1200. (51) Becke, A. Phys. Rev. A 1988, 38, 3098. (52) Perdew, J. P. Phys. Rev. B 1986, 33, 8822. (53) Perdew, J. P. Phys. Rev. B 1986, 34, 7406. (54) Margl, P. M.; Ziegler, T. Organometallics 1996, 15, 5519. (55) Margl, P. M.; Ziegler, T. J. Am. Chem. Soc. 1996, 118, 7337. (56) Willis, B. G.; Jensen, K. F. J. Phys. Chem. A 1998, 102, 2613. (57) Jensen, V. R.; Børve, K. J. J. Comput. Chem. 1998, 19, 947.

Macromolecules, Vol. 33, No. 24, 2000 (58) Nayak, S. K.; Jena, P. J. Am. Chem. Soc. 1999, 121, 644. (59) Jacobsen, H.; Ziegler, T. J. Am. Chem. Soc. 1994, 116, 3667. (60) Ritter, E.; Kvisle, S.; Nirisen, Ø.; Ystenes, M.; Øye, H. A. In Transition Metal Catalyzed Polymerization: Ziegler-Natta and Metathesis Polymerization; Quirk, R. P., Ed.; Cambridge University Press: Cambridge, 1988; p 292. (61) Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B 1988, 37, 785. (62) To further investigate these somewhat different conclusions, we optimized a 5-fold coordinated structure by fixing the TiCl(surface) distance to 3.20 Å (which is the Ti-Cl(surface) distance found for a TiCl4 molecule adsorbed on the (100) MgCl2 face). The optimized 5-fold geometry is higher in energy by only 2.1 kcal/mol relative to 3b, indicating a substantially flat potential energy surface. It is reasonable that different levels of theory can give quite different results in such cases. (63) Novaro, O.; Blaisten-Barojas, E.; Clementi, E.; Giunchi, E.; Ruiz-Vizcaya, M. E. J. Chem. Phys. 1978, 68, 2337. (64) The very short Ti-Ti bond distance we calculated for 6-S0 is in a good agreement with the Ti-Ti distances found by Cotton and Wojtczak65 in the diamagnetic complexes Ti2Cl2[H5C6NC(H)NC6H5]4 and Ti2Cl2[H9C12NC(H)NC12H9]4 (TiTi ) 2.899 and 2.916 Å, respectively). Differently, the Ti-Ti distance we calculated for 6-S2 is too short when compared to the Ti-Ti distances found in the paramagnetic complexes Ti2Cl6[1,2-bis(diisopropylphospino)ethane]266 and Ti2Cl6[P(CH3)3]465 (Ti-Ti ) 3.438 and 3.704 Å). To check the performances of the methodology we used, we optimized the structure of the Ti2Cl6(PH3)4 complex with two unpaired electron spins, and we calculated the Ti-Ti distance to be equal to 3.77 Å, in very good agreement with the experimental value for Ti2Cl6[P(CH3)3]4. It is worthy to note that a geometry optimization of the Ti2Cl6(PH3)4 complex without gradient corrections gives a Ti-Ti distance equal to 3.18 Å. (65) Cotton, F. A.; Wojtczak, W. A. Gazz. Chim. Ital. 1993, 123, 499. (66) Hermes, A. R.; Girolami, G. S. Inorg. Chem. 1990, 29, 313. (67) Clearly, Ti4Cl12 species with two unpaired electron spins could have been considered as well. However, they are not needed for the present general discussion. (68) Corradini, P.; Barone, V.; Fusco, R.; Guerra, G. Eur. Polym. J. 1979, 15, 1133. (69) Corradini, P.; Barone, V.; Fusco, R.; Guerra, G. J. Catal. 1982, 77, 32. (70) Arlman, E. J.; Cossee, P. J. Catal. 1964, 3, 99. (71) Cossee, P. J. Catal. 1964, 3, 80. (72) Allegra, G. Makromol. Chem. 1971, 145, 235. (73) Guerra, G.; Cavallo, L.; Moscardi, G.; Vacatello, M.; Corradini, P. Macromolecules 1996, 26, 4834. (74) Rieger, B.; Jany, G.; Fawzi, R.; Steimann, M. Organometallics 1994, 13, 647. (75) Brintzinger, H. H.; Fischer, D.; Mu¨lhaupt, R.; Rieger, B.; Waymouth, R. M. Angew. Chem., Int. Ed. Engl. 1995, 34, 1143. (76) Corradini, P.; Guerra, G.; Barone, V. Eur. Polym. J. 1984, 20, 1177; see particularly Figures 6 and 9.

MA000988H