Table of contents
1. Aims
2. Introduction
a. to the semi empirical method MNDO
b. to the empirical method UFF
c. to Zeolitic Catalysis
3. Experimental
a. Optimization results of [Abbildung in dieser Leseprobe nicht enthalten]
b. HOMO comparisons
c. Structure prediction of [Abbildung in dieser Leseprobe nicht enthalten]
d. Optimization Results
e. Reaction modeling
f. Tables
4. Appendix
a. Summary, Conclusions & further work
b. Acknowledgements
c. Computational Details
d. References
e. Definition of bond and angle names
Aims
The aims of this research project were:
a) Optimization of the (known) aluminosilsesquioxane structure of
[Abbildung in dieser Leseprobe nicht enthalten] using empirical (UFF, Dreiding) and semiempirical (MNDO) methods and comparison and evaluation of those methods
b) Derivation of the as-small-as-possible but as-big-as-necessary model for modelling of the reaction of this silsesquioxane with an olefin, based on structure and HOMO properties
c) Prediction of the hitherto unknown structure of another aluminosilsesquioxane with the sum formula [Abbildung in dieser Leseprobe nicht enthalten]
d) Reaction modeling of both silsesquioxanes with an olefin
MNDO [Dewar and Thiel 1977]
The variational principle says that the energy for the exact wave function is lower than the energy for the ground state of any other antisymmetric normalized wave function. So the problem of approximating the wave function is a problem of finding coefficients ciµ in [Abbildung in dieser Leseprobe nicht enthalten] so that the energy is at a minimum. (φi are the molecular orbitals and [Abbildung in dieser Leseprobe nicht enthalten] the basis functions).
This leads to a description derived independently and at the same time by Roothan [Roothan 1951] and Hall [Hall 1951], written in matrix form: FC = SCE. This is a nonlinear equation as the coefficients cµi appear as well in the Fock matrix as in the molecular orbitals and has therefore to be solved iteratively, named the SCF or self-consistent field method. F is the Fock matrix and stands for the average effects of the field of all electrons on each orbitals, S is the overlap matrix and E the energy matrix, where each element [Abbildung in dieser Leseprobe nicht enthalten] is the one-electron energy of a molecular orbital [Abbildung in dieser Leseprobe nicht enthalten].
In Hartree-Fock (HF) ab-inito calculations all elements of the Fock matrix are calculated using
[Abbildung in dieser Leseprobe nicht enthalten], irrespective of location of the basis functions [Abbildung in dieser Leseprobe nicht enthalten] and [Abbildung in dieser Leseprobe nicht enthalten]. P is the density matrix and represents the energy of a single electron in the field of bare nuclei and [Abbildung in dieser Leseprobe nicht enthalten] are the two-electron repulsion integrals, each electron sees all other electrons as an average distribution.
Most of the computation time in an HF-calculation is required for calculating integrals, so it is an obvious approach to reduce computational effort by neglecting or approximating parts of the integrals. Therefore semi-empirical methods consider (explicitly) only the valence electrons of the system, the core electrons are subsumed into the nuclear core. In addition, the orthogonality of the employed basis functions simplifies some of the integrals. Common to all semi-empirical methods in that the overlap matrix S is set equal to the identity matrix I. The effect is, that the elements that correspond to an overlap between two atomic orbitals on different atoms are set to zero and the Roothaan-Hall equation FC=SCE becomes FC=CE, the standard matrix form.
What is specific to MNDO? MNDO stands for “Modified Neglect of Diatomic Overlap” and is based on NDDO, the “Neglect of Diatomic Differential Overlap” [Pople et. al. 1965]. The basis of NDDO in turn is a method called “CNDO”. CNDO [Pople et. al. 1965], the “Complete Neglect of Differential Overlap”, was the first method to implement the zero-differential overlap approximation, and neglected all differential overlap, regardless whether the orbitals are centered on the same or on different atoms.
NDDO, as the next level of approximation, only neglects differential overlap between atomic orbitals on different atoms. Thus all two-electron, two-center integrals of the form (µν|λσ) where µ and ν are on the same atom and λ and σ are also on the same atom are retained.
MNDO differs from NDDO in respect of the definition of the Fock matrix elements. I will describe these differences now more detailed.
The major improvement are the new terms [Abbildung in dieser Leseprobe nicht enthalten] and [Abbildung in dieser Leseprobe nicht enthalten] in H core,[Abbildung in dieser Leseprobe nicht enthalten] , which represent two-centre, one-electron attractions between an electronic distribution [Abbildung in dieser Leseprobe nicht enthalten] or [Abbildung in dieser Leseprobe nicht enthalten] respectively on atom A and the core of atom B. These are defined as [Abbildung in dieser Leseprobe nicht enthalten] and [Abbildung in dieser Leseprobe nicht enthalten]. Another difference are the core-core repulsion terms with OH and NH being treated separately and a new form for the two-centre, one electron resonance integrals which depend on the overlap Sµν and parameters βµ and βν in H core,[Abbildung in dieser Leseprobe nicht enthalten]
An advantage of MNDO over an earlier semiempirical method, MINDO/3 [Bingham et al. 1975], is the use solely of monoatomic parameters. Almost any element is defined in MNDO. In the original method only an s, p basis set was employed, but versions that are more recent explicitly include d orbitals for heavier elements, although they are of dubious quality.
A problem in MNDO in the overestimation of repulsion between atoms separated by a distance equal to the sum of their Van-der-Waals radii, so hydrogen bonds are not sufficiently modelled. Another problem in MNDO are conjugated systems and sterically crowded molecules, where MNDO energies are too positive and fourmembered rings, where the energy is too negative.
Force Fields
When applied to large molecules (i.e. with more than a few hundred atoms) even semiempirical methods are too time consuming to be feasible in a sensible time range. Thus, although in principle accessible by those methods, a classical model is introduced, where the forces between atoms are only based on simple functions describing a “ball-and-spring” behaviour, e.g. Hooke’s law. This approach reduces the computational costs by about two orders of magnitude (a factor of 100), so even biological systems are accessible, although they consist of thousands of atoms. This area can still be considered as the domain of the non-quantum mechanical approach of force fields.
A “reference” value for bond lengths, angels and torsion angels is defined (derived from experiments) and an energy function of the general form
Abbildung in dieser Leseprobe nicht enthalten
can be conceived, where we find (in this order) terms how the energy varies with the bond length, bond angles, torsion angles and interaction between non-bonded atoms. The latter in simple force fields often consists of a Coulomb potential for electrostatic and a Lennard-Jones (6,12) potential for van-der-Waals interactions. In more sophisticated force fields also cross terms, as angle-angle terms etc. are included.
Some force fields have been derived to lead satisfactory results on special groups of molecules, e.g. MM2 (originally) for hydrocarbons [Allinger 1977], SHAPES [Allured et al. 1991] especially for transition metal complexes or AMBER [Weiner et al. 1984] for proteins and nucleic acids.
Here we use the UFF “Universal Force Field” [Rappe et al. 1992] and the Dreiding Force Field [Mayo et al. 1990].
The UFF force field is designed to model the entire periodic table and features a novel treatment of angle bending. The commonly used force fields use an harmonic potential for angle bending, but this approach is misleading when angles approach 180°. Therefore, UFF uses (as the SHAPES force field) a cosine Fourier series [Abbildung in dieser Leseprobe nicht enthalten] for each angle ABC.
The special reason for using UFF in this case is that the Si-O-Si bond energy changes only very little over a wide range (~120 - 180 °) [Nicholas et al. 1991] and the Fourier expansion in UFF is designed to deal with this particularity.
The Dreiding force field in contrary uses a harmonic cosine form for angle bending and is therefore expected to yield less satisfactory results. In a geometry optimisation using force fields the energy has to be minimized. A variety of difficulties arise from this question as for example most optimisation algorithms don’t allow the energy to rise and so often yield the nearest (and local) instead of the global minimum, that is the aim of the optimization. Minimization algorithms can be divided into groups that use the derivative of the energy (as e.g. the Newton-Raphson algorithm does) and those who don’t (e.g. the simplex algorithm).
All of the algorithms have advantages and disadvantages, dependent on the actual form of the potential energy surface.
An introduction can be found in [Leach 1996].
Introduction to Zeolitic Catalysis
Acidic aluminum sites in zeolites play an important role in catalysis as for example in isomerization and cracking of hydrocarbons [Meisel et. al. 1976].
The unique environment of the Brønsted acidic sites controls the overall behavior in Zeolites significantly; one of the consequences is the resulting stereochemistry of the catalysts. In addition, transport mechanisms play an important role since for sterical reasons not every reactant can be transported to the catalytic center and not every product out of the cluster. What makes those studies difficult is the heterogeneity of these systems.
A wide range of reactions between Zeolites and other molecules has been explored with theoretical methods, for example with CO [Bates and Dwyer 1993] or CH3CN [Meijer et al. 1996]. Since I performed reactions with C-C bonds, I will restrain my description of the catalytic activity of Zeolites to this particular field.
Zeolites have a structure in common where the H next to Al acts as a Brønsted acid and the adjacent oxygens, because of their lone pairs, as Lewis bases.
Abbildung in dieser Leseprobe nicht enthalten
This particular structure determines the catalytic behaviour to a significant extent. A good review of theoretical aspects of zeolite catalysis can be found in [van Santen and Kramer 1995].
I shall confine my description further to the activation of Olefins; some of the earliest calculations in this area were performed by [Senchenya and Kazansky 1987] and [Pelmenschikov et al. 1987].
The reaction with ethene revealed the following mechanism [Senchenya and Kazansky 1991], based on calculations for a 1T-cluster with the 3-21 G basis set:
Abbildung in dieser Leseprobe nicht enthalten
The ð-bonded complex is formed (a), stabilized by approximately 30 kJ/mol. As an intermediate, we observe (b), reacting to the ó-bonded complex (c). We tried to follow this reaction path on a 3T model and on both of our silsesquioxanes.
Experimental
Our first goal is the comparison and evaluation of the UFF, Dreiding and MNDO method with respect to the structure optimization of [Abbildung in dieser Leseprobe nicht enthalten]
The cation was removed in all models because it forms only a separated ion pair with the anion and is expected to have a negligible effect [Edelmann et al 1991]. Therefore, the anion structure obtained from x-ray data was optimized employing all of the mentioned methods. Afterwards the structure of the optimized models was compared with the experimental data.
The second aim is to derive the as-small-as-possible but as-big-as-necessary model in respect to both geometry around the central Aluminum and shape and expansion of the HOMO, as this frontier orbital is thought to be very important for catalytic activity. Thus, we successively removed the cyclohexanes and substituted them by hydrogen to observe how those properties change when we simplify the molecule.
The X-ray structure we compare our computed results with can be found in [Edelmann et.al.1999].
The full anion of the molecule [Abbildung in dieser Leseprobe nicht enthalten] is as shown below.
Abbildung in dieser Leseprobe nicht enthalten
This anion was reduced by removing the four cyclohexanes furthest away from the Aluminum, yielding the so-called layer 3. The SiMe3 groups were also exchanged by OH yielding a model compound with the formula [Abbildung in dieser Leseprobe nicht enthalten]- and a structure shown below.
[Abbildung in dieser Leseprobe nicht enthalten
From this model we reduced the size by removing the six cyclohexanes separated by one Silicon from the Al in the center, the result is the layer 2 anion with the formula [Abbildung in dieser Leseprobe nicht enthalten]
Abbildung in dieser Leseprobe nicht enthalten
Finally we cut the last cyclohexanes from the molecule, yielding layer 1, [Abbildung in dieser Leseprobe nicht enthalten]
Abbildung in dieser Leseprobe nicht enthalten
All these molecules were optimized using the UFF, Dreiding and MNDO method. For all these layers the HOMO values were computed and isosurfaces were plotted to follow the change in chemical behavior.
Results of optimizations
The bond lengths and angles of the original anion and those of the optimized structures are given in the table 1 (p. 26).
To find a way for measuring the rate of systematic errors we computed an arithmetic mean of the deviations of the bond length and an arithmetic mean of the absolute value of the derivations of the bond lengths. If we then divide the former by the latter, we have a possibility to measure the rate of consistency of the systematic error.
As we can see, the UFF method always overestimates bond lengths by 0.106 to 0.142 Ångstroms with an average of 0.116 Angstroms within Al-O bond lengths. With respect to O-Si bonds the bond lengths are overestimated by an average of 0.138 Ångstroms, the minimal overestimation is 0.105 Ångstroms and the maximal overestimation 0.160 Ångstroms. The overall average deviation is 0.127 Ångstroms.
As all deviations are within a small range, we can say that UFF gives bond lengths with a systematic error of about +0.13 Angstroms (100% consistent systematic error, the bond lengths are always overestimated). Therefore, we can expect to transfer this result to other molecules of this class and introduce an amount correction.
To describe the results of the MNDO method, we see that the bond lengths deviate from the experimental values to a much lesser extent. The greatest negative deviations (i.e. underestimation of bond lengths) can be found between layer 1 and the original structure with a deviation between -0.079 Angstroms (bond a). The greatest overestimation is found between both layer 2 and layer 3 and the original structure with a deviation of +0.027 Angstroms (bond f). To describe the direction of the deviations, we can say that they are to a lesser extent consistent than those of the UFF method. We derive degrees of consistency of the deviations between 27% and 64%, so we can say that the MNDO method underestimates bond lengths, but not in every case and depending on the type of the bond. Al-O bond lengths are rather underestimated where O-Si bond lengths are rather overestimated.
In table 2 we see the results for the bond angles. We have to distinguish between AlO-Si angles and O-Al-O bond angles. Both UFF and MNDO yield 100% consistent systematic errors in Al-O-Si angles, UFF yields angles that are underestimated by an average of -27.3 degree. The minimal underestimation amounts to -19.7 degrees whereas the maximal underestimation amounts to -41.2 degrees.
MNDO overestimates bond angles by an average of +14.0 degrees. The minimal overestimation amounts to 1.0 degrees whereas the maximal overestimation has the value of 20.4 degrees.
Thus, none of the methods describes the geometry around the Aluminum correctly, both methods show consistent systematic errors.
The O-Al-O angles are computed with much smaller deviations from the experimental values. Both UFF and MNDO give non-consistent systematic errors (consistent to a degree of 19% respectively 12 %) in a much smaller range than applied to Al-O-Si angles. Where the absolute arithmetic mean of the deviation within Al-O-Si angles was -27.3 degrees (UFF) and 14.0 degrees (MNDO), it is here as small as 3.6 degrees (UFF) and 4.3 degrees (MNDO).
The Dreiding Force Field did not converge in any of the performed optimizations, regardless whether they were started from the original structure or from a structure that has already been optimized with other methods.
Therefore, it could not be applied to this problem. A possible explanation is high sensitivity with regard to the input structure, resulting from the already mentioned treatment of angle bending with a cosine term that can be misleading.
HOMO comparisons
The next step was to compute HOMO values over a 3D grid and plotting isosurfaces at values of +0.05 and -0.05 e/Å3. From the change in shape and expansion as the cyclohexanes were removed we wanted to predict the change in chemical behavior.
We computed the HOMO of the anion, as no major difference between the HOMO auf the anion and the HOMO of the protonated molecule is expected.
Abbildung in dieser Leseprobe nicht enthalten
In this figure, we see the isosurfaces of the HOMO frontier orbitals plotted at values of +0.05 and -0.05 e/Å3. We notice that the HOMOs are mainly centered on the oxygens next to the central Aluminum, to a smaller extent on the oxygens separated by one silicon atom from the Aluminum.
Abbildung in dieser Leseprobe nicht enthalten
In this figure, we see the HOMO of our layer 3, where the four cyclohexanes furthest away from the Aluminum were removed. The perspective is the same as in the figure before; we notice no major change in size or shape of the orbitals. Therefore, we can assume that this model is still accurate enough for a prediction of the catalytic activity.
Abbildung in dieser Leseprobe nicht enthalten
We again reduced the size of our model, now only the four cyclohexanes closest to the central Aluminum remain. Even in this model, layer 2, we notice no major differences in size and shape of the HOMO orbitals, with regard to the full model.
Abbildung in dieser Leseprobe nicht enthalten
In this picture, we see the HOMO of layer 1, where all cyclohexanes were substituted by hydrogens. Again, we notice no major differences in shape and size of these parts of the orbitals centered on the oxygens. However, we can see an unrealistic size of the HOMO centered on the hydrogen (on the bottom of the picture). Thus, this model can be assumed to be to crude a model for the reaction modeling.
Therefore, because of the HOMO comparison and the geometric comparisons we can state that even layer 2 should be accurate enough to make predictions concerning the catalytic activity. This is an advantage because the computational costs of the later reaction modeling could be efficiently reduced.
Structure prediction of [Abbildung in dieser Leseprobe nicht enthalten]
In this step we tried to predict the structure of the anion with the sum formula [Abbildung in dieser Leseprobe nicht enthalten]. This structure was at the beginning of our calculations completely unknown, so we had to conceive a sensible input structure for our optimization algorithms.
As we compare this sum formula with the one above, we notice that formally spoken the OSiMe3 groups of the silsesquioxane with the known structure are substituted by the smaller OH groups. Thus, a coordination number of six seems to be possible (Yurii Gun’ko, private conversation), leading to the following input structure.
Abbildung in dieser Leseprobe nicht enthalten
In this structure, all the O-Al and O-Si bonds have an initial length of 1.4 Ångstroms. The hydrogens were initially set on one line with the centers of the cubes and the oxygens in the corners of the cubes with an O-H distance of 1.0 Ångstroms. All intracubic angles have the size of 90 degrees. The hydrogens bonded to the oxygens have been set on one line with the centers of the cubes and the oxygens to which they are bonded. The initial distance was again 1.0 Ångstroms.
This starting point seemed to be sensible, knowing that the multidimensional energy
landscape can give misleading results.
This structure was optimized using UFF and MNDO methods, leading to those results:
Abbildung in dieser Leseprobe nicht enthalten
In this figure, we see the UFF optimized structure starting from the cubic input structure.
Abbildung in dieser Leseprobe nicht enthalten
Here we see the MNDO optimized structure. Obvious are the greater Si-O-Si angles if we compare this structure with the UFF optimized structure.
At this stage, we had two optimized structures; one of them employed a semiempirical method where the other one used force fields. If we look at the structure, we see the most obvious difference in the Si-O-Si and Al-O-Si angles. The Al-O-Si angles all have a size of about 150 degrees when using the MNDO method where they have a size of about 115 degrees when using the Universal Force Field. The next step was to compare the computational results with the later resolved x-ray data.
Abbildung in dieser Leseprobe nicht enthalten
In this figure, we see the experimental core-structure, derived from x-ray data. If we compare this molecule graphically with our computed models, we notice a better agreement with the MNDO than with the UFF optimized structure. In the next section follows a quantitative analysis of both methods.
Results of structure prediction employing UFF and MNDO
The UFF optimized structure (the guess was the original structure) leads in general to the same results as we found within the optimization of the first silsesquioxane. The deviations with respect to bond lengths are 100% consistent systematic (positive) within UFF
The Al-O bond lengths of the close coordinated oxygens are overestimated by an average of 0.105 Ångstroms, the greatest overestimation amounts to 0.122 Ångstroms whereas the smallest overestimation has the value of 0.073 Ångstroms.
The Al-O bond lengths of the far coordinated oxygens (distances i and j in table 3) are overestimated by an average of 0.421 Ångstroms within a small range, the maximal deviation amounts to 0.429 Ångstroms and the minimal deviation to 0.413 Ångstroms.
The O-Si bond lengths are overestimated by an average of 0.122 Ångstroms with a minimum value of 0.094 Ångstroms and a maximum value of 0.149 Ångstroms.
MNDO applied to the original structure and to our cubic guess leads to quite similar results, what suggests that the energy minimization algorithm found two local minima close to the global minimum.
Bond lengths are in general slightly underestimated, but within x-ray resolution. The errors are not consistent systematic. In both cases MNDO gives a great overestimation of the bond length from Al to the far coordinated oxygens, applied to the cubic structure by an average of 1.17 Ångstroms, applied to the original structure by an average of 0.912 Ångstroms. All other bond lengths are predicted very well within a deviation of 0.04 Ångstroms.
So with respect to bond lengths MNDO predicted the right geometric arrangement of the center of the molecule.
The angle predictions by both methods are within small ranges of the actual value, UFF deviates by an abs. ar. mean of 5.6 degrees; the MNDO deviations are smaller, 2.9 degrees at the cubic input and 3.4 degrees applied to the original structure.
One could think that even UFF gives satisfying results, but a glimpse at the calculated structure compared to the experimental structure shows us that the angles further apart are predicted totally wrong. MNDO gives a much better result.
We can state that we predicted the structure of the unknown aluminosilsesquioxane right, using the MNDO method.
No bond length is predicted with a deviation greater than 0.04 Ångstroms and no angle is predicted wrong more than 2.9 degrees (concerning the near coordinated oxygens). In addition, we predicted the four-fold coordination of the central Aluminum right. We have an error of about one Ångstrom with regard to the bond lengths of the far away coordinated oxygens, but that does not affect the coordination number of the Aluminum.
Reaction modeling
I tried to model the reaction using various methods (MNDO, HF, B3LYP) starting with ethene about 2.5 Ångstroms apart from the acidic hydrogen of the silsesquioxane, as carried out by Senchenya and Kazansky (1991) on a 1T-cluster, employing a 3-21G basis set. An energy minimization was used, as we assumed this part of the reaction to decrease the energy content of the system.
They found that in a first step the ë-bond of ethene and the acidic hydrogen of the zeolite form a bond. This result could not be reproduced, with none of the employed methods and with none of various orientations. In all calculations the ethene disappeared from the hydrogen it was assumed to react with; possibly there is an energy barrier that could not be overcome with energy optimization methods.
In a next step, we had to figure out whether the reactants do not react or our used methods were not capable of dealing with this reaction. So we tried to perform the calculations on a 3T cluster and ethene, as for this reaction Senchenya and Kazansky (1991) already revealed a possible reaction path applied to a 1T cluster. The catalytic behavior was assumed not to differ between those molecules.
Therefore, both educts and product have been MNDO and B3LYP (with a 3- 21G and 6-31G basis set) optimized and the QST2 transition state search was performed This method is suggested for more difficult transition states and needs both reactants and products as input data.
In the input structure of the educts, ethene was simply about four Ångstroms apart from the acidic hydrogen; the product structure was the following in case of the 3T cluster.
Abbildung in dieser Leseprobe nicht enthalten
Even in the case of employing the QST2 method on a 3T cluster and ethene, we did not get a sensible transition state. Thus, the results by Senchenya and Kazansky (1991) on a 1T cluster could not be reproduced on a 3T cluster, although the catalytic activity of the 3T cluster was observed experimentally.
We have to assume that the employed methods were not able to deal with this case of reaction modeling
Abbildung in dieser Leseprobe nicht enthalten
Summary, Conclusions and further work
We compared the predictive capacities of UFF and Dreiding force fields and MNDO applied to a known and to an unknown aluminosilsesquioxane.
For the optimization of the known structure [Abbildung in dieser Leseprobe nicht enthalten] UFF gave consistent systematic errors of significant size, whereas MNDO gave results with a much smaller error and conserved the geometry around the center.
The Dreiding force field did not converge in energy and could therefore not be applied to this problem.
We derived the smallest suitable model for reaction modeling based on HOMO and geometric properties with the formula [Abbildung in dieser Leseprobe nicht enthalten]. We predicted the structure of the unknown aluminosilsesquioxane [Abbildung in dieser Leseprobe nicht enthalten] right using MNDO, whereas UFF is not applicable. The predicted structure was later confirmed by X-ray data.
None of the employed methods gave good results concerning the modeling of the reaction with ethene.
Further work should be made in the area of reaction modeling, where the reaction mechanism of the 3T cluster with ethene has either to be confirmed or disproved. Thereafter the studied aluminosilsesquioxanes have to be examined again with regard to their catalytic activity.
Acknowledgements
I want to thank Dr. Simon Bates for his motivating support and Dr. Yurii Gun’ko for his suggestions from experimental work. In addition, I want to thank Trinity College Dublin for facilitating this research project.
Computational Details
All calculations were performed with Gaussian98 on Pentium III - 450 MHz computers under Red Hat Linux. The convergence criteria were used as predefined. For visualization and measurement of structural properties, we used gOpenMol.
References
Allinger, N.L. 1977 J. Am. Chem. Soc. (99) 8127
Allured, V.S., C.M. Kelly and C.R. Landis 1991 J. Am. Chem. Soc. (113) Bates, S.P. and J. Dwyer 1993 J. Phys. Chem. (97) 5897
Bingham, R.C., M.J.S. Dewar and D.H. Lo 1975 J. Am. Chem. Soc. (97) 1302-1306 Dewar, M.J.S. and W. Thiel 1977 J. Am. Chem. Soc. (99) 4899-4907 Edelmann, F.T., Y.K.Gun’ko, S. Giessmann and F. Olbrich 1999 Inorg. Chem. (38) 210
Hall, G.G. 1951 Proc. Roy. Soc. (London) (A343) 541-552
Leach 1996, Molecular Modeling, Addison Wesley Longman Ltd., p 211 ff.
Mayo, S.L., B.D. Olafson and W.A. Goddard III 1990 J. Phys. Chem. (94) 8897 Meijer, E.L., R.A. van Santen and A.P.J. Jansen 1996 J. Phys. Chem. (100) 9282 Meisel, S.L., J.P. McCullogh, C.H. Lechthaler and P.B. Weisz 1976 Chem. Tech. (Heidelberg) (6) 86
Nicholas, J.B., A.J. Hopfinger, F.R. Trouw and L.E. Eton 1991 J. Am. Chem. Soc. (113) 4792
Pelmenschikov, A.G., N.U. Zhanpeisov, E.A. Paukshtis, L.V. Malyshave, G.M. Zhdirmov and
K.I.Zamaraev 1987 Dokl. Akad. Nauk SSSR (393) 915
Pople, J.A. and G.A. Segal 1965 J. Chem. Phys. (43) S136-S139
Rappe, A.K., C.J. Casewit, K.S. Colwell, W.A. Goddard III and W.M. Skiff 1992 J. Am. Chem. Soc. (114) 10624
Roothan, C.C.J. 1951 Rev. Mod. Phys. (23) 69-89
Senchenya, I.N. and V.B. Kazansky 1987 Kinet. Catal. (28) 566 Senchenya, I.N. and V.B. Kazansky 1991 Catal. Lett. (8) 317 Van Santen, R.A. and G.J. Kramer 1995 Chem. Rev. (95) 637
Weiner, S.J., P.A. Kollmann, D.A. Case, U.C. Singh, C. Ghio, G. Alagona, S. Profeta and
P. Weiner 1984 J. Am. Chem. Soc. (106) 765
32
Definition of bond and angle names
Abbildung in dieser Leseprobe nicht enthalten
-
-
-
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X. -
-
Téléchargez vos propres textes! Gagnez de l'argent et un iPhone X.