Institut für Physikalische und Theoretische Chemie
Refine
Has Fulltext
- yes (16)
Is part of the Bibliography
- yes (16)
Year of publication
Document Type
- Doctoral Thesis (14)
- Journal article (1)
- Preprint (1)
Keywords
- Exziton (16) (remove)
Institute
EU-Project number / Contract (GA) number
- 614623 (2)
The goal of this thesis was the development and application of higher-order spectroscopic techniques. In contrast to ordinary pump–probe (PP) and two-dimensional (2D) spectroscopy, higher-order coherently detected spectroscopic methods measure a polarization that has an order of nonlinearity higher than three. The key idea of the techniques in this thesis is to isolate the higher-order signals from the lower-order signals either by their excitation frequency or by their excitation intensity dependence. Due to the increased number of interactions in higher-order spectroscopy, highly excited states can be probed. For excitonic systems such as aggregates and polymers, the fifth-order signal allows one to directly measure exciton–exciton annihilation (EEA). In polymers and aggregates, the exciton transport is not connected to a change of the absorption and can therefore not be investigated with conventional third-order techniques. In contrast, EEA can be used as a probe to study exciton diffusion in these isonergetic systems. As a part of this thesis, anisotropy in fifth-order 2D spectroscopy was investigated and was used to study geometric properties in polymers.
In 2D spectroscopy, the multi-quantum signals are separated from each other by their spectral position along the excitation axis. This concept can be extended systematically to higher signals. Another approach to isolate multi-quantum signals in PP spectroscopy utilizes the excitation intensity. The PP signal is measured at specific excitation intensities and linear combinations of these measurements result in different signal contributions. However, these signals do not correspond to clean nonlinear signals because the higher-order signals contaminate the lower-order multi-quantum signals. In this thesis, a correction protocol was derived that uses the isolated multiquantum signals, both from 2D spectroscopy and from PP spectroscopy, to remove the contamination of higher-order signals resulting in clean nonlinear signals. Using the correction on the third-order signal allows one to obtain annihilation-free signals at high excitation intensities, i.e., with high signal-to-noise ratio. Isolation and correction in PP and 2D spectroscopy were directly compared by measuring the clean third-order signals of squaraine oligomers at high excitation intensities. Furthermore, higher-order PP spectroscopy was used to isolate up to the 13th nonlinear order of squaraine polymers.
The demonstrated spectroscopic techniques represent general procedures to isolate clean signals in terms of perturbation theory. The technique of higher-order PP spectroscopy needs only small modifications of ordinary PP setups which opens the field of higher-order spectroscopy to the broad scientific community. The technique to obtain clean nonlinear signals allows one to systematically increase the number of interacting (quasi)particles in a system and to characterize their interaction energies and dynamics.
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.
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.
Zahlreiche theoretische und experimentelle Untersuchungen haben erwiesen, dass in halbleitenden Kohlenstoffnanoröhren durch Absorption von Licht hauptsächlich Exzitonen erzeugt werden. Die photophysikalischen Eigenschaften und insbesondere die Prozesse nach der optischen Anregung sind aber gegenwärtig noch nicht vollständig verstanden. Zeitaufgelöste Spektroskopie bietet die Möglichkeit, diese Prozesse zu verfolgen und somit detaillierten Einblick in das photophysikalische Verhalten von Kohlenstoffnanoröhren zu nehmen. Hierbei scheinen auch extrinsische Faktoren - zu nennen sind die Herstellungsmethode, die Art der Probenpräparation, der Aggregationsgrad sowie der durch das Lösungs- bzw. Dispersionsmittel bedingte Einfluss - eine entscheidende Rolle zu spielen.
In dieser Dissertation wurden die Exzitonengröße sowie die exzitonische Dynamik in einwandigen Kohlenstoffnanoröhren mittels transienter Absorptionsspektroskopie sowie stationärer und zeitaufgelöster Photolumineszenzmessungen untersucht. Alle Experimente fanden dabei an halbleitenden (6,5)-Kohlenstoffnanoröhren statt, deren chirale Anreicherung durch Dichtegradientenultrazentrifugation gelang. Für die temperaturabhängigen Messungen wurde ein Verfahren zur Herstellung von tensidstabilisierten Gelatinefilmen entwickelt. Diese zeichnen sich durch eine hohe Temperaturstabilität bei gleichzeitiger Minimierung von Streulichteffekten aus.
Die Bestimmung der Exzitonengröße erfolgte mit Hilfe des Phasenraumfüllmodells, das die intensitätsabhängige Änderung der Oszillatorstärke eines Übergangs mit der Exzitonengröße verknüpft. Hierfür wurden leistungsabhängige Messungen der transienten Absorption durchgeführt und die Signalintensität des Photobleichens gegen die absorbierte Photonenflussdichte aufgetragen. Da diese beiden Größen nur bei geringer Exzitonendichte in einer linearen Beziehung stehen, aus der sich die Exzitonengröße berechnen lässt, wurde im Experiment besonderer Wert auf niedrige Anregungsfluenzen und deren exakte Bestimmung gelegt. Um den Einfluss der Aggregation quantifizieren zu können und den Vergleich mit der Literatur zu erleichtern, fanden die Untersuchungen sowohl an individualisierten als auch an aggregierten Röhrenproben statt. Die Datenanalyse, bei der erstmalig die stimulierte Emission sowie der spektrale Überlapp von Photoabsorptions- und Photobleichbande Berücksichtigung fanden, ergab für individualisierte (6,5)-Nanoröhren einen Wert von 12.0 nm für die Größe des S1-Exzitons, während diese bei der aggregierten Röhrenprobe nur 5.6 nm beträgt. Die Probenabhängigkeit der Exzitonengröße macht den Vergleich mit anderen experimentell ermittelten Werten schwierig. Diese liegen fast ausschließlich zwischen 1 nm und 4.5 nm, ihre Bestimmung fand aber teilweise an stark aggregierten bzw. polydispersen Proben statt. Theoretische Berechnungen liefern für die Exzitonengröße Werte zwischen 1 nm und 4 nm. Zwar gelten einige der Berechnungen für Vakuum, was verglichen zu einer experimentell in Lösung bzw. im Film bestimmten Exzitonengröße einen kleineren Wert mit sich bringt, jedoch kann allein hierdurch die Diskrepanz zu der in dieser Arbeit ermittelten Exzitonengröße von 12.0 nm nicht erklärt werden. Setzt man experimentell und theoretisch für Vakuum bestimmte Werte für die Exzitonengröße und die Bindungsenergie in einen einfachen Zusammenhang, entspricht eine Exzitonengröße von 12.0 nm einer Bindungsenergie zwischen 0.21 eV und 0.27 eV. Die mittels Zweiphotonenexperimenten ermittelten Werte für die Bindungsenergie von (6,5)-Kohlenstoffnanoröhren befinden sich zwischen 0.37 eV und 0.42 eV; diese wurden allerdings unter Zuhilfenahme eines vereinfachten zylindrischen Modells abgeschätzt. Weitere experimentelle und theoretische Untersuchungen könnten klären, inwieweit eine exzitonische Bindungsenergie zwischen 0.21 eV und 0.27 eV für (6,5)-SWNTs in Betracht kommt.
Strahlender und nichtstrahlender Zerfall in den Grundzustand scheinen in (6,5)-Kohlenstoffnanoröhren durch eine Dynamik zwischen verschiedenen Zuständen sowie durch die Diffusion der Exzitonen beeinflusst zu werden. Um diese für die Rekombination maßgeblichen Prozesse besser zu verstehen, wurden temperaturabhängige Messungen der stationären und zeitaufgelösten Photolumineszenz sowie der transienten Absorption durchgeführt. Die Ergebnisse der stationären PL-Experimente deuten darauf hin, dass die Exzitonen zwischen dem optisch aktiven Singulettzustand mit A2-Symmetrie - im Folgenden mit [B] bezeichnet - und einem energetisch tiefer liegenden dunklen Zustand [D] gestreut werden. Mit einem Wert von 5 meV für die energetische Aufspaltung zwischen [B] und [D] gelingt eine gute Anpassung an die Daten, was mit Blick auf die Bandstruktur von (6,5)-SWNTs vermuten lässt, dass es sich bei [D] um den A1-Singulettzustand handelt. Außerdem scheint eine nichtthermische Verteilung der Exzitonen auf [B] und [D] vorzuliegen, wobei strahlende Rekombination nur vom Zustand [B] aus möglich ist. Mit diesen Annahmen kann das temperaturabhängige Verhalten der stationären Photolumineszenz modelliert werden, die Ergebnisse der zeitaufgelösten PL-Messungen jedoch nicht. Mit einem rein diffusionsdominierten Modell gelingt dies ebenso wenig, so dass zur Interpretation des PL-Zerfalls vermutlich ein Modell entwickelt werden muss, in dem sowohl die Streuung der Exzitonen zwischen [B] und [D] als auch das durch Diffusion bedingte Löschen an Defektstellen oder Röhrenenden Berücksichtigung findet. Die Bedeutung der Diffusion von Exzitonen zu Defektstellen oder Röhrenenden, an denen bevorzugt nichtstrahlender Zerfall stattfindet, kann durch spektral- und zeitaufgelöste PL-Messungen belegt werden. Abhängig von der zur Verfügung stehenden thermischen Energie und der Höhe der Potenzialbarrieren des untersuchten Systems kann die Diffusion niederenergetischer Exzitonen, die sich in Potenzialminima befinden, soweit eingeschränkt werden, dass diese eine fast bis um den Faktor zwei längere PL-Lebensdauer aufweisen als höherenergetische Exzitonen. Das unterschiedliche Verhalten von transienter Absorption und zeitaufgelöster Photolumineszenz bei Temperaturen zwischen 14 K und 35 K zeigt, dass die Repopulation des Grundzustands hauptsächlich von einem anderen Zustand aus erfolgt als die strahlende Rekombination. Ob es sich hierbei aber um den mit [D] bezeichneten A1-Singulettzustand oder einen anderen dunklen Zustand handelt, kann nicht abschließend geklärt werden.
Aufgrund inhomogener Verbreiterung stellt die Halbwertsbreite der Banden im Absorptionsspektrum ein Maß für die Höhe der Potenzialbarrieren bzw. für die energetische Verteilung der Exzitonen im angeregten Zustand dar. In dieser Arbeit wurde anhand vier verschiedener Nanorohrsuspensionen gezeigt, dass Sättigungsverhalten der transienten Absorption von (6,5)-Kohlenstoffnanoröhren und Bandenbreite im Absorptionsspektrum demselben Trend folgen. Begründen kann man dies damit, dass das Sättigungsverhalten der transienten Absorption durch Exziton-Exziton-Annihilation bestimmt wird. Aufgrund ihrer eindimensionalen Struktur unterliegen Kohlenstoffnanoröhren einer starken Beeinflussung durch die Umgebung. Abhängig vom Lösungs- bzw. Dispersionsmittel resultiert eine unterschiedliche inhomogene Verbreiterung der Absorptionsbanden und damit unterschiedlich hohe Potenzialbarrieren im angeregten Zustand. Niedrige Potenzialbarrieren erlauben eine weitreichende Diffusion der Exzitonen, sodass effiziente Exziton-Exziton-Annihilation schon bei einer vergleichsweise geringen Exzitonendichte stattfindet und das Signal der transienten Absorption bei einer niedrigen Impulsfluenz sättigt.
Two thematic complexes were addressed within this work. One part is related to improvements and new implementations into the CAST program package. Thereby the main focus laid on the delivery of a tool which can be used to characterize complex reactions and their mechanisms. But also within the new force field (FF) method (SAPT-FF) within the CAST program, several improvements were made. The second topic is related to the description of dye molecules and their spectral properties. The main focus within these studies was set on the influence of the environment on these properties. In the first topic improvements of the local acting NEB (nudged elastic band) methods were included and the number of available methods was extended. The initial pathway generation was improved by implementing the IDPP (image dependent pair potential) method and a new method was implemented for describing temperature dependent pathways. Additionally, improvements have been made to the optimization routines (global NEB). As a second part the Pathopt (PO) method was considerably improved. In the beginning of the work the original PO idea was used. In this approach one starts with a global optimization on one n-1 dimensional hyperplane which divides the reaction into two sub-areas for obtaining guesses of TSs (transition states). These found TS guesses were used to optimize to the ”true” TS. Starting from the optimized ones a relaxation to the next connected minima is done. This idea has been automatically implemented and extended to several number of hyperplanes. In this manner a group of pathsegments is obtained which needs to be connected, but within this work it was realized that such a procedure might be not very efficient. Therefore, a new strategy was implemented which is founded on the same constrained global optimization scheme (MCM) for which the user defines the number of hyperplanes generated. The number of such generated hyperplanes should be large enough
134
to describe the space between the concerning reactants in a sufficient way. The found minima are directly used to built up the reaction pathway. For this purpose a RMSD (root mean square deviation) criterion is used to walk along ways of minimal change from one to another hyperplane. To prove the implementations various test calculations were carried out and extensions included to prove the capabilities of the new strategy. Related to these tests a new strategy for applying the move steps in MCM (Monte Carlo with minimization) was realized which is also related to the question of the coordinates representation. We were able to show that the hopping steps in MCM can be improved by applying Cartesian steps in combination of random dihedral moves with respect to the constraint. In this way it was possible to show that a large variety of systems can be treated. An additional chapter shows the improvements of the SAPT-FF implementation and related test cases. It was possible to treat benzene dimer and cluster systems of different sizes consistently also in accordance with high level ab initio based approaches. Furthermore, we showed that the SAPT-FF with the right parameters outperforms the standard AMOEBA implementation which is the basis of the SAPT-FF implementation. In the last three chapters deal with the description of perlyene-based dyes. In the first smaller chapter ground state chemistry description of macro cycles of PBI (perylene bisimide) derivatives were investigated. Therefore, AFM (atomic force microscopy) based pictures were explained within our study. The methods to explain aggregation behavior in dependency of the ring size were MD simulations and configuration studies. The last two chapters deal with opto-electronic or photo-physical properties of PBI and PTCDA (perylene-3,4,9,10-tetracarboxylic dianhydride). In detail, we investigated the role of the environment and the aggregate or crystal surrounding by applying different models. In that way implicit and explicit solvation models, the size of aggregates and vibration motions were used. In the case of PBI the recent work is found on preliminary studies related to my bachelor thesis and extends it. It was shown that the direct influence of a polarizable surrounding, as well as explicit inclusion of solvent molecules on the overall description of the excitations and nature of the excited states is weaker as one might expect. However the inclusion of intra-molecular degrees of freedom showed a stronger influence on the state characteristics and can induce a change of the order of states within the dimer picture. For the PTCDA molecule the main focus was set on the description of the absorption spectrum of crystalline thin films. Related to this older works exist which already gave a description and assignment of the absorption band, but are based on different approaches compared to the one used in this work. We used the supermolecule ansatz, whereas the environment and different aggregate sizes were investigated. Within the dimer based approach we were able to show that using continuum solvation (IEFPCM/COSMO) based description for the environment the relative order of states remains unchanged. Similar to the PBI calculations the influence of the vibrational motions /distortions is larger. The simulation of the crystal environment by using QM/MM (quantum mechanics/molecular mechanics) approaches delivered that an asymmetric charge distribution might induce a localization of the excitation and a stronger mixing of states. For obtaining further insights we go beyond the dimer picture and aggregates of different sizes were used, whereas the simulations up to the octadecamer mono- and even dual-layer stack were carried out. Within these calculations it was shown that the H-coupling is dominating over a weaker J-coupling between different stacks. Additionally the calculations based on DFT (density functional theory) and semi-empirics showed that the lowest state in terms of energy are mostly of Frenkel type, whereas the higher lying states are CT ones which mix with embedded Frenkel type states. The first band of the absorption spectrum was explained by inclusion of vibrational motions within the stacks which induce an intensity gain of the first excited state. This intensity was not explainable by using the undistorted stacks. Also relaxations at the crystal surface might play a role, but are experimentally not explainable.
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.
We present a theoretical study on exciton–exciton annihilation (EEA) in a molecular dimer. This process is monitored using a fifth-order coherent two-dimensional (2D) spectroscopy as was recently proposed by Dostál et al. [Nat. Commun. 9, 2466 (2018)]. Using an electronic three-level system for each monomer, we analyze the different paths which contribute to the 2D spectrum. The spectrum is determined by two entangled relaxation processes, namely, the EEA and the direct relaxation of higher lying excited states. It is shown that the change of the spectrum as a function of a pulse delay can be linked directly to the presence of the EEA process.
We present a theoretical study on exciton–exciton annihilation (EEA) in a molecular dimer. This process is monitored using a fifth-order coherent two-dimensional (2D) spectroscopy as was recently proposed by Dostál et al. [Nat. Commun. 9, 2466 (2018)]. Using an electronic three-level system for each monomer, we analyze the different paths which contribute to the 2D spectrum. The spectrum is determined by two entangled relaxation processes, namely, the EEA and the direct relaxation of higher lying excited states. It is shown that the change of the spectrum as a function of a pulse delay can be linked directly to the presence of the EEA process.
The present work presents investigations on energy and charge transport properties in organic crystals. Chapter 4 treats exciton transport in anthracene, which is an example for weakly coupled π-systems. The electronic coupling parameter is evaluated by the monomer transition density approach. With these and the reorganization energy hopping rates are calculated in the framework of the Marcus theory. Together with the knowledge of the crystal structure, these allow us to calculate the experimental accessible exciton diffusion lengths, whose isotropic part fits nicely within the scattering of experimental values found in the literature. Furthermore, the anisotropy of the exciton diffusion lengths is reproduced qualitatively and quantitatively correct. This chapter also contains studies about electron and hole transport in both polymorphs (α and β) of perylene. Reorganization energies as well as diffusion coefficients for both crystal structures and types of charge transport were calculated. The best transport is hole transport in β-perylene, but it is strongly isotropic. The preferred transport direction is along the b-axis of the unit cell with couplings of greater than 100 meV. However, there is no transport along the c-axis. The diffusion constant in b-direction is bigger by two orders of magnitude than in c-direction (62.7•10-6 m2/s vs. 0.4•10-6 m2/s). Charge transport is calculated to be strongly anisotropic for holes as well as electrons in both modifications. To verify these results experimental electron mobilities have been compared to the simulations. Good agreement was found with errors of less than 27%. As it was shown above, the calculation and measurement of transport properties between weakly coupled systems is possible. However, it is difficult to exactly determine the quality of the electronic coupling. For this reason a collaboration about strongly interacting π-systems was started between us and the research group of Prof. Ingo Fischer. There, [2.2]paracyclophanes and its derivates were investigated to show how hydroxyl substitution influences absorption properties. Overall, a combination of SCS-MP2 and SCS-CC2 performs best to address the description of geometric and electronic structures for both ground and excited states of these model systems as well as their parent compounds benzene and phenol. Only [2.2]paracyclophane shows a double minimum potential regarding a twist and shift motion between the benzene/phenol subunits towards each other. All other systems are less flexible due to their substitution pattern. Almost all [2.2]paracyclophanes display minor changes in their geometric structure upon excitation to the S1 state: The inter-ring distance shortens, but qualitatively they keep their shift and twist characteristics, although the extent of these deformations diminishes. The exception is p-DHPC, which turns from a shifted ground state structure into a twisted excited state structure. Consequently, the intensity of the 0-0 transition cannot be observed experimentally due to small Franck-Condon factors and impurities of o-DHPC. In the present thesis, the structures and their changes due to excitation are explained by electrostatic potentials as well as antibonding (bonding) HOMO (LUMO) orbitals. Adiabatic excitation energies have been corrected by ZPEs and result in accuracies with errors smaller than 0.1 eV. Note that corrections on the B3LYP level worsen the results and one has to apply SCS-CC2 to achieve this accuracy. These calculations allow an interpretation of the experimental [1+1]REMPI spectra. Band progressions of the twist, shift and breathing of the [2.2]paracyclophane skeleton vibrations have been identified and show good agreement to the experiment. This work shows that the substitution pattern in [2.2]paracyclophanes can have a significant impact on spectroscopic properties. Because these properties are directly linked to the transport properties of these materials, the hereby gained insight can be used to design materials with customized transport properties. It was shown that the SCS-CC2 method is very appropriate to predict the interaction between the π-systems
Im Rahmen dieser Dissertation wurden optische Eigenschaften von halbleitenden, einwandigen Kohlenstoffnanoröhren (SWNTs) der (6,5)-Chiralität untersucht. Dies gelang durch Ensemblemessungen aber vor allem durch den Aufbau eines Mikroskops zur Messung an einzelnen SWNTs. Dieses Einzel- SWNT-Mikroskop ermöglichte nebst „normaler“ Bildgebung durch Sammlung und Abbildung der nahinfraroten Photolumineszenz (PL) der (6,5)-SWNTs auch die spektral- und zeitaufgelöste Untersuchung der PL. Durch Verwendung von Dichtegradientenultrazentrifugation (DGU) zur chiralen Aufreinigung des SWNT-Rohmaterials konnten alle Messungen unter Minimierung des störenden Einflusses von Aggregaten oder SWNTs anderer Chiralität durchgeführt werden. Untersucht und bestimmt wurde der Absorptionsquerschnitt und die Exzitonengröße, die PL-Eigenschaften aggregierter SWNTs und der Einfluß der Permittivität auf die PL einzelner SWNTs.