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 - 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 -