Refine
Has Fulltext
- yes (4)
Is part of the Bibliography
- yes (4)
Year of publication
- 2016 (4) (remove)
Document Type
- Doctoral Thesis (4) (remove)
Keywords
- Quantendynamik (2)
- Aggregat <Chemie> (1)
- Coherent two-dimensional Spectroscopy (1)
- Dimension 2 (1)
- Exziton (1)
- Exzitonen (1)
- Femtosekundenspektroskopie (1)
- Globale Optimierung (1)
- Kurzzeitphysik (1)
- Ladungstransfer (1)
- Molecular Structure (1)
- Molekularbewegung (1)
- Molekülstruktur (1)
- Physikalische Chemie (1)
- Polychromophores System (1)
- QM/MM (1)
- Quantenmechanik (1)
- Reaktionsmechanismus (1)
- Reaktionspfadsuche (1)
- Relaxation Dynamics (1)
- Spektroskopie (1)
- Transient-Absorption Sectroscopy (1)
- Ultrafast Spectroscopy (1)
- Ultrakurzzeitspektroskopie (1)
- Verbindungen (1)
- Wellenpaket (1)
- Zweidimensionale elektronische Spektroskopie (1)
- organische Halbleiter (1)
- quantum dynamics (1)
Institute
- Institut für Physikalische und Theoretische Chemie (4) (remove)
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 den spektralen Signaturen molekularer Aggregate sowie mit ihrer Wellenpakets- und Populationsdynamik in angeregten Zuständen unter dem Einfluss externer Störungen und photoinduzierter Asymmetrie.
Hierzu werden quantendynamische numerische Berechnungen mit der Multi-Configuration Time-Dependent Hartree-Methode
durchgeführt, um die angesprochenen Prozesse zu charakterisieren.
Durch die Konzentration auf Modellrechnungen sind die qualitativen Ergebnisse dieser Arbeit auf viele Systeme übertragbar.
Zunächst widmet sich die Arbeit den linearen UV/Vis-Absorptions- und Emissionsspektren von Aggregaten.
Hier zeigt sich, dass die Anzahl der Größen, die ein Absorptionsspektrum bestimmen -- etwa die Anzahl der Chromophore, ihre geometrischen Anordnung und die elektronische Kopplung zwischen ihnen -- zu groß ist, um ihre numerischen Werte eindeutig aus den Spektren bestimmen zu können.
Insbesondere können sich die Auswirkungen der Aggregatgröße und der Kopplungsstärke gegenseitig so beeinflussen, dass die Form der Absorptionsbande bei sehr unterschiedlichen Systemen nahezu identisch ist. Daraus ergeben sich Schwierigkeiten bei der Interpretation experimenteller Spektren, insbesondere von selbst-aggregierten Oligomeren, deren Größe unbekannt ist.
Es ist daher notwendig, entweder die elektronische Kopplung oder die Anzahl der Monomere in einem Aggregat
durch andere experimentelle Methoden unabhängig zu bestimmen. Ist die Aggregatgröße jedoch bekannt,
können die Absorptionsspektren sehr wohl zur Bestimmung anderer Eigenschaften des Systems herangezogen werden.
Dies wird durch die Untersuchung der Spektren kovalent gebundener zyklischer Aggregate aus drei und vier cis-Indolenin-Squarain-Molekülen als Beispiel für Systeme mit bekannter Größe dargestellt.
Das zweite Hauptthema der Arbeit ist die Populationsdynamik in angeregten Zuständen molekularer Aggregate.
Dazu werden numerische Rechnungen an Dimeren, Pentameren und Nonameren durchgeführt. Eine Asymmetrie, sei es im System selbst oder am Wellenpaket, das durch die Anregung entsteht, kann dazu führen, dass ein einzelnes Monomer dauerhaft bevorzugt populiert ist. Wenn durch eine externe Störung die Energie des angeregten Zustands bestimmter Monomere für eine gewisse Zeit erhöht ist,
kommt es zu einer Lokalisation der Population in diesem energetisch höheren Zustand.
In einem System mit weiteren internen Freiheitsgraden wird die Population auf benachbarte Monomere übertragen,
wenn der Betrag der Energieverschiebung des gestörten Zustands mit dem Abstand der Schwingungsniveaus zusammenfällt.
Der anfängliche Lokalisierungseffekt ist darüber hinaus zustandsspezifisch:
Er wird durch die Überlappintegrale der Schwingungskomponenten der Wellenfunktion in den diabatischen angeregten elektronischen Zuständen bestimmt.
Durch die Kombination von zwei Laserpulsen kann auch ein Wellenpaket in den angeregten Zuständen erzeugt werden,
dessen Symmetrieachsen nicht mit denen der Potentialflächen des Systems zusammenfallen.
Dadurch, dass hier die Asymmetrie schon im Wellenpaket vorliegt, kann es auch ohne äußere Störung zu einer Lokalisation der Population auf einem Monomer kommen.
Diese Arbeit befasst sich mit zeitaufgelösten Prozessen in molekularen Systemen. Dabei wurde sowohl die Wellenpaketdynamik nach Photoanregung betrachtet als auch spektrale Eigenschaften mittels Absorptions- und zweidimensionaler Spektroskopie untersucht.
Zunächst widmet sich die Arbeit der Wellenpaket- und Populationsdynamik in zwei diabatischen, gekoppelten Zuständen. Nach impulsiver Anregung aus dem zu Beginn besetzten Zustand treten in der Populationsdynamik zwei deutlich verschiedene Oszillationen auf. Der langsamer variierende Populationstransfer besitzt die Periodendauer der Vibrationsbewegung und ist auf einen Wechsel der Zustände beim Durchlaufen des Wellenpakets durch die Kreuzungsregion der diabatischen Potentiale zurückzuführen. Die ultraschnelle Komponente mit einer Periodendauer von etwa 4 fs lässt sich als eine Art Rabi-Oszillation beschreiben, die durch die (zeitunabhängige) Kopplung hervorgerufen wird. Sie wurde mit Hilfe von analytischen Berechnungen ausführlich charakterisiert. Damit dieser Prozess auftreten kann müssen mehrere Bedingungen erfüllt werden: Das Wellenpaket muss über die Dauer der Oszillationen annähernd örtlich lokalisiert bleiben; dies ist an den Umkehrpunkten der Wellenpaketsbewegung der Fall. Die Amplitude der Oszillationen in den Populationen ist proportional zum Verhältnis der Kopplung zum Energieabstand der Zustände. Deshalb muss an den stationären Stellen die Kopplung groß im Vergleich zum Energieabstand sein. Die Amplitude der Oszillationen hängt außerdem von dem Populationsverhältnis und den Phasen der Komponenten des Wellenpakets in den beiden Zuständen ab. Die ultraschnellen Oszillationen bleiben auch in mehrdimensionalen Systemen mit unterschiedlichen Vibrationsfrequenzen je Freiheitsgrad erhalten.
Das gleiche Modell wurde benutzt, um Ladungstransferprozesse mittels linearer und 2D-Spektroskopie zu untersuchen. Eine Kopplung an die Umgebung wurde, aufbauend auf einer Quanten-Master-Gleichung in Markov-Näherung, wellenfunktionsbasiert mittels eines Quantum-Jump-Algorithmus mit expliziter Dephasierung beschrieben. Dabei findet mit vorher definierten Wahrscheinlichkeiten zu jedem Zeitschritt einer von drei stochastischen Prozessen statt. Neben kohärenter Propagation können Sprünge in einen anderen Eigenzustand des Systems und Dephasierungen auftreten. Zwei Dissipationsparameter spielen dabei eine Rolle. Dies ist zum einen die Stärke der System-Bad-Kopplung, welche die Gesamtrate der Energierelaxation beschreibt. Weiterhin beeinflusst die Dephasierungskonstante den Verlust kohärenter Phasen ohne Energieänderung. Fallenzustände wurden identifiziert, die durch sehr geringe Sprungraten in niedrigere Zustände charakterisiert sind. Die Langlebigkeit kann durch die Form der Eigenfunktionen erklärt werden, die eine deutlich andere Wahrscheinlichkeitsverteilung als die der Nicht-Fallenzustände besitzen. Dadurch werden die in die Sprungraten eingehenden Matrixelemente klein. Das Absorptionsspektrum zeigt Peaks an der Stelle der Fallenzustände, da nur die Eigenfunktionen der Fallenzustände große Franck-Condon-Faktoren mit der Anfangswellenfunktion besitzen. Verschiedene Kombinationen der Dissipationsparameter führen zu Änderungen der relativen Peakintensitäten und der Peakbreiten.
Die 2D-Spektren des Ladungstransfersystems werden störungstheoretisch über die Polarisation dritter Ordnung berechnet. Sie zeigen viele eng nebeneinander liegende Peaks in einer schachbrettmusterförmigen Anordnung, die sich auf Übergänge unter Mitwirkung der Fallenzustände zurückführen lassen. Höhere System-Bad-Kopplungen führen aufgrund der effizienten Energiedissipation zu einer Verschiebung zu kleineren Energien. Peaks, die mit schneller zerfallenden Fallenzuständen korrespondieren, bleichen schneller aus. Höhere Dephasierungskonstanten resultieren in verbreiterten Peaks. Um den Einfluss der Dissipation genauer zu charakterisieren, wurden gefilterte 2D-Spektren betrachtet. Dazu wurden Ausschnitte der Polarisation dritter Ordnung zu verschiedenen Zeiten fouriertransformiert. Längere Zeiten führen zu einer effektiveren Energierelaxation entlang der entsprechenden Zeitvariablen. Die Entvölkerung der höher liegenden Zustände lässt sich somit zeit- und energieaufgelöst betrachten. Weiterhin wurde gezeigt, dass sich der Zerfall eines einzelnen Peaks mit dem Populationsabfall des damit korrespondierenden Eigenzustandes in Einklang bringen lässt, obwohl die Zuordnung der Peaks im 2D-Spektrum zu Übergängen zwischen definierten Eigenzuständen nicht eindeutig ist.
Mit dem benutzten eindimensionalen Modell können auch Ladungstransferprozesse in organischen gemischtvalenten Verbindungen beschrieben werden. Es wurde die Frage untersucht, welche Prozesse nach einem optisch induzierten Energietransfer in solchen Systemen ablaufen. Experimentelle Daten (aufgenommen im Arbeitskreis von Prof. Lambert) deuten auf eine schnelle interne Konversion (IC) gefolgt von Thermalisierung hin. Um dies theoretisch zu überprüfen, wurden Absorptionsspektren bei verschiedenen Temperaturen berechnet und mit den gemessenen transienten Spektren verglichen. Es findet sich, abhängig von der Stärke der elektronischen Kopplung, eine sehr gute bis gute Übereinstimmung, was die Annahme eines schnellen ICs stützt.
Im letzten Teil der Arbeit wurden vibronische 2D-Spektren von molekularen Aggregaten betrachtet. Dazu wurde die zeitabhängige Schrödingergleichung für ein Monomer-, Dimer- und Trimersystem mit der Multi-Configuration Time-Dependent Hartree-Methode gelöst und die Polarisation nicht-störungstheoretisch berechnet. Der Hamiltonoperator des Trimers umfasst hierbei sieben gekoppelte elektronische Zustände und drei bzw. sechs Vibrationsfreiheitsgrade. Der betrachtete Photonenecho-Beitrag der Polarisation wurde mittels phasencodierter Laserpulse extrahiert. Die resultierenden Spektren sind geometrieabhängig, ein Winkel zwischen den Übergangsdipolmomenten der Monomere von 0° (180°) resultiert in einem H-Aggregat (J-Aggregat). Die Lage und Intensität der Peaks im rein elektronischen Trimer wurde analytisch erläutert. Die Spektren unter Einbeziehung der Vibration zeigen eine ausgeprägte vibronische Struktur. Es wurde gezeigt, wie die Spektren für höhere Aggregationsgrade durch die höhere Dichte an vibronischen Zuständen komplexer werden. Im J-Aggregat ist mit zunehmender Aggregation eine stärkere Rotverschiebung zu sehen. Das Spektrum des H-Aggregats zeigt eine im Vergleich zum J-Aggregat kompliziertere Struktur. Die Verwendung zweier Vibrationsfreiheitsgrade je Monomer führt zu Spektren mit überlappenden Peaks und einer zusätzlichen vibronischen Progression. Der Vergleich von Spektren verschiedener Mischungen von Monomer, Dimer und Trimer, entsprechend einem von Temperatur und Konzentration abhängigen Aggregationsgrad, zeigt den Einfluss dieser experimentellen Faktoren. Schließlich wurden mögliche Ansätze aufgezeigt, anhand der Spektren auf den Aggregationsgrad zu schließen.
Time-resolved spectroscopy allows for analyzing light-induced energy conversion and
chromophore–chromophore interactions in molecular systems, which is a prerequisite in
the design of new materials and for improving the efficiency of opto-electronic devices.
To elucidate photo-induced dynamics of complex molecular systems, transient absorption
(TA) and coherent two-dimensional (2D) spectroscopy were employed and combined
with additional experimental techniques, theoretical approaches, and simulation models
in this work.
A systematic series of merocyanines, synthetically varied in the number of chromophores
and subsitution pattern, attached to a benzene unit was investigated in cooperation with
the group of Prof. Dr. Frank Würthner at the University of Würzburg. The global analysis
of several TA experiments, and additional coherent 2D spectroscopy experiments, provided
the basis to elaborate a relaxation scheme which was applicable for all merocyanine
systems under investigation. This relaxation scheme is based on a double minimum on the
excited-state potential energy surface. One of these minima is assigned to an intramolecular
charge-transfer state which is stabilized in the bis- and tris-chromophoric dyes by
chromphore–chromophore interactions, resulting in an increase in excited-state lifetime.
Electro-optical absorption and density functional theory (DFT) calculations revealed a
preferential chromophore orientation which compensates most of the dipole moment of
the individual chromophores. Based on this structural assignment the conformationdependent
exciton energy splitting was calculated. The linear absorption spectra of the
multi-chromophoric merocyanines could be described by a combination of monomeric and
excitonic spectra.
Subsequently, a structurally complex polymeric squaraine dye was studied in collaboration
with the research groups of Prof. Dr. Christoph Lambert and Prof. Dr. Roland Mitric
at the University of Würzburg. This polymer consists of a superposition of zigzag and
helix structures depending on the solvent. High-level DFT calculations confirmed the previous
assignment that zigzag and helix structures can be treated as J- and H-aggregates,
respectively. TA experiments revealed that in dependence on the solvent as well as the
excitation energy, ultrafast energy transfer within the squaraine polymer proceeds from
initially excited helix segments to zigzag segments or vice versa. Additionally, 2D spectroscopy
confirmed the observed sub-picosecond dynamics. In contrast to other conjugated
polymers such as MEH-PPV, which is investigated in the last chapter, ultrafast
energy transfer in squaraine polymers is based on the matching of the density of states
between donor and acceptor segments due to the small reorganization energy in cyanine-like
chromophores.
Finally, the photo-induced dynamics of the aggregated phase of the conjugated polymer
MEH-PPV was investigated in cooperation with the group of Prof. Dr. Anna Köhler at the University of Bayreuth. Our collaborators had previously described the aggregation of MEH-PPV upon cooling by the formation of so-called HJ-aggregates based on exciton
theory. By TA measurements and by making use of an affiliated band analysis distinct
relaxation processes in the excited state and to the ground state were discriminated. By
employing 2D spectroscopy the energy transfer between different conjugated segments
within the aggregated polymer was resolved. The initial exciton relaxation within the
aggregated phase indicates a low exciton mobility, in contrast to the subsequent energy
transfer between different chromophores within several picoseconds.
This work contributes by its systematic study of structure-dependent relaxation dynamics
to the basic understanding of the structure-function relationship within complex
molecular systems. The investigated molecular classes display a high potential to increase
efficiencies of opto-electronic devices, e.g., organic solar cells, by the selective choice of
the molecular morphology.