A R T I C L E S
Joubert et al.
has been illustrated by the understanding at a molecular level
of the difference of reactivity of [ZrR4] and [W(tCR′)R3]
complexes with alumina.
on alumina are active polymerization catalysts, while their
corresponding silica supported analogues, neutral, are inactive.7-9
The combination of theory and physical chemistry characteriza-
tions hence appears as a promising tool for the understanding,
at a molecular scale, of well-defined surface catalytic species,
enabling the rational design of highly efficient catalysts.
The first step of the reaction of ZrR4 with an alumina partially
dehydroxylated at 500 °C is an electrophilic cleavage of a
ZrsC bond preferentially on HOsAlIV (∆rE ) -201 kJ/mol)
through σ-bond metathesis with a low activation barrier (Ea )
37 kJ/mol). This reaction gives [(AlsO)ZrR3] along with one
alkane molecule. In the case of W(tCR′)R3, the first step is
also an electrophilic cleavage of the WsR bond preferentially
on HOsAlIV (∆rE ) -150 kJ/mol) to form [(AlsO)W(tCR′)-
R2], but this reaction is associated with a much higher activation
barrier (Ea ) 81 kJ/mol). The lower reactivity of the W complex
is due to the carbyne ligand, which hinders the necessary
deformation of the complex for the approach of the AlsOH
toward the WsC bond. Although preferred on HOsAlIV, the
grafting reaction is also possible on other accessible OH groups
such as HOs(Alv,Alv), especially for the case of the Zr complex
where the reaction barrier is low.
Computational Methods and Systems
The calculations were performed on the framework of density
functional theory (DFT) using a periodic description of the system as
implemented in the VASP code.39,40 The generalized gradient ap-
proximation was used in the formulation of Perdew and Wang PW91.41
Atomic cores were described with the projected augmented wave
method (PAW) which is equivalent to an all electron frozen core
approach.42,43 The one electron wave functions are developed on a basis
set of plane waves. With the selected PAW potentials, a cutoff energy
of 275 eV is adequate and yields a converged total energy.
Brillouin zone integration was converged with a 331 k-point mesh
generated by the Monkhorst-Pack algorithm.44 Vibrational frequencies
were calculated in the harmonic approximation by a numerical
evaluation of the Hessian matrix. An anharmonicity term of 80 cm-1
,
For Zr, the favored second step is the reaction of [(AlsO)-
ZrR3] with an adjacent hydroxyl to give a neutral complex,
[(AlsO)2ZrR2], and the evolution of a second molecule of alkane.
This second step is still exoenergetic (-154 kJ/mol), its barrier
is lower than that of the first step and is not affected by the
size of the ligands (R ) Me or CH2tBu). This surface complex
further reacts with the alumina surface in an almost barrierless
process through the transfer of one of its alkyl ligands onto an
adjacent Als Lewis center, giving a cationic surface complex
{[(AlsO)2ZrR]+[(Als)R-]}. The alkyl transfer to give [(AlIV)R-]
calculated previously on hydroxyl groups on boehmite, has been applied
a posteriori.45
Chemical shifts have been evaluated with the GIAO method46
implemented in the Gaussian03 code47 at a DFT/B3LYP level.48-51 The
IGLO-II basis set52 has been used for carbon and hydrogen. For other
atoms, the Hay and Wadt effective core potential53-55 has been used
with the adapted LANL2DZ basis set. The chemical shift calculations
have been performed on selected clusters. The interested reader will
find more details in the Supporting Information in which a validation
of computational parameters on well-defined molecules is also proposed.
and [(AlVI)R-] has been confirmed by the combined use of 13
C
Experimental Section
CP-MAS solid-state NMR and chemical shift calculations.
Moreover, this complete reaction pathway involves the forma-
tion of two alkane molecules per grafted Zr as observed
experimentally.
General Procedure. All experiments were conducted under
strict inert atmosphere or vacuum conditions using standard Schlenk
techniques. Solvents were purified and dried according to standard
procedures. tBuCH2MgCl was prepared from tBuCH2Cl (98%, Aldrich)
and Mg turnings (99%, Lancaster). [(1-13C 33%) tBu13CH2MgCl]
was prepared according to the literature procedure (It was pre-
pared from a mixture of 1:2 [(1-13C 99%) tBu13CH2MgCl] and un-
labeled tBuCH2MgCl. [(1-13C 99%) tBu13CH2Cl] was prepared from
(1-13C 99%) H13CONMe2 (Cambridge Isotope Laboratories) and
unlabeled tBuLi followed by reduction to the alcohol, chloration with
a Vilsmeier reagent, and finally conversion to the Grignard with Mg).35
[Zr(CH2tBu)4] (1) and [(13C 33%) Zr(CH2tBu)4] (1*) were prepared
by alkylation of ZrCl4 (Clariant, used as received) with tBuCH2MgCl
and [(1-13C 33%) tBu13CH2MgCl], respectively, and were sublimed prior
The situation is completely different for the W carbyne
complex. After the first step, the [(AlsO)W(tCR′)R2] species
is found to be nonreactive. Further reaction with surface
hydroxyls involves a high reaction barrier (∼100 kJ/mol) for
the real CH2tBu ligand. This is in agreement with the evolution
of a single alkane molecule in the reactive adsorption process.
Moreover, transfer of an alkyl ligand to an adjacent Al center
does not occur either as it involves the formation of unstable
cationic species (∆rE ) +50 kJ/mol). This is also confirmed
by the absence of typical 13C NMR signals from surface alkyl
aluminum species. Nonetheless, this neutral complex shows
weak lateral interactions with the remaining adjacent surface
hydroxyls as evidenced by a combination of IR spectroscopy
and frequency calculations. It is noteworthy that the alumina
supported W species, [(AlsO)W(tCR′)R2], is in principle well-
defined but that its 13C NMR spectrum is complicated simply
because this well-defined species interacts through its ligands
(alkyl and carbyne) with different remaining hydroxyls present
on the surface of alumina, thus providing different conformers,
whose calculated NMR spectra are consistent with experimental
data. This reflects the complexity of γ-alumina surfaces. In
contrast, the zirconium complex, which has reacted with its
surrounding hydroxyl groups and an adjacent aluminum, gives
a fairly well-defined complex as evidenced by a rather simple
NMR spectrum. It is a true cationic alkyl zirconium species,
and this explains quite well why zirconium alkyls supported
(39) Kresse, G.; Furthmuller, J. Comput. Mater. Sci. 1996, 6, 15.
(40) Kresse, G.; Furthmueller, J. Phys. ReV. B: Condens. Matter 1996, 54,
11169.
(41) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M.
R.; Singh, D. J.; Fiolhais, C. Phys. ReV. B: Condens. Matter 1992, 46,
6671.
(42) Bloechl, P. E.; Forst, C. J.; Schimpl, J. Bull. Mater. Sci. 2003, 26, 33.
(43) Bloechl, P. E. Phys. ReV. B: Condens. Matter 1994, 50, 17953.
(44) Monkhorst, H. J.; Pack, J. D. Phys. ReV. B: Condens. Matter 1976, 13,
5188.
(45) Raybaud, P.; Digne, M.; Iftimie, R.; Wellens, W.; Euzen, P.; Toulhoat, H.
J. Catal. 2001, 201, 236.
(46) Ditchfield, R. J. Chem. Phys. 1972, 56, 5688.
(47) Pople, J. A. et al. Gaussian 98; Gaussian Inc.: Pittsburgh, PA, 1998.
(48) Becke, A. D. Phys. ReV. A: At., Mol., Opt. Phys. 1988, 38, 3098.
(49) Dickson, R. M.; Becke, A. D. J. Chem. Phys. 1993, 99, 3898.
(50) Becke, A. D. J. Chem. Phys. 1993, 98, 1372.
(51) Lee, C.; Yang, W.; Parr, R. G. Phys. ReV. B: Condens. Matter 1988, 37,
785.
(52) Kutzelnigg, W.; Fleischer, U.; Schindler, M. In NMR basic principles and
progress; Springer: Berlin, 1990; Vol. 23, p 165.
(53) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 82, 270.
(54) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 82, 299.
(55) Wadt, W. R.; Hay, P. J. J. Chem. Phys. 1985, 82, 284.
9
9168 J. AM. CHEM. SOC. VOL. 128, NO. 28, 2006