TY - THES A1 - Schaupp, Thomas T1 - Quantendynamik korrelierter Elektronen- und Kernbewegung T1 - Quantum dynamics of correlated electron and nuclear motion N2 - Im Rahmen dieser Arbeit werden unterschiedliche Aspekte der korrelierten Elektronen-Kerndynamik, anhand verschiedener Modellsysteme untersucht. Dabei wird vor allem auf den Vergleich numerisch exakter und approximativer Methoden zur Beschreibung der Wellenpaketdynamik eingegangen, wobei bei letzterem das Augenmerk auf der Born-Oppenheimer (BO) Näherung liegt. Die verwendeten Modellsysteme erlauben es, die gekoppelte Elektronen-Kern-Dynamik exakt zu beschreiben. Die daraus gewonnenen Ergebnisse dienen als Referenz für den Vergleich mit den Näherungsmethoden. Im ersten Teil der Arbeit wird die Dynamik eines Wellenpakets in der Umgebung einer Konischen-Durchschneidung (CI) untersucht, wobei die Beschreibung des Wellenpakets quantenmechanisch und durch die klassische Mechanik im Phasenraum erfolgt. Im zweiten Teil wird die Wahrscheinlichkeitsflussdichte untersucht. Zuerst wird ein Fall konstruiert, in welchem die Bewegung im elektronischen Grundzustand stattfindet, sodass die Bedingungen der BO Näherung erfüllt sind. Dabei wird vor allem auf das Verschwinden der elektronischen Wahrscheinlichkeitsflussdichte innerhalb der BO Näherung eingegangen. Im weiteren Verlauf werden die Flussdichten in der Umgebung einer CI untersucht, wobei unterschiedliche Situationen modelliert werden. Im dritten Teil wird die Berechnung des elektronischen Impulserwartungswerts innerhalb der BO Näherung untersucht. Dieser verschwindet innerhalb der BO Näherung, wenn man diesen direkt berechnet (Geschwindigkeitsform), während man über das Ehrenfest Theorem (Längenform) sehr gute Werte erhält. Im vierten Teil wird eine neue Flussdichte, die Translationsflussdichte, vorgestellt. Diese ergibt sich aus der Überlegung, dass die Geschwindigkeitsform des Impulserwartungswerts durch die Wahrscheinlichkeitsflussdichte ausgedrückt werden kann. Demnach muss auch die Längenform einer Flussdichte entsprechen und man erhält die Translationsflussdichte. N2 - In the present work different aspects of the correlated electron nuclear dynamics are investigated based on different model systems. The focus is set on the comparison of numerically exact and approximate methods, especially on the Born-Oppenheimer (BO) approximation, for the description of wave packet dynamics. The employed model systems allows a numerically exact solution for the coupled dynamics which serves as a reference for the comparison with the approximate description. In the first part of this thesis, wave packet dynamics are studied in the vicinity of a conical intersection (CI). The wave packet motion is treated quantum mechanically and also using classical mechanics. In the second part of this work the probability flux density is examined. A case is constructed where nearly the complete dynamics take place in the electronic ground state, so that the BO approximation is valid and the disappearance of the electronic probability flux density within the BO approximation is discussed. After that flux densities are studied in the vicinity of a CI whereby different situations are constructed In the third part of this thesis, the calculation of the electronic momentum expectation value within the BO approximation is examined. The electronic momentum expectation value disappears within the BO approximation if it is directly calculated (velocity form). The Ehrenfest theorem offers the possibility to calculate the momentum expectation value via the expectation value of the position operator (length form). The length form can be used to calculated an electronic momentum expectation value within the BO approximation which is in very good agreement with the one obtained from the coupled dynamics. In the fourth part of the thesis a new flux density, the translational flux density, is presented. This flux density emerges from the property that the velocity form of the momentum expectation value can be expressed by the probability flux density. Accordingly, the length form must also result from a flux density which identifies the translation flux density. KW - Quantenmechanik KW - Flussdichte KW - Born-Oppenheimer-Näherung KW - Wellenpaket KW - Konischen-Durschneidung KW - Shin-Metiu-Modell KW - Elektronen-Kernbewegung Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-237743 ER - TY - THES A1 - Süß, Jasmin T1 - Theoretische Untersuchungen an molekularen Aggregaten: 2D-Spektroskopie und Exzitonendynamik T1 - Theoretical studies on molecular aggregates: 2D spectroscopy and exciton dynamics N2 - Diese Dissertation beschäftigt sich mit der Exzitonendynamik molekularer Aggregate, die nach Mehrphotonen-Anregung auf ultrakurzer Zeitskala stattfindet. Hierbei liegt der Fokus auf der Charakterisierung der Exziton-Exziton-Annihilierung (EEA) mithilfe von zweidimensionaler optischer Spektroskopie fünfter Ordnung. Dazu werden zwei verschiedene Modellsysteme implementiert: Das elektronische Homodimer und das elektronische Homotrimer-Modell, wobei Letzteres eine Erweiterung des Dimer-Modells darstellt. Die Kopplung des quantenmechanischen Systems an die Umgebung wird mithilfe des Quantum-Jump-Ansatzes umgesetzt. Besonderes Interesse kommt der Analyse des Signals fünfter Ordnung in Abhängigkeit der Populationszeit T zu. Anhand des Dimer-Modells als kleinstmögliches Aggregat lassen sich bereits gute Vorhersagen auch über das Verhalten größerer molekularer Aggregate treffen. Der Zerfall des oszillierenden Signals für lange Populationszeiten korreliert mit der EEA. Dies zeigt, dass die zweidimensionale optische Spektroskopie genutzt werden kann, um den Annihilierungsprozess zu charakterisieren. Innerhalb des Modells des Dimers wird weiterhin der Einfluss der Intraband-Relaxation untersucht. Zunehmende Intraband-Relaxation verhindert den Austausch zwischen den lokalen Zuständen, der essentiell für den Annihilierungsprozess ist, und die EEA wird blockiert. Das elektronische Trimer-Modell erweitert das Dimer-Modell um eine Monomereinheit. Somit befinden sich die Exzitonen im Anschluss an die Anregung nicht mehr unvermeidlich nebeneinander. Es gibt somit eine Konfiguration, bei der sich die Exzitonen zunächst zueinander bewegen müssen, bevor die Startbedingung des Annihilierungsprozesses gegeben ist. Dieser zusätzliche Schritt wird auch Exzitonendiffusion genannt. Die Ergebnisse dieser Arbeit legen nahe, dass das erwartete Verhalten nur zu sehr kurzen Zeiten im Femtosekundenbereich auftritt und somit die Zeitskala der Exzitonendiffusion im Falle des Trimers nicht sichtbar wird. Es bedarf demnach eines größeren Modellsystems, bei dem sich der Effekt der zeitverzögert eintretenden EEA deutlich in der Zerfallsdynamik manifestieren kann. N2 - This work addresses the exciton dynamics of molecular aggregates which occur after femtosecond multi-photon laser excitation. Thereby, the focus is on the characterization of exciton-exciton annihilation (EEA) via fifth order two dimensional optical spectroscopy. Two model systems are employed: the electronic homodimer model and the electronic homotrimer model, where the latter one is an extension of the dimer system. The systems are coupled to the surrounding. In the numerical calculation, the system-bath interaction is realized via the quantum jump approach. Particular attention is payed to energy-integrated spectra as a function of the population time T. The dimer is the smallest molecular aggregate, but it is a good reference system if larger aggregates are supposed to be understood. The decay of the oscillating fifth-order signal corresponds to the EEA. This indicates that two dimensional optical spectroscopy can be used to monitor the annihilation process. Furthermore, the effect of intraband relaxation is studied within the dimer model. The results display that increasing the intraband relaxation inhibits the population transfer between the localized states of the system. This blocks the EEA. In extending the dimer model system by one monomer unit, one obtains the electronic trimer model system. Within this model, the situation after excitation differs from the one in the dimer model. The excitons do not exclusively reside next to each other so that EEA is immediately possible. In that case, the excitons have to diffuse to each other before they eventually meet and the annihilation process starts. The results suggest that the expected properties are merely correct at very short times around a few femtoseconds. Within the trimer model, the additional time scale for the exciton diffusion doesn't show in the results. In particular, it requires a larger model system for the effect of the delayed EEA to be seen in the regarded signal. KW - Molekulardynamik KW - Quantenmechanik KW - Spektroskopie KW - Exziton KW - Exziton-Exziton-Annihilierung KW - Quantum-Jump-Ansatz KW - Wellenpaketdynamik Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-247136 ER - TY - THES A1 - Sauer, Susanne T1 - Implementation and Application of QM/MM Hybrid Methods T1 - Implementierung und Anwendung von QM/MM-Hybridmethoden N2 - Within this work, an additive and a subtractive QM/MM interface were implemented into CAST. The interactions between QM and MM system are described via electrostatic embedding. Link atoms are used to saturate dangling bonds originating from the separation of QM and MM system. Available energy evaluation methods to be combined include force fields (OPLSAA and AMBER), semi-empirical programs (Mopac and DFTB+), and quantum-chemical methods (from Gaussian, Orca, and Psi4). Both the additive and the subtractive interface can deal with periodic boundary conditions. The subtractive scheme was extended to enable QM/QM, three-layer, and multi-center calculations. Another feature only available within the subtractive interface is the microiteration procedure for local optimizations. The novel QM/MM methods were applied to the investigation of the reaction path for the complex formation between rhodesain and K11777. Benchmark calculations show a very good agreement with results from Gaussian-ONIOM. When comparing the relative energies obtained with different options to a computation where the whole system was treated with the “QM method” DFTB3, the electrostatic embedding scheme with option “delM3” gives the best results. “delM3” means that atoms with up to three bonds distance to the QM region are ignored when creating the external charges. This is done in order to avoid a double counting of Coulomb interactions between QM and MM system. The embedding scheme for the inner system in a three-layer calculation, however, does not have a significant influence on the energies. The same is true for the choice of the coupling scheme: Whether the additive or the subtractive QM/MM interface is applied does not alter the results significantly. The choice of the QM region, though, proved to be an important factor. As can be seen from the comparison of two QM systems of different size, bigger is not always better here. Instead, one has to make sure not to separate important (polar) interactions by the QM/MM border. After this benchmark study with singlepoint calculations, the various possibilities of CAST were used to approximate the solution of a remaining problem: The predicted reaction energy for the formation of the rhodesain-K11777 complex differs significantly depending on the starting point of the reaction path. The reason for this is assumed to be an inadequate adjustment of the environment during the scans, which leads to a better stabilization of the starting structure in comparison to the final structure. The first approach to improve this adjustment was performing the relaxed scan with a bigger QM region instead of the minimal QM system used before. While the paths starting from the covalent complex do not change significantly, those starting from the non-covalent complex become more exothermic, leading to a higher similarity of the two paths. Nevertheless, the difference of the reaction energy is still around 15 kcal/mol, which is far from a perfect agreement. For this reason, Umbrella Samplings were run. Here, the adjustment of the environment is not done by local optimizations like in the scans, but by MD simulations. This has the advantage that the system can cross barriers and reach different local minima. The relative free energies obtained by Umbrella Samplings with suitable QM regions are nearly identical, independently of the starting point of the calculation. Thus, \(\Delta A\) evaluated by these computations can be assumed to reproduce the real energy change best. An MD simulation that was started from the transition state in order to mimic a “real-time” reaction indicates a very fast adjustment of the environment during the formation of the complex. This confirms that Umbrella Sampling is probably better suitable to describe the reaction path than a scan, where the environment can never move strong enough to leave the current local minimum. N2 - In dieser Arbeit wurden ein additives und ein subtraktives QM/MM-Interface in CAST implementiert. Die Wechselwirkungen zwischen QM- und MM-System werden durch elektrostatische Einbettung beschrieben. Link-Atome dienen dazu, lose Bindungen abzusättigen, die durch die Trennung von QM- und MM-System entstehen. Als Methoden zur Energieberechnung, die kombiniert werden können, stehen Kraftfelder (OPLSAA und AM- BER), semiempirische Programme (Mopac und DFTB+) und quantenchemische Verfahren (aus Gaussian, Orca und Psi4) zur Verfügung. Sowohl das additive als auch das subtraktive Interface können mit periodischen Randbedingungen verwendet werden. Erweiterungen des subtraktiven Schemas ermöglichen Berechnungen mit QM/QM, drei Schichten o der mehreren QM-Zentren. Ebenfalls nur im subtraktiven Interface verfügbar ist die lokale Optimierung mittels Mikroiterationsschema. ... KW - Quantenmechanik KW - Molekularmechanik KW - QM/MM KW - Umbrella Sampling KW - Computational Chemistry KW - Theoretical Chemistry Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-243213 ER - TY - THES A1 - Wehner, Johannes T1 - Wellenfunktionsbasierte Analyse zweidimensionaler Spektren: Wellenpaketbewegung in Dimeren und Quantendiffusionsdynamik T1 - Wavefunction based analysis of two-dimensional spectra: vibronic wave-packet motion in dimers and quantum-state-diffusion dynamics N2 - Diese Arbeit befasst sich mit der störungstheoretischen Berechnung von zweidimensionalen Photonen-Echo-Spektren für das elektronische und vibronische Modell eines Homo- und Hetero-Dimers sowie für ein vibronisches Modell eines Monomers unter dem Einfluss einer System-Bad-Wechselwirkung. Bei der Analyse der Dimerspektren steht neben der Orientierungsmittelung der Polarisation dritter Ordnung der Unterschied zwischen elektronischen und vibronischen Spektren sowie der Vergleich der Spektren von Homo- und Hetero-Dimeren im Zentrum des Interesses. Bei der Analyse der Monomer-Spektren steht die Behandlung einer dissipativen Dynamik bzw. des vibrational-coolings innerhalb eines stochastischen Ansatzes im Vordergrund. Der erste Teil dieser Arbeit konzentriert sich auf die störungstheoretische Berechnung der Polarisation dritter Ordnung in Dimeren. Dabei werden alle Aspekte und Ergebnisse für verschiedene Geometrien der Übergangsdipolmomente analysiert und diskutiert. Die Berechnungen berücksichtigen dabei auch die zufällige Anordnung der Moleküle in der Probe. Die Zusammenhänge zwischen den 2D-Spektren und den Eigenschaften der Monomereinheiten, die Abhängigkeit der Intensitäten mancher Peaks von der zeitlichen Abfolge der Pulse sowie der Einfluss der elektronischen Kopplung und verschiedener Übergangsdipolmomente ermöglichen ein grundlegendes Verständnis der elektronischen Photonen-Echo-Spektren. Im elektronischen Dimer wird der Hetero-Dimer-Charakter durch verschiedene Monomeranregungsenergien sowie unterschiedliche Übergangsdipolmomente der Monomereinheiten bestimmt. Der Einfluss dieser Größen auf die Photonen-Echo-Spektren kann durch die Kombination einer detaillierten analytischen Betrachtung und numerischen Rechnungen anschaulich nachvollzogen werden. In der vibronischen Betrachtungsweise zeigt sich, dass die Spektren deutlich an Komplexität gewinnen. Durch die Vibrationsfreiheitsgrade vervielfachen sich die möglichen Übergänge im System und damit die möglichen Peakpositionen im Spektrum. Jeder Peak spaltet in eine Vibrationssubstruktur auf, die je nach ihrer energetischen Position mit anderen überlagern kann. Der Vergleich zwischen Homo- und Hetero-Dimer-Spektren wird durch die Wahl verschiedener Vibrationsfrequenzen und unterschiedlicher Gleichgewichtsabstände entlang der Vibrationskoordinaten erweitert. Die Berechnung des Orientierungsmittels erfolgt mit zwei verschiedenen Ansätzen. Zum einen wird das Mittel durch den numerischen sampling-Ansatz berechnet. Dabei werden Azimutal- und Polarwinkel in kleinen Winkelinkrementen abgetastet und für jede Kombination ein 2D-Spektrum berechnet. Die Einzelspektren werden anschließend gemittelt. Diese Methode erweist sich im Dimer als sehr effektiv. Zum anderen erlaubt die analytische Auswertung der Polarisation dritter Ordnung, das gemittelte Spektrum direkt in einer einzelnen Rechnung durch winkelgemittelte Gewichtungsfaktoren zu bestimmen. Bei der Berechnung der elektronischen 2D-Spektren ist diese Methode sehr leistungsfähig, da alle Ausdrücke analytisch bekannt sind. Für vibronische Systeme ist dieser Ansatz ebenfalls sehr leistungsstark, benötigt aber eine einmalige aufwendige Analyse vor der Berechnung. Trotz der deutlich erhöhten Anzahl an Zustandsvektoren, die propagiert werden müssen, ist diese Methode circa zweimal schneller als die direkte Mittelung mit der sampling-Methode. Im zweiten Teil konzentriert sich die Arbeit auf die Beschreibung eines Monomers, das sich in einer dissipativen Umgebung befindet. Dabei wird auf die Lösung einer stochastischen Schrödingergleichung zurückgegriffen. Speziell wird die sogenannte quantum-state-diffusion-Methode benutzt. Dabei werden nicht nur die Erwartungswerte für die Energie und den Ort, sondern auch die Polarisation dritter Ordnung – eine phasensensitive Größe – bestimmt. In der theoretischen Fragestellung wird dabei, ausgehend von der von-Neumann Gleichung, die Zeitentwicklung der reduzierten Dichtematrix durch die Integration einer stochastischen zeitabhängigen Schrödingergleichung reproduziert. In Rechnungen koppelt die Stochastik über die Erwartungswerte von Ort und Impuls die verschiedenen störungstheoretischen Korrekturen der Wellenfunktion miteinander. Die Spektren, die aus den numerischen Simulationen erhalten werden, spiegeln das dissipative Verhalten des Systems detailliert wider. Eine Analyse der Erwartungswerte von Ort und Energie zeigt, dass sich die einzelnen elektronischen Zustände wie gedämpfte harmonische Oszillatoren verhalten und jeweils einen exponentiellen Zerfall abhängig von der Dissipationskonstante zeigen. Dieser Teil der Arbeit erweitert vorausgehende Untersuchungen, bei denen ein vereinfachter Ansatz zu Einsatz kam, der die korrelierte Stochastik nicht berücksichtigte. N2 - This PhD-thesis is centered around the calculation of two-dimensional photon-echo spectra for different model systems. Two systems are investigated in detail, the electronic and vibronic homo- and hetero-dimer as well as a vibronic monomer unit treated as an open-quantum system. Dimer-spectra are obtained within an perturbative approach, which takes the random orientation of the molecules in a sample into account. The orientationally averaged spectra of electronic homodimers are influenced by the coupling strength between the monomer units and different dipole orientations. By analysing these spectra analytically, a fundamental understandig of the photon-echo spectra is obtained. The prediction of energetic positions and relative intensities of the spectral peaks is possible. Furthermore, it is possible to extract the dipole geometry of the dimer system by comparing different peak-intensities. For an even deeper insight, the oscillatory behaviour of some peaks as a function of the time ordering of laser pulses is analysed. Switching to the electronic heterodimer increases complexity. The different excitation energies of the monomer units and the different transition-dipole strengths influence the two-dimensional spectra. The energetics of the heterodimer can be understood similarly to the homodimer. Peak intensities are difficult to analyse due to the more complex system. In a next step one vibrational degree of freedom per monomer unit is included. This vibronic dimer shows a very dense set of eigenenergies leading to a manifold of peaks in the spectra. The underlying transitions contributing to one peak are in many cases hard to identify: Different vibrational signatures appear at the same position in the spectrum and the underlying vibronic transitions cannot be determined easily. In addition, in a heterodimer different frequencies and equilibrium distances are encountered, which have a huge impact on the spectra. All presented aspects are based on analytical treatment of the third-order polarisation or the numerical calculations of 2d photon-echo spectra. To perform the orientational average two different methods are introduced: a brute-force method (sampling method) and a method using analytically derived specific weighting factors. Within the analytical approach electronic spectra can be calculated very efficiently, though this method becomes more complex if vibrational degrees of freedoms are taken into account. To analytically average the spectra a huge number of statevectors needs to be propagated. Here the brute-force method comes into play. Sampling the orientation of the polarisation vector of the incoming fields over a set of discrete angles yields spectra for fixed orientations. The desired spectrum is obtained by taking the average over spectra for fixed orientations. This method is very effective, because averaging over a small set of spectra with fixed orientations (in the presented example five) is sufficient to yield a reliable result. Because the sampling-method is easy to implement and transferable to other systems, the minor time advantage gained in the analytical approach doesn’t compensate the demanding system-specific analytical treatment. The second main topic, which is addressed in this thesis, is a vibronic monomer coupled to its environment. The system-bath coupling leads to vibrational cooling within the electronic states. The influence of this dissipative dynamics on two-dimensional photon-echo spectra is analysed. Therfore a stochastic wave-function approach based on the von-Neuman equation and the quantum-state-diffusion method is used. Within this ansatz it is shown, that the spectra can be calculated perturbatively. In addition to our formerly published propagation scheme we take statistically correlated dynamics in the electronic states into account. Otherwise the spectra may show peaks at unphysical positions. The dissipative dynamics can be monitored by the expectation values of the spatial coordinate space and energy. It is shown, that the expectation values are not dependent on the propagation scheme employed. This is different for the calculation of the phase-sensitive third-oder polarisation, where the correlated approach leads to far better results. KW - Molekulardynamik KW - Vierwellenmischen KW - Quantenmechanik KW - Wellenpaket KW - Spektroskopie KW - stochastische Schrödingergleichung KW - System-Bad-Ansatz Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-163555 ER - TY - THES A1 - Reichert, Thorsten T1 - Classification and Reduction of Equivariant Star Products on Symplectic Manifolds N2 - This doctoral thesis provides a classification of equivariant star products (star products together with quantum momentum maps) in terms of equivariant de Rham cohomology. This classification result is then used to construct an analogon of the Kirwan map from which one can directly obtain the characteristic class of certain reduced star products on Marsden-Weinstein reduced symplectic manifolds from the equivariant characteristic class of their corresponding unreduced equivariant star product. From the surjectivity of this map one can conclude that every star product on Marsden-Weinstein reduced symplectic manifolds can (up to equivalence) be obtained as a reduced equivariant star product. N2 - Diese Doktorarbeit klassifiziert äquivariante Sternprodukte (Sternprodukte zusammen mit Quantenimpulsabbildungen) über die äquivariante de Rham Kohomologie. Diese Klassifizierung wird im Folgenen genutzt um ein Analogon der Kirwan-Abbildung zu konstruieren, welches ermöglicht die charakteristische Klasse von bestimmten reduzierten Sternprodukten auf Marsden-Weinstein reduzierten symplektischen Mannigfaltigkeiten direkt aus der äquivarianten charakteristischen Klasse des zugehörigen unreduzierten äquivarianten Sternprodukts zu erhalten. Die Surjektivität dieser Abbildung zeigt schließlich, dass jedes Sternprodukt auf einer Marsden-Weinstein reduzierten symplektischen Mannigfaltigkeit (bis auf Äquivalenz) als Reduktion eines äquivarianten Sternprodukts verstanden werden kann. T2 - Klassifizierung und Reduktion äquivarianter Sternprodukte auf symplektischen Mannigfaltigkeiten KW - Homologische Algebra KW - Differentialgeometrie KW - Quantenmechanik KW - symplectic geometry KW - deformation quantization KW - equivariant cohomology Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-153623 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 - Keß, Martin T1 - Wellenfunktionsbasierte Beschreibung der zweidimensionalen vibronischen Spektroskopie von molekularen Aggregaten und Ladungstransfersystemen T1 - Wave-function based description of the two-dimensional vibronic spectroscopy of molecular aggregates and charge-transfer systems N2 - 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. N2 - This work studies time-resolved phenomena in molecular systems. Both, the wave-packet dynamics after photoexcitation and the spectral properties, examined via absorption and two-dimensional spectroscopy, are regarded. First, the wave-packet and population dynamics in two coupled diabatic states are considered. After an impulsive excitation from the initially populated state, two significantly different oscillatory features are visible in the population dynamics. The slower varying population transfer follows the oscillation period of the vibrational motion and results from the diabatic transition when the wave-packet passes through the crossing region of the respective potentials. The ultrafast oscillatory component with an oscillation period of about 4 fs can be described as a Rabi-like oscillation induced by the (time-independent) coupling. It is characterized in detail via analytic calculations. For this contribution to be visible, some conditions have to be met: The wave-packet needs to be spatially localized during the duration of the oscillations. This is the case at the classical turning points of the wave-packet motion. The oscillations' amplitude seen in the populations is proportional to the ratio between the coupling and the energetic gap between the involved states. This means that the coupling needs to be large compared to the energy separation at the points where the wave-packet is stationary. Additionally, the amplitude depends on the relative populations and the phases of the wave-packet components in the two states. The ultrafast oscillations persist in systems of higher dimensionality with different vibrational frequencies in each degree of freedom. The same model is used to examine charge-transfer processes via linear and 2D spectroscopy. A coupling to the environment is described by a quantum-jump algorithm with explicit treatment of dephasing, based on a quantum-master equation in Markov approximation. At each time step, one of three stochastic processes takes place with a pre-defined probability. Besides coherent propagation, jumps into other eigenstates of the system and dephasing occur. Two dissipation parameters are of relevance. The first is the value of the system-bath coupling which influences the overall energy relaxation rate while, additionally, the dephasing constant causes a loss of phase coherence without energy relaxation. Trap states are identified, which are characterized by very low jump rates to lower states. Their slow decay can be explained by the shape of their respective eigenfunctions, which possess a vastly different probability density than eigenstates of the non-trap states. This results in small matrix elements entering in the equations for the jump rates. The absorption spectrum exhibits peaks at the energies of the trap states because only the trap states' eigenfunctions lead to large Franck-Condon factors with the initial wave function. Different values of the dissipation parameters lead to changes in the relative peak intensities and peak widths. The 2D spectra of the charge-transfer system are calculated via the third-order polarization. They show many close lying peaks in a chessboard like distribution. The peaks can be traced back to transitions involving the trap states. Higher values of the system-bath coupling lead to a shift to lower energies because of the more efficient energy dissipation. Peaks corresponding to faster decaying trap states show more substantial loss in intensity as compared to other peaks. Higher values of the dephasing constant result in broader peaks. To better characterize the influence of the dissipation, we consider filtered 2D spectra. Therefore, cuts of the third-order polarization at different times are Fourier-transformed separately. Cuts at later times map the more effective energy relaxation along the respective time-variable. Via this technique the de-population of higher lying states can be monitored both in time and energy. Additionally it is shown that the decay of a specific peak can be related to the population decay of the corresponding eigenstate, even though the assignment of peaks in the 2D spectrum to transitions between eigenstates is not unique. The one-dimensional model can also be used to examine charge-transfer processes in organic mixed-valence compounds. Here, the question is, which processes take place after an optically induced energy transfer. Transient absorption spectra, recorded in the group of Prof. Lambert, hint to a fast internal conversion (IC) followed by thermalisation. To check this theoretically, absorption spectra at different temperatures are calculated and compared to the measured transient spectra. Depending on the value of the electronic coupling element, a very good to good agreement is found, supporting the existence of a fast IC process. The last part of this work considers vibronic 2D spectra of molecular aggregates. Therefore, the time-dependent Schrödinger equation is solved with the Multi-Dimentional Time-Dependent Hartree-method for a monomer, dimer and trimer system, and the polarization is calculated via a non-perturbative scheme. The trimer Hamiltonian consists of seven coupled electronic states and three or six vibrational degrees of freedom, respectively. The photon-echo contribution of the polarization is extracted via phase-coded laser pulses. This results in geometry dependent spectra: An angle between the monomer transition dipole moments of 0° (180°) leads to an H-aggregate (J-aggregate). In the purely electronic system, the location and intensities of the peaks are explained analytically. The spectra including vibrations show a rich vibronic structure. It is shown that spectra for higher degrees of aggregation are more complex because of the higher density of vibronic states. The J-aggregate is stronger red shifted in larger aggregates. The spectrum of the H-aggregate possesses a more complicated structure as compared to the J-aggregate spectrum. The inclusion of a second vibrational degree of freedom into each monomer results in spectra with overlapping peaks and an additional vibrational progression. Spectra of different mixtures of monomer, dimer and trimer are compared. Because the level of aggregation depends on temperature and concentration, this documents the influence of the experimental conditions on the 2D spectra. Finally, possible approaches to infer the degree of aggregation from the spectra are discussed. KW - Quantenmechanik KW - Ladungstransfer KW - Quantendynamik KW - Zweidimensionale elektronische Spektroskopie KW - Aggregat KW - Spektroskopie KW - Dimension 2 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-136458 ER - TY - THES A1 - Falge, Mirjam T1 - Dynamik gekoppelter Elektronen-Kern-Systeme in Laserfeldern T1 - Dynamics of Coupled Electron-Nuclei-Systems in Laser Fields N2 - Die vorliegende Arbeit beschäftigt sich mit der theoretischen Untersuchung zweier Themenkomplexe: der Erzeugung Hoher Harmonischer in Molekülen und dem Einfluss von gekoppelter Elektronen-Kern-Dynamik auf Ultrakurzpuls-Ionisationsprozesse und Quantenkontrolle. Während bei der Untersuchung der Hohen Harmonischen die Auswirkungen der Kernbewegung auf die Spektren im Mittelpunkt des Interesses stehen, wird bei der Analyse der gekoppelter Elektronen-Kern-Dynamik das Hauptaugenmerk auf die nicht-adiabatischen Effekte gerichtet, die auftreten, wenn Kern- und Elektronenbewegung sich nicht, wie es im Rahmen der Born-Oppenheimer-Näherung in der Quantenchemie häufig angenommen wird, voneinander trennen lassen. N2 - This work aims at the theoretical analysis of high harmonic generation in molecules and the influence of coupled electron and nuclear dynamics on ultra-short pulse ionization processes. KW - Nichtadiabatischer Prozess KW - Laserstrahlung KW - Quantenmechanik KW - Molekulardynamik KW - Quantendynamik KW - nicht-adiabatische Effekte KW - Hohe Harmonische KW - Photoelektronenspektroskopie KW - quantum dynamics KW - nonadiabatic effects KW - high harmonic generation KW - photoelectron spectroscopy Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-72889 ER - TY - THES A1 - Gräfe, Stefanie T1 - Laser-control of molecular dynamics T1 - Lasergesteuerte Kontrolle molekularer Dynamik N2 - In this work a new algorithm to determine quantum control fields from the instantaneous response of systems has been developed. The derived fields allow to establish a direct connection between the applied perturbation and the molecular dynamics. The principle is most easily illustrated in regarding a classical forced oscillator. A particle moving inside the respective potential is accelerated if an external field is applied acting in the same direction as its momentum (heating). In contrary, a deceleration is achieved by a field acting in the opposite direction as the momentum (cooling). Furthermore, when the particle reaches a classical turning point and then changes its direction, the sign of the field has to be changed to further drive the system in the desired way. The frequency of the field therefore is in resonance with the oscillator. This intuitively clear picture of a driven classical oscillator can be used for directing (or controlling) quantum mechanical wave packet motion. The efficiency of the instantaneous dynamics algorithm was demonstrated in treating various model problems, the population transfer in double well potentials, excitation and dissociation of selective modes, and the population transfer between electronic states. Although it was not tried to optimize the fields to gain higher yields, the control was found to be very efficient. Driving population transfer in a double well potential could be shown to take place with nearly 100% efficiency. It was shown that selective dissociation within the electronic ground state of HOD can be performed by either maximizing a selected coordinate's differential momentum change or the energy absorption. Concerning the population transfer into excited electronic states, a direct comparison with common control algorithms as optimal control theory and genetic algorithms was accomplished using a one-dimensional representation of methyl iodide. The fields derived from the various control theories were effective in transferring population into the chosen target state but the underlying physical background of the derived optimal fields was not obvious to explain. The instantaneous dynamics algorithm allowed to establish a direct relation between the derived fields and the underlying molecular dynamics. Bound-to-bound transitions could be handled more effectively. This was demonstrated on the sodium dimer in a representation of 3 electronic states being initially in its vibronic ground state. The objective was to transfer population into a predefined excited state. Choosing the first or the second state as a target, the control fields exhibited quite different features. The pulse-structure is related to the excited state wave packet, moving in, and out of the Franck-Condon region. Changing the control objective, the derived control field performed pure electronic transitions on a fast time-scale via a two-step transition. Futhermore, orientational effects have been investigated. The overall-efficiency of the population transfer for differently oriented molecules was about 70 % or more if applying a control field derived for a 45° orientation. Spectroscopic methods to gain information about the outcome of the control process have been investigated. It was shown that pump/probe femtosecond ionization spectroscopy is suited to monitor time-dependent molecular probability distributions. In particular, time-dependent photoelectron spectra are able to monitor the population in the various electronic states. In the last chapter a different possibility of controlling molecules was regarded by investigating molecular iodine with a setup similar to the STIRAP (“Stimulated Raman Adiabatic passage”) scenario. The possibility to extend this technique to a fs-time scale was examined in theory as well as in experiments, the latter being performed by Dr. Torsten Siebert in the Kiefer group, University of Würzburg. It was shown that off-resonant excitation with implementation of the pulses with a higher intensity of the Stokes pulse as compared to the pump pulse - describing a so-called f-STIRAP like configuration - was shown to effectively transfer population into excited ground-state vibrational levels. This was theoretically underlined by comparing the numerically exact coupling case with the adiabatic picture. The process was described to run in the vicinity of adibaticity. A new model explaining the process by the system's vector rotating around the dressed state vector will be adopted in future calculations. Altogether, a new promising algorithm to control dynamical processes based on the instantaneous response has been developed. Because the derived control fields have been shown to be very efficient in selectively influencing molecules, it is to be expected that farther reaching applications can be realized in future investigations. N2 - In dieser Arbeit wurde ein neuer Algorithmus zur Bestimmung von Kontrollfeldern aus der instantanen Respons von Systemen auf die Wirkung von Laserfeldern entwickelt. Die damit berechneten Felder ermöglichen es, eine Verbindung zwischen der Störung durch das Laserfeld und der molekularen Dynamik herzustellen. Das Prinzip lässt sich an einem klassischen Oszillator veranschaulichen: Ein sich innerhalb dieses Oszillatorpotenzials bewegendes Teilchen wird durch ein externes Feld beschleunigt, wenn dieses und der Impuls des Teilchens in die gleiche Richtung weisen. Ein Abbremsen des Teilchens wird durch ein Feld erzielt, welches dem Impuls des Teilchens entgegen gerichtet ist. Wenn das Teilchen in dem Oszillator einen Umkehrpunkt erreicht und dort seine Richtung ändert, wird das Vorzeichen des Feldes an die neue Richtung angepasst: Die Frequenz des Feldes befindet sich in Resonanz mit der Oszillatorfreuqenz. Dieses klassische Bild der erzwungenen Schwingung eines Oszillators kann für die Kontrolle quantenmechanischer Wellenpaketbewegungen angewendet werden. Die Effizienz des Algorithmus' wurde an verschiedenen Problemen, wie dem Populationstransfer (PT) in Doppelminimum-Potenzialen, Anregung und Dissoziation selektiver Moden und den PT in unterschiedliche el. Zuständen aufgezeigt. Obwohl keine Optimierung der Felder bezüglich höherer Ausbeuten durchgeführt wurde, konnte eine hohe Effizienz der Prozesse nachgewiesen werden. Ein PT in Doppelminimum-Potentialen wurde nahezu vollständig erreicht. Selektive Dissoziation innerhalb des el. Grundzustandes des HOD-Moleküls wurde unter Verwendung zweier unterschiedlicher Methoden, der Maximierung der zeitlichen Änderung des Impulses oder der Energieabsorption einer Koordinate, erzielt. Bezüglich des PT in el. angeregte Zustände wurden bekannte Kontrollalgorithmen wie die Theorie der optimalen Kontrolle und genetischer Algorithmen mit dem in dieser Arbeit entwickelten Prinzip der instantanen Respons anhand einer 1D Darstellung des Methyliodids verglichen. Die aus den verschiedenen Theorien konstruierten Felder erzielten einen effektiven PT in den zuvor definierten Zielzustand, jedoch ist der dem zu Grunde liegende, physikalische Hintergrund nicht einfach zu beschreiben. Mit Hilfe des Instantanen-Respons-Algorithmus' konnte eine direkte Relation zwischen den Feldern und der molekularen Dynamik hergestellt werden. Anhand des Na2 in einer Darstellung von 3 elektronischen Zuständen sollte nur ein Zustand selektiv angeregt werden. Je nach Wahl des Zielzustandes zeigten sich deutliche Unterschiede. Selektive Anregung des 1. Zustandes erzeugte ein Feld bestehend aus einer Pulsfolge, die durch ein Wellenpaket im angeregten Zustand, welches sich in und aus dem Franck-Condon Fenster heraus bewegt, erklärt werden konnte. Anregung des 2. Zustandes führte zu einem Feld, welches nicht auf Vibration, sondern rein elektronischer Anregung in einem 2-Stufen-Prozess beruht. Bei der Betrachtung von Orientierungseffekten konnte gezeigt werden, dass PT für alle Orientierungen mit einem Feld, welches aus einer mittleren Orientierung bestimmt wurde, effizient ist. Untersuchungen spektroskopischer Methoden, um Informationen über die Effizienz von Kontrollprozessen zu liefern, zeigten, dass Pump-Probe Ionisationsspektroskopie im Femtosekundenbereich (fs) dazu sehr gut dazu geeignet ist. Im Speziellen konnte mit zeitabh. Photoelektronenspektren die Populationen in den elektronischen Zuständen nach Anlegen des jeweiligen Feldes „beobachtet“ werden. Im letzten Kapitel wurde eine andere Methode der Kontrolle von Molekülen in Anlehnung an einen STIRAP ("Stimulated Raman Adiabatic Passage“) Prozess am Beispiel molekularen Iods vorgestellt. Dabei wurde die Möglichkeit, diese Technik auf die fs-Zeitskala auszudehnen,in Theorie und Experiment untersucht, wobei die Messungen von Dr. Torsten Siebert (Universität Würzburg, Arbeitskreis Prof. Kiefer) durchgeführt worden sind. Nicht-resonante Anregung, mit einer Abfolge der Pulse, in der der Stokes-Puls mit der höheren Intensität im Vergleich zum Pump-Puls in einer f-STIRAP-artigen Anordnung dem Pump-Puls vorausgeht, führte zu einem effizienten PT in einen schwingungsangeregten Zustand im el. Grundzustand. Dies konnte durch einen Vergleich des numerisch exakten Falls mit einer adiabatischen Behandlung theoretisch untermauert werden. Die zu Grunde liegenden Prozesse sind näherungsweise durch adiabatisches Verhalten charakterisiert. Dazu wird gerade ein neues Modell entwickelt, welches den Prozess mit einem um einen dressed-state rotierenden Vektor im Hilbertraum erklärt. Zusammenfassend wurde in dieser Arbeit ein Algorithmus zur Kontrolle von Moleküldynamik entwickelt, der auf der instantanen Antwort eines Systems bei Wechselwirkung mit einem elektrischen Feld beruht. Die daraus berechneten Kontrollfelder sind sehr effizient bezüglich einer selektiven Kontrolle von Molekülen und versprechen noch viele zukünftige Anwendungsmöglichkeiten. KW - Laserstrahlung KW - Molekulardynamik KW - Mehrphotonenprozess KW - Quantenmechanik KW - Quantendynamik KW - Kontrolltheorie KW - STIRAP KW - geformte Laserfelder KW - Multi-Photonen Prozesse KW - quantum dynamics KW - control theory KW - STIRAP KW - pulse shaping KW - multi-photon processes Y1 - 2005 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-13388 ER -