On the structure and reactivity of large complexes of metals and metalloids (PDF)

Summary

In this thesis four challenging problems concerning metal- or metalloid-containing systems were carefully studied in order to explain their properties, reactivities and selectivities. The first project focused on the dialkylzinc additions to unsaturated aldehydes. The reaction was proved to have bimetallic intermediates. The experimentally observed non-linear effect was addressed to the differences in stability of various zinc-ligand aggregates. Furthermore, the differences in regioselectivity between cinnamaldehyde and N-formylbenzylimine were traced back to the differences in the π-conjugation of both compounds. In contrary, the stereoselectivity of the reaction was controlled by the paracyclophane unit as well as by bulky substituents in the side-chain of the ligand. Based on these findings, a set of modifications for the ligand system were proposed and tested in silico. One of the compounds is expected to be superior to any previous catalysts. An important outcome from this study is the survey of methods used in calculations of reaction barrier heights. It is shown that economic DFT calculations, when extended by empirical dispersion corrections allow for the prediction of stereoselectivities. However, the prediction of the regioselectivities is a much more demanding task and only a consistent treatment of electron correlation can yield qualitative agreement with experimental findings. In this case, the computationally more expensive – but even for large systems still applicable – LPNO-CEPA/1 approach seems to be an efficient and reliable choice.

In the second project systems with significantly elongated bonds were studied. Germanium (Ge14[Si(SiMe3)3]5Li3(THF)6) and tin (Sn4Si[Si(SiMe3)3]4(SiMe3)2) cluster compounds appeared to be interesting examples of stable molecules with a number of undercoordinated atoms (three instead of four coordination partners for Ge or Sn) that are too far from each other to form a normal bond but are close enough to expect some kind of interaction. In both systems rhombic M4R6 (M = Ge, Sn; R – ligand) entities were found and the spin situation between two bridgehead atoms is best described as singlet biradicaloid. Moreover, an important contribution of the quantum-chemical studies was the demonstration that the stretching of the bond is not solely responsible for the observed biradicaloid character because the orientation and size of the ligands can also cause transitions between a radicaloid, a simple bounded system and a biradical with a quasi-degenerate singlet and triplet state. Furthermore, these calculations clearly show the limit of the applicability of single reference methods like DFT. Indeed, a more sophisticated treatment of static correlation, e.g. by using CASSCF wave functions is needed in order to describe the electronic state of such systems properly.

Smaller cluster compounds of nickel were considered in the next step. Experimental measurements concerning the dimer and the trimer are rather limited and are inconsistent with many theoretical findings. Within this work a benchmark study was performed on the Ni2 molecule, and the bond length was calculated to be 2.22 Å with a corresponding vibrational frequency of 286 cm-1. The NEVPT2 method was shown to provide data that is in excellent agreement with MRACPF results. The perturbation approach was used to estimate the electronic ground state including spin-orbit interactions. It was found to be 0𝑔+ in agreement with experiment. Moreover, it was pointed out that the TPSSh functional when combined with a broken symmetry approach provides bond lengths similar to those obtained with the NEVPT2 method. According to the literature, the bond is elongated by 0.1 Å upon electron attachment to the nickel dimer. This effect was observed both at the NEVPT2 and the TPSSh level and the bond expansion was addressed to the occupation of a 􀄱* orbital by the excess electron. A similar analysis was performed in the case of the nickel trimer. The NEVPT2 and TPSSh calculations agree that the ground state of the neutral cluster has a triangular geometry and a quintet spin state. However, at the NEVPT2 level the anion has clearly a linear form while the TPSSh functional shows linear and triangular isomers to be nearly equal in energy. The final distinction can be made by measuring the O-H stretching vibration in the systems where the nickel cluster is solvated by two ethanol molecule. According to TPSSh calculations, the difference between the O-H stretching frequencies of the two ethanol molecules will be much smaller (~15 cm-1) for the linear isomer than for the triangular geometry (~90 cm-1). Upcoming experimental IR studies should clearly be able to distinguish between these two forms.

