Institut für Physikalische und Theoretische Chemie
Refine
Has Fulltext
- yes (64)
Is part of the Bibliography
- yes (64)
Year of publication
Document Type
- Doctoral Thesis (61)
- Journal article (3)
Language
- German (64) (remove)
Keywords
- Kohlenstoff-Nanoröhre (9)
- Spektroskopie (8)
- Exziton (7)
- Photodissoziation (6)
- Photoionisation (6)
- Quantenmechanik (6)
- Raman-Spektroskopie (6)
- Molekulardynamik (5)
- Quantendynamik (5)
- REMPI (5)
Institute
EU-Project number / Contract (GA) number
- 646737 (1)
Die Entwicklung maßgeschneiderter Proteinliganden ist ein integraler Bestandteil unterschiedlicher wissenschaftlicher Disziplinen, wie z.B. Wirkstoffentwicklung. Die vorliegende Arbeit befasst sich mit der reversiblen Inhibition in Form von kovalent gebundenen Enzym-Inhibitor-Komplexen der humanen Golgi-alpha-Mannosidase II (GM II) und der Cysteinprotease Rhodesain. Beide Enzyme sind erfolgversprechende Targets in der Bekämpfung von zwei sehr unterschiedlichen Erkrankungen. Einerseits die Golgi-alpha-Mannosidase bei der Behandlung der Tumorprogression und andererseits die Cysteinprotease Rhodesain bei der Behandlung der Afrikanischen Schlafkrankheit. Die Arbeit an den zwei Enzymen unterteilt sich in zwei Teilprojekte.
Die Entwicklung von maßgeschneiderten kovalent-reversiblen Inhibitoren für die genannten Enzyme wurde im Rahmen eines in-house entwickeltes Protokolls zwecks des rationalen Designs kovalenter Inhibitoren, durchgeführt. Dieses Protokoll basiert auf einer sich gegenseitig unterstützenden Zusammenarbeit zwischen Theorie und Experiment. Die vorliegende Arbeit befasst sich mit den theoretischen Untersuchungen mit Hilfe der quantenmechanischen (QM) als auch mit Hilfe der kombinierten quantenmechanisch/molekülmechanischen (QM/MM) Methoden zu den genannten Enzymen. In einem ersten Schritt des Protokolls geht es um die Anwendung von Screeningverfahren. In einem Screening werden Leitstrukturen, zunächst in Lösung (Schritt I), für eine weitere Untersuchung im Enzym (Schritt II) evaluiert. So können die Inhibitoren, für die experimentelle Mess- oder theoretische Dockingdaten vorhanden sind, als eine Leitstruktur betrachtet werden. Durch das Screening unter Verwendung der quantenmechanischen (QM-Modell) Methode kann eine Reihe von Inhibitoren nach einem sich konsistent veränderndem Muster erstellt werden und auf Bindungsparameter hin untersucht werden (Schritt I). Diese Parameter sind Reaktionsenergien und Höhen der Reaktionsbarriere einer Inhibitionsreaktion. Reaktionsenergien werden in dieser Betrachtung quantenmechanisch innerhalb der Born-Oppenheimer-(BO)-Näherung und im Rahmen des Konzeptes der Potentialhyperflächen (PES) als relative Energien zwischen den optimierten Geometrien der Produkte und der Edukte auf einer Potentialhyperfläche für die Inhibitionsreaktion ermittelt. Die Höhen der Reaktionsbarrieren werden durch die relativen Energien zwischen den Geometrien der Edukte und der Zwischenstufen oder Übergangszustände abgeschätzt. Unter Inhibitionsreaktion wird eine chemische Reaktion verstanden, bei der eine kovalente Bindung zwischen dem Inhibitormolekül und den Aminosäuren in der aktiven Tasche eines Enzyms ausgebildet wird. Für den Schritt I werden die Aminosäuren der aktiven Tasche durch kleine Moleküle, wie Essigsäure und Methanthiol, angenähert. Die kovalent-reversiblen Inhibitoren sollten in dieser Betrachtung nur leicht exotherme Reaktionen mit den relativen Energien im Bereich -5 bis -10 kcal/mol aufweisen. Der experimentelle Teil liefert währenddessen die Synthese der neuen Inhibitoren und die Nachweise zur kovalenten Bindung mit Hilfe massenspektrometrischer Messungen (Schritt I). Die passenden Kandidaten aus dem ersten Schritt des Protokolls, d.h. Inhibitoren mit gewünschten Bindungsparametern, werden durch die QM/MM-Berechnungen im Enzym (Schritt II) und durch die experimentellen Messungen an den Enzym-Inhibitor-Komplexen in Assays (Schritt II) analysiert. Die Untersuchungen für die Stufe II des Protokolls umfassen die Berechnungen der Reaktionsprofile und Minimumenergiereaktionspfade für die chemischen Reaktionen von Inhibitoren im Zielenzym. Ein Pfad minimaler potentieller Energie, der zwei Minima (hier Edukt und Produkt) verbindet, stellt ein Reaktionsprofil für eine chemische Reaktion dar. In der vorliegenden Arbeit wird dies auch als Minimumenergiepfad (MEP) bezeichnet. Der Letztere lässt sich durch die Nudged Elastic Band (NEB)-Methode und mittels Potentialhyperflächen darstellen. Die Reversibilität der Inhibitoren wurde anhand der berechneten chemischen Reaktionen in Form von erstellten Reaktionsprofilen analysiert und diskutiert. Durch Protein-Ligand Docking (Schritt III) wird ein Screening von variierbaren Erkennungseinheiten der neuen Inhibitoren durchgeführt. Die Ergebnisse der Untersuchungen aus dem dritten Schritt liefern Hinweise zur Weiterentwicklung der ausgewählten Inhibitoren. Die letzte Stufe des in-house Protokolls besteht in der erneuten Untersuchung der optimierten Inhibitoren mit Hilfe von Theorie und Experiment (Schritt IV). Die theoretische Untersuchung anhand von QM/MM-Berechnungen überprüft, ob die Inhibitionsreaktion der reaktiven Kopfgruppe nach der Änderung der Erkennungseinheit des Inhibitors weiterhin effektiv und nach dem gleichen Mechanismus mit der aktiven Seite des Enzyms ablaufen kann. Die experimentelle Untersuchung liefert, ähnlich wie im Schritt II, die messbaren Ergebnisse der Inhibition in Hinblick auf die Bindungseigenschaften und die Entstehung der Nebenprodukte.
Die Untersuchungen am System Mannosidase GM II wurden in Zusammenarbeit mit den Arbeitskreisen von Prof.Dr. J. Seibel und Prof.Dr. T. Schirmeister durchgeführt. Die Leitstruktur zur Entwicklung des kovalent-reversiblen Inhibitors stellt der cyclische O,O-Acetal-Inhibitor (bestimmt anhand von Dockingexperimenten an der beta-L-Anhydrogulose durch Arbeitskreis Prof.Dr. J. Seibel) dar.
Die Ergebnisse der theoretischen Studie liefern für den ersten Schritt im Rahmen des Protokolls den geeigneten Kandidaten aus einer Menge von insgesamt 22 modellierten Inhibitoren für die reversible Inhibition der Mannosidase GM II durch die Ausbildung einer kovalenten Bindung. Hierzu zählen zunächst die thermodynamischen Modellberechnungen der Inhibitionsreaktion, welche die Reaktionsenergien für alle Kandidaten des Screenings liefern. Die Inhibitionsreaktion wird in diesem Schritt als Additionsreaktion von Essigsäure an den Inhibitor-Kandidaten modelliert. Für die Leitstruktur resultiert eine thermoneutrale Beschreibung der Reaktion mit Essigsäure und dient im Weiteren als Referenz. Der Inhibitor Nr.7 (der cyclische N,O-Acetal-Inhibitor) zeigt mit -7,7 kcal/mol eine leicht exotherme Reaktion und somit eine bessere Triebkraft der untersuchten Reaktion im Vergleich zur Referenz. Die beiden Inhibitoren wurden dann für Stufe 2 des Protokolls untersucht, in der eine Analyse der Reaktionsprofile im Enzym mit Hilfe der QM/MM-Methodik durchgeführt wurde. Die Ergebnisse des zweiten Teils der Studie zeigen, dass der cyclische N,O-Acetal-Inhibitor eine deutlich bessere Affinität zur aktiven Seite der GM II im Vergleich zu seiner Leitstruktur aufweisen sollte. Dies zeigt sich auch in der deutlich höheren Triebkraft der Inhibitionsreaktion von ca. -13 kcal/mol. Dieser Energiebeitrag ist klein genug, um eine Reversibilität der Inhibitionsreaktion gewährleisten zu können. Das bedeutet auch, dass der N,O-Acetal-Inhibitor im Vergleich zur Referenzstruktur eine deutlich stärkere Inhibition bedingen sollte. Berücksichtigt man dann noch, dass die Reaktion laut Berechnungen nur leicht exotherm sein sollte, erhält man die Möglichkeit einer reversibel stattfinden kovalenten Hemmung. Zusammenfassend liefert dieser Teil der Arbeit, der mit Hilfe der QM- und QM/MM-Berechnungen durchgeführt wurde, ein reaktives molekulares Gerüst mit den gewünschten Eigenschaften.
Durch die theoretischen Untersuchungen (MD-Simulationen am Enzym-Inhibitor-Komplex) konnte außerdem eine zur Komplexbildung geeignete Konformation der Leitstruktur sowie des neuen Inhibitors gefunden werden. Die reversibel agierenden Acetal-Inhibitoren befinden sich in der aktiven Tasche in einer energetisch höher liegenden Twist-Boot-Konformation und begünstigen mit zwei entstehenden Bindungen zum Zn2+-Ion die oktaedrische Koordination im Enzym. Als Teil dieser Arbeit wurden NEB-Berechnungen zur Bestimmung von Minimumenergiepfadendurchgeführt. Dies lieferte erweiterte Einblicke in der Berechnung von Reaktionsmechanismen jeweils auch in Kombination von 2- bzw. 3-dimensionalen Scans. Auch in der Beschreibung von Protonenübertragungsreaktionen nach Grotthus, die einem Umklappen der kovalenten Bindungen entsprechen, erhält man hier Geometrien für die Teilschritte und somit eine detaillierte Beschreibung des Vorgangs.
Der Mechanismus der Inhibition von GM II durch die Leitstruktur beinhaltet einen Wasser-katalysierten (oder auch Wasser-vermittelten) Ringöffnungsschritt in der Tasche des Enzyms. Die Testrechnungen zum Protontransfer haben gezeigt, dass der Protontransfer über ein oder mehrere Wassermoleküle unter Verwendung von Standard-PES-Berechnungen nicht spontan stattfindet. Die Berechnung des MEP durch das Erstellen einer 3-dimensionalen Potentialhyperfläche kann nur dann sinnvolle Ergebnisse liefern, wenn der Protontransfer vom Aspartat Asp341 zum Inhibitor über zwei Wassermoleküle explizit berücksichtigt wird. In diesem Fall ist die Berechnung der PES kein Standard und erfordert eine zusätzliche Variation der Bindungsabstände O-H der beteiligten Moleküle des Protontransfers. Die Details für die zusätzliche Variation der Bindungsabstände O-H bei der Berechnung der 3-dimensionalen PES haben die NEB-Berechnungen geliefert. Der NEB-Formalismus hat sich in der Beschreibung dieser komplexen Reaktionskoordinaten als besser geeignet erwiesen und wurde in dieser Arbeit aus diesem Grund hauptsächlich verwendet. Die Berechnung des Protonentransfers während einer Hemmungsreaktion durch zwei Wassermoleküle mit der NEB-Methode hat den MEP ermittelt, welcher zunächst nicht auf der Grundlage eines 3-dimensionalen Scans ermittelt werden konnte. Solche QM/MM-Rechnungen wurden im Rahmen des in-house Protokolls zum ersten Mal durchgeführt. Dieser Protontransfer ist mit dem Grotthus-Mechanismus konform und kann plausibel anhand einer Klapp-Mechanismus-Betrachtung nachvollzogen werden.
Mit Hilfe der NEB-Methode ist es möglich MEPs effektiv und relativ schnell zu ermitteln. Es werden sowohl die Geometrien entlang des Pfades wie auch die einzelnen relativen Energien erhalten. Zur Überprüfung der gefundenen Übergangszustände wurden die einzelnen Strukturen mit Hilfe der Normalmodenanalyse weiter untersucht und konnten verifiziert werden.
Die MEP-Berechnungen für den Inhibitor Nr.1 ermöglichen die Etablierung eines Protokolls zur Berechnung eines Reaktionspfades über mehrere Moleküle, welches anschließend zur Berechnung des MEP für den Inhibitor Nr.7 angewendet wird. Das Protokoll beinhaltet in seiner einfachen Form die Ermittlung der Two-End-Komponenten einer chemischen Reaktion - Geometrien von Reaktant und Produkt. Betrifft dies eine Reaktion, die über mehrere Moleküle, z.B. Wassermoleküle oder deren Netzwerk, stattfindet, wird die Aufgabe komplexer. In diesem Fall ist eine Berechnung mit Hilfe des NEB-Moduls wesentlich produktiver als die Charakterisierung mit Hilfe der 3-dimensionalen PES. Der Vorteil liegt in der kollektiven Beschreibung der Reaktionskoordinaten, sodass die entscheidenden Reaktionskoordinaten und Variablen für die Durchführung von Scans nicht einzeln bestimmt werden müssen. Dennoch kann es hier bei komplexen Reaktionskoordinaten auch zu Konvergenzproblemen bzw. zu langwierigen Optimierungszyklen kommen.
Als weiteres Resultat liefern die durchgeführten MEP-Berechnungen Einblicke in die katalytischen Eigenschaften der Wassermoleküle für den Protonübertragungsmechanismus nach Grotthus. Die Daten zeigen, dass die Barriere am niedrigsten wird, wenn zwei Wassermoleküle beim Protontransfer beteiligt sind. Wenn nur ein oder gar kein Wassermolekül die Ringöffnung katalysiert, steigt die Barriere auf 12 und 17 kcal/mol.
Die Untersuchung in diesem Teil der Arbeit lässt zudem Einblicke in die nukleophile Substitution der Vollacetale in der Enzym-Tasche der GM II erlangen. Die Rechnungen deuten darauf hin, dass die Vollacetal-Inhibitoren durch Wassermoleküle in der Tasche aktiviert werden. Die ausgebildeten Wasserstoffbrückenbindungen begünstigen die Geometrie des Enzym-Inhibitor-Komplexes. Dies befördert die Ringöffnungreaktion gleichzeitig mit dem nucleophilen Angriff des Aspartatrestes an dem C1-Atom des Inhibitors. Im Falle des gemischten Acetal-Inhibitors hingegen wird die Treibkraft bereits durch die Einführung des Stickstoffatoms deutlich erhöht. Durch die richtig angeordneten Grotthus-Wassermoleküle ist in diesem Fall die Barriere der Protonübertragung durch das Aspartat-Aspartat-System der GM II (Asp341/Asp240) sekundär. Betrachtet man die Schwingungsbewegung entlang der imaginären Moden der Übergangszustände, sind diese in beiden E-I-Komplexen ähnlich. Hierbei wird eine synergistische Bewegung der Bindungsabstände OD2-C1-O6 (Inhibitor Nr.1) bzw. OD2-C1-N (Inhibitor Nr.7) beobachtet.
Die Entwicklung der kovalent-reversiblen Inhibitoren für das Enzym Rhodesain wurde in Zusammenarbeit mit dem Arbeitskreis Prof.Dr. T. Schirmeister durchgeführt. Als Leitstruktur zur Entwicklung des neuen kovalent-reversiblen Vinylsulfon-Inhibitors 4-Pyridyl-Phenylalanyl-Homophenylalanyl-alpha-Fluor-Phenylvinylsulfon dient in diesem Projekt der kovalent-irreversibel bindende Inhibitor K777, für den kristallographische Daten bekannt sind. Im Rahmen des Protokolls wurde eine Reihe von Inhibitoren untersucht, in denen ein alpha-H-Atom der Vinylsulfon-Einheit (im Weiteren VS für Vinylsulfon) durch verschiedene Gruppen X substituiert wurde. Für den zunächst vorgeschlagenen Cyano-Substituent (CN) in einem VS-Inhibitor ergab sich bei Berechnungen in einem polaren Lösungsmittel eine relativ niedrige Reaktionsenergie, d.h. es wurde eine reversible Reaktion vorhergesagt. Dies wurde experimentell bestätigt. Die theoretischen und experimentellen Ergebnisse von Schritt II widersprechen sich aber. Während experimentell eine schwache reversible Hemmung gefunden wurde, sagten die Berechnungen keine Hemmung voraus. Tatsächlich zeigte sich im Nachhinein, dass die experimentell gefundene Hemmung nicht-kompetitiv ist, d.h. nicht in der aktiven Tasche stattfindet. Im Rahmen des Protokolls werden dagegen nur die kompetitiv interagierenden Inhibitoren ausgewertet. An dieser Stelle lassen sich die anhand theoretischen Methoden erhaltenen Daten über die Reversibilität der Hemmung (Reaktion im aktiven Zentrum) mit den Ergebnissen aus den experimentellen Untersuchungen (Reaktion außerhalb des aktiven Zentrums) nicht vergleichen. Durch den Wechsel von CN zu Halogenen wurde schließlich eine neue Reihe von Inhibitoren auf VS-Basis entwickelt. Die Berechnungen von Reaktionsenergien in Lösung und im Enzym haben für diese Inhibitoren eine reversible Hemmung vorhergesagt. Allerdings findet man eine einfache Additionsreaktion an der Doppelbindung der Vinylsulfon-Gruppe. Für X=CN wurde von einer SN2-Reaktion ausgegangen. Für X=Br fand man, dass sich nach der Addition ein HBr-Molekül abspaltet, sodass die Hemmung insgesamt irreversibel ist. Da die Substitutionsreaktion ein irreversibler Prozess ist und die Freisetzung von Bromwasserstoff durch die experimentellen Untersuchungen bestätigt werden konnte, scheint Fluor ein geeigneter Substituent zu sein (X=F). Hier konnte man auch experimentell eine deutlich bessere Hemmung messen. Es wurden daher die Berechnungen im Enzym für Systeme mit den Inhibitoren K777-X mit X=F und X=H (K777-F- und K777-H-Inhibitor) durchgeführt und analysiert. In der vorliegenden Arbeit wurde versucht, die Reversibilität des K777-F-Inhibitors gegen die Irreversibilität des K777-H-Inhibitors durch die quantenmechanischen Berechnungen im Rahmen des Protokolls darzulegen.
Die QM/MM-Berechnungen unterteilen sich in zwei Bereiche. Zunächst wurde das Reaktionsprofil (auch Reaktionspfad) der Additionsreaktion des K777-X-Inhibitors an die aktive Tasche von Rhodesain ausgehend von der vorhandenen Kristallstruktur (PDB-Datei) berechnet. Im Folgenden wird dieses Teilergebnis als XP-Berechnung (im Weiteren XP für X-ray-Pfad) bezeichnet. Alle vier PES (X=H, F, Br und Cl) weisen prinzipiell die gleiche Form auf. Es ergeben sich aber Unterschiede in den berechneten Reaktionsenergien der Additionsreaktion (-20, -16, -10 und -11 kcal/mol). Die Reaktionsenergien der Substituenten Brom und Chlor entsprechen dem Bereich für reversible Reaktionen (ca. -10 kcal/mol), wobei Fluor mit -16 kcal/mol einen Grenzfall darstellt. Die Konturen der beiden PES (X=H vs. X=F) sind allerdings sehr ähnlich: In beiden Fällen findet sich für das anionische Intermediat kein Minimum. In der Potentialhyperfläche für X=F steigt die Barriere der Rückreaktion zwischen dem Intermediat und dem nicht-kovalenten Komplex auf etwa 5 kcal/mol an, die Rückreaktion ist im Vergleich zu dem X=H mit ca. 1,5 kcal/mol leicht exotherm. Das veränderte Verhältnis zwischen der Höhe der Reaktionsbarriere und dem Betrag der Reaktionsenergie (der Übergang von der endothermen zur exothermen Reaktion) auf diesem Abschnitt der PES könnte dazu beitragen, dass die Gesamtreaktion insgesamt reversibel ablaufen kann. Die Reversibilität des Inhibitors mit dem Substituenten Fluor lässt sich auf diesem Schritt der Untersuchung durch die Absenkung der Reaktionsenergie der Additionsreaktion auf etwa 16 kcal/mol erklären, denn die irreversible Reaktionen wurden bisher mit deutlich höheren Reaktionsenergien assoziiert. Die erhaltenen nicht-kovalenten Enzym-Inhibitor-Komplexe der XP-Berechnungen wurden in einem zweiten Teilergebnis weiter verwendet, indem der Reaktionspfad der Additionsreaktion des K777-X-Inhibitors vom nicht-kovalenten Enzym-Inhibitor-Komplex zum kovalenten Enzym-Inhibitor-Komplex hin berechnet wurde. Im Folgenden wird dieses Teilergebnis als NP-Berechnung (NP für Nicht-kovalente-Pfad) bezeichnet.
Die Berechnung der Reaktionsprofile der Additionsreaktion des VS-Inhibitors für X=H und X=F am alpha-Kohlenstoffatom der VS-Kopfgruppe lieferte konsistente Ergebnisse in Bezug auf die Reaktionsenergien. Ähnlich den XP-Berechnungen, wurde ebenfalls die Tendenz der Absenkung der Reaktionsenergie von -7 kcal/mol (X=H) auf -4,3 (X=F) und -0,9 kcal/mol (X=Br) beobachtet. Die Thermodynamik der Additionsreaktion wurde durch einen Wechsel des Substituenten X von H nach F in der VS-Kopfgruppe des K777-X Inhibitors beeinflusst, indem die niedrigere Energiedifferenz zwischen den Edukten und den Produkten erzielt werden konnte. Für beide Teile der Arbeit (XP- und NP-Berechnungen) implizierte dies einen Wechsel von einem irreversiblen zum einem reversiblen Verlauf in der Beschreibung der Reaktionsprofile. Die Ergebnisse des zweiten Teils der Arbeit (NP-Berechnungen) liefern nicht nur die konsistent geringeren Reaktionsenergien (Thermodynamik) sondern auch die höheren Reaktionsbarrieren der Additionsreaktion im Vergleich zu den Ergebnissen der XP-Berechnungen. Die Änderungen der Reaktionsbarrieren im NP-Ansatz weisen zusätzliche Diskrepanzen auf, wenn diese jeweils mittels der PES-Scan- und CI-NEB-Dimer-Methodik berechnet werden. Während die Barriere des irreversiblen Inhibitors K11777 mit dem NEB-Ansatz ca. 11 kcal/mol beträgt und durch den PES-Scan nur um 4 kcal/mol höher liegt, ergibt sich eine umgekehrte Situation beim Übergang zu Fluor als Substituent: Durch die NEB-Berechnung liegt die Barriere bei ca. 18 kcal/mol und durch den PES-Scan ergibt sich eine um 4 kcal/mol niedrig liegende Barriere. Um die Ergebnisse der NP-Berechnungen zu überprüfen, wurden diese QM/MM-Rechnungen wiederholt durchgeführt. In den beiden neu durchgeführten Berechnungen für die Inhibitoren K777-X mit X=H und X=F wurden nur sehr kleine Abweichungen gefunden, die kleiner als die Fehler der Berechnung sind. Die Startstrukturen für die Berechnung des MEP stammten aus der erneut durchgeführten MD-Simulation an der Geometrie des nicht-kovalenten Enzym-Inhibitor-Komplexes, welche die XP-Berechnung resultierte. Die Gesamtdauer der MD-Simulation wurde zu einem Wert von 9 ns gewählt, welche insgesamt 900 Startstrukturen entlang der Simulation lieferte. Die Berechnung ergab die Reaktionsenergie von -8,4 kcal/mol (-7,0 kcal/mol als erstes Ergebnis) und die relative Energie des Int-Komplexes von 13,2 kcal/mol. Somit beträgt die Barriere der Rückreaktion zur Freisetzung des Inhibitors K11777 (X=H) in Form von einem nicht-kovalenten Enzym-Inhibitor-Komplex einen Wert von 21,6 kcal/mol. In analoger Vorgehensweise wurde die Evaluierung der NP-Berechnung für den Inhibitor K777-X mit X=F durchgeführt. Die Reaktionsenergien in den beiden Berechnungen unterscheiden sich in einem marginalen Abstieg zu den Werten von -2,9 kcal/mol und -3,3 kcal/mol (-4,3 kcal/mol als erstes Ergebnis). Beide Berechnungen liefern zudem die relativ kleinen Anstiege der Reaktionsbarriere zu den Werten von 19,8 kcal/mol und 20,9 kcal/mol.
Für die Inhibitoren K777-X mit X=H und X=F entsprechen die gefundenen Barrieren einer verzerrten Konformation des nicht-kovalenten Enzym-Inhibitor-Komplexes, die als eine bioaktive Konformation bezeichnet werden kann. Der anionische Übergangszustand Int*, der oft in der Literatur als ein anionisches Intermediat der Additionsreaktion bezeichnet wird, wurde nur für den Inhibitor mit dem Substituenten Brom (K777-X mit X=Br) identifiziert. Da der Übergangszustand (Int* mit der relativen Energie von 11,1 kcal/mol) nur 1,5 kcal/mol über der bioaktiven Konformation (Int mit der relativen Energie von 12,6 kcal/mol) liegt und die NEB-Reaktionspfade alleine die Barrieren überschätzen, besitzen die anionischen Übergangszustände der Additionsreaktion der Inhibitoren mit X=H und X=F eine geschätzte relative Energie mit vergleichbaren Abweichungen von ca. 2 kcal/mol zu den identifizierten Int-Geometrien. Die durchgeführten Berechnungen zeigen, dass die Substituenten X=H und X=F im Vergleich zum Brom die anionischen Geometrien der nicht-kovalenten Enzym-Inhibitor-Komplexe jedoch mangelhaft bis ausreichend stabilisieren können. Zusätzlich liegt die geschätzte Energiedifferenz zwischen den Geometrien Int* und Int unter der möglichen Fehlergrenze der Berechnungen (ca. 3-4 kcal/mol). Aus diesem Grund misslang die Optimierung in Richtung der metastabilen anionischen Geometrien Int* mit Hilfe der CI-NEB-Dimer-Methodik im Fall der VS-Inhibitoren K777-X mit X=H und X=F.
Der direkte Vergleich der geometrischen Parameter der nicht-kovalenten Enzym-Inhibitor-Komplexe für den Inhibitor K777-F aus den XP-Berechnungen mit solchen aus den NP-Berechnungen lässt darauf schließen, dass die Geometrien der Enzym-Inhibitor-Komplexe der XP-Berechnung nur die lokalen Minima mit der verzerrten Geometrie des Inhibitors auf der PES darstellen und die Gesamtinformation über die Barrieren der Reaktion durch die Ergebnisse aus der NP-Berechnung ergänzt werden sollten.
Zusammenfassend sagen die Berechnungen für die reaktiven Kopfgruppen der Substanzklasse der halogenierten Vinylsulfone K777-X (X=Br, Cl und F) im Vergleich zur Leitstruktur des Vinylsulfon-Inhibitors K11777 deutlich geringere exotherme Additionsreaktionen im aktiven Zentrum von Rhodesain voraus. Darüberhinaus konnte anhand der QM/MM-Berechnungen ein experimentell gemessenen verlangsamten Verlauf der reversiblen Inhibition im Falle von X=F (Inhibitor K777-X) durch die relativ erhöhte Reaktionsbarriere im Vergleich zur Leitstruktur erklärt werden. Dieser Inhibitor dient zunächst als ein erfolgreich selektiertes reaktives Gerüst des neuen Inhibitors K777-X-S3 mit X=F und S3=4-Pyridyl (K777-F-Pyr), welcher mit Hilfe des Docking-Experiments (Schritt III durch die Arbeitsgruppe Prof.Dr. T. Schirmeister) deutlich verbessert werden konnte. Die Affinität des durch Docking verbesserten VS-Inhibitors mit Fluor als Substituent durch die eingeführte Seitenkette S3=4-Pyridyl (4-Pyridyl-Phenylalanyl-Homophenylalanyl-(Phenyl)-alpha-F-Vinylsulfon) stieg im Rhodesain von 190 nM zu 32 nM (Schritt IV, experimenteller Teil). Gleichzeitig konnte durch die QM/MM-Berechnungen in Schritt IV gezeigt werden, dass die Reaktion der reaktiven Kopfgruppe im neuen Inhibitor immer noch eine kovalent-reversible Hemmung von Rhodesain darstellt, auch wenn die Erkennungseinheit geändert wurde. Hierfür kann man die beiden Reaktionsprofile der NP-Berechnungen vergleichen. Die beiden fluorierten VS-Inhibitoren weisen eine Ähnlichkeit bezüglich der Barrierenhöhe und der Reaktionsenergie auf. Der fluorierte Vinylsulfon-Inhibitor K777-F wurde somit als ein neuer kovalent-reversibler Vinylsulfon-Inhibitor der Cysteinprotease Rhodesain erfolgreich eingefügt.
Simulationen zur transienten Absorptionsspektroskopie an Energie- und Ladungstransfersystemen
(2022)
Anregungsinduzierte Ladungstransferprozesse gemischtvalenter Verbindungen in einem, bzw. zwei Vibrationsfreiheitsgraden werden mithilfe vibronischer Modellsysteme untersucht. Anhand transienter und linearer Absorptionsspektren werden die berechneten mit experimentell bestimmten Daten verglichen. Eine detailliertere theoretische Analyse erfolgt unter den Gesichtspunkten der Populations- und Wellenpaketdynamik.
Darüber hinaus wird der Prozess der Exziton-Exziton-Annihilierung mithilfe eines elektronischen Modellsystems untersucht. Zu diesem Zweck werden, zusätzlich zu den oben genannten Methoden, spektroskopische Signale unterschiedlicher Emissionsrichtungen zum Vergleich herangezogen.
Um den jahrtausendealten Weg der Menschheit vom Papyrus über Buchdruck und siliziumbasierte Halbleiter in Richtung noch leistungsfähigerer Technologien zu gehen und weiterhin Heureka-Momente zu schaffen, bieten Kohlenstoffnanoröhren ein weites Forschungsfeld. Besonders die halbleitenden Charakteristika von SWNTs sowie die Manipulation dieser durch Dotierung bergen viele Möglichkeiten für zukünftige Anwendungen in moderner Elektrotechnologie. Der Weg zu einer industriellen Implementierung von SWNTs in neuartigen optoelektronischen Bauteilen ließe sich durch eine Ausweitung des Wissens bezüglich SWNTs und der dotierungsbasierten Anpassung ihrer Eigenschaften ebnen.
Mit dieser Erkenntniserweiterung als Zielsetzung wurden im Rahmen dieser Dissertation halbleitende, einwandige (6,5)-Kohlenstoffnanoröhren als chiralitätsreine, polymerstabilisierte Proben untersucht. Die ultrakurzzeitaufgelöste Spektroskopie der SWNTs erfolgte an organischen Suspensionen wie auch Dünnschichtfilmen, die je mittels eines gewissen Quantums an Gold(III)-chlorid dotiert worden waren. So konnten die ablaufenden Dynamiken auf einer ps-Zeitskala untersucht werden.
In Kapitel 4 konnte mittels transienter Absorptionsexperimente an redoxchemisch p-dotierter SWNT-Suspensionen zunächst gezeigt werden, dass sich die bei optischer Anregung gebildeten Trionen nicht analog zu Exzitonen diffusiv entlang der Nanoröhre bewegen, sondern lokalisiert vorliegen. Die längere trionischen Zerfallsdauer nach X$_1$- verglichen mit X$_1^+$-resonanter Anregung zeugt außerdem davon, dass das Trion aus dem Exziton gespeist wird. Der Einfluss der Dotierung auf die Zerfallsdynamiken von X$_1$ und X$_1^+$ wurde an SWNT-Dünnschichtfilmen untersucht. Das Photobleichsignal des Exzitons verschiebt hypsochrom und zerfällt schneller mit zunehmender Ladungsträgerdichte durch höherer Gold(III)-chloridkonzentrationen. Dies resultiert aus dem verringerten Abstand zwischen den Ladungsträgern, welche als nichtstrahlende Löschstellen fungieren. Für das X$_1^+$-PB ist ein ähnliches Verhalten zu beobachten. Dabei wird dieses Signal mit weiter steigender Dotierung von einer der H-Bande zuzuordnenden Photoabsorption überlagert. Diese lässt sich in einer starken Sättigung der Dotierung wie auch einer hohen Bandkantenverschiebung begründen.
In Kapitel 5 wurde die Größe der Exzitonen und Trionen in dotierten SWNT-Dünnschichtfilmen mittels des Phasenraumfüllmodells bestimmt. Dabei lag besonderes Augenmerk auf der Kompensation des PB/PA-Überlapps, dem schnellen Zerfall, einem Ausgleich von Differenzen zwischen Anrege- und Absorptionsspektrum sowie dem Anteil intrinsischer/dotierter Nanorohrsegmente, um korrigierte Größen $\xi_\mathrm{k}$ zu erhalten. Für die Trionengröße wurde zusätzlich der Überlapp der Absorptionsbanden einbezogen, um korrigierte Werte $\xi_{\mathrm{T,k}}$ zu bestimmen. $\xi_\mathrm{k}$ beträgt in der intrinsischen Form 6$\pm$2\,nm und bleibt bis zu einer Ladungsträgerdichte $n_{\mathrm{LT}}<0.10$\,nm$^{-1}$ etwa gleich, anschließend ist ein Absinken bis auf etwa 4\,nm bei $n_{\mathrm{LT}}\approx0.20$\,nm$^{-1}$ zu beobachten. Für diesen Trend ist die Überlagerung von Exziton- und H-Bande verantwortlich, da so der Faktor zur Bestimmung des Anteils intrinsischer Nanorohrsegmente an der SWNT verfälscht wird. Die Abweichung der intrinsischen Größe von den in der Literatur berichteten 13$\pm$3\,nm ist möglicherweise auf Unterschiede in der Probenpräparation zurückzuführen. Für die Trionengröße ergibt sich bei steigender Dotierung ein ähnliches Verhalten: Sie beträgt für $n_{\mathrm{LT}}<0.20$\,nm$^{-1}$ 1.83$\pm$0.47\,nm, was in der Größenordnung in guter Übereinstimmung mit der Literatur ist. Für höhere Dotierungen sinkt $\xi_{\mathrm{T,k}}$ bis auf 0.92$\pm$0.26nm ab. Dies erklärt sich dadurch, dass bei höherer $n_{\mathrm{LT}}$ die H-Bande das Spektrum dominiert, sodass der Einfluss der Absorptionsbandenüberlagerung nicht mehr vollständig durch den entsprechenden Korrekturfaktor kompensiert werden kann.
Kapitel 6 beschäftigte sich anstelle redoxchemischer Dotierung der nanoskaligen Halbleiter mit der (spektro-)elektrochemischen Untersuchung von Vorläufern molekularer Radikale. SWV-Messungen weisen dabei darauf hin, dass die Pyrene Pyr1-Pyr3 entsprechend der Anzahl ihrer Substituenten bei Reduktion Mono-, Bi- beziehungsweise Tetraradikale bilden. Die strukturelle Ähnlichkeit der Moleküle äußert sich in gleichen Reduktionspotentialen wie auch ähnlichen potentialabhängigen Absorptionsspektren. Während nur marginale Unterschiede in den PL-Spektren der neutralen und reduzierten Spezies festgestellt werden konnte, lieferte das zeitkorrelierte Einzelphotonenzählen aufschlussreichere Ergebnisse: So wird die Fluoreszenzlebensdauer stark von der Polarität der Umgegbung beeinflusst - bereits die Zugabe des Leitsalzes führt hier zu Änderungen. Die durchschnittliche Fluoreszenzlebensdauer $\tau_{\mathrm{av}}$ sinkt außerdem mit Reduktion und Radikalbildung; für höhere Emissionswellenlängen ist $\tau_{\mathrm{av}}$ außerdem höher. Insgesamt verdeutlichten die Experimente die gute Abschirmung zwischen Pyrenkern und Naphthalimidsubstituenten der Moleküle sowie die Sensibilität gegenüber dem Medium durch TICT, das Vorhandensein von Bi- und Tetraradikalen kann allerdings nicht vollständig belegt werden, wofür EPR-Messugen notwendig wären.
Nanoröhren, die auf dem Element Kohlenstoff basieren, besitzen ein großes Potential für ihre
Anwendung als neuartige und nachhaltige Materialien im Bereich der Optoelektronik und weiteren
zukunftsweisenden Technologiefeldern. Um jedoch hierfür genutzt werden zu können, ist
ein tiefgreifendes Kenntnis über ihre außergewöhnlichen photophysikalischen Eigenschaften notwendig.
Kohlenstoffnanoröhren sind als eindimensionale Halbleiter sehr vielseitige Materialien.
Jedoch ist der Zusammenhang zwischen ihrer Eignung als Halbleiter und der dafür notwendigen
Dotierung nur sehr unzureichend verstanden.
Die Ziele der vorliegenden Dissertation waren deshalb, ein grundlegendes Verständnis der photophysikalischen
Energietransferprozesse in Nanoröhren zu erlangen und den Einfluss von gezielten
Dotierungen auf diese Prozesse im Hinblick auf ihre Eigenschaften als eindimensionale Halbleiter
detailliert zu untersuchen. Die Grundlage für die Experimente bildeten verschiedene Filme
aus einwandigen (6,5)-Kohlenstoffnanoröhren, die durch ein Polyfluoren-Copolymer in einer
organischen Lösungsmittelumgebung isoliert wurden. Mit Hilfe der Ultrakurzzeitspektroskopie
wurden die auf einer schnellen (ps-ns) Zeitskala ablaufenden photophysikalischen Prozesse an
diesen Filmen unter verschiedenen Bedingungen untersucht und analysiert.
In Kapitel 4 wurde der generelle Energietransfer der Kohlenstoffnanoröhren in Polymermatrizen
im Detail studiert. Hierbei wurden durch Simulationen theoretische dreidimensionale
Verteilungen von Kohlenstoffnanoröhren erzeugt und die nach einem Energietransfer vorliegenden
Polarisationsanisotropien berechnet. Verschiedene Berechnungsansätze ergaben, dass die
Nanorohrdichte ϱSWCNT für ein Massenüberschuss X der Matrix nahezu unabhängig von dem
Röhrenvolumen war und durch ϱSWCNT = X−1 · 40 000 μm−1 angenähert werden konnte. Die
Simulationen lieferten von der Röhrendichte abhängige Gaußverteilungen der zwischen den
Nanoröhren vorliegenden Abständen. Aus den Verteilungen konnte weiterhin der Anteil an Röhren
bestimmt werden, die für einen Energietransfer zur Verfügung stehen. Weitere Simulationen
von Nanorohrverteilungen lieferten die Polarisationsanisotropie in Abhängigkeit von der Anzahl
an durchgeführten Energietransferschritten. Die Ergebnisse aus den Simulationen wurden zur
Interpretation der Ultrakurzzeitmessungen angewandt. Hierbei wurden durch die Variation der
Polymermatrix die zwischen den Nanoröhren vorliegenden Abstände verändert und damit die
Art und Intensität des Energietransfers kontrolliert. In Messungen der transienten Anisotropie
zeigte sich, dass ein Exziton nach seiner Erzeugung zwei depolarisierende Energietransferschritte
durchführte. Die Zerfallsdynamiken des Exzitons gaben auch klare Hinweise auf weitere nicht
depolarisierende Energietransferprozesse, die durch parallel zueinander stehende Übergangsdipolmomente
ermöglicht wurden. Eine Erklärung für dieses Verhalten lieferte die faserige
Struktur der Filme, die sich in Aufnahmen durch das Elektronenmikroskop zeigte.
Das Kapitel 5 beschäftigte sich mit dem Aufbau eines transienten Nahinfrarotspektrometers
und den nötigen experimentellen Umbauten zur Messung der transienten Absorption für energiearme
Signale im Spektralbereich unterhalb von 1.4 eV. Hierzu wurde die Weißlichterzeugung
für die Verwendung von Calciumfluorid umgebaut. Das erzeugte Weißlicht wurde in das aufgebaute
Prismenspektrometer eingekoppelt, um es weitestgehend linear auf einer Energieskala zu dispergieren. Auf diese Weise wurden energiearme Spektralkomponenten nicht auf unverhältnismäßig
viele Pixel verteilt und konnten mit ausreichender Intensität detektiert werden. Die
Lichtdetektion erfolgte mittels zweier Detektorzeilen aus Indiumgalliumarsenid, die das transiente
Signal durch eine direkte Referenzierung stabilisierten. Weiterhin wurde in diesem Kapitel die
Justage und die programmierte Ansteuerung des Systems detailliert beschrieben. Hierbei wurde
auf die Justage der Einkopplung per Freistrahl, die Kalibrierung mittels Bandpassspektren
sowie auf die Aufnahme von Weißlichtspektren und transienten Karten detailliert eingegangen.
An Nanorohrdispersionen durchgeführten Testmessungen zeigten, dass das transiente Nahinfrarotspektrometer
mit direkter Signalreferenzierung einwandfrei funktionierte und daher den
beobachtbaren Spektralbereich auf den Bereich von Energien bis unterhalb von 1 eV erweiterte.
Damit ermöglichte der Aufbau einen Zugang zu der Beobachtung größerer Nanorohrchiralitäten
sowie zu der Untersuchung von energiearmen, spektralen Signaturen von Nanorohrdefekten.
In Kapitel 6 wurde das transiente Nahinfrarotspektrometer genutzt, um das zeitabhängige
Verhalten von redoxchemisch p-dotierten Nanoröhren zu charakterisieren und quantitativ zu
beschreiben. Hierzu wurden die spektralen Eigenschaften von SWCNT-Dünnfilmen als Funktion
eines steigenden Dotierungsgrades durch die Messungen der transienten und linearen Absorption
studiert. In der linearen Absorption im Bereich von 0.9 - 2.5 eV vereinfachte sich das Spektrum
mit ansteigender Dotierung stark und verlor vor allem im Bereich des ersten Subbandes
deutlich an Oszillatorstärke. Bei starker Dotierung verschwanden die Signalbeiträge von X1
und der Phononenseitenbande. Weiterhin bleichte auch die bei mittleren Dotierungsgraden
auftauchende Trionenabsorption aus und ging in die breite Absorptionsbande der H-Bande über.
Das Erscheinen und Verschwinden der trionischen sowie exzitonischen Absorption war ebenfalls
in der transienten Absorption durch zeitgleich auftretende/verschwindende Photobleichsignale
zu erkennen. Sowohl der Zerfall des exzitonischen PB-Signals wie auch des Trions beschleunigte
sich mit einer steigenden Dotierung. Die Zerfallszeit des Exzitons im undotierten Film betrug
6.87 ps und verkürzte sich auf 0.732 ps bei höheren Dotierungsgraden. Die Zerfallszeit des Photobleichens
des Trions reduzierte sich von 2.02 ps auf 0.440 ps. Auffallend war hierbei, dass das
Trion im Vergleich zu dem Exziton exponentiell zerfiel und damit auf eine Lokalisierung dieses
Zustandes hinweist. Bei höheren Dotierungsmittelkonzentrationen tauchte in der transienten
Absorption ein neuer Signalbeitrag auf. Die Existenz dieses Signals konnte auf die H-Bande
zurückgeführt werden und könnte auf einer Verschiebung des linearen Absorptionsspektrums
aufgrund einer Renormalisierung der Bandlücke oder der Sättigung von Ladungsträgern beruhen.
Das Signal zeigte eine klare Abhängigkeit vom Dotierungsgrad des Nanorohrfilmes. So wies es
eine hypsochrome Verschiebung auf, wurde spektral breiter und seine Zerfallsdauer reduzierte
sich von 1.62 ps auf 0.520 ps mit steigendem Dotierungsgrad.
Quantenchemische Untersuchungen von Umgebungseinflüssen bei offen- und geschlossenschaligen Systemen
(2022)
In dieser Dissertation werden die Umgebungseinflüsse auf die strukturellen und elektronischen Eigenschaften von verschiedenen offen- und geschlossenschaligen Systemen mittels quantenchemischer Methoden berechnet.
Ein Kernpunkt umfasst die Untersuchung von verdreht angeordneten, biradikalischen Diborylalkenen, welche eine ungesättigte C2R2-Brücke (R = Et, Me) besitzen und durch cyclische (Alkyl)(amino)carbene (CAACs) stabilisiert werden. Quantenchemische Berechnungen zeigen, dass hauptsächlich sterische Effekte für die Ausbildung einer verdrehten Molekülanordnung verantwortlich sind, während bei geringen sterischen Wechselwirkungen (R = H) die Delokalisationseffekte überwiegen, wodurch eine planare Struktur begünstigt wird. Die Bevorzugung einer offenschaligen Singulettkonfiguration anstelle eines Tripletts ist auf den großen Energieunterschied der beiden einfach besetzten Molekülorbitale zurückzuführen. Durch die Berechnung der Lösungsmitteleffekte mithilfe von polarisierbaren Kontinuumsmodellen kann gefolgert werden, dass mit zunehmender statischer Dielektrizitätskonstante eine planare und geschlossenschalige Struktur stärker stabilisiert wird als eine verdrehte Anordnung.
Ein weiteres Thema dieser Dissertation befasst sich mit der quantenchemischen Analyse eines makrozyklischen Perylenbisimid-Trimersystems, welches eingebettet in einer Polymethylmethacrylat-Matrix bei Temperaturen nahe dem absoluten Nullpunkt eine Lokalisierung der ersten drei angeregten Zustände zeigt. Quantenchemische Vakuumberechnungen ergeben, dass unabhängig von der gegenseitigen geometrischen Orientierung der drei Perylenbisimid-Chromophore der Übergang vom Grundzustand in den S1-Zustand verboten ist und dass die ersten drei angeregten Zustände delokalisiert vorliegen. Mithilfe von expliziten Lösungsmittelmodellen kann jedoch gezeigt werden, dass das Auftreten dieser Lokalisierungen auf eine inhomogene Polymethylmethacrylat-Umgebung zurückzuführen ist, die zu einem Symmetriebruch und somit zu einer Zunahme der Oszillatorstärke für S1 und der Lokalisierungsgrade für S1, S2 und S3 führt.
Darüber hinaus wird der Lösungsmitteleinfluss auf die angeregten Zustände des Azulens mittels impliziter und expliziter Lösungsmittelmodelle berechnet. Bei einer Erhöhung der dynamischen Dielektrizitätskonstante im impliziten Modell nehmen die Anregungsenergien der vertikalen Singulettzustände ab, wobei der Effekt mit steigender Oszillatorstärke zunimmt. Die Auswirkung der statischen Dielektrizitätskonstante auf die Anregungsenergien ist dagegen deutlich schwächer ausgeprägt. Im expliziten Modell bewirkt das Lösungsmittel ebenfalls eine Abnahme der Anregungsenergie des hellen Singulettzustands, wenn auch in geringerem Umfang als im impliziten Modell.
Als letztes Thema wird der Inhibitionsmechanismus der Cysteinprotease Rhodesain durch zwei modifizierte 1,4-Naphthoquinone untersucht. Während beide Naphthoquinone an der 2-Position eine Dipeptideinheit aufweisen, besitzen sie an der 3-Position entweder einen Nitril- oder Chloridsubstituenten. Zwar erfolgt bei beiden Derivaten die Inhibition über einen kovalent-reversiblen Mechanismus, jedoch verläuft die Hemmung im Falle des Nitrilderivats erheblich effektiver. Die quantenchemischen Berechnungen eines vereinfachten Modells zeigen, dass die Cysteineinheit (HS-R) bevorzugt in einer exothermen und reversiblen Additionsreaktion an die elektronenarme C-C-Doppelbindung der Naphthoquinone anlagert. Dabei werden kleinere Reaktionsenergien für die Reaktion des Chlorderivats als für die Reaktion des Nitrilderivats erhalten. Durch die Berücksichtigung von Wasser in einem impliziten Lösungsmittelmodell kommt es bei fast allen Reaktionsprodukten zu einer Energiezunahme, die bei der Reaktion des Nitrilderivats stärker ausfällt als bei der Reaktion des Chlorderivats.
Im Rahmen dieser Arbeit werden unterschiedliche Aspekte der korrelierten Elektronen-Kerndynamik, anhand verschiedener Modellsysteme untersucht. Dabei wird vor allem auf den Vergleich numerisch exakter und approximativer Methoden zur Beschreibung der Wellenpaketdynamik eingegangen, wobei bei letzterem das Augenmerk auf der Born-Oppenheimer (BO) Näherung liegt. Die verwendeten Modellsysteme erlauben es, die gekoppelte Elektronen-Kern-Dynamik exakt zu beschreiben. Die daraus gewonnenen Ergebnisse dienen als Referenz für den Vergleich mit den Näherungsmethoden.
Im ersten Teil der Arbeit wird die Dynamik eines Wellenpakets in der Umgebung einer Konischen-Durchschneidung (CI) untersucht, wobei die Beschreibung des Wellenpakets quantenmechanisch und durch die klassische Mechanik im Phasenraum erfolgt.
Im zweiten Teil wird die Wahrscheinlichkeitsflussdichte untersucht. Zuerst wird ein Fall konstruiert, in welchem die Bewegung im elektronischen Grundzustand stattfindet, sodass die Bedingungen der BO Näherung erfüllt sind. Dabei wird vor allem auf das Verschwinden der elektronischen Wahrscheinlichkeitsflussdichte innerhalb der BO Näherung eingegangen. Im weiteren Verlauf werden die Flussdichten in der Umgebung einer CI untersucht, wobei unterschiedliche Situationen modelliert werden.
Im dritten Teil wird die Berechnung des elektronischen Impulserwartungswerts innerhalb der BO Näherung untersucht. Dieser verschwindet innerhalb der BO Näherung, wenn man diesen direkt berechnet (Geschwindigkeitsform), während man über das Ehrenfest Theorem (Längenform) sehr gute Werte erhält.
Im vierten Teil wird eine neue Flussdichte, die Translationsflussdichte, vorgestellt. Diese ergibt sich aus der Überlegung, dass die Geschwindigkeitsform des Impulserwartungswerts durch die Wahrscheinlichkeitsflussdichte ausgedrückt werden kann. Demnach muss auch die Längenform einer Flussdichte entsprechen und man erhält die Translationsflussdichte.
Die ersten Beispiele für Lewis-Basen-Addukte des Stammboraphosphaketens H\(_{2}\)B-PCO und ihre cyclischen Dimere wurden hergestellt. Eines dieser Addukte zeigt unter milden Bedingungen eine Decarbonylierung und anschließende Insertion des Phosphinidens in die B-C-Bindung eines Borols, was in der Bildung sehr seltener Beispiele für 1,2-Phosphaborinine, B,P-Isostere von Benzol, resultiert. Die starken Donoreigenschaften dieser 1,2-Phosphaborinine wurden durch die Synthese ihrer π-Komplexe mit Metallen der Gruppe 6 bestätigt.
Bei der Einelektronenreduktion eines durch eine cyclisches (Alkyl)(amino)carben (CAAC) stabilisierten Arylborylen-Carbonylkomplexes erfolgt die Bildung eines dimeren Borylketyl-Radikalanions, bedingt durch eine intramolekulare Arylmigration zum CO Kohlenstoffatom. Computergestützte Analyse liefert Hinweise auf eine radikalanionische [(CAAC)B(CO)Ar]\(^{.-}\) Zwischenstufe. Weiterführende Reduktion des entstandenen Komplexes liefert ein hoch nukleophiles (Boranyliden)methanolat.
Diese Dissertation beschäftigt sich mit der Exzitonendynamik molekularer Aggregate, die nach Mehrphotonen-Anregung auf ultrakurzer Zeitskala stattfindet. Hierbei liegt der Fokus auf der Charakterisierung der Exziton-Exziton-Annihilierung (EEA) mithilfe von zweidimensionaler optischer Spektroskopie fünfter Ordnung. Dazu werden zwei verschiedene Modellsysteme implementiert: Das elektronische Homodimer und das elektronische Homotrimer-Modell, wobei Letzteres eine Erweiterung des Dimer-Modells darstellt. Die Kopplung des quantenmechanischen Systems an die Umgebung wird mithilfe des Quantum-Jump-Ansatzes umgesetzt. Besonderes Interesse kommt der Analyse des Signals fünfter Ordnung in Abhängigkeit der Populationszeit T zu.
Anhand des Dimer-Modells als kleinstmögliches Aggregat lassen sich bereits gute Vorhersagen auch über das Verhalten größerer molekularer Aggregate treffen. Der Zerfall des oszillierenden Signals für lange Populationszeiten korreliert mit der EEA. Dies zeigt, dass die zweidimensionale optische Spektroskopie genutzt werden kann, um den Annihilierungsprozess zu charakterisieren. Innerhalb des Modells des Dimers wird weiterhin der Einfluss der Intraband-Relaxation untersucht. Zunehmende Intraband-Relaxation verhindert den Austausch zwischen den lokalen Zuständen, der essentiell für den Annihilierungsprozess ist, und die EEA wird blockiert.
Das elektronische Trimer-Modell erweitert das Dimer-Modell um eine Monomereinheit. Somit befinden sich die Exzitonen im Anschluss an die Anregung nicht mehr unvermeidlich nebeneinander. Es gibt somit eine Konfiguration, bei der sich die Exzitonen zunächst zueinander bewegen müssen, bevor die Startbedingung des Annihilierungsprozesses gegeben ist. Dieser zusätzliche Schritt wird auch Exzitonendiffusion genannt. Die Ergebnisse dieser Arbeit legen nahe, dass das erwartete Verhalten nur zu sehr kurzen Zeiten im Femtosekundenbereich auftritt und somit die Zeitskala der Exzitonendiffusion im Falle des Trimers nicht sichtbar wird. Es bedarf demnach eines größeren Modellsystems, bei dem sich der Effekt der zeitverzögert eintretenden EEA deutlich in der Zerfallsdynamik manifestieren kann.
Das erste Ziel der vorliegenden Dissertation bestand darin, ein bereits bestehendes TOF-MS-Setup dahingehend zu erweitern, um damit Velocity Map Imaging-Experimente durchführen zu können. Dies erforderte zunächst die Konzipierung und Programmierung einiger für die Datenaufnahme, -verarbeitung und -analyse benötigter LabView-Anwendungen. Anschließend konnten erste Kalibrierexperimente an Methyliodid, in denen wichtige experimentelle Parameter identifiziert und optimiert wurden, durchgeführt werden. Außerdem gelang es dadurch, die Messgenauigkeit des Setups auf 0.7 % und dessen Auflösungsvermögen auf 4.4 % zu bestimmen, was im Bereich für VMI-Apparaturen typischer Werte liegt. Zur weiteren Überprüfung der Funktionstüchtigkeit des Setups wurde in ersten zeitaufgelösten Experimenten im Folgenden die Desaktivierung des S1-Zustands von Pyridin untersucht. Neben der Reproduktion einiger bereits literaturbekannter Resultate konnten dabei zusätzlich die im Multiphotonen-Ionisationsschritt populierten Rydberg-Zustände identifiziert werden. Anschließend wurde mit Experimenten an bisher weniger gut untersuchten organischen Aromaten und Heteroaromaten fortgefahren. Das Ziel dieser Studien lag in der Aufklärung der photoinduzierten Dynamiken der Verbindungen, wobei das zur Verfügung stehende ps-Lasersystem die Möglichkeit bot, die Desaktivierung elektronisch angeregter Zustände gezielt in Abhängigkeit von deren Schwingungsenergie zu untersuchen. Der darin bestehende Vorteil zeigte sich vor allem in Studien an Tolan und Phenanthridin, deren erste angeregte, optisch aktive Zustände am Origin Lebensdauern im ns-Bereich aufweisen, die sich mit zunehmender vibronischer Anregung jedoch auf bis zu 10 ps verringern. Als Grund dafür konnten nichtstrahlende Desaktivierungsprozesse, für deren Eintreten eine energetische Barriere überwunden werden muss, identifiziert werden. Während in Tolan nach Photoanregung ein Übergang in einen (πσ∗)-Zustand, der zur Ausbildung einer trans-bent-Struktur führt, erfolgt, ist im Falle von Phenanthridin vermutlich ein El-Sayed-erlaubter ISC-Übergang in einen 3(nπ∗)-Zustand für die drastische Verkürzung der S1-Lebensdauer verantwortlich. Ein solcher konnte weder im zu Phenanthridin isomerischen Benzo[h]quinolin, noch in dessen PAH-Muttermolekül Phenanthren beobachtet werden, was auf die höhere energetische Lage bzw. die Abwesenheit des mittels ISC populierten 3(nπ∗)-Zustands in diesen Molekülen zurückgeführt werden kann. In weiteren im Rahmen der vorliegenden Arbeit durchgeführten Experimente wurden zudem die aromatischen Moleküle Acenaphthylen und 4-(Dimethylamino)benzethin (DMABE) untersucht. Zeitaufgelöste Studien zeigten dabei, dass die Desaktivierung der S2-Zustände beider Moleküle auf der sub-ps-Zeitskala stattfindet und mit dem vorhandenen Lasersystem daher nicht aufgelöst werden kann. In Acenaphthylen erfolgt die S2-Relaxation größtenteils über einen sequentiellen IC-Mechanismus, innerhalb dem der S1-Zustand des Moleküls intermediär besetzt wird. Dessen Lebensdauer konnte am Origin auf 380 ps bestimmt werden, fällt mit steigender Schwingungsanregung jedoch auf bis zu 55 ps ab. Für die Desaktivierung des S2-Zustands von DMABE konnte hingegen ein paralleles Relaxationsmodell, in dem neben dem S1-Zustand ein weiterer elektronisch angeregter Zustand populiert wird, nachgewiesen werden. Bei diesem könnte es sich möglicherweise um einen (πσ∗)-Zustand, dessen Besetzung die Ausbildung einer trans-bent-Geometrie innerhalb der Acetylen-Einheit des Moleküls zur Folge hat, handeln. Einen weiteren großen Teil der vorliegenden Dissertation nahmen Experimente an van-der-Waals-gebundenen Clustersystemen ein. Im Fokus der Studien standen dabei Moleküle mit ausgedehnten aromatischen π-Systemen, da solche eine hohe Relevanz für verschiedene materialwissenschaftliche Forschungsgebiete besitzen. Ein Beispiel hierfür ist Tetracen, welches als Modellsystem für die Untersuchung von Singlet Fission-Prozessen angesehen wird. In Kombination mit nichtadiabatischen Surface-Hopping-Simulationen zeigten Experimente an Tetracen-Dimeren, dass nach deren S2-Anregung zunächst ein schneller S1←S2-Übergang (τ < 1 ps), gefolgt von der Ausbildung einer Excimerstruktur, stattfindet. Letztere erfolgt mit einer Zeitkonstante von 62 ps und führt zu einem Anstieg des transienten Ionensignals, wohingegen die Desaktivierung des Excimer-Zustands von einem abklingenden Signalbeitrag mit τ = 123 ps repräsentiert wird. Wenngleich über die weitere Relaxation der Excimerspezies zum gegenwärtigen Zeitpunkt keine Aussage getroffen werden kann, besteht damit die Möglichkeit, dass Excimer-Zustände als Zwischenstufe im SF-Mechanismus isolierter Tetracen-Dimere auftreten. In zeitaufgelösten Experimenten an Phenanthren-Dimeren konnte ebenfalls ein Anstieg des transienten Signals mit einer vergleichbaren Zeitkonstante von τ = 86 ps, der jedoch auf einem konstanten Signaloffset endet, gefunden werden. Dies deutet darauf hin, dass auch Phenanthren-Dimere in der Lage sind, Excimerstrukturen, die im Gegensatz zu denen des Tetracens jedoch deutlich langlebiger sind, auszubilden. Studien an den Dimerspezies der Azaphenanthrene Benzo[h]quinolin und Phenanthridin offenbarten hingegen etwas schnellere Relaxationen mit Zeitkonstanten von 15 bzw. 40 ps. Zudem zeigten beide Spezies eine stark ausgeprägte Fragmentation, sodass für deren Untersuchung auf die VMI-Detektionsmethode zurückgegriffen werden musste. Dadurch wurde deutlich, dass sich Photoionen-Imaging-Experimente hervorragend für Studien an schwach gebundenen Clustersystemen eignen, da diese die Separation verschiedener Signalbeiträge innerhalb eines betrachteten Massenkanals ermöglichen.