Back to Search Start Over

Implementierung neuer Reaktionspfad bestimmender Methoden und Untersuchung von Lösungsmitteleinflüssen auf die Natur der angeregten Zustände perylen-basierter Farbstoffe

Authors :
Bellinger, Daniel
Publication Year :
2016

Abstract

Two thematic complexes were addressed within this work. One part is related to improvements and new implementations into the CAST program package. Thereby the main focus laid on the delivery of a tool which can be used to characterize complex reactions and their mechanisms. But also within the new force field (FF) method (SAPT-FF) within the CAST program, several improvements were made. The second topic is related to the description of dye molecules and their spectral properties. The main focus within these studies was set on the influence of the environment on these properties. In the first topic improvements of the local acting NEB (nudged elastic band) methods were included and the number of available methods was extended. The initial pathway generation was improved by implementing the IDPP (image dependent pair potential) method and a new method was implemented for describing temperature dependent pathways. Additionally, improvements have been made to the optimization routines (global NEB). As a second part the Pathopt (PO) method was considerably improved. In the beginning of the work the original PO idea was used. In this approach one starts with a global optimization on one n-1 dimensional hyperplane which divides the reaction into two sub-areas for obtaining guesses of TSs (transition states). These found TS guesses were used to optimize to the ”true” TS. Starting from the optimized ones a relaxation to the next connected minima is done. This idea has been automatically implemented and extended to several number of hyperplanes. In this manner a group of pathsegments is obtained which needs to be connected, but within this work it was realized that such a procedure might be not very efficient. Therefore, a new strategy was implemented which is founded on the same constrained global optimization scheme (MCM) for which the user defines the number of hyperplanes generated. The number of such generated hyperplanes should be large enough 134 to describe the space between the concerning reactants in a sufficient way. The found minima are directly used to built up the reaction pathway. For this purpose a RMSD (root mean square deviation) criterion is used to walk along ways of minimal change from one to another hyperplane. To prove the implementations various test calculations were carried out and extensions included to prove the capabilities of the new strategy. Related to these tests a new strategy for applying the move steps in MCM (Monte Carlo with minimization) was realized which is also related to the question of the coordinates representation. We were able to show that the hopping steps in MCM can be improved by applying Cartesian steps in combination of random dihedral moves with respect to the constraint. In this way it was possible to show that a large variety of systems can be treated. An additional chapter shows the improvements of the SAPT-FF implementation and related test cases. It was possible to treat benzene dimer and cluster systems of different sizes consistently also in accordance with high level ab initio based approaches. Furthermore, we showed that the SAPT-FF with the right parameters outperforms the standard AMOEBA implementation which is the basis of the SAPT-FF implementation. In the last three chapters deal with the description of perlyene-based dyes. In the first smaller chapter ground state chemistry description of macro cycles of PBI (perylene bisimide) derivatives were investigated. Therefore, AFM (atomic force microscopy) based pictures were explained within our study. The methods to explain aggregation behavior in dependency of the ring size were MD simulations and configuration studies. The last two chapters deal with opto-electronic or photo-physical properties of PBI and PTCDA (perylene-3,4,9,10-tetracarboxylic dianhydride). In detail, we investigated the role of the environment and the aggregate or crystal surrounding by applying different models. In that way implicit and explicit solvation models, the size of aggregates and vibration motions were used. In the case of PBI the recent work is found on preliminary studies related to my bachelor thesis and extends it. It was shown that the direct influence of a polarizable surrounding, as well as explicit inclusion of solvent molecules on the overall description of the excitations and nature of the excited states is weaker as one might expect. However the inclusion of intra-molecular degrees of freedom showed a stronger influence on the state characteristics and can induce a change of the order of states within the dimer picture. For the PTCDA molecule the main focus was set on the description of the absorption spectrum of crystalline thin films. Related to this older works exist which already gave a description and assignment of the absorption band, but are based on different approaches compared to the one used in this work. We used the supermolecule ansatz, whereas the environment and different aggregate sizes were investigated. Within the dimer based approach we were able to show that using continuum solvation (IEFPCM/COSMO) based description for the environment the relative order of states remains unchanged. Similar to the PBI calculations the influence of the vibrational motions /distortions is larger. The simulation of the crystal environment by using QM/MM (quantum mechanics/molecular mechanics) approaches delivered that an asymmetric charge distribution might induce a localization of the excitation and a stronger mixing of states. For obtaining further insights we go beyond the dimer picture and aggregates of different sizes were used, whereas the simulations up to the octadecamer mono- and even dual-layer stack were carried out. Within these calculations it was shown that the H-coupling is dominating over a weaker J-coupling between different stacks. Additionally the calculations based on DFT (density functional theory) and semi-empirics showed that the lowest state in terms of energy are mostly of Frenkel type, whereas the higher lying states are CT ones which mix with embedded Frenkel type states. The first band of the absorption spectrum was explained by inclusion of vibrational motions within the stacks which induce an intensity gain of the first excited state. This intensity was not explainable by using the undistorted stacks. Also relaxations at the crystal surface might play a role, but are experimentally not explainable.<br />In der Arbeit wurden zwei große Themenkomplexe bearbeitet. Zum einen wurden Verbesserungen und neue Methoden in CAST, unserem Entwicklungstool, implementiert. Hierbei geht es vor allem darum ein Werkzeug bereit zu stellen, mit dem es möglich ist Reaktionen genauer zu charakterisieren. Aber auch neue Beschreibungen innerhalb der Kraftfeldmethoden (SAPT-FF) wurden bereitgestellt bzw. erweitert. Der zweite Themenkomplex behandelt die Beschreibung von Farbstoffmolekülen und ihrer spektralen Eigenschaften. Insbesondere liegt in dieser Studie der Fokus auf Umgebungseinflüsse. Im ersten Abschnitt wurden Erweiterungen in den lokalen Methoden, die auf NEB (nudged elastic band) basieren, implementiert. Hier wurde zum einen das Spektrum an Methoden erweitert. So wurde der initiale Vorschlag für den Startpfad durch Implementierung der IDPP (image dependent pair potential) Methode verbessert. Des Weiteren wurde eine Methode zur temperaturabhängigen NEB Beschreibung integriert, die auf Maximierung des Fluxes beruht. Auch wurden Verbesserungen hinsichtlich der Optimierungsroutinen durchgeführt. Der wesentliche Teil im ersten Themenbereich beschäftigt sich mit der Verbesserung und Automatisierung von Pathopt (PO). Zu Beginn der Arbeit wurde die ursprüngliche Idee aufgegriffen. Hierbei ermittelt man Vorschläge für Übergangszustände (¨UZ) durch eine globale Optimierung mit Nebendbedingungen auf einer n-1 dimensionalen Hyperfläche, die den Reaktionsraum teilt. Diese ¨UZ bilden den Startpunkt, um mittels einer ”uphill” Optimierung hin zum ”wirklichen ÜZ” zu gelangen. Ausgehend von diesen wurde in die nächst verknüpften Minima relaxiert. Diese Idee wurde automatisiert und auf mehrere Hyperflächen ausgeweitet. So erhält man eine Schar an Pfadsegmenten, die verknüpft werden müssen. Im Laufe der Arbeit, stellte sich jedoch heraus, dass diese Vorgehensweise nicht sehr effizient ist und daher wurde eine neue Idee verwirklicht. Diese beruht wiederum auf der globalen Optimierung mittels Monte Carlo mit Minimierung und Nebenbedingungen auf einer vom Nutzer bestimmten Anzahl an n-1 dimensionalen Hyperflächen. Nun wählt man diese Anzahl entsprechend groß genug aus, um den Raum zwischen den Reaktanden zu beschreiben. Die so gefundenen Mininima auf den n-1 Hyperflächen werden für die direkte Pfaderzeugung genutzt. Dies geschieht mittels eines RMSD (root mean square deviation) Kriteriums, um so den Weg der geringsten Änderungen anhand der Hyperflächen zu wählen. Im Zuge der Implementierung der Methode wurden zahlreiche Testrechnungen und Methodenerweiterungen durchgeführt, um die Funktionalität zu überpüfen und zu verbessern. Diese Verbesserungen liegen zum Bsp. in den Sprungstrategien bzw. der Wahl des Koordinatensystems. Hier konnte gezeigt werden, dass eine Verbindung unterschiedlicher Strategien für die Durchführung des ”Hüpfens” in Monte Carlo zu entscheidenden Verbesserungen führt. Diese Verbesserung besteht in der Verknüpfung von Kartesischen Schritten und zufälliger Veränderungen der Diederwinkel im Rahmen der Nebenbedingungen. Mit Hilfe dieser Verbesserungen konnte eine Vielzahl von Systemen behandelt werden. Ein weiteres Kapitel beschreibt Verbesserungen zum SAPT-FF (FF=Kraftfeld). Testrechnungen zu strukturellen Eigenschaften von Benzol Clustern belegen die Genauigkeit der Ansätze. Auch wurde aufgezeigt, dass das SAPT-verbesserte AMOEBA Kraftfeld der Standard Parametrisierung überlegen ist. Die letzten drei Abschnitte dieser Arbeit behandeln Perylen-basierte Farbstoffe. In einem ersten kleinen Kapitel geht es um die Grundzustandseigenschaften von PBI (Perylenbisimide) Makrozyklen und Erklärung von AFM (Atomic Force Mycroscopy) Messungen. Hier konnten wir mittels MD-Simluation (Molekular Dynamik) und deren Analyse, sowie Beschreibungen unterschiedlicher Konfigurationen, das Aggregationsverhalten in Abhängigkeit der Ringgröße genauer beleuchten. In den beiden letzten Kapiteln geht es um die optoelektronischen Eigenschaften bzw. die photophysikalische Beschreibung von PBI und PTCDA (Perylen-3,4,9,10-Tetracarboxyl Dianhydrid). Im Genaueren wurde die Rolle der Umgebung in Aggregat und Kristall durch unterschiedliche methodische Ansätze untersucht. So wurden implizite Solvensmodelle und explizite Solvatation, Aggregatgröße und vibronische Freiheitsgrade untersucht. In den Arbeiten zum PBI konnte gezeigt werden, dass ein direkter Einfluss durch die Beschreibung mittels impliziter Solvatation, als auch expliziter Solvensmoleküle, auf die Lage der Zustände auch in Hinsicht auf deren Charakterisik nicht auftritt. Berücksichtigt man intra-molekulare Freiheitsgrade, so wird die Lage der Zustände deutlich stärker beeinflusst und sogar ein Wechsel der Zustände wird induziert. Im Fall von PTCDA lag vor allem die Beschreibung und Erklärung der Absorptionsspektren von kristallinem PTCDA im Fokus. Hierzu gibt es ¨altere Arbeiten, die bestimmte Zuordnungen der Banden und ihrer Übergänge postuliert haben. In dieser Arbeit sollte diese Beschreibung im Rahmen eines Supermolekül Ansatzes geklärt und weiter beschrieben werden, wobei Umgebungseinflüsse und auch Eigenschaften verschiedener Aggregate untersucht wurden. Im Dimer Bild konnten wir zeigen, dass die Umgebung, beschrieben durch Continuums Ansätze (IEFPCM/COSMO) die Lage der Zustände nicht beeinflusst und im Wesentlichen nur Zustände mit großer Oszillatorstärke stabilisiert werden. Ähnlich wie im Falle des PBIs hat die Berücksichtigung vibronischer Freiheitsgrade einen wesentlich größeren Einguss. Die Simulation der Kristallumgebung durch QM/MM-Ansätze (Quantenmechanik/Molekularmechanik) ergab, dass eine asymmetrische Ladungsverteilung zu einer Lokalisierung der Anregung und einem stärkeren Durchmischen der Zustände fuhrt. Für eine noch weitergehende Beschreibung wurde das Dimer Bild verlassen und unterschiedliche Aggregate, bis hin zum Oktadekamer im mono- bzw. zweifach-Lagen-Aggregat untersucht. Hier konnte gezeigt werden, dass die Kopplung im H-Aggregat die dominierende Rolle einnimmt und die J-Aggregat Kopplung vernachlässigt werden kann. Zudem zeigen die Rechnungen, die mittels DFT (Dichtefunktionaltheorie) und semi-empirischen Ansätzen durchgeführt wurden, dass die energetisch niedriger liegenden Zustände im wesentlichen Frenkel Charakter aufweisen während die energetisch höher liegenden Zustände CT (Charge Transfer) Charakter haben. Das Auftreten der ersten Bande im Absorptionsspektrum wurde zudem auf das Vorhandensein von möglichen Schwingungsanregungen (mehrere Moden) zurückgeführt, da diese zu einer Zunahme an Intensität des ersten angeregten Zustandes führen, die ohne Berücksichtigung nicht in diesem Maße erhalten wird. Auch könnten Oberflächenrelaxationen eine Rolle spielen, wobei diese experimentell nicht beobachtbar sind.

Details

Language :
English
Database :
OpenAIRE
Accession number :
edsair.od.......713..28e4544cce3a764b64702df4dee900bd