Similar problems with DFT methods were encountered in the last application covered by this work. The aim was to analyse the molecular oxygen activation process by cobalt(II) diketonate complexes like bis[trifluoroacetylacetonato(−1)]cobalt(II) [Co(tfa)2]. This reaction is the first step in a mild and cheap aerobic oxidation of alkenols into functionalized tetrahydrofurans (molecules of high biological activity). However, various mechanisms can be proposed depending on the electronic structure of the adduct. We found, that the isolated catalyst with one water molecule [Co(tfa)2(H2O)] has a high spin quartet state. As dioxygen approaches the cobalt atom, the quartet state couples with a triplet dioxygen molecule and forms a sextet, a quartet, and a doublet spin state with the high-spin state being the lowest in energy. At the equilibrium Co–O2 distance of 1.9 Å, Co(tfa)2(H2O)(O2) has a doublet superoxo Co(III) ground state with the unpaired electron residing on the oxygen moiety, in a nearly unchanged O2 π* orbital. Here, it should be emphasised that even if DFT properly describes the electronic ground states of the catalyst and adduct, the relative spin state energies are, however, highly dependent on the functional used. A clear distinction would not be possible without NEVPT2 multireference calculations.

The calculations presented in this work established a solid base for further investigations concerning large metal- and metalloid-containing systems. The results evidently show that there is a great need to go beyond DFT in such computations. Elongated bonds or neardegeneracy problems within 3d transition metals are certainly the limit of single reference methods. In such cases, the CASSCF approach together with an inexpensive estimation of the remaining correlation energy (NEVPT2) was proven to be as reliable as the CI-based methods. Moreover, the regioselectivity prediction in reactions involving only closed shell species requires a consistent treatment of electron correlation in the transition states which can have very different structures. The progress made with efficient and accurate approximations permits the use of high-level wave function based methods like CEPA for molecules with over 100 atoms, even with a decent basis set. Through the careful evaluation of our calculations, the stereoselectivity was shown to be properly described with most of the modern DFT functionals when augmented with dispersion corrections. In this work also a number of predictions concerning the chemical structures and their reactivities have been made. The provided quantitative data can be directly compared to the perspective experimental studies what will certainly assist the deep understanding of the metal-based catalytic processes.


Zusammenfassung

In der vorliegenden Arbeit wurden verschiedene metallhaltige Komplexe und metalloide Clusterverbindungen mit Hilfe quantenchemischer Methoden untersucht. Dadurch war es unter anderem möglich, die Eigenschaften und die Reaktivität dieser Systeme zu verstehen sowie die Selektivität bei Reaktionen zu erklären.

Die erste Fragestellung befasste sich mit der Addition von Dialkylzink an ungesättigte Aldehyde. Hier konnte gezeigt werden, dass die Reaktion über bimetallische Intermediate verläuft. Der im Experiment beobachtete nicht-lineare Effekt kommt dabei durch die Stabilitätsunterschiede verschiedener Zink-Ligand-Aggregate zustande. Darüberhinaus konnte gezeigt werden, dass die Unterschiede in der Regioselektivität von Zimtaldehyd und Nformylbenzylimin auf der unterschiedlichen π-Konjugation beruhen. Die Stereoselektivität der Reaktion selbst wird durch die Paracyclophan-Einheit sowie durch sterisch anspruchsvolle Substituenten in der Seitenkette des Liganden kontrolliert. Basierend auf diesen Erkenntnissen wurden mehrere Modifikationen für den Liganden vorgeschlagen und in silico überprüft. Eine dieser Verbindungen sollte allen bisher verwendeten Katalysatoren überlegen sein. Ein weiteres wichtiges Ergebnis ist der erhaltene Überblick über die Qualität verschiedener quantenchemischer Methoden für die Berechnung von Reaktionsbarrieren. Dichtefunktionalberechnungen sind beispielsweise zur Vorhersage der Stereoselektivität geeignet, sofern eine empirische Korrektur für Dispersionswechselwirkungen berücksichtigt wird. Im Gegensatz dazu ist die Vorhersage von Regioselektivitäten deutlich schwieriger und für die qualitative Übereinstimmung mit experimentellen Daten ist eine konsistente Beschreibung der Elektronenkorrelation erforderlich. Eine geeignete Wahl hierfür könnte die in dieser Arbeit verwendete LPNO-CEPA/1-Methode darstellen, welche zwar deutlich rechenintensiver, aber dennoch für große Systeme anwendbar ist.

