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. / 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.
| Identifer | oai:union.ndltd.org:uni-wuerzburg.de/oai:opus.bibliothek.uni-wuerzburg.de:12349 |
| Date | January 2015 |
| Creators | Walter, Christof |
| Source Sets | University of Würzburg |
| Language | English |
| Detected Language | English |
| Type | doctoralthesis, doc-type:doctoralThesis |
| Format | application/pdf |
| Rights | https://opus.bibliothek.uni-wuerzburg.de/doku/lic_ohne_pod.php, info:eu-repo/semantics/openAccess |
Page generated in 0.0051 seconds