TY - THES A1 - Schmid, Paul T1 - Quantenchemische Untersuchungen von Umgebungseinflüssen bei offen- und geschlossenschaligen Systemen T1 - Quantum chemical studies of environmental effects in open-shell and closed-shell systems N2 - In dieser Dissertation werden die Umgebungseinflüsse auf die strukturellen und elektronischen Eigenschaften von verschiedenen offen- und geschlossenschaligen Systemen mittels quantenchemischer Methoden berechnet. Ein Kernpunkt umfasst die Untersuchung von verdreht angeordneten, biradikalischen Diborylalkenen, welche eine ungesättigte C2R2-Brücke (R = Et, Me) besitzen und durch cyclische (Alkyl)(amino)carbene (CAACs) stabilisiert werden. Quantenchemische Berechnungen zeigen, dass hauptsächlich sterische Effekte für die Ausbildung einer verdrehten Molekülanordnung verantwortlich sind, während bei geringen sterischen Wechselwirkungen (R = H) die Delokalisationseffekte überwiegen, wodurch eine planare Struktur begünstigt wird. Die Bevorzugung einer offenschaligen Singulettkonfiguration anstelle eines Tripletts ist auf den großen Energieunterschied der beiden einfach besetzten Molekülorbitale zurückzuführen. Durch die Berechnung der Lösungsmitteleffekte mithilfe von polarisierbaren Kontinuumsmodellen kann gefolgert werden, dass mit zunehmender statischer Dielektrizitätskonstante eine planare und geschlossenschalige Struktur stärker stabilisiert wird als eine verdrehte Anordnung. Ein weiteres Thema dieser Dissertation befasst sich mit der quantenchemischen Analyse eines makrozyklischen Perylenbisimid-Trimersystems, welches eingebettet in einer Polymethylmethacrylat-Matrix bei Temperaturen nahe dem absoluten Nullpunkt eine Lokalisierung der ersten drei angeregten Zustände zeigt. Quantenchemische Vakuumberechnungen ergeben, dass unabhängig von der gegenseitigen geometrischen Orientierung der drei Perylenbisimid-Chromophore der Übergang vom Grundzustand in den S1-Zustand verboten ist und dass die ersten drei angeregten Zustände delokalisiert vorliegen. Mithilfe von expliziten Lösungsmittelmodellen kann jedoch gezeigt werden, dass das Auftreten dieser Lokalisierungen auf eine inhomogene Polymethylmethacrylat-Umgebung zurückzuführen ist, die zu einem Symmetriebruch und somit zu einer Zunahme der Oszillatorstärke für S1 und der Lokalisierungsgrade für S1, S2 und S3 führt. Darüber hinaus wird der Lösungsmitteleinfluss auf die angeregten Zustände des Azulens mittels impliziter und expliziter Lösungsmittelmodelle berechnet. Bei einer Erhöhung der dynamischen Dielektrizitätskonstante im impliziten Modell nehmen die Anregungsenergien der vertikalen Singulettzustände ab, wobei der Effekt mit steigender Oszillatorstärke zunimmt. Die Auswirkung der statischen Dielektrizitätskonstante auf die Anregungsenergien ist dagegen deutlich schwächer ausgeprägt. Im expliziten Modell bewirkt das Lösungsmittel ebenfalls eine Abnahme der Anregungsenergie des hellen Singulettzustands, wenn auch in geringerem Umfang als im impliziten Modell. Als letztes Thema wird der Inhibitionsmechanismus der Cysteinprotease Rhodesain durch zwei modifizierte 1,4-Naphthoquinone untersucht. Während beide Naphthoquinone an der 2-Position eine Dipeptideinheit aufweisen, besitzen sie an der 3-Position entweder einen Nitril- oder Chloridsubstituenten. Zwar erfolgt bei beiden Derivaten die Inhibition über einen kovalent-reversiblen Mechanismus, jedoch verläuft die Hemmung im Falle des Nitrilderivats erheblich effektiver. Die quantenchemischen Berechnungen eines vereinfachten Modells zeigen, dass die Cysteineinheit (HS-R) bevorzugt in einer exothermen und reversiblen Additionsreaktion an die elektronenarme C-C-Doppelbindung der Naphthoquinone anlagert. Dabei werden kleinere Reaktionsenergien für die Reaktion des Chlorderivats als für die Reaktion des Nitrilderivats erhalten. Durch die Berücksichtigung von Wasser in einem impliziten Lösungsmittelmodell kommt es bei fast allen Reaktionsprodukten zu einer Energiezunahme, die bei der Reaktion des Nitrilderivats stärker ausfällt als bei der Reaktion des Chlorderivats. N2 - In this thesis, the environmental effects on the structural and electronic properties of various open-shell and closed-shell systems are calculated using quantum chemical methods. A key issue involves the investigation of twisted, biradical diborylalkenes, which contain an unsaturated C2R2 bridge (R = Et, Me) and are stabilized by cyclic (alkyl)(amino)carbenes (CAACs). Quantum chemical calculations show that mainly steric effects are responsible for the formation of a twisted molecular arrangement, while delocalization effects predominate at low steric interactions (R = H), favoring a planar structure. The preference for an open-shell singlet configuration over a triplet is due to the large energy difference between the two singly occupied molecular orbitals. By calculating the solvent effects using polarizable continuum models, it can be concluded that with an increasing static dielectric constant, a planar and closed-shell structure is more stabilized compared to a twisted conformation. Another topic of this thesis deals with the quantum chemical analysis of a macrocyclic perylene bisimide trimer system, which is embedded in a poly(methyl methacrylate) matrix and shows a localization of the first three excited states at temperatures near absolute zero. Quantum chemical vacuum calculations reveal that regardless of the mutual geometric orientation of the three perylene bisimide chromophores, the transition from the ground state to the S1 state is forbidden and that the first three excited states are delocalized. However, with the help of explicit solvent models, it can be shown that the occurrence of these localizations is caused by an inhomogeneous poly(methyl methacrylate) environment, which leads to symmetry breaking and thus to an increase in the oscillator strength for S1 and the degrees of localization for S1, S2, and S3. In addition, the influence of the solvent on the excited states of azulene is calculated using implicit and explicit solvent models. With an increase in the dynamic dielectric constant in the implicit model, the excitation energies of the vertical singlet states decrease. This effect becomes stronger with rising oscillator strength. In contrast, the effect of the static dielectric constant on the excitation energies is much weaker. In the explicit model, the solvent also causes a decrease in the excitation energy of the bright singlet state, although the extent is smaller than in the implicit model. As a final topic, the inhibition mechanism of the cysteine protease rhodesain by two modified 1,4-naphthoquinones is investigated. While both naphthoquinones have a dipeptide unit at the 2-position, they have either a nitrile or chlorine substituent at the 3-position. Although for both derivatives the inhibition takes place via a covalent reversible mechanism, the inhibition in the case of the nitrile derivative is considerably more effective. The quantum chemical calculations of a simplified model show that the cysteine moiety (HS-R) preferentially attaches to the electron-deficient C-C double bond of the naphthoquinones in an exothermic and reversible addition reaction. Smaller reaction energies are obtained for the reaction of the chlorine derivative than for the reaction of the nitrile derivative. By considering water in an implicit solvent model, the energies of almost all reaction products rise, whereby the energetic increase for the reaction of the nitril derivative is greater than for the reaction of the chlorine derivative. KW - Umgebungseinfluss KW - Lösungsmitteleffekt KW - Biradikal KW - Berechnung KW - CAAC KW - Diborylalkene KW - Quantenchemie Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-265106 ER - TY - THES A1 - Kaiser, Dustin T1 - Non-standard computational approaches applied to molecular systems T1 - Unkonventionelle Berechnungsansätze für molekulare Systeme N2 - In this thesis, several contributions to the understanding and modeling of chemical phenomena using computational approaches are presented. These investigations are characterized by the usage of non-standard computational modeling techniques, which is necessitated by the complex nature of the electronic structure or atomic fluctuations of the target molecules. Multiple biradical-type molecules and their spectroscopic properties were modeled. In the course of the investigation, it is found that especially the impact of correct molecular geometries on the computationally predicted absorption properties may be critical. In order to find the correct minimum geometries, Multi-Reference methods may have to be invoked. The impact of geometry relaxation on the excitonic properties of Perylene Bisimide dimers were investigated. Oftentimes, these geometry factors are neglected in Organic Semiconductor modeling as an approximation. This present investigation suggests that this approximation is not always valid, as certain regimes are identified where geometrical parameters have critical impact on the localization and energetic properties of excitons. The mechanism of the Triazolinedione (TAD) tyrosine bioconjugation reaction is investigated using quantum-chemical methods. By comparison of different conceivable mechanisms and their energetic ordering, the TAD tyrosine bioconjugation is found to proceed by means of a base-mediated electrophilic aromatic substitution reaction. The kth nearest neighbor entropy estimation protocol is investigated. This estimator promises accurate entropy estimates even for flexible molecules with multiple structural minima. Our granular investigation of formal and practical properties of the estimator suggests that the uneven variance of a molecule’s vibrational modes is the cause of the observed slow convergence of the estimator. A rescaling procedure to reestablish fast convergence is suggested and benchmarks are performed. N2 - Im Rahmen dieser Arbeit wurden Berechnungsansätze für die Modellierung mit Standardverfahren näherungsweise berechenbarer molekularer und spektroskopischer Probleme diskutiert. Zunächst wurden die spektroskopischen Eigenschaften der biradikalischen Moleküle Methylbismut, Diphenylacetylen, Pentadiynylidin (sowie dessen methylsubstituierte Derivate) und Diphenylpropynyliden untersucht. Diese Forschungsvorhaben wurden in Zusammenarbeit mit experimentell arbeitenden physi kalischen und synthetischen Chemikern durchgeführt. Der Vergleich von sowohl spek troskopisch als auch durch theoretische Rechnungen erlagter Erkenntnisse erlaubt eine genaue Charakterisierung der physiko-chemischen Eigenschaften der Moleküle. Mit Hilfe von computergestützter Modellierung wurden die Ionisierungsenergien und die Absorptionsspektren der Ionisierung in der Franck-Condon-Approximation berechnet. Für Methylbismut wurden skalar-relativistische und Pseudopotenzial-basierte Ansätze zur Berechnung verwendet. Durch Analyse der angeregten elektronischen Zustände und ihrer geometrieabhängigen energetischen Eigenschaften, konnte das transiente Absorptionsspektrum von Diphenylpropynyliden genau charakterisiert und rational isiert werden. ... KW - Quantenchemie KW - Theoretische Chemie KW - Entropie KW - Computational Chemistry KW - Quantum Chemistry Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-276641 ER - TY - THES A1 - Grüne, Marvin T1 - Solid-state NMR Spectroscopic, X-Ray Diffraction and Quantum Chemical Investigations of the Crystalline Cancer Drug Paclitaxel and Paclitaxel incorporated into Polymer Micelles T1 - Festkörper-NMR-, Röntgendiffraktometrie- und quantenchemische Untersuchungen des kristallinen Krebs-Wirkstoffs Paclitaxel und Paclitaxel eingebettet in Polymermizellen N2 - Paclitaxel (PTX) is one of the leading drugs against breast and ovarian cancer. Due to its low solubility, treatment of the patients with this drug requires a very well-suited combination with a soluble pharmaceutical excipient to increase the bioavailability and reduce the strong side ef-fects. One efficient way to achieve this in the future could be the incorporation of PTX into pol-ymeric micelles composed of poly(2-oxazoline) based triblock copolymers (POL) which ena-bles PTX loadings of up to 50 wt.%. However, structural information at an atomic level and thus the knowledge of interaction sites within these promising but complex PTX-POL formula-tions were not yet available. Such results could support the future development of improved excipients for PTX and suitable excipients for other pharmaceutical drugs. Therefore, a solid-state MAS NMR investigation of these amorphous formulations with different POL-PTX com-positions was performed in this thesis as this gives insights of the local structure at an atomic level in its solid state. NMR in solution showed very broad 13C signals of PTX for this system due to the reduced mobility of the incorporated drug which exclude this as an analytical meth-od. In a first study, crystalline PTX was structurally characterized by solid-state NMR as no com-plete 13C spectrum assignment and no 1H NMR data existed for the solid state. In addition, the asymmetric unit of the PTX crystal structure consists of two molecules (Z'=2) that can only be investigated in its solid state. As crystalline PTX in total has about 100 different 13C and 1H chemical shifts with very small differences due to Z’=2, and furthermore, its unit cell consisting of more than 900 atoms, accompanying GIPAW (CASTEP) calculations were required for NMR signal assignments. These calculations were performed using the first three available purely hydrous and anhydrous PTX structures, which were determined by XRD and published by Vel-la-Zarb et al. in 2013. Within this thesis, is was discovered that two investigated batches of commercially available PTX from the same supplier both contained an identical and so far un-known PTX phase that was elucidated by PXRD as well as solid-state NMR data. One of the two batches consists of an additional phase that was shown to be very similar to a known hy-drated phase published in 2013.[1] By heating the batch with the mixture of the two phases un-der vacuum, it is transformed completely to the new dry phase occurring in both PTX batches. Since the drying conditions to obtain anhydrous PTX in-situ on the PXRD setup described by Vella-Zarb et. al.[1] were much softer than ours, we identify our dry phase as a relaxed version of their published anhydrate structure. The PXRD data of the new anhydrate phase was trans-ferred into a new structural model, which currently undergoes geometry optimization. Based on solid-state NMR data at MAS spinning frequencies up to 100 kHz, a 13C and a partial 1H signal assignment for the new anhydrous structure were achieved. These results provided sufficient structural information for further investigations of the micellar POL-PTX system. In a second study, the applicability and benefit of two-dimensional solid-state 14N-1H HMQC MAS NMR spectra for the characterization of amorphous POL-PTX formulations was investi-gated. The mentioned technique has never been applied to a system of similar complexity be-fore and was chosen because around 84% of the small-molecule drugs contain at least one nitrogen atom. In addition, the number of nitrogen atoms in both POL and PTX is much smaller than the number of carbons or hydrogens, which significantly reduces the spectral complexity. 14N has a natural abundance of 99.6% but leads to quadrupolar broadening due to its nuclear spin quantum number I = 1. While this is usually undesirable due to broadening in the resulting 1D 14N NMR spectra, this effect is explicitly used in the 2D 14N-1H HMQC MAS experiment. The indirect 14N measurement can avoid the broadening while maintaining the advantage of the high natural abundance and making use of the much more dispersed signals due to the additional quadrupolar shifts as compared to 15N. This measurement method could be successfully applied to the complex amorphous POL-PTX mixtures. With increasing PTX loading of the formulations, additional peaks arise as spatial proximities of the amide nitrogens of POL to NH or OH groups of PTX. In addition, the 14N quadrupolar shift of these amide nitrogens decreases with increasing PTX content indicating a more symmetric nitrogen environment. The latter can be explained by a transformation of the trigonal planar coordination of the tertiary amide nitrogen atoms in pure POL towards a more tetrahedral environment upon PTX loading induced by the formation of hydrogen bonds with NH/OH groups of PTX. In the third and last project, the results of the two abovementioned studies were used and ex-tended by solid state 13C and two-dimensional 1H-13C as well as 1H-1H MAS NMR data with the aim to derive a structural model of the POL-PTX formulations at an atomic level. The knowledge of the NMR signal assignments for crystalline PTX was transferred to amorphous PTX (present in the micelles of the formulations). The 13C solid-state NMR signals were evalu-ated concerning changes in chemical shifts and full widths of half maximum (FWHM) for the different PTX loadings. In this way, the required information about possible interaction sites at an atomic level becomes available. Due to the complexity of these systems, such proximities often cannot be assigned to special atoms, but more to groups of atoms, as the individual de-velopments of line widths and line shifts are mutually dependent. An advantageous aspect for this analysis was that pure POL already forms unloaded micelles. The evaluation of the data showed that the terminal phenyl groups of PTX seem to be most involved in the interaction by the establishment of the micelle for lowest drug loading and that they are likely to react to the change in the amount of PTX molecules as well. For the incorporation of PTX in the micelles, the following model could be obtained: For lowest drug loading, PTX is mainly located in the inner part of the micelles. Upon further increasing of the loading, it progressively extends to-ward the micellar shell. This could be well shown by the increasing interactions of the hydro-phobic butyl chain of POL and PTX, proceeding in the direction of the polymer backbone with rising drug load. Furthermore, due to the size of PTX and the hydrodynamic radius of the mi-celles, even at the lowest loading, the PTX molecules partially reach the core-shell interface of the micelle. Upon increasing the drug loading, the surface coverage with PTX clusters increas-es based on the obtained model approach. The latter result is supported by DLS and SANS data of this system. The abovementioned results of the 14N-1H HMQC MAS investigation of the POL-PTX formulations support the outlined model. As an outlook, the currently running geometry optimization and subsequently scheduled calcu-lation of the chemical shieldings of the newly obtained anhydrous PTX crystal structure can further improve the solid-state NMR characterization through determination of further spatial proximities among protons using the existing 2D 1H(DQ)-1H(SQ) solid-state MAS NMR spec-trum at 100 kHz rotor spinning frequency. The 2D 14N-1H HMQC MAS NMR experiments were shown to have great potential as a technique for the analysis of other disordered and amor-phous drug delivery systems as well. The results of this thesis should be subsequently applied to other micellar systems with varying pharmaceutical excipients or active ingredients with the goal of systematically achieving higher drug loadings (e.g., for the investigated PTX, the similar drug docetaxel or even different natural products). Additionally, it is planned to transfer the knowledge to another complex polymer system containing poly(amino acids) which offers hy-drogen bonding donor sites for additional intermolecular interactions. Currently, the POL-PTX system is investigated by further SANS studies that may provide another puzzle piece to the model as complementary measurement method in the future. In addition, the use of MD simu-lations might be considered in the future. This would allow a computerized linking of the differ-ent pieces of information with the aim to determine the most likely model. N2 - Paclitaxel (PTX) ist eines der führenden Medikamente gegen Brust-und Eierstockkrebs. Aufgrund seiner geringen Löslichkeit erfordert die Behandlung der Patienten mit diesem Medikament eine sehr gut geeignete Kombination mit einem löslichen pharmazeutischenHilfsstoff, um die Bioverfügbarkeit zu erhöhen und die starken Nebenwirkungen zu reduzieren. Ein effizienter Weg, dies in Zukunft zu erreichen, könnte der Einbau von PTX in polymere Mizellen sein, die aus Poly(2-oxazolin)-basierten Triblock-Copolymeren (POL) bestehen und PTX-Beladungen von bis zu 50 Gew.-% ermöglichen. Strukturelle Informationen auf atomarer Ebene und damit die Kenntnis von Wechselwirkungeninnerhalb dieser vielversprechenden, aber komplexen PTX-POL-Formulierungen waren jedoch bisher nichtverfügbar. Solche Ergebnisse könnten die zukünftige Entwicklung von verbesserten Hilfsstoffen für PTX und von geeigneten Hilfsstoffen für andere pharmazeutische Wirkstoffe unterstützen. Aus diesem Grund wurdenin der vorliegenden DissertationFestkörper-NMR-Untersuchungen andiesenamorphen Formulierungen mit unterschiedlichen POL-PTX Zusammensetzungen durchgeführt, weil damit Einblickein die lokale Struktur auf atomarer Ebene im festen Zustand erhalten werden können. Aufgrund der verringerten Mobilität des eingebrachten Wirkstoffs in diesem System ergeben NMR-Messungen in Lösung sehr breite 13C-PTX-Signale, was diese Technikals Analysemethode ausschließt. ... KW - Wirkstoff-Träger-System KW - NMR-Spektroskopie KW - Röntgendiffraktometrie KW - Taxol KW - Quantenchemie KW - Solid-State NMR Spectroscopy KW - X-Ray Diffraction KW - Quantum Chemical Calculations KW - Drug Delivery System KW - Taxol KW - Festkörper-NMR KW - quantenchemische Berechnungen Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-237199 ER - TY - THES A1 - Hader, Kilian T1 - Lokalisierungsdynamik unter Berücksichtigung von Molekül-Feld-Wechselwirkung, Kern-Elektron-Kopplung und Exziton-Exziton-Annihilierung T1 - Localization dynamics considering molecule-field interaction, nuclear-electron coupling and exciton-exciton annihilation N2 - Diese Arbeit befasst sich mit verschiedenen Aspekten der Dynamik von Kernen, Elektronen und gekoppelten Kern-Elektron-Systemen, wobei je nach System unterschiedliche Herangehensweisen gewählt wurden. Zentrale Punkte sind bei allen drei Kapiteln einerseits die Lokalisierung von Teilchen und Energie und andererseits eine hohe Sensitivität in Bezug auf die Wahl der Anfangsbedingungen. Im ersten Teil wurden von der Carrier-Envelope-Phase (CEP) abhängende, laser-induzierte Lokalisierungen betrachtet. Das zentrale Element ist dabei das entwickelte Doppelpulsschema, mit welchem eine CEP-Abhängigkeit in beobachtbaren Größen erzeugt wird. Als Beispielsysteme wurden die Fragmentation im D₂⁺-Modellsystem und eine Isomerisierung im Doppelminimumpotential (DMP) untersucht. Als Observable wird die Asymmetrie betrachtet Im DMP kann die Asymmetrie mit dem Entantiomeren/Isomerenüberschuss gleich gesetzt werden kann und im D₂⁺-Modellsystem mit der Lokalisierung des Elektrons auf einem der beiden dissoziierenden Kerne. Eine Phasenabhängigkeit der Asymmetrien besteht nur für die CEP des zweiten Pulses φ₂, für welchen keine Begrenzungen für die Anzahl an Laserzyklen auftreten. Im DMP wurde die CEP-Abhängigkeit der Asymmetrien auch bei unterschiedlichen Startkonfigurationen untersucht. Für alle untersuchten Startkonfigurationen konnte ein Laserparametersatz gefunden werden, der für zumindest eine der beiden Asymmetrien eine CEP-Abhängigkeit liefert. Aufgrund der aufgehobenen energetischen Entartung der Paare gerader und ungerader Symmetrie ist die resultierende Lokalisierung zeitabhängig. Zur Messung der vorhergesagten Dynamiken ist z.B. die Aufnahme eines Photoelektronen-Spektrums denkbar. In nächsten Kapitel wurden unterschiedliche Dynamiken innerhalb eines 4d Kern-Elektron-Modells in der Nähe einer konischen Durchschneidung (CI) zweier Potentiale betrachtet. Hierbei ist hervorzuheben, dass eine solche gleichzeitige Untersuchung von Kern- und Elektron-Dynamik in Systemen mit CIs in der Literatur, nach Wissen des Autors, bisher nicht veröffentlicht ist. Das 4d-Potential wurde mit Hilfe des sogenannten Potfit-Algorithmus gefittet. Dieser Fit wurde anschließend verwendet, um die Dynamik des gekoppelten Systems mit Hilfe der ”Multi-Configuration Time-Dependent Hartree”(MCTDH)-Methode zu berechnen. Aus der Analyse der gekoppelten Kern-Elektron-Wellenfunktion ergaben sich zwei grundlegend unterschiedliche Klassen von Dynamiken: • Diabatisch: Kern- und Elektrondynamik sind nahezu entkoppelt. Der Kern bewegt sich und das Elektron bleibt statisch. • Adiabatisch: Kern- und Elektrondynamik sind stark gekoppelt. Die Kerndynamik findet auf Kreisbahnen statt. Mit der Rotation der Kerndichte um den Winkel φ geht eine Rotation der Elektron-Dichte einher. Die diabatische Bewegung entspricht der Dynamik durch die konische Durchschneidung und die adiabatische Bewegung der Dynamik auf der unteren Potentialfläche. Welche der beiden Dynamiken stattfindet, wird durch die Wahl der Anfangsbedingung bestimmt. Der wesentliche Unterschied zwischen den beiden Startzuständen ist dabei die Lage des Knotens im elektronischen Anteil der Wellenfunktion. In den diabatischen Bewegungen bleibt z.B. der pₓ -artige Charakter der elektronischen Wellenfunktion konstant, wohingegen sich bei der adiabatischen Dynamik der Charakter mit der Kernbewegung ändert. Die Zeitersparnis durch die Verwendung des MCTDH-Ansatzes im Vergleich zur Split-Operator-Methode liegt etwa bei einem Faktor 5. Das letzte Kapitel widmet sich der mikroskopischen Beschreibung von Exziton-Exziton- Annihilierung (EEA). Dabei werden numerische Lösungen der aus einem mikro- skopischen Modell hergeleiteten Ratengleichungen mit Messungen ( transienter Absorption) verglichen. Es wurden zwei Systeme untersucht: ein Squarain-basiertes Heteropolymer (SQA-SQB)ₙ und ein [2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenvinylen]-Polymer, auch bekannt als MEH-PPV. In beiden Fällen gelang die systematische Parameterbestimmung mit Hilfe einer Aufteilung in lokalisierte Subsysteme. Diese Subsysteme werden einzeln gewichtet und anschließend aufsummiert, wobei die Gewichte optimiert werden können. Aus den so erhaltenen Parametern ergibt sich für beide Systeme ein ähnliches Bild: • Durch ultraschnelle Lokalisierung der Anregung im fs-Bereich auf kleinere Aggregateinheiten bilden sich voneinander getrennte Subsysteme. • Die in den Subsystemen lokalisierten Exzitonen können sich nur innerhalb dieser Bereiche frei bewegen. Es ist ausreichend, direkt benachbarte Mono-, Bi-, Tri- und Tetra-Exzitonen in bis zu zwei Dimensionen zu berücksichtigen. • Auf einer fs-Zeitskala annihilieren direkt benachbarte Exzitonen. • Im MEH-PPV ergibt sich der Signalzerfall im fs-Bereich als Mittelwert aus einer schnellen (zwischen Ketten) und einer langsamen (innerhalb von Ketten) Annihilierung. • Im ps- bis ns-Bereich wird sowohl durch Diffusion vermittelte Annihilierung, also auch der Zerfall der ersten angeregten Zustände bedeutsam. N2 - In the present work the dynamics of nuclei, electrons, and coupled nuclei-electron systems are examined in different ways. Items that are central in all three chapters are, on the one hand localization of particles and energy and, on the other, a high sensitivity to the choice of initial condition. In the first chapter carrier-envelope-phase (CEP) dependent, laser induced localization is examined. The main element of the considerations is a double pulse scheme, which creates a CEP-dependence in the monitored observables. As example systems the fragmentation of a D₂⁺-model and the isomerization in a double well potential (DWP) are investigated. As an observable the asymmetry is chosen. In the DWP this entity can be related to enantiomeric or isomeric excess and in the D₂⁺-model it describes the localization of the electron on a fragment. The phase dependent part of the asymmetries only relies on the CEP φ₂ of the second pulse which does not have any restrictions on the amount of laser cycles. In the DWP a CEP-dependence of the asymmetries could be examined starting from different initial configurations. For all different initial conditions a set of laser parameters could be found which produces at least one CEP-dependent asymmetry. Due to the removed degeneracy between states of even and odd parity, the resulting localization in the left or right potential well is time-dependent. - fied such that the first pulse does not populate all states equally. A possible way to test the predicted behavior experimentally is the measurement of photo-electron spectra. In the next chapter coupled nucleus-electron-dynamics in the vicinity of a conical intersections (CI) of two potentials are investigated in a 4-d model system.Such examinations of coupled nucleus-electron-dynamics on equal footing in systems containing CIs is, to the author’s knowledge, not published in literature. The 4-d potential has been fitted by use of the so-called Potfit-algorithm which subsequently could be used to calculate the dynamics of the coupled system in the ”Multi-Configuration Time-Dependent Hartree”(MCTDH)-framework. The analysis of the coupled nucleus-electron-wavefunction yielded two fundamentally different classes of dynamics: • Diabatic: Nucleus- and electron dynamics are nearly uncoupled. The nucleus moves while the electron remains static. • Adiabatic: Nucleus- and electron dynamics are strongly coupled. The dynamic of the nucleus takes place on orbits. The rotation of the nuclear density by the angle φ is accompanied by a rotation of the electron-density at the same angle. The diabatic dynamics are present if the wave packet is passing through the conical intersection and the adiabatic dynamics can be attributed to a wave packet moving on the lower potential surface. Which of the two classes of dynamics takes place can be controlled by choice of the initial wavefunction. The most significant difference between the two initial wavefunctions is the plane in which the node of the electronic wavefunction is located. In case of a diabatic motion the pₓ -like character of the wavefunctions remains constant, while in case of a adiabatic motion the character changes with the motion of the nucleus.The time saving by usage of the MCTDH-method in comparison to the Split-Operator-method is about a factor of 5. The last chapter is dedicated to the microscopic description of exciton-exciton-annihilation (EEA). The numerical solution of the rate equations which are derived starting from a microscopic Hamiltonian, are compared with measurements. The experimental data are time-dependent traces of transient absorption measurements at different laser fluences which were available for two different systems:a squaraine-based copolymer (SQA-SQB)ₙ and a [2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene] polymer also know as MEH-PPV. In both cases a systematic parameter determination could be achieved by introduction of localized subsystems. These subsystems are weighted independently and are summed up whereby the weighs can be optimized. The resulting interpretation of the obtained parameters is similar for both systems: • Ultrafast localization of the excitation energy takes place in the fs-regime which leads to excitons residing on smaller subsystems. • Excitons in these subsystems can only move inside of these domains. A re- construction of experimental data is feasible by inclusion of mono-, bi-, tri- and tetra-excitons in up to two dimensions. • In the fs-regime neighbouring excitons annihilate • In the MEH-PPV polymer the signal decay in the fs-regime can be described as the average of a fast annihilation (between chains) and a slow annihilation (inside chains). •On a longer time-scale (ps to ns) diffusion-meditated annihilation and decay of the first excited states take place KW - Quantenmechanik KW - Quantenchemie KW - Laserstrahlung KW - Nichtadiabatischer Prozess KW - Exziton KW - multicycle CEP control KW - exciton exciton annihilation KW - exact conical intersection dynamics KW - coupled nuclear-electron MCTDH KW - quantum dynamics KW - Quantentheoretische Chemie Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-146735 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 - Settels, Volker T1 - Quantum chemical description of ultrafast exciton self-trapping in perylene based materials T1 - Quanten-chemische Beschreibung von ultraschnellem Self-trapping von Exzitonen in Perylen-basierten Materialien N2 - Im Rahmen dieser Dissertation wurden sehr lange Exzitonen-Diffusionslängen (LD) unter idealen Bedingungen für Perylen-basierte Materialien simuliert. Dies ist ein Indiz dafür, dass die sehr kurzen LD in realen Materialien aus einer extrinsischen sowie einer intrinsischen Immobilisierung resultieren. Letztere basiert auf einer Relaxation in sogenannten „Self-Trapping“-Zustände. Ein tieferes Verständnis der dem Self-Trapping zugrunde liegenden atomistischen Prozesse ist notwendig, um zukünftig Materialien mit langen LD entwickeln zu können, bei denen eine intrinsische Exzitonen-Immobilisierung verhindert wird. Für die Entwicklung eines solchen mechanistischen Verständnisses ist das Vorliegen einer eindeutigen Korrelation zwischen der molekularen Anordnung und der LD unabdingbar. Diese weisen Einkristalle von Diindenoperylen (DIP) und α-Perylen-tetracarboxyl-anhydrid (α-PTCDA) auf. Bei ersteren wurde eine außergewöhnlich lange LD von 90 nm und bei letzteren nur 22 nm gemessen. Teil dieser Arbeit war es, Gründe für diesen Unterschied in der LD zu finden. Nur Self-Trapping kommt als Ursache in Frage. Aus diesem Grund eignen sich diese Materialien, um ein atomistisches Verständnis des Self-Trappings exemplarisch an ihnen zu erarbeiten. Mutmaßlich könnten Differenzen in der elektronischen Struktur in DIP und α-PTCDA für das unterschiedliche Self-Trapping verantwortlich sein. Allerdings konnte gezeigt werden, dass es für viele Perylen-basierte Materialien keine signifikanten Unterschiede in der elektronischen Struktur gibt, wodurch diese für die Aufklärung von Immobilisierungsmechanismen zu vernachlässigen sind. Eine weitere mögliche Begründung wäre in Polarisationseffekten im Kristall zu suchen, welche die elektronische Struktur in Perylen-basierten Materialien unterschiedlich beeinflussen. Vor allem ihr Einfluss auf Ladungstrennungs-Zustände (CT), die oberhalb des optisch hellen Frenkel-Zustandes liegen, war fraglich, weil sie energetisch abgesenkt werden könnten. Ein signifikanter Einfluss von Polarisationseffekten konnte aber für alle Zustände mittels eines polarisierbaren Kontinuum-Modells ausgeschlossen werden. Die geringe LD im α-PTCDA ist folglich ein Indiz für ein Self-Trapping, das durch die Kristallstruktur aus π-Stapeln evoziert wird, welche in DIP fischgrätenartig ist. Da Polarisationseffekte auszuschließen sind, übt der Kristall lediglich durch sterische Restriktionen einen Einfluss auf das Dimer aus. Daher muss die Methode für die Beschreibung von Self-Trapping nur diese Effekte berücksichtigen, so dass sich für den Einsatz des mechanical embedding QM/MM-Ansatzes entschieden wurde. Nun konnten Potentialflächen berechnet werden, auf denen anschließend eine Wellenpaketdynamik durchgeführt wurde. Diese Methode erlaubt es erstmals, Mechanismen der Exzitonen-Immobilisierung in organischen Materialien auf einer atomistischen Ebene zu beschreiben. Als Erklärung für Self-Trapping in α-PTCDA dienten Potentialflächen, die eine intermolekulare Verschiebung des Dimers im Kristall abbilden. So wurde eine Exzitonen-Immobilisierung innerhalb von 500 fs gefunden, die aus einem irreversiblem Energieverlust und einer lokalen Verzerrung der Kristallstruktur resultiert und auf diese Weise den weiteren Transport des Exzitons verhindert. Im Fall von DIP kann diese Immobilisierung aufgrund hoher Energiebarrieren nicht stattfinden. Diese Barrieren resultieren aus der fischgrätenartigen Kristallstruktur des DIP. Diese Diskrepanzen in der Dynamik erklären die unterschiedlichen LD-Werte für DIP und α-PTCDA. In einem weiteren Fall wurde eine Exzitonen-Immobilisierung in helikalen π Aggregaten von Perylen-tetracarboxyl-bisimid (PBI) Molekülen festgestellt. Hier wird Self-Trapping durch einen Relaxationsmechanismus verursacht, in dem das Exziton durch geringe asymmetrische Schwingungen des Aggregats innerhalb von 200 fs von dem hellen Frenkel- in den dunklen Frenkel-Zustand transferiert wird, wobei dieser Übergang von einem CT-Zustand vermittelt wird. Der gesamte Vorgang ist nur bei helikalen Aggregaten möglich, weil nur hier CT-Zustände sehr dicht bei dem hellen Frenkel-Zustand vorhanden sind. Im finalen Frenkel-Zustand tritt eine Torsionsbewegung um die π-Stapelachse ein, so dass ein Energieverlust und eine lokale Änderung der Aggregatstruktur erfolgt – also ein Self-Trapping des Exzitons. Dieser modellierte Mechanismus steht im Einklang zu allen vorliegenden experimentellen Daten. Diese Erkenntnisse lassen die Schlussfolgerung zu, dass in künftigen Materialen für organische Solarzellen eine irreversible und ultraschnelle Deformation des Aggregats nach der Photoanregung vermieden werden muss - will man lange LD erreichen. Nur so kann Self-Trapping von Exzitonen verhindert werden. N2 - In the context of this dissertation very long ranged exciton diffusion lengths (LD) were simulated for perylene-based materials under ideal conditions. This leads to the conclusion that the short LD values in existing materials result from an extrinsic and intrinsic immobilization. The latter, which is a specific material property, is based on a relaxation of the exciton into self-trapping states. An in-depth understanding of the atomistic processes defining self-trapping is essential to developing materials with long LD in the future, in which intrinsic immobilization is prevented. For the development of such a mechanistic understanding it is crucial that a clear relationship between molecular structure and LD is available. This is given by single crystals of diindeno perylene (DIP) and α-perylene tetracarboxylic anhydride (α-PTCDA). An extraordinary large LD of 90 nm was measured for the first one, while the latter possesses only 22 nm. Part of this thesis was to deliver reasons for this discrepancy. Only self-trapping comes into question to explain the different LD values. One reason for the different self-trapping in DIP and α-PTCDA could lie in the electronic structure. However, it was possible to demonstrate that a wide range of perylene-based materials possess no significant differences in their electronic structures. Consequently, such differences can be neglected for the explanation of immobilization mechanisms for the exciton. A further possible explanation could be polarization effects in the crystal, which influences the electronic structure of perylene based materials differently. Especially their influence on charge transfer (CT) states, which are located above the optically bright Frenkel state, was in question because such states could be stabilized by a polarizable surrounding. A significant influence of polarization effects on all considered states were excluded by using a polarizable continuum model. Hence, the small LD values in α-PTCDA are an evidence for self-trapping, which produces a crystal structure built up by π-stacks, while the one of DIP is of herringbone type. Since polarization effects can be neglected, is the dimer only via steric restrictions influenced by the crystal. Hence, a method describing self-trapping has to consider such effects, so that a mechanical embedding QM/MM approach is sufficient. Now, potential energy surfaces were calculated, on which wave packet dynamics were subsequently performed. In this way, atomistic mechanisms for the immobilization of excitons were described for the first time in organic materials. Self-trapping was studied in crystals of α-PTCDA by potential energy surfaces, which map an intermolecular shift motion of the dimer in the crystal. An immobilization of excitons occurs within 500 fs, which results from an irreversible energy loss together with a local deformation of the crystal lattice. This prevents a further transport of the exciton. In the case of DIP, this immobilization does not proceed due to high barriers. These barriers result from the herringbone type packing motif in the DIP crystal. This discrepancy in the dynamics explains the different LD values in DIP and α-PTCDA. In a further example, an exciton immobilization was found in helical π-aggregates of perylene tetracarboxylic bisimide (PBI) molecules. Self-trapping is caused by a relaxation mechanism, in which the exciton is transferred by asymmetric vibrations of the aggregate from the bright to a dark Frenkel state within 200 fs, whereby the transition is mediated by a CT state. However, the CT state is almost non-populated during the whole mechanism so that its participation could not yet be proven experimentally. This entire procedure is solely possible in helical aggregates, because only for such structures is there a CT state located next to the bright Frenkel state. At the final Frenkel state a torsional motion around the π-stacking axis is possible so that the loss in energy and the local rearrangement of the aggregate structure occurs, which means a self-trapping of the exciton. This mechanism is in perfect agreement with all available experimental data. These insights allow the conclusion that in future materials for organic solar cells an irreversible and ultrafast deformation of aggregates after photo-absorption must be avoided. Only in this way long LD values can be achieved and exciton self-trapping can be prevented. However, small LD values are always predicted in helical aggregates of perylene-based materials, because exciton immobilization occurs already due to small molecular motions. For this reason such aggregates are inappropriate for the use in organic solar cells. Long LD values are expected for aggregate structures with long intermolecular shifts or molecules with bulky substituents. KW - Exziton KW - Quantenchemie KW - Angeregter Zustand KW - Self-Trapping KW - CC2 KW - exciton KW - self-trapping KW - quantum chemistry KW - excited state KW - Perylenderivate Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-69861 ER - TY - THES A1 - Liu, Wenlan T1 - Exciton Coupling in Valence and Core Excited Aggregates of pi-Conjugated Molecules T1 - Exzitonen-Kopplung in valenz- und rumpfangeregten Aggregaten pi-konjugierter Moleküle N2 - Im Rahmen dieser Arbeit werden theoretische Modelle zur Beschreibung von Valenz- und Rumpf-angeregten elektronischen Zuständen diskutiert. Im Fall der Valenz-Anregungen wurden time-dependend Hartree-Fock (TD-HF) und timedependent Dichtefunktionaltheorie (TD-DFT)Methoden mit verschiedenen Funktionalen für ein Perylenbisimid (PBI) System validiert. Eine einfache Analyse der Charaktäre der angeregten Zustände wurde vorgeschlagen, die auf den berechneten Übergangsdipolmomenten basiert. Dieser Ansatz ist allerdings auf Zustände beschränkt, die ein signifikantes Übergangsdipolmoment aufweisen. Deshalb wurde eine allgemeinere und fundiertere Methode entwickelt, die auf einer Analyse der berechneten CISWellenfunktion basiert. Darüberhinaus wurde ein literaturbekannter Model-Hamiltonoperator Ansatz von einem lokalisierten Molekülorbitalbild (MO) abgeleitet, das aus der generelleren Analyse-Methode resultiert. Auf diesem Weg ist ein Zugang zu diabatischen angeregten Zuständen und korrespondierenden Kopplungsparametern auf der Basis von ab initio Rechnungen gegeben. Für rumpfangeregte elektronische Zustände wurden drei Methoden für C 1s-angeregte und ionisierte Zustände verschiedener kleiner Moleküle validiert. Darüberhinaus wurde die Basissatzabhängigkeit dieser Zustände untersucht. Anhand der Resultate wurde die frozen core Näherung ausgewählt um rumpfangeregte Zustände von Naphthalintetracarbonsäuredianhydrid (NTCDA) zu berechnen. Um experimentelle Ergebnisse zu erklären, wurde ein Algorithmus entwicklet, der die Exzitonenkopplungsparameter im Fall von nicht-orthogonalen MOs berechnet. N2 - This work focuses on theoretical approaches for predicting the valence and core excited states of aggregate systems. For the valence excitations, TD-HF and TD-DFT with different functionals have been tested at the Perylene bisimide (PBI) system. A simple character analysis method based on the calculated transition dipole moments is proposed. However, this method does not work for excited states without any transition dipole moment. Thus, we proposed a more general and more valid method based on a calculated CIS type wavefunction for the character analysis. Furthermore, a model Hamiltonian method is derived from a localized picture. The energies of the diabatic states and the corresponding coupling parameters were also determined on the basis of ab initio calculations. For the core excitation, three different methods were validated for C 1s-excited and ionized states if several small molecules. Also we tested the basis sets dependence of these core excited states. Based on those results, we chose the frozen core approximation method to evaluate the core excited states of NTCDA molecules. In order to explain the findings in the experiments, we developed an algorithm to evaluate the exciton coupling parameter where non-orthogonal MOs are used. KW - Exziton KW - Dichtefunktionalformalismus KW - Hartree-Fock-Methode KW - Aggregat KW - Angeregter Zustand KW - Quantenchemie KW - Förster-Kopplung KW - zeitabhängige Dichtefunktionaltheorie KW - TD-DFT KW - angeregte Zustände in Aggregaten KW - Quamtum chemistry KW - Förster coupling KW - Exciton KW - time-dependent density functional theory KW - TD-DFT KW - excited states in aggregates Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-56169 ER - TY - THES A1 - Maksimenka, Katsiaryna T1 - Absolute Configuration by Circular Dichroism: Quantum Chemical CD Calculations T1 - Absolute Konfiguration durch Circular-Dichroismus: Quantenchemische Rechnungen des Circular-Dichroismus (CD) N2 - Quantum chemical calculations of circular dichroism (CD) spectra in combination with experimental CD studies are one of the most efficient analytical tools for the elucidation of the three-dimensional structure of a chiral molecule. In the present work 18 chiral compounds of most different molecular structures and origins were investigated using various theoretical methods (the semiempirical CIS methods, the time-dependent DFT and DFT/MRCI approaches). The advantages and limitations of the applied methods were discussed in the context of the studied compounds. Furthermore, the last part of this work deals with the CD investigations of a chiral compound in the crystalline state. A well-known natural product with a specific conformation/CD spectrum behavior was used as a model compound to examine a novel solid-state CD method and to investigate the possibility of its improvement to provide a higher reliability for the assignment of the absolute configuration. N2 - Quantenchemische Rechnungen des Circular-Dichroismus (CD) in Kombination mit experimentellen CD-Studien sind eines der besten analytischen Werkzeuge zur Aufklärung der dreidimensionalen Struktur eines chiralen Moleküls. In der vorliegenden Arbeit wurden 18 chirale Verbindungen mit unterschiedlichsten Strukturen und von unterschiedlichster Herkunft unter Verwendung von verschiedenen theoretischen Methoden (semiempirische CIS-Methoden, zeitabhängige DFT- und DFT/MRCI-Ansätze) untersucht. Die Vorteile und Grenzen der angewandten Methoden wurden im Rahmen der untersuchten Verbindungen diskutiert. Desweiteren befasst sich der letzte Teil der vorliegenden Arbeit mit den CD-Untersuchungen einer chiralen Verbindung im kristallinen Zustand. Ein bekannter Naturstoff mit einem spezifischen Konformation/CD-Spektrum-Verhalten wurde als Modellverbindung verwendet, um einen neuartige Festkörper-CD-Ansatz zu testen, und um die Möglichkeit seiner Verbesserung zu untersuchen, um eine höhere Zuverlässlichkeit für die Zuordnung der absoluten Konfiguration zu bieten. KW - Circular-Dichroismus KW - Quantenchemie KW - Konfiguration KW - Quantenchemische Rechnungen KW - Absolute Konfiguration KW - Absolute Configuration KW - Circular Dichroism KW - Quantum Chemical CD Calculations Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-56552 ER - TY - THES A1 - Bahmann, Hilke T1 - Implementation, Development and Assessment of Local Hybrid Density Functionals T1 - Implementierung, Entwicklung und Validierung lokaler Hybriddichtefunktionale N2 - In order to describe complex molecular systems theoretically, an efficient and reliable solution to the underlying quantum mechanical equations of motion is required. Density functional theory (DFT) represents in most cases the best compromise between accuracy and efficiency for the treatment of electronic interactions. In Kohn-Sham DFT, the non-classical contribution to electron-electron interactions is gathered in the exchange-correlation functional, which has to be approximated in practice. While a large number of exchange-correlation functionals are of semi-empirical nature, some have been derived from physical considerations exclusively. In so-called global hybrid functionals a constant amount of the integrated DFT exchange-energy density is replaced by the exact-exchange energy from Hartree-Fock theory. The most popular functional, B3LYP, contains 20 % exact exchange and several empirical parameters. It has been discovered that the optimal amount of exact exchange depends to a large extent on the molecular property to be computed. A possible solution to this problem is to use local hybrid functionals. Therein, the admixture of exact exchange is controlled by a position-dependent local mixing function (LMF), leading to molecule-specific amounts of exact exchange. In this work a semi-empirical approach is pursued for the development of new local hybrid functionals. Parameterized LMFs are introduced in the exchange-energy density integrals, for which the DFT contributions are taken from established approximations to the exchange-correlation functional. The LMFs developed here contain at least one empirical parameter and a variable that depends on the ratio of the von-Weizsäcker single-particle kinetic energy density to the correlated kinetic-energy density (the so-called t-LMFs), or on the reduced density gradient (referred to as s-LMF). Additional LMFs are obtained by inclusion of the spin polarization. All parameters are fitted to atomization energies and reaction barriers of well-established test sets. Visualization of the LMFs provides an additional tool for analyzing their physical and chemical behavior, potentially leading to further developments. As a general trend, an increasing exact-exchange admixture is observed upon bond stretching for all LMFs, with a more pronounced effect for t-LMFs. This observation correlates with a better performance for reaction barriers of t -LMF-based local hybrid functionals. Most of the local hybrid functionals discussed in this work are based on the exchange and correlation functional from the local spin density approximation (LSDA) and contain therefore no gradient correction such as in the generalized gradient approximation (GGA). The new functionals were initially implemented non-self-consistently into a development version of the quantum chemical Turbomole program package. That is, only the total energy is calculated for a given set of molecular orbitals or electron density, respectively. This is a reliable approximation that allows for significant time savings especially during parameter optimizations. In order to calculate orbital-dependent molecular properties, the local hybrid potential corresponding to the local hybrid energy is required as well. It is obtained as a functional derivative of the exchange-correlation energy with respect to the orbitals. Some of the resulting integrals contain the LMF-weighted non-local exact-exchange potential. These terms as well as the exact-exchange energy density itself cannot be calculated analytically. Following a well-established approach, they have been approximated using a basis set expansion of the exact-exchange potential. For simplicity, the underlying atomic basis set is employed in this resolution of the identity (RI) approximation. For comparison and in view of the optimization of auxiliary basis sets, the optional calculation of the potential by numerical integration has also been implemented in this work. The computational cost of local hybrid calculations for a given basis set, using the RI approximation is comparable to the one of gobal hybrid functionals: a slightly larger prefactor applies to a calculation with a local hybrid functional as compared to a meta-GGA global hybrid, while the scaling of computational effort as a function of system size is the same. Several molecular test sets including atomization energies, barrier heights, dissociation energies and equilibrium distances have been considered for the assessment. Some of them represent particular challenges for current density functional approximations. All of the discussed local hybrid functionals yield significantly better results for the 223 atomization energies of the G3 test set than the B3LYP functional. Especially local hybrid functionals with spin-polarized t -LMFs gives impressively small mean absolute errors for the G3 set. Most of our functionals are in addition significantly superior to B3LYP for the calculation of barrier heights. Some other global hybrid functionals perform even better than our functionals for barriers, but their intrinsic amount of exact exchange is inappropriately high for thermochemical property calculations. For the first time, LSDA-based local hybrid functional have thus been presented that gives accurate results for thermochemistry and reaction barriers simultaneously. The dissociation behavior of symmetric radical cations remains a challenge for the local hybrid functionals presented here. Dissociation energies are significantly overestimated, and the equilibrium distances are too short. The results are overall only slightly better than those obtained using the B3LYP functional. A larger amount of exact exchange is most likely needed for these systems to reduce self-interaction errors. Additionally, the performance of local hybrid functionals for 3d transition metal dimers and monohydrides has been studied. An accurate description of dynamical and nondynamical correlation is essential for the former. The poor performance of most exchange-correlation functionals for transition metal monohydrides can be attributed to self-interaction errors. Our local hybrid functionals perform similarly to B3LYP for the dimers and marginally better for the monohydrides. They do not provide any improvement for the atomic s-d transfer energies of 3d metals. The most suitable local hybrid functional for this particular property uses a s-LMF in the exchange functional and the LYP correlation functional. It yields, however, only average-quality results for thermochemistry and kinetics. Satisfactory results similar to B3LYP are obtained for the isotropic hyperfine coupling constants (HFCCs) of small main group compounds with a t-LMF-based local hybrid functional. The RI approximation to the local hybrid potential has been validated by comparing it to the numerically exact potential for the calculation of total energies, isotropic HFCCs and orbital energies. The error in total energies due to the RI approximation is comparatively small considering the rather large deviations from experimental values. Comparison of mean absolute errors from experimental values of the 26 isotropic HFCCs reveals only small differences between the RI and the numerically exact local hybrid potential. Further analysis shows that inaccuracies in the RI potential may have a larger impact on the isotropic HFCCs or the orbital energies of a particular molecule, especially if only small or medium-sized basis sets are employed. Several of the local hybrid functionals are suitable for the calculation of thermochemical and kinetic properties. Different functionals yield also results similar to other commonly used functionals for isotropic HFCCs of small main group compounds, as well as for the dissociation energies and equilibrium distances of 3d transition metal dimers and monohydrides. The local hybrid functionals studied in this work represent therefore an important step towards the development of universal approximations to the exchange-correlation functional. For a more accurate description of certain transition metal properties and the dissociation behavior of symmetric radical cations while maintaining a good performance for thermochemistry and kinetics, more complex LMFs will have to be considered. Ultimately a local hybrid functional with meta-GGA exchange and correlation energy densities that fulfills more exact constraints is desirable. Therefore further studies on the different gauges of the exchange energy densities are necessary. Another possibility would be the development of a specifically designed correlation functional to be combined with a local hybrid exchange functional based on the LSDA. More detailed studies on the quality of the RI approximation are recommended. Possible properties for this purpose include, e.g., ionization energies and electron affinities. Auxiliary basis sets should be implemented and optimized for the expansion of the exact-exchange potential in order to avoid additional deviations due to the RI-approximation or even fortuitously good results in the assessment of local hybrid functionals with normally contracted basis sets. Since density functional methods are applied extensively for structure optimizations, the gradient of the local hybrid energy with respect to the nuclear coordinates should be implemented to enable this feature in future versions of the code. N2 - Für die Lösung der quantenmechanischen Bewegungsgleichungen, die komplexe, molekulare Systeme beschreiben, sind effiziente und verlässliche Näherungsverfahren erforderlich. Die Dichtefunktionaltheorie (DFT) stellt für die Behandlung der Elektronenwechselwirkung in vielen Fällen den besten Kompromiss zwischen Effizienz und Genauigkeit dar. Im Rahmen der DFT wird die gesamte nicht-klassische Elektron-Elektron-Wechselwirkung im so genannten Austausch-Korrelationsfunktional angenähert. Viele solcher Näherungen sind semi-empirischer Natur, andere wurden ausschließlich von physikalischen Überlegungen abgeleitet. In globalen Hybridfunktionale wird ein konstanter Anteil der integrierten DFT-Austauschenergiedichte durch exakten Austausch aus der Hartree-Fock Näherung ersetzt. Das populärste Funktional B3LYP enthält 20 % exakten Austausch und mehrere empirische Parameter. Der optimale Prozentsatz hängt allerdings sehr stark von den zu berechnenden Systemen und molekularen Eigenschaften ab. Eine Lösung dieses Problems sollten lokale Hybridfunktionale liefern, in denen die Beimischung der exakten Austauschenergiedichte über eine lokale Mischfunktion (LMF) gesteuert wird und daher positions- und molekülabhängig ist. In dieser Arbeit wird ein semi-empirischer Ansatz für die Entwicklung neuer lokaler Hybridfunktionale verfolgt: während die Energiedichten unverändert aus etablierten Näherungen zum Austauschkorrelationsfunktional übernommen werden, stehen parametrisierte LMFs im Zentrum der Untersuchungen. Die verschiedenen LMFs beinhalten neben mindestens einem empirischen Parameter eine Variable die vom Quotienten der von-Weizsäcker kinetischen Energiedichte und der korrelierten kinetischen Energiedichte (sogenannte t-LMFs) bzw. dem reduzierten Dichtegradienten (bezeichnet als t-LMFs) abhängt. Weitere LMFs werden durch zusätzliche Berücksichtigung der Spinpolarisation erhalten. Alle Parameter werden an Atomisierungsenergien bzw. Reaktionsbarrieren bekannter molekularer Testsätze gefittet. Durch Visualisierung der LMFs können zusätzlich Einblicke in den physikalischen Hintergrund und in Möglichkeiten der Weiterentwicklung gewonnen werden. Es wurde beispielsweise beobachtet, dass entlang einer gedehnten Bindung höhere Werte der LMF und damit größere Beimischungen exakter Austauschenergie in Übergangszuständen einhergehen. Dieser Effekt ist für t-LMFs am ausgeprägtesten und korreliert mit besseren Ergebnissen für Reaktionsbarrieren mit lokalen Hybridfunktionalen, die auf einer t-LMF basieren. Bis auf wenige Ausnahmen leiten sich die lokalen Hybridfunktionale in dieser Arbeit aus dem Austausch- und Korrelationsfunktional der lokalen Dichtenäherung (LSDA) ab und enthalten keine Gradientenkorrektur im Sinne der GGA (generalized gradient approximation). Die neuen Funktionale wurden zunächst nicht-selbstkonsistent in eine Entwicklerversion des quantenchemischen Programmpaketes Turbomole implementiert. Das bedeutet, für gegebene Molekülorbitale bzw. eine gegeben Elektronendichte kann lediglich die Gesamtenergie berechnet werden. Dies ist eine anerkannte Näherung, die vor allem für die Optimierung der Parameter eine große Zeitersparnis darstellt. Um letztlich orbitalabhängige, molekulare Eigenschaften berechnen zu können wird neben der Gesamtenergie auch noch das zugehörige lokale Hybridpotential benötigt. Für die selbstkonsistente Implementierung wird die funktionale Ableitung der Austauschkorrelationsenergie nach den Orbitalen bestimmt. Daraus resultierend müssen neben den üblichen lokalen Austauschkorrelationspotentialtermen auch Integrale berechnet werden, die das mit der LMF gewichtete nicht-lokale exakte Austauschpotential enthalten. Die entsprechenden Terme kann man, genauso wie die exakte Austauschenergiedichte an sich, nicht analytisch berechnen. Früheren Ansätzen folgend wurden sie in der vorliegenden Arbeit in einer Basissatzentwicklung angenähert, wobei der Einfachheit halber die atomaren Basisfunktionen verwendet wurden. Um die Genauigkeit dieser sogenannten RI (resolution of the identity)-Näherung validieren zu können und auch schon im Hinblick auf die Anpassung einer Hilfsbasis, wurde darüber hinaus die numerische Berechnung aller Integrale, die das exakte Austauschpotential und die entsprechende Energiedichte enthalten, implementiert. Unter Verwendung der RI-Näherung ist der Rechenaufwand lokaler Hybride vergleichbar mit dem globaler Hybridfunktionale: Während die formale Skalierung in Abhängigkeit der Systemgröße gleich ist, ergab sich ein etwas höherer Vorfaktor für die lokalen Hybride. Verschiedene Literaturbekannte Testsätze mit Atomisierungsenergien, Reaktionsbarrieren, Dissoziationsenergien oder Gleichgewichtsabständen, die teilweise einige Schwächen bisheriger Dichtefunktionalnäherungen aufdecken, wurden berücksichtigt. Für die 223 Atomisierungsenergien des G3 Testsatzes stellen alle unsere Funktionale eine signifikante Verbesserung gegenüber B3LYP dar. Atomisierungsenergien sind insofern ein sensibler Test, da alle Bindungen gebrochen werden und Fehlerkompensation eine untergeordnete Rolle spielt. Vor allem lokale Hybridfunktionale, deren LMFs neben der kinetischen Energiedichte explizit von der Spinpolarisation abhängen, lieferten hervorragende Resultate. Obwohl im Vergleich zu Atomisierungsenergien für die korrekte Berechnung von Reaktionsbarrieren im Allgemeinen mehr exakter Austausch benötigt wird, sind unsere Funktionale auch für zwei Testsätze mit jeweils 38 Reaktionsbarrieren besser als B3LYP. Zwar kann mit einem globalen Hybrid mit 50 % exaktem Austausch eine geringere Abweichung von den Richtwerten erzielt werden, aber ein solches Funktional ist für thermochemische Daten unzureichend. Hier wurde erstmals gezeigt, dass lokale Hybridfunktionale ohne Gradientenkorrektur sowohl für Thermochemie als auch für Kinetik zufrieden stellende Ergebnisse liefern können. Das Dissoziationsverhalten symmetrischer Radikalkationen stellt für die hier diskutierten Dichtefunktionale nach wie vor eine Herausforderung dar: Die Dissoziationsenergien von sieben Modellsystemen werden mit unseren Funktionalen stark überschätzt und Gleichgewichtsabstände unterschätzt. Insgesamt sind die Werte nur marginal besser als mit B3LYP. Neben Eigenschaften von Hauptgruppenverbindungen wurden zudem Übergangsmetalldimere und -monohydride untersucht. Für erstere ist eine gute Beschreibung dynamischer sowie statischer Elektronenkorrelation ausschlaggebend. In den Hydriden andererseits dominiert mit gängigen Dichtefunktionalen die unphysikalische Selbstwechselwirkung eines Elektrons mit sich selbst. Für die 3d-Übergangsmetalldimere sind die getesteten Funktionale genauso gut wie B3LYP und für die Hydride etwas besser. Atomare s-d Transferenergien von 3d Übergangsmetallen verbleiben auch für unsere lokalen Hybridfunktionale, die insgesamt schlechtere Ergebnisse erzielen als B3LYP, noch problematisch. Das hierfür geeignetste lokale Hybridfunktional basiert auf einer s-LMF und beinhaltet LYP Korrelation. Für die isotropen Hyperfeinkopplungskonstanten (HFCCs) kleiner Hauptgruppenverbindungen wurden zufriedenstellende Ergebnisse (ähnlich wie B3LYP) mit einem t-LMF basierten lokalen Hybrid erzielt. Die RI Näherung zum lokalen Hybridpotential wurde dem numerisch exakten Potential für die Berechnung von Gesamtenergien, isotrope HFCCs und Orbitalenergien für verschiedene Basissätze gegenübergestellt. Wie erwartet ist der Fehler für Gesamtenergien mit der RI-Näherungen vergleichsweise gering, vor allem relativ zu den verbleibenden Abweichungen von experimentellen Energien. Der Vergleich der mittleren absoluten Abweichung von experimentellen Werten für 26 isotrope HFCCs zeigt sogar für mittelgroße und kontrahierte IGLO Basissätze nur geringe Unterschiede zwischen dem RI-Potential und dem numerisch exakten lokalen Hybridpotential. Die Analyse der HFCCs einzelner Moleküle und der Orbitalenergien des CN Moleküls offenbart allerdings, dass Ungenauigkeiten aufgrund der RI-Näherung hier eine größere Rolle spielen, vor allem wenn zu kleine atomare Basissätze verwendet werden. Von den untersuchten lokalen Hybriden stellen sich einige als hervorragende Kandidaten für die Berechnung thermochemischer und kinetischer Eigenschaften heraus. Jeweils unterschiedliche Funktionale erzielen darüber hinaus mit den besten bekannten Funktionalen vergleichbare Ergebnisse für isotrope Hyperfeinkopplungskonstanten und ausgewählte Eigenschaften kleiner Übergangsmetallverbindungen. Die in dieser Arbeit präsentierten lokalen Hybridfunktionale stellen daher einen wichtigen Schritt in der Entwicklung universeller Näherungen zum Austauschkorrelationsfunktional dar. Zur akkuraten Beschreibung molekularer Eigenschaften von Übergangsmetallkomplexen und dem Dissoziationsverhalten von Radikal-Kation-Dimeren neben Thermochemie und Kinetik, werden in Zukunft wohl komplexere LMFs benötigt. Um konkurrenzfähige lokale Hybride mit gradientenkorrigierter Austausch- und Korrelationsenergiedichte zu entwickeln, müssen darüber hinaus weitere Studien zum Einfluss des abweichenden Eichursprungs der miteinander kombinierten Austauschenergiedichten durchgeführt werden. Eine andere Möglichkeit ist die Entwicklung speziell abgestimmter Korrelationsfunktionale für lokale Hybride. Außerdem sollte die Qualität der RI-Näherung zum lokalen Hybridpotential detaillierter untersucht werden. Hierfür könnten zum Beispiel Ionisierungsenergien und Elektronenaffinitäten herangezogen werden. Um zusätzliche Abweichungen oder sogar fälschlicherweise "zu gute" Ergebnisse bei Validierungsrechnungen zu vermeiden, sollten Hilfsbasen für die Entwicklung des nicht-lokalen exakten Austauschpotentials implementiert und optimiert werden. Einer der nächsten Implementierungsschritte sollte auch Gradienten bezüglich der Kernkoordinaten beinhalten, um die Validierung der neuen lokalen Hybridfunktionale auf Strukturoptimierungen auszuweiten. KW - Validierung KW - Quantenchemie KW - Dichtefunktional KW - Atomisierungsenergien KW - Reaktionsbarrieren Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-55641 ER - TY - THES A1 - Seibt, Joachim T1 - Theoretical investigations on the spectroscopy of molecular aggregates T1 - Theoretische Untersuchungen zur Spektroskopie von Molekülaggregaten N2 - Die spektroskopischen Eigenschaften von Molekülaggregaten wurden mittels quantendynamischer Berechnungen untersucht. Hierbei wurden sowohl lineare als auch nichtlineare Spektroskopietechniken einbezogen. Zur Simulation von Absorptions- und CD-Spektroskopie wurden Kopplungseffekte sowie die relative Orientierung der Monomer-Einheiten in den Modellen berücksichtigt, um gemessene Spektren reproduzieren und so die entsprechenden Parameter zu bestimmen. Zur genaueren Beschreibung wurden auch Ergebnisse quantenchemischer Rechnungen verwendet. Darüber hinaus wurden Untersuchungen zur nichtlinearen optischen Spektroskopie an Dimeren durchgeführt. N2 - The spectroscopic properties of molecular aggregates have been investigated by means of quantum dynamical calculations. Thereby both linear and nonlinear spectroscopic techniques have been taken into account. For the simulation of absorption and CD-spectra, coupling effects were regarded as well as the relative orientation of the monomer units in order to determine the parameters by reproducing measured spectra. For a more detailled description, results from quantum chemical calculations have also been included. Furthermore, investigations on nonlinear spectroscopy of molecular dimers have been performed. KW - Theoretische Chemie KW - Aggregat KW - Nichtlineare Spektroskopie KW - CD-Spektroskopie KW - Quantenchemie KW - Quantendynamik KW - quantum dynamics Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-37218 ER - TY - THES A1 - Mladenovic, Milena T1 - Theoretical Investigation into the Inhibition of Cystein Proteases T1 - Theoretische Untersuchung der Inhibierung von Cysteinproteasen N2 - Although known about and investigated since the late 1970’s, the picture of the basic principles governing inhibitor strengths and the structure-activity relationships of the cysteine protease inhibition mechanism is still very incomplete. Computational approaches can be a very useful tool for investigating such questions, as they allow the inspection of single, specific effects in isolation from all others, in a manner very difficult to achieve experimentally. The ab initio treatments of such large systems like proteins are still not feasible. However, there is a vast number of computational approaches capable of dealing with protein structures with reasonable accuracy. This work presents a summary of theoretical investigations into cysteine protease cathepsin B using a range of methods. We have concentrated on the investigation of cysteine protease inhibition by epoxide- and aziridine-based inhibitors in order to obtain better insight into these important topics. Various model systems are simulated by means of pure quantum mechanical methods and by hybrid (QM/MM) methods. Both approaches provide a static picture. Dynamical effects are then accounted for by additional molecular dynamics (MD) simulations, using both classical and QM/MM MD approaches. The quantum mechanical approach was used to study very small model systems consisting only of the electrophilic warhead of the inhibitor (both substitituted and not) and molecular moieties simulating a very simplified protein active site (methylthiolate instead of Cys29 and methylimidazolium instead of His199 residue) and solvent surroundings (two waters or two ammonium ions, in combination with a continuum solvent model). Although simple, such a system provides a good description of the most important interactions involved in the inhibition reaction. It also allows investigation of the influence of the properties of the electrophilic warhead on the reaction rate. Beside the properties of the electrophilic warhead, the protein and solvent environment is also an important factor in the irreversible deactivation of the enzyme active site by the inhibitor. The non-covalent interactions of the inhibitor with the oxyanion hole and other subsites of the enzyme, as well as its interaction with the solvent molecules, need to be explicitly taken into account in the calculations, because of their possible impact on the reaction profile. As molecular modeling methods allow the treatment of such large systems, but lack the possibility of describing covalent interactions, our method of choice was the combined quantum mechanics/molecular modeling approach. By splitting the system into a smaller part that undergoes the bond cleavage/formation process and must be treated quantum mechanically, and a larger part, comprised of the rest of the protein, which could be treated using force fields, we managed to simulate the system at the desired precision. Our investigations concentrated on the role of His199 in the inhibition mechanism as well as on the structure-reactivity relationships between cysteine protease and various inhibitors, yielding new insight into the kinetics, regio- and stereospecificity of the inhibition. In particular, our calculations provide the following insights: i.) an explanation for the regioselectivity of the reaction, and original insight into which interactions affect the stereoselectivity; ii.) a clear model which explains the known structure-activity relationships and connects these effects with the pH-dependency of the inhibition; iii.) our computations question the generally accepted two-step model by showing that substituent effects accelerate the irreversible step to such an extent that the achievement of an equilibrium in the first step is doubtful; iv.) by way of theoretical characterizations of aziridine models, the reasons for similarities and differences in the mode of action of epoxide- and aziridine-based inhibitors are elucidated; and finally, v.) combining our results with experimental knowledge will allow rational design of new inhibitors. To account for dynamical effects as well, molecular dynamics (MD) computations were also performed. In these calculations the potential energy was computed at the force field level. The results not only supported and clarified the QM/MM results, but comparison with previous X-ray structures helped correct existing errors in the available geometrical models and resolved inconsistencies in the weighting of various factors governing the inhibition. In the work the first QM/MM MD calculations on the active site of the cysteine proteases are presented. In contrast to the MD simulations, these calculations used potential energies computed at the QM/MM-level. With the help of these computations we sought to address strongly disputed questions about the reasons for the existence of the active site ion pair and its role in the high activity of the enzyme. N2 - Obwohl bereits seit den späten 70ern bekannt und untersucht, ist das Bild über die grundsätzlichen Prinzipien, welche die Wirksamkeit der Inhibition und die Struktur-Aktivitäts-Beziehungen (SAB) der Cysteinprotease-Hemmmechanismen beeinflussen, immer noch sehr unvollständig. In dieser Arbeit wurden quantenmechanische (QM), molekulardynamische (MD) und gemischte quantenmechanische/molekularmechanische (QM/MM) Berechnungen durchgeführt, um die irreversible Inaktivierung des Enzyms durch den Inhibitor, zu untersuchen. Die Stärke von computergestützten Verfahren liegt in der Möglichkeit, den Einfluss einzelner spezifischer Effekte durch das Ausblenden von Umgebungseffekten zu untersuchen. Diese Herangehensweise ist nur sehr schwer im Experiment zu erreichen. Die quantenmechanische Methode wurde benutzt, um sehr kleine Modellsysteme zu untersuchen, welche lediglich aus dem elektrophilen Kerngerüst des Inhibitors und einigen, das aktive Zentrum im Protein vereinfachend darstellenden, Molekülen bestanden. Die Solvensumgebung wurde durch zwei explizite Wassermoleküle und zwei Ammoniumionen in Kombination mit einem Kontinuum-Solvensmodell berücksichtigt. Ein solches vereinfachtes System liefert bereits eine gute Beschreibung der meisten wichtigen Wechselwirkungen während der Inhibierungsreaktion. Es erlaubt die Untersuchung des Einflusses der Eigenschaften des elektrophilen „Warheads“ auf den Reaktionsverlauf. Neben den Eigenschaften des elektrophilen „Warheads“ ist sowohl die Protein- als auch die Solvensumgebung von großer Bedeutung für die irreversible Deaktivierung des aktiven Zentrums des Enzyms durch den Inhibitor. Aufgrund eines möglichen Einflusses auf das Reaktionsprofil müssen Solvensmoleküle sowie die nicht-kovalenten Wechselwirkungen des Inhibitors mit dem Oxyanionloch und anderen Nebenbindungstellen des Enzyms explizit behandelt werden. Berücksichtigt man, dass nur molekularmechanische Methoden eine Behandlung von solch großen Systeme erlauben, im Gegenzug aber nicht in der Lage sind, kovalente Wechselwirkungen zu beschreiben, so wurde als Methode der Wahl ein kombinierter QM/MM Ansatz gewählt. Durch das Aufteilen des Gesamtsystems in einen kleinen Bereich, der die Bindungsspaltung- und Bindungsbildungsreaktionen beinhaltet (mit QM behandelt), und in einen großen Teil, welcher den Rest des Proteins umfasst (mit MM behandelt), waren wir in der Lage, das System mit gewünschter Genauigkeit zu simulieren. In den Kapiteln 3.2.2 bis 3.2.4 werden die Struktur-Reaktivitäts-Beziehungen (SRB) zwischen der Cysteinprotease und verschiedenen Inhibitoren vorgestellt. Es wird gezeigt, dass die SRB eine entscheidende Rolle in Kinetik, Regio- und Stereoselektivität der Inhibierung spielen. Um auch die dynamischen Effekte im Anspruch zu nehmen, wurden klassische molekulardynamische (MD) Simulationen durchgeführt. Die Ergebnisse dieser Berechnungen haben nicht nur unterstützt und die QM/MM Ergebnisse näher erklärt, sondern auch durch der Vergleich mit bereits bestehenden Röntgenstrukturen verschiedener Cysteineprotease-Inhibitor-Komplexen geholfen, Fehler in vorhandenen Proteinkristallstrukturen zu korrigieren. Somit wurden Widersprüche in der Bewertung verschiedener Funktionen, die die Inhibierungsreaktion bestimmen, aufgelöst (Kapitel 3.2.5). In der Arbeit werden auch die ersten QM/MM MD Berechnungen vom aktiven Zentrum der Cysteineprotease vorgestellt. Im Kontrast zu den klassischen MD Simulationen, wird in diesen Untersuchungen die potentielle Energie auf dem QM/MM Theorieniveau berechnet. Hier haben wir versucht die stark umstrittene Frage über die Gründe für das Vorhandensein des aktiven Zentrums als His199+/Cys29- Ion-paar und seine Rolle für die hohe Aktivität des Enzyms zu beantworten. KW - Quantenchemie KW - Quantum mechanics / molecular modeling KW - Molecular dynamics Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-25763 ER - TY - THES A1 - Auer, Dominik T1 - Lithiosilane mit stereogenen Silicium-Zentren : Synthese, Struktur und Reaktionsverhalten T1 - Lithiosilanes with stereogenic silicon centres : synthesis, structure and reactivity N2 - Die vorliegende Arbeit liefert einen Beitrag zur Synthese und zum Reaktionsverhalten enantiomerenangereicherter Silyllithium- und Silylmagnesiumverbindungen. Dabei sind mehrere Punkte von Bedeutung: die Aufklärung der Struktur im Festkörper wie in Lösung, die Bestimmung der Selektivitäten und des stereochemischen Verlaufs hoch enantiomerenangereicherter polarer Silylmetallverbindungen in Umsetzungen mit Elektrophilen, die konfigurative Stabilität hoch enantiomerenangereicherter Silyllithium- und Silylmagnesiumverbindungen sowie die Aufklärung ungewöhlicher NMR-Verschiebungen. Im Rahmen dieser Arbeit durchgeführte quantenchemische Berechnungen haben sich vor allem auf die beiden letzten Punkte konzentriert. Insgesamt wurden folgende Themenkomplexe bearbeitet: 1.) Synthese einer hoch enantiomerenangereicherten Silyllithiumverbindung 2.) Reaktionsverhalten der hoch enantiomerenangereicherten Silyllithiumverbindung 2 gegenüber Elektrophilen 3.) Konfigurative Stabilität der hoch enantiomerenangereicherten Silyllithiumverbindung 2 4.) Weitere Untersuchungen zur Darstellung von Silyllithiumverbindungen 5.) Untersuchungen zur Struktur von Silyllithiumverbindungen im Festkörper und in Lösung 6.) Quantenchemische Analysen der Ursachen der NMR-Verschiebungen von Disilenen und heteroatom-substituierten Silanen N2 - This work introduces syntheses and reactivity studies of enantiomerically enriched silyllithium and silylmagnesium compounds. The following main topics are of importance here: structure determination in the solid state and solution, evaluation of selectivities and stereochemical path in reactions of highly enantiomerically enriched polar silylmetallic compounds with electrophiles, the configurational stability of highly enantiomerically enriched silyllithium and silylmagnesium compounds and rationalisation of unexpected NMR-shifts. The quantum chemical calculations presented here focus mainly on the latter two topics. Alltogether this work focusses on the following fields: 1.) Synthesis of a highly enantiomerically enriched silyllithium compound 2.) Reactivity of the highly enantiomerically enriched silyllithium compound 2 with respect to electrophiles 3.) Configurational stability of the highly enantiomerically enriched silyllithium compound 2 4.) Further investigations concerning the preparation of silyllithium compounds 5.) Investigations regarding the structure of silyllithium compounds in the solid state and in solution 6.) Quantum chemical analyses of the cause of NMR-shifts of disilenes and heteroatom-substituted silanes KW - Lithiumorganische Verbindungen KW - Siliciumorganische Verbindungen KW - Stereochemie KW - Lithium KW - Silyllithium KW - stereogen KW - NMR KW - Quantenchemie KW - lithium KW - silyllithium KW - stereogenic KW - NMR KW - quantum chemistry Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-11300 ER -