Im Rahmen des zweiten Projekts wurden metallhaltige Systeme mit deutlich aufgeweiteten Bindungsabständen untersucht. Ausgewählte Clusterverbindungen mit Germanium (Ge14[Si(SiMe3)3]5Li3(THF)6) und Zinn (Sn4Si[Si(SiMe3)3]4(SiMe3)2) sind interessante Beispiele für stabile Verbindungen, bei denen einige Atome eine Koordinationszahl von drei anstatt von vier aufweisen. Der Abstand zwischen den dreifach koordinierten Atomen ist einerseits zu groß für die Ausbildung einer gewöhnlichen Bindung, aber andererseits klein genug um eine Wechselwirkung zu erwarten. Bei beiden genannten Systemen liegen rhombische M4R6-Einheiten (M = Ge, Sn; R – Ligand) vor und die elektronische Situation zwischen den beiden subvalenten Atomen lässt sich am besten als Singulett-Biradikaloid beschreiben. Mit Hilfe der quantenchemischen Untersuchungen konnte gezeigt werden, dass die Aufweitung der Bindung nicht allein für den beobachteten biradikaloiden Charakter verantwortlich ist. Vielmehr können durch die Orientierung und die Größe der Liganden Übergänge zwischen radikaloiden, gebundenen Zuständen und echten Biradikalen mit quasientarteten Singulett- und Triplett-Zuständen verursacht werden. Des Weiteren wurde anhand dieser Berechnungen die Grenze der Anwendbarkeit von Eindeterminanten-Ansätzen wie DFT deutlich. Für die korrekte Beschreibung des elektronischen Zustands solcher Systeme ist eine angemessene Berücksichtigung der statischen Korrelation, beispielsweise durch CASSCF-Verfahren, notwendig.

Im dritten Teil der Arbeit wurden kleinere Nickelcluster untersucht. Für das Dimer und das Trimer liegen bisher nur wenige experimentelle Daten vor, welche jedoch nicht im Einklang mit den meisten Berechnungen sind. Im Rahmen der vorliegenden Arbeit wurde daher eine Benchmark-Studie für das Ni2-Molekül durchgeführt. Für die Bindungslänge ergab sich ein Wert von 2.22 Å und die Schwingungsfrequenz beträgt 286 cm-1. Die NEVPT2-Methode liefert hierbei Daten, welche in sehr guter Übereinstimmung mit den MRACPF-Ergebnissen sind. Der störungstheoretische Ansatz wurde verwendet, um den elektronischen Grundzustand unter Berücksichtigung von Spin-Bahn-Wechselwirkungen zu bestimmen. In Übereinstimmung mit dem Experiment konnte ein 0𝑔+-Zustand zugeordnet werden. Darüberhinaus wurde gezeigt, dass das TPSSh-Funktional im Rahmen des Broken-Symmetrie-Ansatzes ähnliche Bindungslängen wie die NEVPT2-Methode ergibt. Nach Literaturangaben verlängert sich die Bindung im Nickel-Dimer um 0.1 Å durch Anlagerung eines Elektrons. Dieser Effekt konnte sowohl mit NEVPT2 als auch mit TPSSh reproduziert werden und beruht auf der Besetzung eines 􀄱*-Orbitals durch das zusätzliche Elektron. Ähnliche Untersuchungen wurden für das Nickel-Trimer durchgeführt. Hier ergaben die NEVPT2- und TPSSh-Berechnungen, dass der Grundzustand des neutralen Clusters eine dreieckige Geometrie und einen Quintett-Zustand besitzt. Für das entsprechende Anion konnte mit NEVPT2 eindeutig eine lineare Geometrie zugeordnet werden, während im Fall von TPSSh das lineare und das dreieckige Isomer ähnliche Energien aufweisen. Eine Möglichkeit zur Unterscheidung beider Spezies bietet die Messung der O-H-Streckschwingungsfrequenz im System [Ni3(EtOH)2]-. Gemäß den TPSSh-Berechnungen unterscheiden sich die Frequenzen der beiden O-H-Streckschwingungen im Fall des linearen Isomers weniger (~15cm-1) als bei der dreieckigen Geometrie (~90cm-1). Zukünftige IRExperimente sollten daher in der Lage sein, die Geometrie eindeutig zu bestimmen.

