TY - THES A1 - Vogt, Matthias Guido T1 - Elektronische Eigenschaften von Wabengittern mit starker Spin-Bahn-Kopplung T1 - Electronic Properties of honeycomb lattices with strong spin-orbit coupling N2 - Im Rahmen dieser Arbeit wurden die elektronischen Eigenschaften von Graphen auf Metalloberflächen mittels Rastertunnelmikroskopie und Quasiteilcheninterferenz (englisch quasiparticle interference, QPI)-Messungen untersucht. Durch das Verwenden schwerer Substrate sollte die Spin-Bahn-Wechselwirkung des Graphen verstärkt werden und damit eine Bandlücke am K-Punkt der Bandstruktur mittels QPI beobachtet werden. Um das Messen von QPI auf Graphen zu testen, wurde auf der Oberfläche eines SiC(0001)-Kristalls durch Erhitzen Graphen erzeugt und mit dem Rastertunnelmikroskop untersucht. Dieses System wurde schon ausführlich in der Literatur beschrieben und bereits bekannte QPI-Messungen von Streuringen, die auf den Dirac-Kegeln des Graphen am K-Punkt basieren, konnte ich auf gr/SiC(0001) in guter Qualität erfolgreich reproduzieren. Anschließend wurde Graphen nach einem wohlbekannten Verfahren durch Aufbringen von Ethylen auf ein erhitztes Ir(111)-Substrat erzeugt. Dieses gr/Ir(111)-System diente auch als Grundlage für Interkalationsversuche von Bismut (gr/Bi/Ir(111)) und Gadolinium (gr/Gd/Ir(111)) zwischen das Graphen und das Substrat. Auf gr/Bi/Ir(111) wurde ein schon aus der Literatur bekanntes Netzwerk aus Versetzungslinien beobachtet, dem zusätzlich eine Temperaturabhängigkeit nachgewiesen werden konnte. Beim Versuch, Gadolinium zu interkalieren, wurden zwei verschieden Oberflächenstrukturen beobachtet, die auf eine unterschiedlich Anordnung bzw. Menge des interkalierten Gadoliniums zurückzuführen sein könnten. Auf keinem dieser drei Systeme konnten allerdings Streuringe mittels QPI beobachtet werden. Als Vorbereitung der Interkalation von Gadolinium wurden dessen Wachstum und magnetische Eigenschaften auf einem W(110)-Kristall untersucht. Dabei konnte eine aus der Literatur bekannte temperaturabhängige Austauschaufspaltung reproduziert werden. Darüber hinaus konnten sechs verschieden magnetische Domänen beobachtet werden. Zusätzlich sind auf der Oberfläche magnetische Streifen auszumachen, die möglicherweise auf einer Spinspirale basieren. Als Grundlage für die mögliche zukünftige Erzeugung Graphen-artiger Molekülgitter wurde das Wachstum von H-TBTQ und Me-TBTQ auf Ag(111) untersucht. Die Moleküle richten sich dabei nach der Oberflächenstruktur des Silber aus und bilden längliche Inseln, deren Kanten in drei Vorzugsrichtungen verlaufen. Auf H-TBTQ wurde zudem eine zweite, Windmühlen-artige Ausrichtung der Moleküle auf der Oberfläche beobachtet. Auf den mit den Molekülen bedeckten Stellen der Oberfläche wurde eine Verschiebung des Ag-Oberflächenzustands beobachtet, die mit einem Ladungstransfer vom Ag(111)-Substrat auf die TBTQ-Moleküle zu erklären sein könnte. N2 - In this thesis, the electronic properties of graphene on metal surfaces were investigated by scanning tunneling microscopy and quasiparticle interference (QPI) measurements. In order to enhance the spin orbital interaction of the graphene and possibly observe a band gap at the K-point of the band structure via QPI, substrates with heavy atoms were used. To test the ability to measure QPI on graphene, graphene was produced on the surface of a SiC(0001) crystal by heating and examined with a scanning tunneling microscope. This system has already been described in detail in the literature and I was able to successfully reproduce QPI measurements of clearly recognizable scattering rings, which are due to the Dirac cones of the graphene at the K-point Afterwards, graphene was produced by a well-known process by applying ethylene to a heated Ir(111) substrate. This gr/Ir(111) system also served as a basis for intercalation experiments of bismuth (gr/Bi/Ir(111)) and gadolinium (gr/Gd/Ir(111)) between the graphene and the substrate. On gr/Bi/Ir(111), a network of dislocation lines known from literature was observed, which also showed a temperature dependence. In the attempt to intercalate gadolinium, two different surface structures were observed which could be due to a different arrangement or quantity of the intercalated gadolinium. However, on none of these three systems scattering rings were observed by QPI. In preparation for the intercalation of gadolinium, its growth and magnetic properties were investigated on a W(110) substrate. A temperature-dependent exchange splitting of the surface density of states known from the literature could be reproduced. In addition, six different magnetic domains and magnetic stripes were observed on the surface, which may be based on a spin spiral. The growth of H-TBTQ and Me-TBTQ on Ag(111) was investigated as a basis for a possible subsequent generation of graphene-like molecular lattices in the future. The molecules are aligned to the surface structure of the silver and form elongated islands with edges in three preferred directions. H-TBTQ also appeared in a second, windmilllike orientation of the molecules on the surface. A shift of the Ag surface state was observed on the surface areas covered by the molecules, which might be explained by a charge transfer from the Ag(111) substrate to the TBTQ molecules. KW - Spin-Bahn-Wechselwirkung KW - Graphen KW - Rastertunnelmikroskopie KW - Wabengitter KW - Tribenzotriquinacen KW - Quasiteilcheninterferenz Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-207506 ER - TY - THES A1 - Brückner, Charlotte T1 - The Electronic Structure and Optoelectronic Processes at the Interfaces in Organic Solar Cells Composed of Small Organic Molecules - A Computational Analysis of Molecular, Intermolecular, and Aggregate Aspects T1 - Die elektronische Struktur und die optoelektronischen Prozesse an den Grenzflächen in organischen Solarzellen aus kleinen organischen Molekülen - eine theoretische Analyse auf molekularer, intermolekularer und Aggregatebene N2 - Describing the light-to-energy conversion in OSCs requires a multiscale understanding of the involved optoelectronic processes, i.e., an understanding from the molecular, intermolecular, and aggregate perspective. This thesis presents such a multiscale description to provide insight into the processes in the vicinity of the organic::organic interface, which are crucial for the overall performance of OSCs. Light absorption, exciton diffusion, photoinduced charge transfer at the donor-acceptor interface, and charge separation are included. In order to establish structure-property relationships, a variety of different molecular p-type semiconductors are combined at the organic donor-acceptor heterojunction with fullerene C60, one of the most common acceptors in OSCs. Starting with a comprehensive analysis of the accuracy of diverse ab initio, DFT, and semiempiric methods for the properties of the individual molecules, the intermolecular, and aggregate/device stage are subsequently addressed. At all stages, both methodological concepts and physical aspects in OSCs are discussed to extend the microscopic understanding of the charge generation processes. N2 - Um die Umwandlung von Licht zu Strom in organischen Solarzellen zu verstehen, müssen die beteiligten optoelektronischen Prozesse sowohl auf molekularem als auch auf intermolekularem und auf dem Aggregatniveau beschrieben werden. Diese Arbeit stellt eine solche mehrstufige Beschreibung dar, um zum grundlegenden Verständnis derjenigen Prozesse am organisch::organischen Interface beizutragen, die für die finale Gesamtleistung der Zelle ausschlaggebend sind. Dabei werden die wesentlichen Schritte von der Lichtabsorption und Exzitonendiffusion über den photoinduzierte Charge-Transfer-Schritt am Donor-Akzeptor-Interface bis hin zur endgültigen Ladungstrennung berücksichtigt. Um auf Struktur-Eigenschafts-Beziehungen rückschließen zu können, wurden verschiedene molekulare p-Halbleiter in der heterojunction mit Fulleren C60 kombiniert, einem der gängigsten Akzeptoren in organischen Solarzellen. Nach einer umfangreichen Bewertung der Eignung verschiedener ab initio und semiempirischer Methoden sowie diverser DFT-Funktionale für die Beschreibung der molekularen Eigenschaften wurden intermolekulare und Aggregataspekte diskutiert. Auf allen Ebenen, d.h. auf der molekularen, intermolekularen und auf der Aggregatebene, stehen sowohl methodische Ansätze als auch grundlegendende photophysikalische Überlegungen im Mittelpunkt, um das mikroskopische Verständnis der Ladungsträgererzeugung in organischen Solarzellen zu erweitern. KW - Benchmark KW - Solarzelle KW - organic interfaces KW - benchmark KW - charge carrier generation KW - organische Grenzflächen KW - Benchmark KW - Ladungsträgererzeugung KW - solar cell Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-141652 ER - TY - THES A1 - Heilos, Anna T1 - Mechanistic Insights into the Inhibition of Cathepsin B and Rhodesain with Low-Molecular Inhibitors T1 - Mechanistische Untersuchungen zur Inhibition von Cathepsin B und Rhodesain mit niedermolekularen Inhibitoren N2 - Cysteine proteases play a crucial role in medical chemistry concerning various fields reaching from more common ailments like cancer and hepatitis to less noted tropical diseases, namely the so-called African Sleeping Sickness (Human Arfican Trypanosomiasis). Detailed knowledge about the catalytic function of these systems is highly desirable for drug research in the respective areas. In this work, the inhibition mechanisms of the two cysteine proteases cathepsin B and rhodesain with respectively one low-molecular inhibitor class were investigated in detail, using computational methods. In order to sufficiently describe macromolecular systems, molecular mechanics based methods (MM) and quantum mechanical based method (QM), as well as hybrid methods (QM/MM) combining those two approaches, were applied. For Cathespin B, carbamate-based molecules were investigated as potential inhibitors for the cysteine protease. The results indicate, that water-bridged proton-transfer reactions play a crucial role for the inhibition. The energetically most favoured pathway (according to the calculations) includes an elimination reaction following an E1cB mechanism with a subsequent carbamylation of the active site amino acid cysteine. Nitroalkene derivatives were investigated as inhibitors for rhodesain. The investigation of structurally similar inhibitors showed, that even small steric differences can crucially influence the inhibition potential of the components. Furthermore, the impact of a fluorination of the nitroalkene inhibitors on the inhibition mechanism was investigated. According to experimental data measured from the working group of professor Schirmeister in Mainz, fluorinated nitroalkenes show – in contrast to the unfluorinated compounds – a time dependent inhibition efficiency. The calculations of the systems indicate, that the fluorination impacts the non-covalent interactions of the inhibitors with the enzymatic environment of the enzyme which results in a different inhibition behaviour. N2 - Cysteinproteasen spielen eine wichtige Rolle in der medizinischen Chemie. Nicht nur im Bereich bekannterer Krankheiten wie Krebs oder Hepatitis, sondern auch bezüglich weniger verbreiteter, tropischer Krankheiten wie der sogenannten afrikanischen Schlafkrankheit (Afrikanische Trypanosomiasis) haben diese Enzyme eine große Bedeutung. Im Bereich der Wirkstofffindung ist ein detailliertes Wissen über die katalytische Funktion der an einer Krankheit beteiligten Enzyme unabdingbar .In der vorliegenden Arbeit wurden die Inhibitionsmechanismen der beiden Cysteinproteasen Cathepsin B und Rhodesain in Verbindung mit zwei niedermolekularen Inhibitorklassen anhand theoretischer Berechnungen untersucht. Um die makromolekularen Systeme ausreichend genau beschreiben zu können, wurden neben molekularmechanischen (MM) und quantenmechanischen (QM) Ansätzen auch Hybridmethoden verwendet, welche beide Ansätze (QM/MM) verbinden. Für Cathepsin B wurden Carbamat-basierte Moleküle als potenzielle Inhibitoren der Cysteinprotease untersucht. Die Ergebnisse weisen darauf hin, dass wasser-verbrückte Protonentransferreaktionen eine entscheidende Rolle für die Inhibition spielen. Der laut den Rechnungen energetisch günstigste Mechanismus beinhaltet eine Eliminierungsreaktion nach einem E1cB Mechanismus gefolgt von der Carbamylierung der Aminosäure Cystein in der aktiven Tasche des Enzyms. Nitroalken-Derivate wurden als potenzielle Rhodesain Inhibitoren untersucht. Der Vergleich strukturell ähnlicher Verbindungen weist darauf hin, dass schon kleine sterische Veränderungen einen großen Einfluss auf das Inhibitionspotenzial der Nitroalkene haben können. Außerdem wurde der Einfluss einer Fluorierung der Inhibitoren anhand von Berechnungen untersucht. Messungen der Arbeitsgruppe von Prof. Schirmeister in Mainz zu fluorierten und unfluorierten Nitroalkenen zeigen, dass die fluorierten Verbindungen ein zeitabhängiges Inhibitionspotenzial in Rhodesain aufweisen. Die Berechnungen der Systeme deuten darauf hin, dass die Fluorierung die nicht-kovalenten Wechselwirkungen der Inhibitoren mit der enzymatischen Umgebung des Systems beeinflussen, was zu einem unterschiedlichen Inhibitionsverhalten führt. KW - Cysteinproteasen KW - Inhibitor KW - Mechanismus KW - Berechnung KW - Inhibition Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-178228 ER - TY - THES A1 - Bellinger, Daniel T1 - Implementation of new reaction pathway determining methods and study of solvent effects on the excited state nature of perylene based dyes T1 - Implementierung neuer Reaktionspfad bestimmender Methoden und Untersuchung von Lösungsmitteleinflüssen auf die Natur der angeregten Zustände perylen-basierter Farbstoffe N2 - 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. N2 - 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. KW - Globale Optimierung KW - Reaktionsmechanismus KW - Exziton KW - Reaktionspfadsuche KW - QM/MM KW - organische Halbleiter KW - Exzitonen Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-144435 ER - TY - THES A1 - Le, Thien Anh T1 - Theoretical investigations of proton transfer and interactions or reactions of covalent and non-covalent inhibitors in different proteins T1 - Theoretische Untersuchungen des Protontransfers und Interaktion oder Reaktion von kovalenten und nicht-kovalenten Inhibitoren in verschiedenen Proteinen N2 - Nowadays, computational-aided investigations become an essential part in the chemical, biochemical or pharmaceutical research. With increasing computing power, the calculation of larger biological systems becomes feasible. In this work molecular mechanical (MM) and quantum mechanical approaches (QM) and the combination of both (QM/MM) have been applied to study several questions which arose from different working groups. Thus, this work comprises eight different subjects which deals with chemical reactions or proton transfer in enzymes, conformational changes of ligands or proteins and verification of experimental data. This work firstly deals with reaction mechanisms of aromatic inhibitors of cysteine proteases which can be found in many organisms. These enzymes are responsible for various cancer or diseases as for example Human African Trypanosomiasis (HAT) or the Chagas disease. Aromatic SNAr-type electrophiles might offer a new possibility to covalently modify these proteases. Quantum mechanical calculations have been performed to gain insights into the energetics and possible mechanisms. The next chapter also deals with Trypanosomiasis but the focus was set on a different enzyme. The particularity of Trypanosomiasis is the thiol metabolism which can also be modified by covalent inhibitors. In this context, the wild type and point mutations of the enzyme tryparedoxin have been investigated via molecular dynamic (MD) simulations to examine the influence of specific amino acids in regard to the inhibitor. Experimental data showed that a dimerization of the enzyme occurs if the inhibitor is present. Simulations revealed that the stability of the dimer decreases in absence of the inhibitor and thus confirms these experiments. Further investigations concerning cysteine proteases such as cruzain and rhodesain have been conducted with respect to experimental kinetic data of covalent vinylsulfone inhibitors. Several approaches such as QM or QM/MM calculations and docking, MD or MMPBSA/MMGBSA simulations have been applied to reproduce these data. The utilization of force field approaches resulted in a qualitatively accurate prediction. The kinase AKT is involved in a range of diseases and plays an important role in the formation of cancer. Novel covalent-allosteric inhibitors have been developed and crystallized in complex with AKT. It was shown that depending on the inhibitor a different cysteine residue is modified. To investigate these differences in covalent modification computational simulations have been applied. Enoyl-(acyl carrier) (ENR) proteins are essential in the last step of the fatty acid biosynthesis II (FAS) and represent a good target for inhibition. The diphenylether inhibitor SKTS1 which was originally designed to target the ENR’s of Staphylococcus aureus was also crystallized in InhA, the ENR of Mycobacterium tuberculosis (TB). Crystal structures indicate a change of the inhibitor's tautomeric form. This subject was investigated via MD simulations. Results of these simulations confirmed the tautomerization of the inhibitor. This work also deals with the development of a covalent inhibitor originating from a non-covalent ligand. The target FadA5 is an essential enzyme for the degradation of steroids in TB and is responsible for chronic tuberculosis. This enzyme was crystallized in complex with a non-covalent ligand which served as starting point for this study. Computations on QM or QM/MM level and docking and MD simulations have been applied to evaluate potential candidates. The next chapter focuses on the modification of the product spectrum of Bacillus megaterium levansucrase, a polymerase which catalyzes the biosynthesis of fructans. The covalent modification of the wild type or mutants of the enzyme lead to an accumulation of oligosaccharides but also to polymers with higher polymerization degree. To understand these changes in product spectra MD simulations have been performed. Finally, the proton transfer in catalytic cysteine histidine dyads was investigated. The focus was set on the influence of the relaxation of the protein environment to the reaction. Calculations of the enzymes FadA5 and rhodesain revealed that the preferred protonation state of the dyade depends on the protein environment and has an impact on the reaction barrier. Furthermore, the adaptation of the environment to a fixed protonation state was analyzed via MD simulations. N2 - Heutzutage sind computergestützte Untersuchungen ein essentieller Teil in der chemischen, biochemischen oder pharmazeutischen Forschung. Durch die in den Jahren gestiegene Rechenleistung ist die Berechnung biologischer Systeme möglich. Im Rahmen dieser Arbeit wurden molekularmechanische (MM) und quantenmechanische (QM) Methoden sowie die Kombination beider (QM/MM) für verschiedene Studien eingesetzt, die teilweise aus Fragestellungen verschiedener Arbeitsgruppen hervorgegangen sind. Dadurch umfasst diese Arbeit acht verschiedene Themenkomplexe, bei denen chemische Reaktionen, aber auch der Protonentransfer in Enzymen, Konformationsänderungen von Liganden oder Proteinen und die Verifizierung experimenteller Daten im Fokus standen. Die Arbeit befasst sich anfangs mit Reaktionsmechansimen aromatischer Inhibitoren für Cysteinproteasen, Enzyme, welche in vielen Organismen enthalten sind. Diese Enzyme sind für verschiedene Karzinome oder Krankheiten wie der Afrikanischen Trypanosomiasis oder der Chagas-Krankheit verantwortlich. Aromatische SNAr-Elektrophile bieten hierbei eine neue Möglichkeit der kovalenten Modifikation dieser Proteasen. Quantenmechanische wurden durchgeführt, um Einblicke in die Energetik und mögliche Mechanismen zu erhalten. Das nächste Kapitel befasst sich ebenfalls mit Trypanosomiasis, setzt aber den Fokus auf ein anderes Enzym. Die Besonderheit von Trypanosomiasis ist der Thiol Metabolismus, welcher durch kovalente Inhibitoren modifiziert werden kann. In diesem Kontext wurden der Wildtyp und Punktmutationen des Enzyms Tryparedoxin mittels Molekulardynamik Simulationen untersucht, um Interaktionen einzelner Aminosäuren mit dem kovalenten Inhibitor zu evaluieren. Experimentelle Daten zeigten, dass eine Dimerisierung des Enzyms in Anwesenheit des Inhibitors stattfindet. Durch MD-Simulationen konnte gezeigt werden, dass die Stabilität des Dimers in Abwesenheit des Inhibitors sinkt, wodurch experimentellen Daten bestätigt wurden. Weitere Untersuchungen zu Cysteinproteasen wie Cruzain und Rhodeasin wurden durchgeführt, um experimentelle kinetische Daten von kovalenten Vinylsulfon Inhibitoren zu reproduzieren. Hierbei wurden Methoden wie QM oder QM/MM Rechnungen aber auch Docking, MD und MMPBSA/MMGBSA Simulationen angewandt, um diese Daten zu reproduzieren. In den Untersuchungen zeigte sich, dass die Verwendung der Kraftfeld-basierten Methoden zu qualitativ richtigen Vorhersagen führte. Die Kinase AKT ist in einer Reihe von Krankheiten involviert und spielt eine wichtige Rolle bei der Entstehung von Krebs. Neue kovalent-allosterische Inhibitoren wurden entwickelt und im kovalenten Komplex mit AKT kristallisiert. Die Kristallstrukturen zeigten, dass je nach Inhibitor ein anderes Cystein adressiert wurde. Um diese Unterschiede zu untersuchen, wurden computergestützte Simulationen verwendet. Enoyl-(acyl carrier) (ENR) Proteine sind essentiell für den letzten Schritt in der Fettsäurebiosynthese II (FAS) und bilden ein gutes Target zur Inhibition. Der Diphenylether Inhibitor SKTS1, welchen man ursprünglich als Target für den ENR von Staphylococcus aureus entwarf, wurde auch in InhA, dem ENR von Mycobacterium Tuberculosis (TB), kristallisiert. Die Kristallstrukturen weisen je nach Protein auf einen Wechsel der tautomeren Form des Inhibitors hin. Dieser Sachverhalt wurde mittels MD Simulationen untersucht. Hierbei zeigten die Ergebnisse eine Übereinstimmung mit den experimentellen Daten. Diese Arbeit befasst sich ebenfalls mit der Entwicklung eines kovalenten Inhibitors ausgehend von einem nicht-kovalenten Liganden. Das Target FadA5 ist ein integrales Enzym zur Degradation von Steroiden in TB und ist für die chronische Tuberkulose verantwortlich. Dieses Enzym wurde im Komplex mit einem nicht-kovalenten Liganden kristallisiert, welches als Startpunkt dieser Untersuchungen diente. QM, QM/MM, Docking und MD Simulationen wurden hierbei verwandt, um potentielle Kandidaten zu evaluieren. Das nächste Kapitel befasst sich mit der Modifikation des Produktspektrums von Bacillus megaterium Levansucrase, eine Polymerase, welche die Biosynthese von Fruktanen katalysiert. Durch kovalente Modifikatoren im Wildtyp oder bei Mutanten des Enzyms konnte sowohl eine Anreicherung von Oligosacchariden, aber auch von Polymeren mit höherem Polymerisationsgrad erzielt werden. Um diese Änderungen im Produktspektrum zu verstehen, wurden MD Simulationen durchgeführt. Schließlich wurde die Untersuchung des Protonentransfers in katalytischen Cystein Histidin Dyaden durchgeführt. Hierbei stand der Einfluss der Relaxation der Proteinumgebung auf diese Reaktion im Fokus. Berechnungen in den Enzymen FadA5 und Rhodesain zeigten, dass der präferierte Protonierungszustand der Diade von der Proteinumgebung abhängt und einen großen Einfluss auf die Reaktionsbarriere hat. Um dynamische Effekte einzubeziehen, wurde die Adaption der Umgebung auf einen fixierten Protonierungszustand mittels MD Simulationen analysiert. KW - Computational chemistry KW - Molekularbewegung KW - QM/MM KW - proteins KW - covalent and non-covalent inhibitors KW - Protonentransfer KW - Enzyminhibitor KW - molecular dynamics Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-170511 ER - TY - THES A1 - Welz, Eileen T1 - Theoretical Investigations on Inorganic, Boron-containing Biradicals with a unique Structure T1 - Theoretische Untersuchungen von anorganischen, biradikalischen Borverbindungen mit einer einzigartigen Struktur N2 - In this work, biradical boron containing systems with various structures are investigated to reveal the dependency of the biradical character on the ligated carbene (NHC, CAAC) and the related steric demands of the substituents. N2 - In dieser Arbeit werden biradikalische Bor-Systeme mit unterschiedlichen Strukturen untersucht, um die Abhängigkeit des biradikalischen Charakter von der Carben-Einheit, die am Bor (NHC, CAAC) gebunden ist, und den sterischen Anspruch der Substituenten zu analysieren. KW - Biradical KW - Bor KW - theoretical chemistry KW - SF-DFT KW - Benchmark KW - CAAC KW - NEVPT2 KW - NHC Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-202598 ER - TY - THES A1 - Deißenberger, Andrea T1 - Dibortetrahalogenide für die Darstellung neuer borhaltiger Verbindungen in niedrigen Oxidationsstufen T1 - Diborontetrahalides for the synthesis of novel boron-containing compounds in low oxidation-states N2 - Im Rahmen dieser Arbeit konnten nasschemische Synthesen für Dibortetrafluorid und chlorid ausgehend von Dibortetrabromid entwickelt werden, die durch einfachen Halogenaustausch mit SbF3 bzw. GaCl3 realisiert wurden. In Verbindung mit Arbeiten von Dr. Jonas Müssig zur Synthese von B2I4 gelang die Darstellung aller vier Dibortetrahalogenide mittels einfacher Schlenktechnik basierend auf der Synthese von B2Br4 durch Nöth und Pommerening im Jahr 1981. Dibortetrachlorid konnte mit Phosphanen (PMe3, PCy3 und PPh3) und Singulett-Carbenen (IDipp und MeCAAC) zu den klassischen Bisaddukten 44−46 bzw. 54 und 55 umgesetzt werden. Die Addition eines Isonitrils (CNtBu) an B2Cl4 führte zunächst zur Ausbildung des Bisadduktes 53, allerdings konnte in Lösung eine Umlagerung beobachtet werden, deren Verlauf 11B-NMR-spektroskopisch verfolgt wurde, jedoch nicht final aufgeklärt werden konnte. Durch die Zugabe eines Unterschusses der Lewis-Basen IDipp bzw. PCy3 sollten zunächst Monoaddukte von B2Cl4 dargestellt werden, deren Umsetzung mit einer weiteren Lewis-Base die Synthese asymmetrischer Lewis-Basen-Addukte von B2Cl4 ermöglichen sollte. Die sp2-sp3-Diborane 56 und 57 konnten bei tiefen Temperaturen 11B-NMR-spektroskopisch nachgewiesen werden, allerdings führte eine Abfangreaktion mit diversen Lewis-Basen nicht zu den gewünschten asymmetrischen Addukten. Bei Raumtemperatur konnte eine Folgereaktion von 56 zur Chlorid-verbrückten kationischen Spezies 58 mit einem Tetrachloroborat-Anion beobachtet werden. Im Fall von Dibortetrafluorid konnten keine Lewis-Basen-Addukte (LB = PMe3 und MeCAAC) isoliert werden. Die Reaktivität von B2Cl4 gegenüber ungesättigten Substraten wurde anhand mehrerer literaturbekannter Beispiele (Acetylen, 2-Butin, 3-Hexin, Diphenylacetylen und Bis(trimethylsilyl)acetylen) nachvollzogen und um die terminalen Alkine Propin und 1 Hexin erweitert. Eine selektive Addition von B2Br4 an Dreifachbindungen gelang nicht. Die so erhaltenen Diborylalkene sollten zur Darstellung von 1,2-Diboreten genutzt werden, wobei zunächst über eine von Siebert et al. entwickelte Route die Bis(N,N-dialkylaminochlorboryl)alkene 67g, h, j und k dargestellt wurden. Ein nachfolgender Ringschluss unter reduktiven Bedingungen verlief nur für die Diisopropyl¬amino-substituierten Diborylalkene 67g und j selektiv und lieferte das 1,2-Dihydro-1,2-diboret 71g und das umgelagerte 1,3-Dihydro-1,3-diboret 68j. Der Austausch der Aminosubstituenten gegen Halogenide, der für eine weitere Reduktion zur B-B-Doppelbindung nötig wäre, gelang nicht. Die Umsetzung der Diborylalkene 61 (R = Me), 62 (R = Et) und 65 (R = Ph) mit Singulett-Carbenen (LB = IMe, IiPr, IDipp und MeCAAC) führte zu den chloridverbrückten Monoaddukten 74−76 und 79−81. Alle Verbindungen dieses Typs zeigten in NMR-spektroskopischen Untersuchungen ein sp2- und ein sp3-koordiniertes Borzentrum, welche für die CAAC-stabilisierten Verbindungen auch röntgenkristallografisch nachgewiesen werden konnten. Theoretische Untersuchungen bestätigten die Relevanz des verbrückenden Chloratoms zur Stabilisierung dieser Verbindungen. Für die Stammverbindung der Diborylalkene (59 (R = H)) konnte bei der Umsetzung mit MeCAAC eine unlösliche Verbindung erhalten werden, deren Struktur als Bisaddukt 82 mittels NMR-spektroskopischen Untersuchungen im Festkörper und durch Verbrennungsanalyse bestätigt werden konnte. Die Reduktion der CAAC-stabilisierten Diborylalkene 79 und 80 in Gegenwart von MeCAAC führte zu den captodativ-stabilisierten Diborylradikalen 83 und 84, deren Strukturanalyse eine orthogonale Anordnung der C2-Brücke zur B(CAAC)-Einheit offenlegt. Ausführliche EPR-spektroskopische Untersuchungen bei variabler Temperatur und theoretische Berechnungen bestätigen eine schwache Wechselwirkung der beiden Radikalzentren und einen offenschaligen Singulett-Grundzustand mit einem energetisch tiefliegenden Triplett-Zustand (ΔES T = 0.017 kcal mol−1). Der experimentell bestimmte Spin-Spin-Abstand und die Analyse der einfach besetzten Molekülorbitale (SOMO) bestätigen eine Delokalisierung der Spindichte über die NCAAC-CCAAC-B-Einheit. Der Austausch der verbrückenden Einheit und die somit einhergehende Verringerung der Sterik führt zu einer Planarisierung des Moleküls im Festkörper (87). Theoretische Untersuchungen und die Auswertung der strukturellen Parameter ergeben eine Delokalisierung der Elektronendichte über das gesamte planare System. EPR- und NMR-spektroskopische Untersuchungen ergaben dennoch Hinweise auf das Vorliegen einer paramagnetischen Verbindung. Untersuchungen zum Reduktionsverhalten von zweifach CAAC-stabilisiertem 1,4-Bis-(dibromboryl)benzol (97) ergaben die vollständige Enthalogenierung der Borzentren. Im Zuge dessen entstand ein hochreaktives, lineares Borylen, welches eine CH-Aktivierung mit dem Isopropylsubstituenten des CAAC-Liganden eingeht (98). Zur Stabilisierung des Borylens wurde die Reduktion in Gegenwart weiterer Lewis-Basen (Pyridin (Pyr), IiPr, IMeMe, PMe3, CNtBu und CO) durchgeführt, die in der Ausbildung der Diborylene 99−104 resultierten. Die Darstellung einer para-Phenylen-verbrückten Donor-Akzeptor-Verbindung (D: Borylen, A: BMes2) gelang nicht. N2 - Within the scope of this work, solution-phase syntheses of diborontetrafluoride and chloride were developed using simple halogen exchange reactions, by reacting the precursor B2Br4 with SbF3 or GaCl3, respectively. In combination with the work of Dr. Jonas Müssig on the synthesis of B2I4, the preparation of all four diborontetrahalides was achieved by usual Schlenk technique, based on the synthesis of diborontetrabromide by Nöth and Pommerening in 1981. Diborontetrachloride was reacted with phosphines (PMe3, PCy3 and PPh3) and singlet carbenes (IDipp and MeCAAC) to yield bisadducts 44, 45, 46, 54 and 55, respectively. Reaction with an isonitrile (CNtBu) initially resulted in the formation of bisadduct 53, which underwent rearrangement in solution. The process was investigated by 11B NMR spectroscopy but was not finally resolved. Addition of a substoichiometric amount of the Lewis bases IDipp and PCy3 was supposed to lead to a monoadduct which should be subsequently reacted with another Lewis base to yield an asymmetric Lewis base adduct of B2Cl4. The sp2-sp3-diboranes 56 and 57 were indeed detected by 11B NMR spectroscopy at low temperatures but could not be trapped by another Lewis base to form the desired asymmetric adducts. At room temperature 56 underwent a follow-up reaction to give the chloride-bridged cationic species 58, supported by a tetrachloroborate anion. In the case of diborontetrafluoride no Lewis base adducts (LB = PMe3 and MeCAAC) could be isolated. The reactivity of B2Cl4 towards unsaturated substrates was confirmed for several literature known examples (acetylene, 2-butyne, 3-hexyne, diphenylacetylene and bis(trimethylsilyl)acetylene) and expanded to the terminal alkynes propyne and 1-hexyne. The addition of B2Br4 to triple bonds was not selective. Those diborylalkenes should be utilized for the syntheses of 1,2-diboretes. Bis(N,N-dialkylaminochloroboryl)alkenes 67g, h, j and k were synthesized according to a route developed by Siebert et al. but the subsequent reductive ring closure was only selective for diborylalkenes bearing diisopropylamino substituents 67g and j, yielding the 1,2-dihydro-1,2-diborete 71g and the rearranged 1,3-dihydro-1,3-diborete 68j.[155] An exchange of the amino substituents for halides, which would have been necessary for a further reduction to yield a B-B double bond, was not possible. Reaction of the diborylalkenes 61 (R = Me), 62 (R = Et) and 65 (R = Ph) with singlet carbenes (LB = IMe, IiPr, IDipp and MeCAAC) led to the chloride-bridged monoadducts 74−76 and 79−81. NMR spectroscopic investigations of those examples showed an sp2- as well as an sp3-hybridized boron atom, which could be structurally determined for the CAAC-stabilized compounds by X-ray crystallography. Theoretical analyses showed the importance of the bridging chloride for the stability of these compounds. Addition of MeCAAC to the parent diborylalkene 59 (R = H) resulted in an insoluble precipitate, whose bisadduct structure 82 was confirmed by solid-state NMR spectroscopy and elemental analysis. Reduction of CAAC-stabilized diborylalkenes 79 and 80 in the presence of MeCAAC led to the formation of captodatively stabilized diborylradicals 83 and 84, whose olefinic bridging unit lies perpendicular to the B(CAAC) unit, as shown by X-ray crystallography. Detailed EPR spectroscopic investigations at variable temperature alongside with theoretical studies confirmed a weak interaction of both radical centers and an open shell singlet ground state with a triplet state negligibly higher in energy (ΔES T = 0.017 kcal mol−1). The experimentally determined interspin distance and the analysis of the singly occupied molecular orbitals (SOMO) revealed a delocalization of the spin density over the NCAAC-CCAAC-B π system. Exchange of the bridging unit accompanied by the reduction of steric crowding led to a planarization of the molecule in the solid state (87). Theoretical studies and the evaluation of the structural parameters showed the delocalization of electron density over the entire planar system. EPR and NMR spectroscopic investigations indicated a paramagnetic nature of 87. Reduction of doubly CAAC-stabilized 1,4-bis(dibromoboryl)benzene (97) resulted in complete dehalogenation of the boron centers. As a result, a highly reactive, linear borylene was formed, which underwent CH activation with the isopropyl substituent of the CAAC ligand (98). In order to stabilize the borylene, reductions were performed in the presence of other Lewis bases, such as pyridine (Pyr), IiPr, IMeMe, PMe3, CNtBu and CO, resulting in the formation of diborylenes 99−104. The attempted synthesis of a para-phenylene-bridged donor-acceptor compound (D: borylene, A: BMes2) was not successful. KW - Dibortetrahalogenide KW - Biradikal KW - Borylene KW - Borylradikale Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-187758 ER - TY - THES A1 - Walter, Christof T1 - Excitonic States and Optoelectronic Properties of Organic Semiconductors - A Quantum-Chemical Study Focusing on Merocyanines and Perylene-Based Dyes Including the Influence of the Environment T1 - Exzitonische Zustände und optoelektronische Eigenschaften organischer Halbleiter – Eine quantenchemische Untersuchung mit Fokus auf Merocyaninen und perylenbasierten Farbstoffen unter Berücksichtigung der Umgebung N2 - The scope of computational chemistry can be broadened by developing new methods and more efficient algorithms. However, the evaluation of the applicability of the methods for the different fields of chemistry is equally important. In this thesis systems with an unusual and complex electronic structure, such as excitonic states in organic semiconductors, a boron-containing bipolaron and the excited states of pyracene were studied and the applicability of the toolkit of computational chemistry was investigated. Concerning the organic semiconductors the focus was laid on organic solar cells, which are one of the most promising technologies with regard to satisfying the world's need for cheap and environmentally sustainable energy. This is due to the low production and material costs and the possibility of using flexible and transparent devices. However, their efficiency does still not live up to the expectations. Especially the exciton diffusion lengths seem to be significantly too short. In order to arrive at improved modules, a fundamental understanding of the elementary processes occurring in the cell on the molecular and supramolecular level is needed. Computational chemistry can provide insight by separating the different effects and providing models for predictions and prescreenings. In this thesis, the focus was laid on the description of excitonic states in merocyanines and perylene-based dyes taking the influence of the environment into account. At first, the photochemical isomerization between two configurations of 6-nitro BIPS observed experimentally was studied by first benchmarking several functionals against SCS-ADC(2) in the gas phase and subsequently calculating the excited-state potential energy surface. The geometries obtained from a relaxed scan in the ground state as well as from a scan in the excited state were used. The environment was included using different polarizable continuum models. It was shown that the choice of the model and especially the question of the state specificity of the approach is of vital importance. Using the results of the calculations, a two-dimensional potential energy surface could be constructed that could be used to explain the experimental findings. Furthermore, the importance of the excited-state isomerization as a potential deactivation channel in the exciton transport was pointed out. Then the assessment of the suitability of different merocyanines for optoelectronic applications with quantum-chemical methods was discussed. At first, the effect of the environment on the geometry, especially on the bond length alternation pattern, was investigated. It was shown that the environment changes the character of the ground-state wave function of several merocyanines qualitatively, which means that the results of gas-phase calculations are meaningless - at least when a comparison with solution or device data is desired. It was demonstrated that using a polarizable continuum model with an effective epsilon, a qualitative agreement between the calculated geometry and the geometry in the crystal structure can be obtained. Therefore, by comparing the bond length alternation in solution and in the crystal, a rough estimate of the effect of the crystal environment can be made. It was further shown that the connection between the HOMO energy and the open-circuit voltage is not as simple as it is often implied in the literature. It was discussed that it is not clear whether the HOMO of a single molecule or a $\pi$-stack containing several monomers should be used and if the environmental charges of the bulk phase or the interface should be included. Investigating the dependence of the HOMO energy on the stack size yielded no definitive trend. Furthermore, it was discussed that the effect due the optimization of the modules (solvent, bulk heterojunction) during the production masks any potential correlation between the HOMO energy and measured open-circuit values. Therefore, a trend can only be expected for unoptimized bilayer cells. It was concluded that ultimately, the importance of the HOMO energy should not be overestimated. The correlation between the exciton reorganization energy and the so-called cyanine limit, which is predicted by a simple two-state model, was also discussed. By referring to the results of VB calculations, it was discussed that the correlation indeed exists and is non-negligible, although the effect is not as strong as one might have expected. In this context, a potential application of a VB/MM approach was covered briefly. The importance of the molecular reorganization energy and the device morphology was also discussed. It was concluded that the optimization of merocyanines for organic optoelectronic devices is inherently a multiparameter problem and one cannot expect to find one particular parameter, which solely controls the efficiency. The perylene-based dyes were studied with a focus on the description of a potential trapping mechanism involving an intermolecular motion in a dimer. The aim was to find methods which can be applied to larger model systems than a dimer and take the effect of the environment into account. As a test coordinate the longitudinal shift of two monomers against each other was used. At first, it was demonstrated how the character of an excited state in a dimer can be defined and how it can be extracted from a standard quantum-chemical calculation. Then several functionals were benchmarked and their applicability or failure was rationalized using the character analysis. Two recipes could be proposed, which were applied to a constraint optimization (only intermolecular degrees of freedom) in the excited states of the PBI dimer and to the description of the potential energy surfaces of ground and excited states along a longitudinal displacement in the perylene tetramer, respectively. It was further demonstrated that the semi-empirical OMx methods fail to give an accurate description of the excited-state potential energy surfaces as well as the ground-state surface along the test coordinate. This failure could be attributed to an underestimation of overlap-dependent terms. Consequently, it could be shown that the methods are applicable to large intermolecular distances, where the overlap is negligible. The results of DFT calculations with differently composed basis sets suggested that adding an additional single p-function for each atom should significantly improve the performance. QM/MM methods are ideally suited to take the effect of the environment on a a dimer model system into account. However, it was shown that standard force fields also give an incorrect description of the interaction between the monomers along the intermolecular coordinate. This failure was attributed to the isotropic atom-atom interaction in the repulsion term of the Lennard-Jones potential. This was corroborated using two simple proof-of-principle anisotropy models. Therefore, a novel force field called OPLS-AA_O was presented that is based on OPLS-AA, but uses an anisotropic model for the repulsion. The model involves the overlap integral between the molecular densities, which are modeled as a sum of atom-centered p-type Gaussian functions. It was shown that using this force field an excellent agreement with the DFT results can be obtained when the correct parameters are used. These parameters, however, are not very generalizable, which was attributed to the simplicity of the model in its current state (using the same exponential parameter for all atoms). As a short excursion, the applicability of an MO-based overlap model was discussed. It was demonstrated that the repulsion term based on the density overlap can be used to correct the failure of the OMx methods for the ground states. This is in accord with the assumption that an underestimation of the overlap terms is responsible for the failure. It was shown that OPLS-AA_O also gives an excellent description of the longitudinal shift in a PBI tetramer. Using the tetramer as a test system and applying the recipe obtained in the TDDFT benchmark for the QM-part and OPLS-AA_O for the MM-part in conjunction with an electrostatic embedding scheme, a QM/MM description of the excited states of the PBI dimer including the effect of the environment could be obtained. In the last chapter the theoretical description of the Bis(borolyl)thiophene dianion and the excited states of pyracene were discussed. The electronic structure of the Bis(borolyl)thiophene dianion - a negative bipolaron - was elucidated using DFT and CASPT2 methods. Furthermore, an estimation of the extent of triplet admixture to the ground state due to spin-orbit coupling was given. In the second project the S1 and S2 states of pyracene were computed using SCS-CC2 and SCS-ADC(2) and an estimation for the balance between aromaticity and ring strain was given. This also involved computing the vibrational frequencies in the excited states. In both studies the results of the computations were able to rationalize and complete experimental results. N2 - Die Anwendungsmöglichkeiten der Methoden der theoretischen Chemie können erweitert werden, indem neue Methoden und effizientere Algorithmen entwickelt werden. Es ist jedoch ebenso wichtig die Anwendbarkeit der Methoden für die verschiedenen Felder der Chemie zu evaluieren. In dieser Arbeit wurden Systeme mit einer komplexen und ungewöhnlichen Struktur, wie exzitonische Zustände in organischen Halbleitern, ein bor-basiertes Bipolaron und die angeregten Zustände von Pyracen untersucht und die Anwendbarkeit der verschiedenen Methoden evaluiert. Im Bezug auf die organischen Halbleiter wurde der Fokus auf organische Solarzellen gelegt, welche zu den vielversprechendsten Technologien gehören, um dem weltweiten Bedarf an billiger und ökologisch nachhaltiger Energie zu begegnen. Dies liegt an den niedrigen Produktionskosten und der Möglichkeit flexible und transparente Module zu verwenden. Ihre Wirkungsgrade werden den Erwartungen jedoch noch nicht gerecht. Vor allem die Exzitonendiffusionslängen scheinen deutlich zu gering zu sein. Um verbesserte Module zu erhalten ist ein fundamentales Verständnis der Elementarprozesse in der Zelle auf molekularem und supramolekularem Level vonnöten. Die theoretische Chemie kann dabei helfen dies zu erreichen, indem sie die verschiedenen Effekte separiert und Modelle für Vorhersagen und zur Vorauswahl geeigneter Verbindungen bereitstellt. In dieser Arbeit wurde der Fokus auf die Beschreibung von exzitonischen Zuständen in Merocyaninen und perylenbasierten Farbstoffen unter Berücksichtigung von Umgebungseinflüssen gelegt. Zunächst wurde die experimentell beobachtete photochemische Isomerisierung zwischen zwei Konfigurationen von 6-nitro BIPS untersucht, indem zuerst die Anwendbarkeit verschiedener Funktionale im Vergleich zu SCS-ADC(2) in der Gasphase überprüft wurde und anschließend die Potentialfläche des angeregten Zustands berechnet wurde. Es wurden sowohl die Geometrien aus einem relaxed scan im Grundzustand als auch von einem scan im angeregten Zustand verwendet. Umgebungseffekte wurden unter Verwendung verschiedener Kontinuumsansätze (polarizable continuum models) berücksichtigt. Es konnte gezeigt werden, dass die Wahl des Ansatzes und vor allem die Frage nach der Zustandsspezifizität des Kontinuumsansätze sehr entscheidend ist. Mit den Ergebnissen der Berechnungen konnte eine zweidimensionale Potenzialfläche konstruiert werden, mittels welcher die experimentellen Beobachtungen erklärt werden konnten. Außerdem wurde auf die Bedeutung der Isomerisierung im angeregten Zustand als einem potenziellen Deaktivierungskanal für den Exzitonentransport hingewiesen. Anschließend wurde die Möglichkeit einer Bewertung der Eignung verschiedener Merocyanine für optoelektonische Fragestellungen mit quantenchemischen Methoden diskutiert. Zunächst wurde der Einfluss der Umgebung auf die Geometrie und insbesondere auf die Bindungslängenalternanz untersucht. Es wurde gezeigt, dass die Umgebung die Wellenfunktion mehrerer Merocyanine qualitativ verändert, was bedeutet, dass Berechnungen in der Gasphase keinen Sinn machen - zumindest nicht, wenn die Ergebnisse mit Daten, die in Lösung oder in der Zelle erhalten wurden, verglichen werden sollen. Es konnte gezeigt werden, dass unter Verwendung eines Kontinuumsansatzes mit einer effektiven Dielektrizitätskonstante epsilon eine qualitative Übereinstimmung zwischen der berechneten Geometrie und der Geometrie in der Kristallstruktur erzielt werden kann. Dies ermöglicht es, durch einen Vergleich der Bindungslängenalternanz in Lösung und im Kristall eine grobe Abschätzung für den Einfluss der Kristallumgebung zu erhalten. Es wurde außerdem dargelegt, dass der Zusammenhang zwischen der Energie des HOMOs und der Leerlaufspannung nicht so eindeutig ist, wie es oft in der Literatur suggeriert wird. Es stellte sich die Frage, ob die HOMO-Energie eines einzelnen Moleküls oder eines Stapels bestehend aus mehreren Monomeren verwendet werden sollte und ob der Umgebungseffekt der Ladungen der Bulkphase oder der Grenzfläche berücksichtigt werden sollte. Die Untersuchung der Abhängigkeit der HOMO-Energie von der Anzahl der Monomere ergab keinen klaren Trend. Die Tatsache, dass die Optimierung des Moduls während des Produktionsprozesses (Solvent, Bulk-Hereojunction-Konzept) eine potenzielle Korrelation zwischen der HOMO-Energie und der Leerlaufspannung maskiert, wurde ebenfalls diskutiert. Deshalb kann eine Korrelation nur für nicht optimierte Zweischichtzellen erwartet werden. Es wurde der Schluss gezogen, dass die Bedeutung der HOMO-Energie letztendlich nicht überbewertet werden sollte. Der Zusammenhang zwischen der Exzitonenreorganisationsenergie und dem sogenannten Cyaninlimit, welcher von einem einfachen Zwei-Zustands-Model vorhergesagt wird wurde diskutiert. Unter Verweis auf die Ergebnisse von VB-Berechnungen konnte diskutiert werden, dass der Zusammenhang in der Tat existiert und nicht vernachlässigbar, aber auch nicht so groß ist, wie man vermutet haben könnte. In diesem Kontext wurde die potenzielle Anwendbarkeit eines VB/MM-Ansatzes kurz besprochen. Die Bedeutung der molekularen Reorganisationsenergie und der Morphologie der Zelle wurden ebenfalls diskutiert. Es wurde das Fazit gezogen, dass die Optimierung der Merocyanine für die Anwendung in organischen Halbleitern inhärent ein Multiparameterproblem ist und man nicht erwarten kann, einen einzelnen Parameter zu finden, der allein die Effizienz kontrolliert. Die perylenbasierten Farbstoffe wurden mit dem Fokus auf der Beschreibung eines potenziellen Exzitoneneinfangmechanismus, untersucht, welcher auf der intermolekularen Bewegung in einem Dimer basiert. Das Ziel war es Methoden zu finden, die auf größere Systeme anwendbar sind und den Umgebungseinfluss berücksichtigen können. Als Testkoordinate wurde die longitudinale Verschiebung der Monomere gegeneinander verwendet. Zunächst wurde gezeigt, wie der Charakter eines angeregten Zustandes in einem Dimer definiert werden kann und wie ein Maß für den Charakter ausgehend von einer normalen quantenchemischen Berechnung erhalten werden kann. Anschließend wurden verschiedene Funktionale evaluiert und ihre Anwendbarkeit beziehungsweise ihr Versagen mittels der Charakteranalyse rationalisiert. Zwei Ansätze konnten vorgeschlagen werden, welche auf eine Optimierung in den angeregten Zustände des Dimers mit Nebenbedingung (nur intermolekulare Freiheitsgrade) beziehungsweise auf eine Beschreibung der Potenzialflächen des Grundzustandes und der angeregten Zustände für die longitudinale Verschiebung in einem Perylentetramer angewendet wurden. Es wurde außerdem gezeigt, dass die semiempirischen OMx Methoden keine akkurate Beschreibung der Potenzialflächen der angeregten Zustände sowie des Grundzustandes für die Testkoordinate liefern. Dies konnte mit der Unterschätzung der intermolekularen Überlappterme begründet werden. Folglich war es möglich zu zeigen, dass die Methoden für intermolekulare Abstände, bei denen der Überlapp vernachlässigbar ist, anwendbar sind. Die Ergebnisse von DFT-Rechnungen mit unterschiedlich zusammengesetzten Basissätzen ließen ferner den Schluss zu, dass das Hinzufügen einer einzelnen p-Funktion an jedem Atom eine deutliche Verbesserung bringen sollte. QM/MM-Methoden sind ideal geeignet, um den Einfluss der Umgebung auf ein Dimer-Modellsystem zu berücksichtigen. Es wurde jedoch gezeigt, dass gängige Kraftfelder ebenfalls eine inkorrekte Beschreibung der Wechselwirkung zwischen den Monomeren entlang der intermolekularen Koordinate liefern. Dies wurde mit der isotropen Beschreibung der Atom-Atom-Wechselwirkung im Repulsionsterm des Lennard-Jones-Potenzials begründet. Diese Annahme wurde durch die Anwendung zweier Proof-of-Principle-Ansätze untermauert. Folglich wurde ein neues Kraftfeld, genannt OPLS-AA_O, eingeführt, welches auf OPLS-AA basiert, aber eine anisotrope Modellierung der Repulsion verwendet. Diese anisotrope Repulsion basiert auf dem Überlappintegral der molekularen Elektronendichten, welche als Summe aus atomzentrierten p-artigen Gaußfunktionen modelliert wird. Es wurde gezeigt, dass mit diesem Kraftfeld eine hervorragende Übereinstimmung mit den DFT-Ergebnissen erhalten werden kann, wenn die richtigen Parameter verwendet werden. Diese Parameter sind jedoch nicht sehr generalisierbar, was mit der Einfachheit des Models zu seinem momentanen Stand begründet wurde (Verwendung desselben Parameters im Exponenten bei allen Atomen). Als kurzer Exkurs wurde die Anwendbarkeit eines MO-basierten Überlappmodells diskutiert. Es konnte nachgewiesen werden, dass der Repulsionsterm, der auf der Dichteüberlappung basiert, auch als Korrekturterm für die Anwendbarkeit der OMx-Methoden bezüglich des Grundzustandes verwendet werden kann. Dies deckt sich mit der Annahme, dass eine Unterschätzung von Überlapptermen für das Versagen der semiempirischen Methoden verantwortlich ist. Es wurde gezeigt, dass OPLS-AA_O die Potenzialfläche für die longitudinale Verschiebung in einem PBI Tetramer exzellent beschriebt. Unter Verwendung des Tetramers als Testsytem und unter Anwendung eines der vorgeschlagenen TDDFT-Ansätze für den QM-Teil und OPLS-AA_O für den MM-Teil in Verbindung mit einem electrostatic embedding-Ansatz konnte eine QM/MM-Beschreibung der angeregten Zustände des PBI Dimers unter Berücksichtigung des Umgebungseinfluss erhalten werden. Im letzten Kapitel wurde die theoretische Beschreibung des Bis(borolyl)thiophendianions und von Pyracen diskutiert. Die elektronische Struktur des Bis(borolyl)thiophendianions wurde beschrieben unter Verwendung von DFT- und CASPT2-Methoden. Außerdem wurde eine Abschätzung des Ausmaßes der Triplettbeimischung zum Grundzustand durch die Spin-Bahn-Kopplung gegeben. Im zweiten Projekt wurden der S1- und S2- Zustand des Pyracens unter Verwendung von SCS-CC2 und SCS-ADC(2) berechnet und eine Abschätzung des Verhältnisses von Aromatizität und Ringspannung gegeben. Dies beinhaltete auch die Berechnung der Schwingungsfrequenzen im angeregten Zustand. In beiden Studien konnten die Ergebnisse der Berechnungen die experimentellen Daten vervollständigen und rationalisieren. KW - Exziton KW - Angeregter Zustand KW - Quantenchemie KW - organic semiconductors KW - organische Halbleiter Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-123494 ER - TY - THES A1 - Schmidt, Thomas Christian T1 - Theoretical Investigations on the Interactions of Small Compounds with their Molecular Environments T1 - Theoretische Untersuchungen der Wechselwirkungen Kleiner Moleküle mit deren Molekularen Umgebungen N2 - Im ersten Teil dieser Arbeit wird eine Kombination theoretischer Methoden für die strukturbasierte Entwicklung neuer Wirkstoffe präsentiert. Ausgehend von der Kristallstruktur eines kovalenten Komplexes einer Modellverbindung mit dem Zielprotein wurde mit Hilfe von quantenmechanischen und QM/MM Rechnungen die genaue Geometrie des vorausgehenden nicht-kovalenten Komplexes betimmt. Letztere ist der bestimmende Faktor für die Reaktivität des Inhibitors gegenüber der katalytisch aktiven Aminosäure und damit für die Ausbildung einer kovalenten Bindung. Aus diesem Grund wurde diese Geometrie auch für die Optimierung der Substitutionsmusters des Ihnibitors verwendet, um dessen Affinität zum Zielenzyme zu verbessern ohne dass dieser seine Fähigkeit kovalent an das aktive Zentrum zu binden verliert. Die Optimierung des Substitutionsmuster wurde doch Methode des Molekularen Dockings unterstützt, das diese optimal dazu geeignet sind, Bindungsaffinitäten vorherzusagen, die durch eine Modifikation der chemischen Struktur entstehen. Eine Auswahl der besten Strukturen wurde anschließend verwendet, um zu überprüfen, ob die veränderten Moleküle noch genügen Reaktivität gegenüber dem Zielprotein aufweisen. Moleküldynamik Simulationen der neuen Verbindungen haben jedoch gezeigt, dass die veränderten Verbindungen nur so and das Protein binden, dass die Bilung eine kovalenten Bindung zum Enzym nicht mehr möglich ist. Daher wurden in einem weiteren Schritt die Modellverbindungen weiter modifiziert. Neben Änderungen im Substitutionsmuster wurde auch die chemische Struktur im Kern verändert. Die Bindungsaffinitäten wurde wieder mittels Docking überprüft. Für die besten Bindungsposen wurden wieder Simulationen zur Moleküldynamik durchgeführt, wobei diesmal die Ausbildung einer kovalenten Bindung zum Enzyme möglich erscheint. In einer abschließenden Serie von QM/MM Rechnungen unter Berücksichtigung verschiedener Protonierungszustände des Inhibitors und des Proteins konnten Reaktionspfade und zugehörige Reaktionsenergien bestimmt werden. Die Ergebnisse lassen darauf schließen, dass eines der neu entwickelten Moleküle sowohl eine stark verbesserte Bindungsaffinität wie auch die Möglichkeit der kovalenten Bindung an Enzyme aufweist. Der zweite Teil der Arbeit konzentriert sich auf die Umgebungseinflüsse auf die Elektronenverteilung eines Inhibitormodells. Als Grundlage dient ein vinylsulfon-basiertes Moekül, für das eine experimentell bestimmte Kristallstruktur sowie ein theoretisch berechneter Protein Komplex verfügbar sind. Ein Referendatensatz für diese Systeme wurde erstellt, indem der Konformationsraum des Inhibitors nach möglichen Minimumsstrukturen abgesucht wurde, welche später mit den Geometrien des Moleküls im Kristall und im Protein verglichen werden konnten. The Geometrie in der Kristallumgebung konnte direkt aus den experimentellen Daten übernommen werden. Rechnungen zum nicht-kovalenten Protein Komplex hingegen haben gezeigt, dass für das Modellsystem mehrere Geometrien des Inhibiors sowie zwei Protonierungszustände für die katalytisch aktiven Aminosäuren möglich sind. Für die Analyse wurden daher alle möglichen Proteinkomplexe mit der Kristallstruktur verglichen. Ebenso wurden Vergleiche mit der Geometrie des isolierten Moleküls im Vakuum sowie der Geometrie in wässriger Lösung angestellt. Für die Geometrie des Moleküls an sich ergab sich eine gute Übereinstimmung für alle Modellsysteme, für die Wechselwirkungen mit der Umgebung jedoch nicht. Die Ausbildung von Dimeren in der Kristallumgebung hat einen stark stablisierenden Effekt und ist einer der Gründe, warum dieser Kristall so gut wie keine Fehlordungen aufweist. In den Proteinkomplexen hingegen ergibt sich eine Abstoßung zwischen dem Inhibitor und einer der katalytisch aktiven Aminosäuren. Als Ursache für diese Abstoßung konnte die Einführung der Methylaminfunktion ausgemacht werden. Vermutlicherweise führt diese strukturelle Änderung auch dazu, dass der Modellinhibitor nicht in der Lage ist, so wie die Leitstruktur K11777 an das aktive Zentrum des Enzyms zu binden. N2 - In the first part of this work, a combination of theoretical methods for the rational design of covalent inhibitor is presented. Starting from the crystal structure of the covalent complex of a lead compound, quantum mechanical and QM/MM calculations were used to derive the exact geometry of the preceeding non-covalent enzyme inhibitor complex. The geometry of the latter mainly determines the reactivity of the inhibitor against its target enzyme concerning the formation of the covalent bond towards an active site residue. Therefore, this geometry was used as starting point for the optimization of the substitution pattern of the inhibitor such as to increase its binding affinity without loosing its ability to covalently bind to the target protein. The optimization of the chemical structure was supported by using docking procedures, which are best suited to estimate binding affinities that arise from the introduced changes. A screening of the novel substitution patterns resulted in a first generation of model compounds which were further tested for their reactivity against the target. Dynamic simulations on the novel compounds revealed that the orientation that compounds adopt within the active site are such that a covalent interaction with the enzyme is no longer possible. Hence, the chemical structure was further modified, including not only changes in the substituents but also within the core of the molecule. Docking experiments have been conducted to assure sufficiently high binding affinities and to obtain the most favored binding poses. Those have then again been used for dynamic simulations which resulted in structures, for which the bond formation process appeared feasible. A final series of QM/MM calculations considering various protonation states was computed to estimate the reaction energies for the covalent attachment of the inhibitor to the enzyme. The theoretical results indicate a reasonable high inhibition potency of the novel compounds. The second part concentrates on the environmental influences on the electron density of an inhibitor molecule. Therefore, a vinylsulfone-based model compound was selected for which an experimental crystal structure for the pure compound as well as a theoretically determined enzyme-inhibitor complex have been available. To provide reference data for the larger systems, the conformational space of the isolated molecule was screened for favorable geometries which were later compared to those within the crystal and protein surrounding. The geometry of the crystal structure could readily be taken from the experimental data whereas calculations on the protein complex revealed four potential non-covalent complexes exhibiting different arrangements of the molecule within the active site of the protein as well as two possible protonation states of the catalytic dyad. Hence, all four protein complexes have been compared to the crystal structure of the molecule as well as against the more favorable geometries of the isolated molecule being determined within vacuum or aqueous surrounding. Whereas the molecule itself was found to adopt comparable geometries within all investigated environments, the interactions pattern between the crystal surrounding and the protein differed largely from each other. The favorable formation of dimers within the crystal has a strong stabilizing effect and explains the extraordinarily good quality of the crystal. Within the protein however, repulsive forces have been found between the protein and the inhibitor. The origin of the repulsion could be traced back to effect of on of the substituents to the vinyl scaffold. The difference in the chemical structure in comparison to a well known inhibitor might also explain the experimentally found loss of activity for the model compound in comparison to K11777. KW - Theoretische Chemie KW - theoretical chemistry KW - electron density KW - inhibition KW - Elektronendichte KW - Inhibitor Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-127860 ER - TY - THES A1 - Becker, Johannes T1 - Development and implementation of new simulation possibilities in the CAST program package T1 - Entwicklung und Implementierung neuer Simulationsmöglichkeiten in das CAST Programmpaket N2 - The aim of the present work is the development and implementation of new simulation possibilities for the CAST program package. Development included, among other things, the partial parallelization of the already existing force fields, extension of the treatment of electrostatic interactions and implementation of molecular dynamics and free energy algorithms. The most time consuming part of force field calculations is the evaluation of the nonbonded interactions. The calculation of these interactions has been parallelized and it could be shown to yield a significant speed up for multi-core calculations compared to the serial execution on only one CPU. For both, simple energy/gradient as well as molecular dynamics simulations the computational time could be significantly reduced. To further increase the performance of calculations employing a cutoff radius, a linkedcell algorithm was implemented which is able to build up the non-bonded interaction list up to 7 times faster than the original algorithm. To provide access to dynamic properties based on the natural time evolution of a system, a molecular dynamics code has been implemented. The MD implementation features two integration schemes for the equations of motion which are able to generate stable trajectories. The basic MD algorithm as described in Section 1.2 leads to the sampling in the microcanonical (NVE) ensemble. The practical use of NVE simulations is limited though because it does not correspond to any experimentally realistic situation. More realistic simulation conditions are found in the isothermal (NVT) and isothermalisobaric (NPT) ensembles. To generate those ensembles, temperature and pressure control has been implemented. The temperature can be controlled in two ways: by direct velocity scaling and by a Nose-Hoover thermostat which produces a real canonical ensemble. The pressure coupling is realized by implementation of a Berendsen barostat. The pressure coupling can be used for isotropic or anisotropic box dimensions with the restriction that the angles of the box need to be 90� . A crucial simulation parameter in MD simulations is the length of the timestep. The timestep is usually in the rang of 1fs. Increasing the timestep beyond 1fs can lead to unstable trajectories since the fastest motion in the system, usually the H-X stretch vibration can not be sampled anymore. A way to allow for bigger timesteps is the use of a constraint algorithm which constrains the H-X bonds to the equilibrium distance. For this the RATTLE algorithm has been implemented in the CAST program. The velocity Verlet algorithm in combination with the RATTLE algorithm has been shown to yield stable trajectories for an arbitrary length of simulation time. In a first application the MD implementation is used in conjunction with the MOPAC interface for the investigation of PBI sidechains and their rigidity. The theoretical investigations show a nice agreement with experimentally obtained results. Based on the MD techniques two algorithms for the determination of free energy differences have been implemented. The umbrella sampling algorithm can be used to determine the free energy change along a reaction coordinate based on distances or dihedral angles. The implementation was tested on the stretching of a deca-L-alanine and the rotation barrier of butane in vacuum. The results are in nearly perfect agreement with literature values. For the FEP implementation calculations were performed for a zero-sum transformation of ethane in explicit solvent, the charging of a sodium ion in explicit solvent and the transformations of a tripeptide in explicit solvent. All results are in agreement with benchmark calculations of the NAMD program as well as literature values. The FEP formalism was then applied to determine the relative binding free energies between two inhibitors in an inhibitor-protein complex. Next to force fields, ab-initio methods can be used for simulations and global optimizations. Since the performance of such methods is usually significantly poorer than force field applications, the use for global optimizations is limited. Nevertheless significant progress has been made by porting these codes to GPUs. In order to make use of these developments a MPI interface has been implemented into CAST for communication with the DFT code TeraChem. The CAST/TeraChem combination has been tested on the $H_2 O_{10}$ cluster as well as the polypeptide met-Enkephalin. The pure ab-initio calculations showed a superior behavior compared to the standard procedure where the force field results are usually refined using quantum chemical methods. N2 - Das Ziel der hier vorliegenden Arbeit ist die Entwicklung und Implementierung neu- er Simulationsalgorithmen in das CAST Programmpaket. Neben der teilweisen Para- llelisierung der bereits impelentierten Kraftfelder wurde das Programm um einen Mole- kulardynamikcode erweitert. Aufbauend auf diesem Code wurden Algorithmen zur Be- rechnung der freien Energie entlang einer Reaktionskooridnate, sowie eine Erweiter-ung der Behandlung elektrostatischer Wechselwirkungen auf Basis einer Ewald Summation implementiert. Der zeitaufwändigste Teil einer Kraftfeldrechnung stellt die Evaluierung der nichtbin- denden Wechselwirkungen dar. Die Berechnung dieser Wechselwirkungen wurde für die Nutzung von Mehrkernprozessoren optimiert und parallelisiert. Die Parallelisie- rung zeigte eine signifikante Reduktion der benötigten Rechenzeit auf mehreren Re- chenkernen im Vergleich zur seriellen Berechnung auf nur einem Rechenkern für einfa- che Energie- bzw. Gradientenrechnungen sowie für Molekulardynamikrechnungen. Um Rechnungen, die einen cutoff Radius benutzen, weiter zu beschleunigen, wurde der Auf- bau der Verlet-Liste modifiziert. Statt des Standardalgorithmus, der eine Doppelschleife über alle Atome verwendet, wurde ein linked-cell Algorithmus implementiert. Der Auf- bau der Verlet-Liste konnte damit um den Faktor 7 beschleunigt werden. Der Molekulardynamikcode enthält mehrere Algorithmen zur Berechnung von Syste- men in verschiedenen Ensembles. Zur numerischen Integration der Bewegungsgleichun- gen wurden der Velocity-Verlet sowie eine modifizierte Version von Beemans Algorith- mus implementiert. Da der minimale Code, wie er in Kapitel 1.2 beschrieben wird, ein mikrokanonisches Ensemble produziert, und dieses keiner realistischen experimentel- len Situation entspricht, wurden Methoden zur Berechnung und Aufrechterhaltung von Temperatur und Druck implementiert. Die Temperatur kann mittels zweier verschiede- ner Möglichkeiten kontrolliert werden. Die erste Möglichkeit ist die direkte Skalierung der Geschwindigkeiten der Partikel, die zweite Möglichkeit besteht in der Nutzung ei- nes Nòse-Hoover Thermostaten, der ein echtes kanonisches Ensemble generiert. Für die Kontrolle des Drucks wurde ein Berendsen Barostat implementiert. Da die Kontrolle des Drucks die Nutzung von periodischen Randbedingungen voraussetzt, ist die Form der Simulationszelle wichtig. CAST unterstützt aktuell isotrope und anisotrope Simulationszellen, mit der Einschränkung, dass alle Winkel 90◦betragen. Ein kritischer Punkt bei einer MD Simulation ist die Länge des Zeitschritts, der in der Regel bei 1fs liegt. Sollen größere Zeitschritte verwendet werden, müssen die schnell- sten Bewegungnen im System eingeschränkt werden. Dies sind im Normalfall die H-X Streckschwingungen. Zur Einschränkung dieser wurde der RATTLE Algorithmus imple- mentiert der die H-X Bindung mit Hilfe von Lagrange-Multiplikatoren auf den Gleich- gewichtsabstand fixiert. Als erste Anwendung des MD Codes wurde in Kombination mit dem MOPAC Programm die Rigidität und Flexibilität von PBI Seitenketten erfolgreich untersucht. Basierend auf dem MD Code wurden zwei Möglichkeiten zur Bestimmung der freien Energie eingebaut, Umbrella Sampling und Free Energy Perturbation. Umbrella Samp- ling erlaubt die Bestimmung der freien Energie entlang einer Reaktionskoordinate, hier Abstände oder Diederwinkel. Der Algorithmus wurde erfolgreich an zwei Literatur- beispielen, der Streckung von Deca-L-Alanin sowie der Rotation von Butan um den zentralen Diederwinkel getetstet. Die FEP Implementierung wurde an drei Beispielen getestet, einer Nullsummen-Transformation von Ethan in Ethan in explizitem Solvent, dem Lösen eines Natriumions in Wasser und der Transformation von Tyrosin in Alanin in einem Tripeptid. Die Ergebnisse dieser Testrechnungen stimmen hervorragend mit Vergleichsrechnungen mit NAMD sowie Literaturwerten überein. Die FEP Methode wurde schließlich benutzt um die relative freie Bindungsenergie zweier Inhibitoren in einem Inhibitor-Protein-Komplex zu bestimmen. Neben Kraftfeldern können auch ab-initio Methoden für Simulationen benutzt werden. Da die Rechenzeit dieser Methoden um ein vielfaches höher ausfällt als die für Kraftfel- der, ist die Benutzung für die globale Optimierung jedoch limitiert. In den letzten Jah- ren wurde im Bereich der Leistungsfähigkeit dieser Methoden jedoch große Fortschritte erzielt, indem diese Methoden auf Grafikkarten portiert wurden. Um diese Entwick- lung nutzbar zu machen wurde eine MPI-Schnittstelle in CAST implementiert, die mit dem DFT Programm TeraChem kommuniziert. Die Kombination aus beiden Program- men, sowie die Funktionsfähigkeit der Schnittstelle, wurde an H2O10 Clustern sowie dem Polypeptid Met-Enkephalin getestet. Die reinen ab-initio Rechnungen zeigten ein besseres Verhalten gegenüber dem Normalen Protokoll, welches Kraftfeldrechungen mit nachfolgender Optimierung mit qunatenchemischen Methoden vorsieht. KW - Molekulardynamik KW - Molecular dynamics KW - Molecular mechanics KW - Molecular Simulation KW - Free Energy Perturbation (1987 : Princeton, NJ) Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-132032 ER -