Beim letzten Anwendungsprojekt traten im Zusammenhang mit DFT-Methoden ähnliche Schwierigkeiten auf. Das Ziel war es hier, den Aktivierungsprozess von molekularem Sauerstoff durch Kobalt(II)diketonat-Komplexe wie [Co(tfa)2] zu analysieren. Diese Reaktion ist der erste Schritt in einer milden und kostengünstigen aerobischen Oxidation von Alkenolen zu funktionalisierten Tetrahydrofuranen (Moleküle mit hoher biologischer Aktivität). Jedoch sind abhängig von der elektronischen Struktur des Addukts verschiedene Mechanismen denkbar. Es konnte gezeigt werden, dass ein Komplex aus einem Katalysatormolekül mit einem Wassermolekül [Co(tfa)2(H2O)] einen Quartett-Grundzustand besitzt (High-Spin-Zustand). Durch Annäherung des Sauerstoffmoleküls (Triplett) an das Kobaltatom können der Quartett- und der Triplett-Zustand zu einem Sextett-, einem Quartett- und einem Dublett-Zustand koppeln, wobei der High-Spin-Zustand energetisch am günstigsten ist. Beim Gleichgewichtsabstand (Co–O2 1.9 Å) weist der Co(tfa)2(H2O)(O2)-Komplex einen Dublett-Grundzustand auf, wobei sich das ungepaarte Elektron in einem nahezu unveränderten π*-Orbital des Sauerstoffmoleküls befindet. An dieser Stelle muss betont werden, dass die DFTMethode zwar den Grundzustand des Katalysators und des Addukts richtig wiedergibt, aber dass die relativen Energien der einzelnen Spinzustände stark vom verwendeten Funktional abhängen. Für zuverlässige Ergebnisse sind daher NEVPT2-Berechnungen erforderlich.

Die in Rahmen dieser Arbeit durgeführten quantenchemischen Berechnungen bilden eine zuverlässige Grundlage für weitere Untersuchungen an großen metall- und metalloidhaltigen Systemen. Die Ergebnisse verdeutlichen, dass standardmäßig verwendete Dichtenfunktionalmethoden für die gegebenen Fragestellungen oftmals nicht ausreichen. Aufgeweitete Bindungen und nahezu entartete Zustände bei 3d-Übergangsmetallen stellen dabei die Grenze für die Anwendbarkeit von Eindeterminanten-Anzätzen dar. Für diese Fälle wurde gezeigt, dass die NEVPT2-Methode (CASSCF-Verfahren mit Abschätzung der fehlenden Korrelationsenergie) vergleichbar gute Ergebnisse liefert wie CI-basierte Methoden. Zudem erfordert die Vorhersage der Regioselektivität für Reaktionen, selbst wenn nur closed-shell Spezies auftreten, eine konsistente Beschreibung der Elektronenkorrelation in den Übergangszustände, da hier sehr unterschiedliche Strukturen vorliegen können. Durch Fortschritte hinsichtlich der Effizienz können korrelierte wellen-funktionbasierte Methoden wie CEPA für Moleküle mit mehr als 100 Atomen auch in Kombination mit sinnvollen Basissätzen angewendet werden. Anhand der sorgfältigen Auswertung der Ergebnisse dieser Arbeit konnte gezeigt werden, dass die Stereoselektivität bereits mit den meisten Dichtenfunktionalen richtig wiedergegeben werden kann, sofern Dispersionswechselwirkung berücksichtig werden. Im Rahmen der vorliegenden Arbeit wurden zudem mehrere Vorbesagen bezüglich der Struktur und der Reaktivität ausgewählter Verbindungen gemacht. Diese Daten werden für die Auswertung zukünftiger Experimente hilfreich sein und somit zum tieferen Verständnis metallbasierter Katalyseprozesse beitragen.