@phdthesis{Vogt2020, author = {Vogt, Matthias Guido}, title = {Elektronische Eigenschaften von Wabengittern mit starker Spin-Bahn-Kopplung}, doi = {10.25972/OPUS-20750}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-207506}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2020}, abstract = {Im Rahmen dieser Arbeit wurden die elektronischen Eigenschaften von Graphen auf Metalloberfl{\"a}chen mittels Rastertunnelmikroskopie und Quasiteilcheninterferenz (englisch quasiparticle interference, QPI)-Messungen untersucht. Durch das Verwenden schwerer Substrate sollte die Spin-Bahn-Wechselwirkung des Graphen verst{\"a}rkt werden und damit eine Bandl{\"u}cke am K-Punkt der Bandstruktur mittels QPI beobachtet werden. Um das Messen von QPI auf Graphen zu testen, wurde auf der Oberfl{\"a}che eines SiC(0001)-Kristalls durch Erhitzen Graphen erzeugt und mit dem Rastertunnelmikroskop untersucht. Dieses System wurde schon ausf{\"u}hrlich in der Literatur beschrieben und bereits bekannte QPI-Messungen von Streuringen, die auf den Dirac-Kegeln des Graphen am K-Punkt basieren, konnte ich auf gr/SiC(0001) in guter Qualit{\"a}t erfolgreich reproduzieren. Anschließend wurde Graphen nach einem wohlbekannten Verfahren durch Aufbringen von Ethylen auf ein erhitztes Ir(111)-Substrat erzeugt. Dieses gr/Ir(111)-System diente auch als Grundlage f{\"u}r Interkalationsversuche von Bismut (gr/Bi/Ir(111)) und Gadolinium (gr/Gd/Ir(111)) zwischen das Graphen und das Substrat. Auf gr/Bi/Ir(111) wurde ein schon aus der Literatur bekanntes Netzwerk aus Versetzungslinien beobachtet, dem zus{\"a}tzlich eine Temperaturabh{\"a}ngigkeit nachgewiesen werden konnte. Beim Versuch, Gadolinium zu interkalieren, wurden zwei verschieden Oberfl{\"a}chenstrukturen beobachtet, die auf eine unterschiedlich Anordnung bzw. Menge des interkalierten Gadoliniums zur{\"u}ckzuf{\"u}hren sein k{\"o}nnten. Auf keinem dieser drei Systeme konnten allerdings Streuringe mittels QPI beobachtet werden. Als Vorbereitung der Interkalation von Gadolinium wurden dessen Wachstum und magnetische Eigenschaften auf einem W(110)-Kristall untersucht. Dabei konnte eine aus der Literatur bekannte temperaturabh{\"a}ngige Austauschaufspaltung reproduziert werden. Dar{\"u}ber hinaus konnten sechs verschieden magnetische Dom{\"a}nen beobachtet werden. Zus{\"a}tzlich sind auf der Oberfl{\"a}che magnetische Streifen auszumachen, die m{\"o}glicherweise auf einer Spinspirale basieren. Als Grundlage f{\"u}r die m{\"o}gliche zuk{\"u}nftige Erzeugung Graphen-artiger Molek{\"u}lgitter wurde das Wachstum von H-TBTQ und Me-TBTQ auf Ag(111) untersucht. Die Molek{\"u}le richten sich dabei nach der Oberfl{\"a}chenstruktur des Silber aus und bilden l{\"a}ngliche Inseln, deren Kanten in drei Vorzugsrichtungen verlaufen. Auf H-TBTQ wurde zudem eine zweite, Windm{\"u}hlen-artige Ausrichtung der Molek{\"u}le auf der Oberfl{\"a}che beobachtet. Auf den mit den Molek{\"u}len bedeckten Stellen der Oberfl{\"a}che wurde eine Verschiebung des Ag-Oberfl{\"a}chenzustands beobachtet, die mit einem Ladungstransfer vom Ag(111)-Substrat auf die TBTQ-Molek{\"u}le zu erkl{\"a}ren sein k{\"o}nnte.}, subject = {Spin-Bahn-Wechselwirkung}, language = {de} } @phdthesis{Brueckner2017, author = {Br{\"u}ckner, Charlotte}, title = {The Electronic Structure and Optoelectronic Processes at the Interfaces in Organic Solar Cells Composed of Small Organic Molecules - A Computational Analysis of Molecular, Intermolecular, and Aggregate Aspects}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-141652}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2017}, abstract = {Describing the light-to-energy conversion in OSCs requires a multiscale understanding of the involved optoelectronic processes, i.e., an understanding from the molecular, intermolecular, and aggregate perspective. This thesis presents such a multiscale description to provide insight into the processes in the vicinity of the organic::organic interface, which are crucial for the overall performance of OSCs. Light absorption, exciton diffusion, photoinduced charge transfer at the donor-acceptor interface, and charge separation are included. In order to establish structure-property relationships, a variety of different molecular p-type semiconductors are combined at the organic donor-acceptor heterojunction with fullerene C60, one of the most common acceptors in OSCs. Starting with a comprehensive analysis of the accuracy of diverse ab initio, DFT, and semiempiric methods for the properties of the individual molecules, the intermolecular, and aggregate/device stage are subsequently addressed. At all stages, both methodological concepts and physical aspects in OSCs are discussed to extend the microscopic understanding of the charge generation processes.}, subject = {Benchmark}, language = {en} } @phdthesis{Heilos2019, author = {Heilos, Anna}, title = {Mechanistic Insights into the Inhibition of Cathepsin B and Rhodesain with Low-Molecular Inhibitors}, doi = {10.25972/OPUS-17822}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-178228}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2019}, abstract = {Cysteine proteases play a crucial role in medical chemistry concerning various fields reaching from more common ailments like cancer and hepatitis to less noted tropical diseases, namely the so-called African Sleeping Sickness (Human Arfican Trypanosomiasis). Detailed knowledge about the catalytic function of these systems is highly desirable for drug research in the respective areas. In this work, the inhibition mechanisms of the two cysteine proteases cathepsin B and rhodesain with respectively one low-molecular inhibitor class were investigated in detail, using computational methods. In order to sufficiently describe macromolecular systems, molecular mechanics based methods (MM) and quantum mechanical based method (QM), as well as hybrid methods (QM/MM) combining those two approaches, were applied. For Cathespin B, carbamate-based molecules were investigated as potential inhibitors for the cysteine protease. The results indicate, that water-bridged proton-transfer reactions play a crucial role for the inhibition. The energetically most favoured pathway (according to the calculations) includes an elimination reaction following an E1cB mechanism with a subsequent carbamylation of the active site amino acid cysteine. Nitroalkene derivatives were investigated as inhibitors for rhodesain. The investigation of structurally similar inhibitors showed, that even small steric differences can crucially influence the inhibition potential of the components. Furthermore, the impact of a fluorination of the nitroalkene inhibitors on the inhibition mechanism was investigated. According to experimental data measured from the working group of professor Schirmeister in Mainz, fluorinated nitroalkenes show - in contrast to the unfluorinated compounds - a time dependent inhibition efficiency. The calculations of the systems indicate, that the fluorination impacts the non-covalent interactions of the inhibitors with the enzymatic environment of the enzyme which results in a different inhibition behaviour.}, subject = {Cysteinproteasen}, language = {en} } @phdthesis{Bellinger2016, author = {Bellinger, Daniel}, title = {Implementation of new reaction pathway determining methods and study of solvent effects on the excited state nature of perylene based dyes}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-144435}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2016}, abstract = {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.}, subject = {Globale Optimierung}, language = {en} } @phdthesis{Le2020, author = {Le, Thien Anh}, title = {Theoretical investigations of proton transfer and interactions or reactions of covalent and non-covalent inhibitors in different proteins}, doi = {10.25972/OPUS-17051}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-170511}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2020}, abstract = {Nowadays, computational-aided investigations become an essential part in the chemical, biochemical or pharmaceutical research. With increasing computing power, the calculation of larger biological systems becomes feasible. In this work molecular mechanical (MM) and quantum mechanical approaches (QM) and the combination of both (QM/MM) have been applied to study several questions which arose from different working groups. Thus, this work comprises eight different subjects which deals with chemical reactions or proton transfer in enzymes, conformational changes of ligands or proteins and verification of experimental data. This work firstly deals with reaction mechanisms of aromatic inhibitors of cysteine proteases which can be found in many organisms. These enzymes are responsible for various cancer or diseases as for example Human African Trypanosomiasis (HAT) or the Chagas disease. Aromatic SNAr-type electrophiles might offer a new possibility to covalently modify these proteases. Quantum mechanical calculations have been performed to gain insights into the energetics and possible mechanisms. The next chapter also deals with Trypanosomiasis but the focus was set on a different enzyme. The particularity of Trypanosomiasis is the thiol metabolism which can also be modified by covalent inhibitors. In this context, the wild type and point mutations of the enzyme tryparedoxin have been investigated via molecular dynamic (MD) simulations to examine the influence of specific amino acids in regard to the inhibitor. Experimental data showed that a dimerization of the enzyme occurs if the inhibitor is present. Simulations revealed that the stability of the dimer decreases in absence of the inhibitor and thus confirms these experiments. Further investigations concerning cysteine proteases such as cruzain and rhodesain have been conducted with respect to experimental kinetic data of covalent vinylsulfone inhibitors. Several approaches such as QM or QM/MM calculations and docking, MD or MMPBSA/MMGBSA simulations have been applied to reproduce these data. The utilization of force field approaches resulted in a qualitatively accurate prediction. The kinase AKT is involved in a range of diseases and plays an important role in the formation of cancer. Novel covalent-allosteric inhibitors have been developed and crystallized in complex with AKT. It was shown that depending on the inhibitor a different cysteine residue is modified. To investigate these differences in covalent modification computational simulations have been applied. Enoyl-(acyl carrier) (ENR) proteins are essential in the last step of the fatty acid biosynthesis II (FAS) and represent a good target for inhibition. The diphenylether inhibitor SKTS1 which was originally designed to target the ENR's of Staphylococcus aureus was also crystallized in InhA, the ENR of Mycobacterium tuberculosis (TB). Crystal structures indicate a change of the inhibitor's tautomeric form. This subject was investigated via MD simulations. Results of these simulations confirmed the tautomerization of the inhibitor. This work also deals with the development of a covalent inhibitor originating from a non-covalent ligand. The target FadA5 is an essential enzyme for the degradation of steroids in TB and is responsible for chronic tuberculosis. This enzyme was crystallized in complex with a non-covalent ligand which served as starting point for this study. Computations on QM or QM/MM level and docking and MD simulations have been applied to evaluate potential candidates. The next chapter focuses on the modification of the product spectrum of Bacillus megaterium levansucrase, a polymerase which catalyzes the biosynthesis of fructans. The covalent modification of the wild type or mutants of the enzyme lead to an accumulation of oligosaccharides but also to polymers with higher polymerization degree. To understand these changes in product spectra MD simulations have been performed. Finally, the proton transfer in catalytic cysteine histidine dyads was investigated. The focus was set on the influence of the relaxation of the protein environment to the reaction. Calculations of the enzymes FadA5 and rhodesain revealed that the preferred protonation state of the dyade depends on the protein environment and has an impact on the reaction barrier. Furthermore, the adaptation of the environment to a fixed protonation state was analyzed via MD simulations.}, subject = {Computational chemistry}, language = {en} } @phdthesis{Welz2020, author = {Welz, Eileen}, title = {Theoretical Investigations on Inorganic, Boron-containing Biradicals with a unique Structure}, doi = {10.25972/OPUS-20259}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-202598}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2020}, abstract = {In this work, biradical boron containing systems with various structures are investigated to reveal the dependency of the biradical character on the ligated carbene (NHC, CAAC) and the related steric demands of the substituents.}, subject = {Biradical}, language = {en} } @phdthesis{Deissenberger2020, author = {Deißenberger, Andrea}, title = {Dibortetrahalogenide f{\"u}r die Darstellung neuer borhaltiger Verbindungen in niedrigen Oxidationsstufen}, doi = {10.25972/OPUS-18775}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-187758}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2020}, abstract = {Im Rahmen dieser Arbeit konnten nasschemische Synthesen f{\"u}r Dibortetrafluorid und chlorid ausgehend von Dibortetrabromid entwickelt werden, die durch einfachen Halogenaustausch mit SbF3 bzw. GaCl3 realisiert wurden. In Verbindung mit Arbeiten von Dr. Jonas M{\"u}ssig zur Synthese von B2I4 gelang die Darstellung aller vier Dibortetrahalogenide mittels einfacher Schlenktechnik basierend auf der Synthese von B2Br4 durch N{\"o}th und Pommerening im Jahr 1981. Dibortetrachlorid konnte mit Phosphanen (PMe3, PCy3 und PPh3) und Singulett-Carbenen (IDipp und MeCAAC) zu den klassischen Bisaddukten 44-46 bzw. 54 und 55 umgesetzt werden. Die Addition eines Isonitrils (CNtBu) an B2Cl4 f{\"u}hrte zun{\"a}chst zur Ausbildung des Bisadduktes 53, allerdings konnte in L{\"o}sung eine Umlagerung beobachtet werden, deren Verlauf 11B-NMR-spektroskopisch verfolgt wurde, jedoch nicht final aufgekl{\"a}rt werden konnte. Durch die Zugabe eines Unterschusses der Lewis-Basen IDipp bzw. PCy3 sollten zun{\"a}chst Monoaddukte von B2Cl4 dargestellt werden, deren Umsetzung mit einer weiteren Lewis-Base die Synthese asymmetrischer Lewis-Basen-Addukte von B2Cl4 erm{\"o}glichen sollte. Die sp2-sp3-Diborane 56 und 57 konnten bei tiefen Temperaturen 11B-NMR-spektroskopisch nachgewiesen werden, allerdings f{\"u}hrte eine Abfangreaktion mit diversen Lewis-Basen nicht zu den gew{\"u}nschten asymmetrischen Addukten. Bei Raumtemperatur konnte eine Folgereaktion von 56 zur Chlorid-verbr{\"u}ckten kationischen Spezies 58 mit einem Tetrachloroborat-Anion beobachtet werden. Im Fall von Dibortetrafluorid konnten keine Lewis-Basen-Addukte (LB = PMe3 und MeCAAC) isoliert werden. Die Reaktivit{\"a}t von B2Cl4 gegen{\"u}ber unges{\"a}ttigten Substraten wurde anhand mehrerer literaturbekannter Beispiele (Acetylen, 2-Butin, 3-Hexin, Diphenylacetylen und Bis(trimethylsilyl)acetylen) nachvollzogen und um die terminalen Alkine Propin und 1 Hexin erweitert. Eine selektive Addition von B2Br4 an Dreifachbindungen gelang nicht. Die so erhaltenen Diborylalkene sollten zur Darstellung von 1,2-Diboreten genutzt werden, wobei zun{\"a}chst {\"u}ber eine von Siebert et al. entwickelte Route die Bis(N,N-dialkylaminochlorboryl)alkene 67g, h, j und k dargestellt wurden. Ein nachfolgender Ringschluss unter reduktiven Bedingungen verlief nur f{\"u}r die Diisopropyl¬amino-substituierten Diborylalkene 67g und j selektiv und lieferte das 1,2-Dihydro-1,2-diboret 71g und das umgelagerte 1,3-Dihydro-1,3-diboret 68j. Der Austausch der Aminosubstituenten gegen Halogenide, der f{\"u}r eine weitere Reduktion zur B-B-Doppelbindung n{\"o}tig w{\"a}re, gelang nicht. Die Umsetzung der Diborylalkene 61 (R = Me), 62 (R = Et) und 65 (R = Ph) mit Singulett-Carbenen (LB = IMe, IiPr, IDipp und MeCAAC) f{\"u}hrte zu den chloridverbr{\"u}ckten Monoaddukten 74-76 und 79-81. Alle Verbindungen dieses Typs zeigten in NMR-spektroskopischen Untersuchungen ein sp2- und ein sp3-koordiniertes Borzentrum, welche f{\"u}r die CAAC-stabilisierten Verbindungen auch r{\"o}ntgenkristallografisch nachgewiesen werden konnten. Theoretische Untersuchungen best{\"a}tigten die Relevanz des verbr{\"u}ckenden Chloratoms zur Stabilisierung dieser Verbindungen. F{\"u}r die Stammverbindung der Diborylalkene (59 (R = H)) konnte bei der Umsetzung mit MeCAAC eine unl{\"o}sliche Verbindung erhalten werden, deren Struktur als Bisaddukt 82 mittels NMR-spektroskopischen Untersuchungen im Festk{\"o}rper und durch Verbrennungsanalyse best{\"a}tigt werden konnte. Die Reduktion der CAAC-stabilisierten Diborylalkene 79 und 80 in Gegenwart von MeCAAC f{\"u}hrte zu den captodativ-stabilisierten Diborylradikalen 83 und 84, deren Strukturanalyse eine orthogonale Anordnung der C2-Br{\"u}cke zur B(CAAC)-Einheit offenlegt. Ausf{\"u}hrliche EPR-spektroskopische Untersuchungen bei variabler Temperatur und theoretische Berechnungen best{\"a}tigen eine schwache Wechselwirkung der beiden Radikalzentren und einen offenschaligen Singulett-Grundzustand mit einem energetisch tiefliegenden Triplett-Zustand (ΔES T = 0.017 kcal mol-1). Der experimentell bestimmte Spin-Spin-Abstand und die Analyse der einfach besetzten Molek{\"u}lorbitale (SOMO) best{\"a}tigen eine Delokalisierung der Spindichte {\"u}ber die NCAAC-CCAAC-B-Einheit. Der Austausch der verbr{\"u}ckenden Einheit und die somit einhergehende Verringerung der Sterik f{\"u}hrt zu einer Planarisierung des Molek{\"u}ls im Festk{\"o}rper (87). Theoretische Untersuchungen und die Auswertung der strukturellen Parameter ergeben eine Delokalisierung der Elektronendichte {\"u}ber das gesamte planare System. EPR- und NMR-spektroskopische Untersuchungen ergaben dennoch Hinweise auf das Vorliegen einer paramagnetischen Verbindung. Untersuchungen zum Reduktionsverhalten von zweifach CAAC-stabilisiertem 1,4-Bis-(dibromboryl)benzol (97) ergaben die vollst{\"a}ndige Enthalogenierung der Borzentren. Im Zuge dessen entstand ein hochreaktives, lineares Borylen, welches eine CH-Aktivierung mit dem Isopropylsubstituenten des CAAC-Liganden eingeht (98). Zur Stabilisierung des Borylens wurde die Reduktion in Gegenwart weiterer Lewis-Basen (Pyridin (Pyr), IiPr, IMeMe, PMe3, CNtBu und CO) durchgef{\"u}hrt, die in der Ausbildung der Diborylene 99-104 resultierten. Die Darstellung einer para-Phenylen-verbr{\"u}ckten Donor-Akzeptor-Verbindung (D: Borylen, A: BMes2) gelang nicht.}, subject = {Dibortetrahalogenide}, language = {de} } @phdthesis{Walter2015, author = {Walter, Christof}, title = {Excitonic States and Optoelectronic Properties of Organic Semiconductors - A Quantum-Chemical Study Focusing on Merocyanines and Perylene-Based Dyes Including the Influence of the Environment}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-123494}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2015}, abstract = {The scope of computational chemistry can be broadened by developing new methods and more efficient algorithms. However, the evaluation of the applicability of the methods for the different fields of chemistry is equally important. In this thesis systems with an unusual and complex electronic structure, such as excitonic states in organic semiconductors, a boron-containing bipolaron and the excited states of pyracene were studied and the applicability of the toolkit of computational chemistry was investigated. Concerning the organic semiconductors the focus was laid on organic solar cells, which are one of the most promising technologies with regard to satisfying the world's need for cheap and environmentally sustainable energy. This is due to the low production and material costs and the possibility of using flexible and transparent devices. However, their efficiency does still not live up to the expectations. Especially the exciton diffusion lengths seem to be significantly too short. In order to arrive at improved modules, a fundamental understanding of the elementary processes occurring in the cell on the molecular and supramolecular level is needed. Computational chemistry can provide insight by separating the different effects and providing models for predictions and prescreenings. In this thesis, the focus was laid on the description of excitonic states in merocyanines and perylene-based dyes taking the influence of the environment into account. At first, the photochemical isomerization between two configurations of 6-nitro BIPS observed experimentally was studied by first benchmarking several functionals against SCS-ADC(2) in the gas phase and subsequently calculating the excited-state potential energy surface. The geometries obtained from a relaxed scan in the ground state as well as from a scan in the excited state were used. The environment was included using different polarizable continuum models. It was shown that the choice of the model and especially the question of the state specificity of the approach is of vital importance. Using the results of the calculations, a two-dimensional potential energy surface could be constructed that could be used to explain the experimental findings. Furthermore, the importance of the excited-state isomerization as a potential deactivation channel in the exciton transport was pointed out. Then the assessment of the suitability of different merocyanines for optoelectronic applications with quantum-chemical methods was discussed. At first, the effect of the environment on the geometry, especially on the bond length alternation pattern, was investigated. It was shown that the environment changes the character of the ground-state wave function of several merocyanines qualitatively, which means that the results of gas-phase calculations are meaningless - at least when a comparison with solution or device data is desired. It was demonstrated that using a polarizable continuum model with an effective epsilon, a qualitative agreement between the calculated geometry and the geometry in the crystal structure can be obtained. Therefore, by comparing the bond length alternation in solution and in the crystal, a rough estimate of the effect of the crystal environment can be made. It was further shown that the connection between the HOMO energy and the open-circuit voltage is not as simple as it is often implied in the literature. It was discussed that it is not clear whether the HOMO of a single molecule or a \$\pi\$-stack containing several monomers should be used and if the environmental charges of the bulk phase or the interface should be included. Investigating the dependence of the HOMO energy on the stack size yielded no definitive trend. Furthermore, it was discussed that the effect due the optimization of the modules (solvent, bulk heterojunction) during the production masks any potential correlation between the HOMO energy and measured open-circuit values. Therefore, a trend can only be expected for unoptimized bilayer cells. It was concluded that ultimately, the importance of the HOMO energy should not be overestimated. The correlation between the exciton reorganization energy and the so-called cyanine limit, which is predicted by a simple two-state model, was also discussed. By referring to the results of VB calculations, it was discussed that the correlation indeed exists and is non-negligible, although the effect is not as strong as one might have expected. In this context, a potential application of a VB/MM approach was covered briefly. The importance of the molecular reorganization energy and the device morphology was also discussed. It was concluded that the optimization of merocyanines for organic optoelectronic devices is inherently a multiparameter problem and one cannot expect to find one particular parameter, which solely controls the efficiency. The perylene-based dyes were studied with a focus on the description of a potential trapping mechanism involving an intermolecular motion in a dimer. The aim was to find methods which can be applied to larger model systems than a dimer and take the effect of the environment into account. As a test coordinate the longitudinal shift of two monomers against each other was used. At first, it was demonstrated how the character of an excited state in a dimer can be defined and how it can be extracted from a standard quantum-chemical calculation. Then several functionals were benchmarked and their applicability or failure was rationalized using the character analysis. Two recipes could be proposed, which were applied to a constraint optimization (only intermolecular degrees of freedom) in the excited states of the PBI dimer and to the description of the potential energy surfaces of ground and excited states along a longitudinal displacement in the perylene tetramer, respectively. It was further demonstrated that the semi-empirical OMx methods fail to give an accurate description of the excited-state potential energy surfaces as well as the ground-state surface along the test coordinate. This failure could be attributed to an underestimation of overlap-dependent terms. Consequently, it could be shown that the methods are applicable to large intermolecular distances, where the overlap is negligible. The results of DFT calculations with differently composed basis sets suggested that adding an additional single p-function for each atom should significantly improve the performance. QM/MM methods are ideally suited to take the effect of the environment on a a dimer model system into account. However, it was shown that standard force fields also give an incorrect description of the interaction between the monomers along the intermolecular coordinate. This failure was attributed to the isotropic atom-atom interaction in the repulsion term of the Lennard-Jones potential. This was corroborated using two simple proof-of-principle anisotropy models. Therefore, a novel force field called OPLS-AA_O was presented that is based on OPLS-AA, but uses an anisotropic model for the repulsion. The model involves the overlap integral between the molecular densities, which are modeled as a sum of atom-centered p-type Gaussian functions. It was shown that using this force field an excellent agreement with the DFT results can be obtained when the correct parameters are used. These parameters, however, are not very generalizable, which was attributed to the simplicity of the model in its current state (using the same exponential parameter for all atoms). As a short excursion, the applicability of an MO-based overlap model was discussed. It was demonstrated that the repulsion term based on the density overlap can be used to correct the failure of the OMx methods for the ground states. This is in accord with the assumption that an underestimation of the overlap terms is responsible for the failure. It was shown that OPLS-AA_O also gives an excellent description of the longitudinal shift in a PBI tetramer. Using the tetramer as a test system and applying the recipe obtained in the TDDFT benchmark for the QM-part and OPLS-AA_O for the MM-part in conjunction with an electrostatic embedding scheme, a QM/MM description of the excited states of the PBI dimer including the effect of the environment could be obtained. In the last chapter the theoretical description of the Bis(borolyl)thiophene dianion and the excited states of pyracene were discussed. The electronic structure of the Bis(borolyl)thiophene dianion - a negative bipolaron - was elucidated using DFT and CASPT2 methods. Furthermore, an estimation of the extent of triplet admixture to the ground state due to spin-orbit coupling was given. In the second project the S1 and S2 states of pyracene were computed using SCS-CC2 and SCS-ADC(2) and an estimation for the balance between aromaticity and ring strain was given. This also involved computing the vibrational frequencies in the excited states. In both studies the results of the computations were able to rationalize and complete experimental results.}, subject = {Exziton}, language = {en} } @phdthesis{Schmidt2015, author = {Schmidt, Thomas Christian}, title = {Theoretical Investigations on the Interactions of Small Compounds with their Molecular Environments}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-127860}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2015}, abstract = {Im ersten Teil dieser Arbeit wird eine Kombination theoretischer Methoden f{\"u}r die strukturbasierte Entwicklung neuer Wirkstoffe pr{\"a}sentiert. Ausgehend von der Kristallstruktur eines kovalenten Komplexes einer Modellverbindung mit dem Zielprotein wurde mit Hilfe von quantenmechanischen und QM/MM Rechnungen die genaue Geometrie des vorausgehenden nicht-kovalenten Komplexes betimmt. Letztere ist der bestimmende Faktor f{\"u}r die Reaktivit{\"a}t des Inhibitors gegen{\"u}ber der katalytisch aktiven Aminos{\"a}ure und damit f{\"u}r die Ausbildung einer kovalenten Bindung. Aus diesem Grund wurde diese Geometrie auch f{\"u}r die Optimierung der Substitutionsmusters des Ihnibitors verwendet, um dessen Affinit{\"a}t zum Zielenzyme zu verbessern ohne dass dieser seine F{\"a}higkeit kovalent an das aktive Zentrum zu binden verliert. Die Optimierung des Substitutionsmuster wurde doch Methode des Molekularen Dockings unterst{\"u}tzt, das diese optimal dazu geeignet sind, Bindungsaffinit{\"a}ten vorherzusagen, die durch eine Modifikation der chemischen Struktur entstehen. Eine Auswahl der besten Strukturen wurde anschließend verwendet, um zu {\"u}berpr{\"u}fen, ob die ver{\"a}nderten Molek{\"u}le noch gen{\"u}gen Reaktivit{\"a}t gegen{\"u}ber dem Zielprotein aufweisen. Molek{\"u}ldynamik Simulationen der neuen Verbindungen haben jedoch gezeigt, dass die ver{\"a}nderten Verbindungen nur so and das Protein binden, dass die Bilung eine kovalenten Bindung zum Enzym nicht mehr m{\"o}glich ist. Daher wurden in einem weiteren Schritt die Modellverbindungen weiter modifiziert. Neben {\"A}nderungen im Substitutionsmuster wurde auch die chemische Struktur im Kern ver{\"a}ndert. Die Bindungsaffinit{\"a}ten wurde wieder mittels Docking {\"u}berpr{\"u}ft. F{\"u}r die besten Bindungsposen wurden wieder Simulationen zur Molek{\"u}ldynamik durchgef{\"u}hrt, wobei diesmal die Ausbildung einer kovalenten Bindung zum Enzyme m{\"o}glich erscheint. In einer abschließenden Serie von QM/MM Rechnungen unter Ber{\"u}cksichtigung verschiedener Protonierungszust{\"a}nde des Inhibitors und des Proteins konnten Reaktionspfade und zugeh{\"o}rige Reaktionsenergien bestimmt werden. Die Ergebnisse lassen darauf schließen, dass eines der neu entwickelten Molek{\"u}le sowohl eine stark verbesserte Bindungsaffinit{\"a}t wie auch die M{\"o}glichkeit der kovalenten Bindung an Enzyme aufweist. Der zweite Teil der Arbeit konzentriert sich auf die Umgebungseinfl{\"u}sse auf die Elektronenverteilung eines Inhibitormodells. Als Grundlage dient ein vinylsulfon-basiertes Moek{\"u}l, f{\"u}r das eine experimentell bestimmte Kristallstruktur sowie ein theoretisch berechneter Protein Komplex verf{\"u}gbar sind. Ein Referendatensatz f{\"u}r diese Systeme wurde erstellt, indem der Konformationsraum des Inhibitors nach m{\"o}glichen Minimumsstrukturen abgesucht wurde, welche sp{\"a}ter mit den Geometrien des Molek{\"u}ls im Kristall und im Protein verglichen werden konnten. The Geometrie in der Kristallumgebung konnte direkt aus den experimentellen Daten {\"u}bernommen werden. Rechnungen zum nicht-kovalenten Protein Komplex hingegen haben gezeigt, dass f{\"u}r das Modellsystem mehrere Geometrien des Inhibiors sowie zwei Protonierungszust{\"a}nde f{\"u}r die katalytisch aktiven Aminos{\"a}uren m{\"o}glich sind. F{\"u}r die Analyse wurden daher alle m{\"o}glichen Proteinkomplexe mit der Kristallstruktur verglichen. Ebenso wurden Vergleiche mit der Geometrie des isolierten Molek{\"u}ls im Vakuum sowie der Geometrie in w{\"a}ssriger L{\"o}sung angestellt. F{\"u}r die Geometrie des Molek{\"u}ls an sich ergab sich eine gute {\"U}bereinstimmung f{\"u}r alle Modellsysteme, f{\"u}r die Wechselwirkungen mit der Umgebung jedoch nicht. Die Ausbildung von Dimeren in der Kristallumgebung hat einen stark stablisierenden Effekt und ist einer der Gr{\"u}nde, warum dieser Kristall so gut wie keine Fehlordungen aufweist. In den Proteinkomplexen hingegen ergibt sich eine Abstoßung zwischen dem Inhibitor und einer der katalytisch aktiven Aminos{\"a}uren. Als Ursache f{\"u}r diese Abstoßung konnte die Einf{\"u}hrung der Methylaminfunktion ausgemacht werden. Vermutlicherweise f{\"u}hrt diese strukturelle {\"A}nderung auch dazu, dass der Modellinhibitor nicht in der Lage ist, so wie die Leitstruktur K11777 an das aktive Zentrum des Enzyms zu binden.}, subject = {Theoretische Chemie}, language = {en} } @phdthesis{Becker2015, author = {Becker, Johannes}, title = {Development and implementation of new simulation possibilities in the CAST program package}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-132032}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2015}, abstract = {The aim of the present work is the development and implementation of new simulation possibilities for the CAST program package. Development included, among other things, the partial parallelization of the already existing force fields, extension of the treatment of electrostatic interactions and implementation of molecular dynamics and free energy algorithms. The most time consuming part of force field calculations is the evaluation of the nonbonded interactions. The calculation of these interactions has been parallelized and it could be shown to yield a significant speed up for multi-core calculations compared to the serial execution on only one CPU. For both, simple energy/gradient as well as molecular dynamics simulations the computational time could be significantly reduced. To further increase the performance of calculations employing a cutoff radius, a linkedcell algorithm was implemented which is able to build up the non-bonded interaction list up to 7 times faster than the original algorithm. To provide access to dynamic properties based on the natural time evolution of a system, a molecular dynamics code has been implemented. The MD implementation features two integration schemes for the equations of motion which are able to generate stable trajectories. The basic MD algorithm as described in Section 1.2 leads to the sampling in the microcanonical (NVE) ensemble. The practical use of NVE simulations is limited though because it does not correspond to any experimentally realistic situation. More realistic simulation conditions are found in the isothermal (NVT) and isothermalisobaric (NPT) ensembles. To generate those ensembles, temperature and pressure control has been implemented. The temperature can be controlled in two ways: by direct velocity scaling and by a Nose-Hoover thermostat which produces a real canonical ensemble. The pressure coupling is realized by implementation of a Berendsen barostat. The pressure coupling can be used for isotropic or anisotropic box dimensions with the restriction that the angles of the box need to be 90� . A crucial simulation parameter in MD simulations is the length of the timestep. The timestep is usually in the rang of 1fs. Increasing the timestep beyond 1fs can lead to unstable trajectories since the fastest motion in the system, usually the H-X stretch vibration can not be sampled anymore. A way to allow for bigger timesteps is the use of a constraint algorithm which constrains the H-X bonds to the equilibrium distance. For this the RATTLE algorithm has been implemented in the CAST program. The velocity Verlet algorithm in combination with the RATTLE algorithm has been shown to yield stable trajectories for an arbitrary length of simulation time. In a first application the MD implementation is used in conjunction with the MOPAC interface for the investigation of PBI sidechains and their rigidity. The theoretical investigations show a nice agreement with experimentally obtained results. Based on the MD techniques two algorithms for the determination of free energy differences have been implemented. The umbrella sampling algorithm can be used to determine the free energy change along a reaction coordinate based on distances or dihedral angles. The implementation was tested on the stretching of a deca-L-alanine and the rotation barrier of butane in vacuum. The results are in nearly perfect agreement with literature values. For the FEP implementation calculations were performed for a zero-sum transformation of ethane in explicit solvent, the charging of a sodium ion in explicit solvent and the transformations of a tripeptide in explicit solvent. All results are in agreement with benchmark calculations of the NAMD program as well as literature values. The FEP formalism was then applied to determine the relative binding free energies between two inhibitors in an inhibitor-protein complex. Next to force fields, ab-initio methods can be used for simulations and global optimizations. Since the performance of such methods is usually significantly poorer than force field applications, the use for global optimizations is limited. Nevertheless significant progress has been made by porting these codes to GPUs. In order to make use of these developments a MPI interface has been implemented into CAST for communication with the DFT code TeraChem. The CAST/TeraChem combination has been tested on the \$H_2 O_{10}\$ cluster as well as the polypeptide met-Enkephalin. The pure ab-initio calculations showed a superior behavior compared to the standard procedure where the force field results are usually refined using quantum chemical methods.}, subject = {Molekulardynamik}, language = {en} } @phdthesis{Ansorg2015, author = {Ansorg, Kay}, title = {Development of Accurate Physically Grounded Force Fields for Intermolecular Cation-\$\pi\$ Interactions based on SAPT Energy Decomposition Analysis and Computational Investigation of Covalent Irreversible Vinyl Sulfone-based Protease Inhibitors}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-131084}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2015}, abstract = {Part 1 of this work describes the development of accurate physically grounded force fields for intermolecular Cation-π interactions based on SAPT energy decomposition analysis. The presented results demonstrate the benefits of the used DFT-SAPT method to describe non-bonding interactions. First of all, this method is able to reproduce the high level CCSD(T) energy values but using much less computational time. Second it provides the possibility to separate the total intermolecular interaction energy into several physically meaningful contributions. The relative contributions of the dimers investigated can be seen in Fig. 6.16. In Tab. 6.3 the percentage contribution of the attractive energy parts to the stabilization energy is shown. The polarization energy is important for the NH+...C6H6 interaction, whereas it becomes less crucial considering other dimers. The dispersion energy contribution is large in the case of the C6H6...H2O dimers, whereas it is relatively less important for the NH+...C6H6 interaction. The electrostatic energy contributes a large amount of stabilizing energy in all considered dimer interactions. ...}, subject = {Kraftfeld}, language = {en} } @phdthesis{Stehr2015, author = {Stehr, Vera}, title = {Prediction of charge and energy transport in organic crystals with quantum chemical protocols employing the hopping model}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-114940}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2015}, abstract = {As organic semiconductors gain more importance for application, research into their properties has become necessary. This work investigated the exciton and charge transport properties of organic semiconducting crystals. Based on a hopping approach, protocols have been developed for the calculation of Charge mobilities and singlet exciton diffusion coefficients. The protocols do not require any input from experimental data except for the x-ray crystal structure, since all needed quantities can be taken from high-level quantum chemical calculations. Hence, they allow to predict the transport properties of yet unknown compounds for given packings, which is important for a rational design of new materials. Different thermally activated hopping models based on time-dependent perturbation theory were studied for the charge and exciton transport; i. e. the spectral overlap approach, the Marcus theory, and the Levich-Jortner theory. Their derivations were presented coherently in order to emphasize the different levels of approximations and their respective prerequisites. A short reference was made to the empirical Miller-Abrahams hopping rate. Rate equation approaches to calculate the stationary charge carrier mobilities and exciton diffusion coefficients have been developed, which are based on the master equation. The rate equation approach is faster and more efficient than the frequently used Monte Carlo method and, therefore, provides the possibility to study the anisotropy of the transport parameters and their three-dimensional representation in the crystal. The Marcus theory, originally derived for outer sphere electron transfer in solvents, had already been well established for charge transport in organic solids. It was shown that this theory fits even better for excitons than for charges compared with the experiment. The Levich-Jortner theory strongly overestimates the charge carrier mobilities and the results deviate even stronger from the experiment than those obtained with the Marcus theory. The latter contains larger approximations by treating all vibrational modes classically. The spectral overlap approach in combination with the developed rate equations leads to even quantitatively very good results for exciton diffusion lengths compared to experiment. This approach and the appendant rate equations have also been adapted to charge transport. The Einstein relation, which relates the diffusion coefficient with the mobility, is important for the rate equations, which have been developed here for transport in organic crystals. It has been argued that this relation does not hold in disordered organic materials. This was analyzed within the Framework of the Gaussian disorder model and the Miller-Abrahams hopping rate.}, subject = {Exziton}, language = {en} } @phdthesis{Ruetzel2014, author = {R{\"u}tzel, Stefan}, title = {Pulse-Sequence Approaches for Multidimensional Electronic Spectroscopy of Ultrafast Photochemistry}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-98993}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2014}, abstract = {Observing chemical reactions in real time with femtosecond laser pulses has evolved into a very popular � field of research since it provides fascinating insights into the nature of photochemical transformations. Nevertheless, many photochemical reactions are still too complex for which reason the underlying mechanisms and all engaged species cannot be identi� fied thoroughly. In these cases, conventional time-resolved spectroscopy techniques reach their technical limits and advanced approaches are required to follow the conversion of reactants to their products including all reaction intermediates. The aim of this work was therefore the development of novel methods for ultrafast spectroscopy of photoreactive systems. Though the concept of coherent multidimensional spectroscopy has so far exclusively been used to explore photophysical phenomena, it also offers great potential for the study of photochemical processes due to its capability of extracting spectroscopic information along several frequency dimensions. This allows resolving the photochemical connectivity between various interconvertible molecular species with ultrafast temporal resolution on the basis of their absorption and emission properties as the spectral correlations are explicitly visualized in the detected spectra. The ring-open merocyanine form of the photochromic compound 6-nitro BIPS was studied in Chap. 4 of this work. Merocyanines and their associated ring-closed spiropyrans are promising candidates for future applications as, for instance, molecular electronics or optical data storage due to their unique property of being switchable between two stable con� gurations via light illumination. Transient absorption with sub-50 fs temporal resolution and broadband probing was employed to characterize the photodynamics of this system with variable excitation wavelengths. Using global data analysis, it could be inferred that two different merocyanine isomers with differing excited-state lifetimes exist in solution. These isomers differ in the cis/trans con� guration in the last bond of the methine bridge. The minority of isomers exist in the all-trans con� guration (TTT) while the isomer with a cis con� guration of the third dihedral angle (TTC) is dominant. A characteristic band, detected after long pump-probe delays, was attributed to the unidirectional cis->trans photoisomerization reaction of the TTC to the TTT form. The quantum yield of the reaction was estimated to be (18� +-4) \%. In addition, pronounced coherent vibrational wave-packet oscillations were observed and it was concluded that these signatures are related to the product formation. Coherent two-dimensional electronic spectroscopy was successfully implemented using a partially collinear pump-probe beam geometry in combination with a femtosecond pulse shaper. The use of a whitelight probe continuum enabled us to probe contributions far-off the diagonal over the complete visible range. By properly adjusting the relative phase between the � first two laser pulses with the pulse shaper, the principle of phase-cycling was explained and it was demonstrated that the measurement can be carried out in the so-called "rotating frame" in which the observed frequencies detected during the coherence time are shifted to lower values. It was shown that these concepts allow the extraction of the desired background-free photon echo while the amount of necessary data points is highly reduced. In order to put our proposal of multidimensional spectroscopy of photoreactive systems into practice, third-order two- and three-dimensional spectroscopy was then employed for an in-depth analysis of a photoreactive process, in which the photoisomerization of 6-nitro BIPS served as a model system. The measured two-dimensional spectra revealed the cis->trans photoisomerization after long population times. By collecting a large data set of two-dimensional spectra for short population times and by applying a Fourier transform along the population time axis, the third-order three-dimensional spectrum was obtained. The novelty of this approach compared to coherent two-dimensional spectroscopy is the introduction of a third axis associated with the vibrational frequencies of the molecular system. In this way, the formation of the reaction product was evidenced and it was shown that the product is formed in its fi� rst excited singlet state within 200 fs after excitation. This method hence visualizes the photochemical connections between different reactive molecular species in an intuitive manner and further exposes the normal modes connecting reactant and product. Such conclusions cannot be drawn with conventional third-order techniques such as transient absorption since they are not capable of capturing the full third-order response, but only a subset of it. The reaction mechanism and the role of the observed vibrational modes were uncovered by comparing the experimental data with the results of high-level quantum-chemical calculations performed by our collaborators in the group of Prof. B. Engels from the theoretical chemistry department at the University of W{\"u}rzburg. Specifi� c calculated molecular normal modes could be assigned to the experimentally observed vibrational frequencies and potential energy surfaces of the electronic ground state and of the � first excited state were computed. The technique implemented in this chapter is general and is applicable for the time-resolved analysis of a wide range of chemical reaction networks. In the fi� rst part of Chap. 5, coherent two-dimensional spectroscopy was employed to track the reaction paths of the related 6,8-dinitro BIPS after S1 excitation. Several differences to the photochemical properties of 6-nitro BIPS were found. From the 2D spectra, the cis-trans isomerization between the two merocyanine isomers could be excluded as a major reaction path for this compound. To explore the dynamics after reexcitation to higher-lying electronic states, pump-repump-probe spectroscopy was implemented and the formation of a new species, a radical cation, was observed. To identify the precursor isomer, triggered-exchange two-dimensional spectroscopy, a � fifth-order technique previously only available in the infrared regime for vibrational transitions, was implemented for the fi� rst time for electronic excitations in the visible. This approach combines the properties of the pump-repump-probe technique with the potential of coherent two-dimensional spectroscopy. It correlates the absorption frequency of a reactive molecular species with the emission signatures of the product formed from this species after an additional absorption of a photon. Using this method, it was unambiguously proven that only the TTC isomer reacts to the radical cation thus forming the precursor species of the reaction. Electronic triggered-exchange two-dimensional spectroscopy is hence another improved technology for time-resolved spectroscopy with applications in the study of multistep photoreactions and higher-lying electronic states. While in the two preceding chapters third- and � fifth-order experiments were discussed that neglect the vectorial character of light-matter interactions, Chap. 6 focused on a novel theoretical formalism enabling the description of light fi� elds optimized for polarization-sensitive higher-order nonlinearities. This formalism is based on the von Neumann time-frequency representation of shaped femtosecond laser pulses which permits the defi� nition of multipulse sequences on a discrete time-frequency lattice. Hence, not only the temporal spacing between subpulses is adjustable, but also the center frequencies may be adapted such that they � fit the experimental requirements. This method was generalized to the description of pulse sequences with time-varying polarization states. It was shown that by using this description, the polarization ellipticity, orientation angle, relative phase and intensity, and the time-frequency location of each subpulse is explicitly controllable. The accuracy of the transformations from Fourier space to von Neumann domain and vice versa was demonstrated. Moreover, a strict accordance between the von Neumann polarization parameters with the conventional parameters in time domain was found for well separated subpulses. A potential future application of this approach is polarization-sensitive multidimensional spectroscopy in which hidden cross peaks may be isolated by de� fining the pulses in the von Neumann picture with suitable polarization sequences. This method could also be used in quantum control experiments in which the polarization of the light fi� eld is used as a major control knob. This thesis summarizes our efforts to open the � field of femtochemistry to the concept of coherent multidimensional electronic spectroscopy. Making use of femtosecond pulse shaping, sub-50 fs temporal resolution, broadband spectral probing, higher-order nonlinearities, and new types of laser pulse descriptions, the presented methods might stimulate further future advancements in this research area.}, subject = {Ultrakurzzeitspektroskopie}, language = {en} } @phdthesis{Schmid2022, author = {Schmid, Paul}, title = {Quantenchemische Untersuchungen von Umgebungseinfl{\"u}ssen bei offen- und geschlossenschaligen Systemen}, doi = {10.25972/OPUS-26510}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-265106}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2022}, abstract = {In dieser Dissertation werden die Umgebungseinfl{\"u}sse auf die strukturellen und elektronischen Eigenschaften von verschiedenen offen- und geschlossenschaligen Systemen mittels quantenchemischer Methoden berechnet. Ein Kernpunkt umfasst die Untersuchung von verdreht angeordneten, biradikalischen Diborylalkenen, welche eine unges{\"a}ttigte C2R2-Br{\"u}cke (R = Et, Me) besitzen und durch cyclische (Alkyl)(amino)carbene (CAACs) stabilisiert werden. Quantenchemische Berechnungen zeigen, dass haupts{\"a}chlich sterische Effekte f{\"u}r die Ausbildung einer verdrehten Molek{\"u}lanordnung verantwortlich sind, w{\"a}hrend bei geringen sterischen Wechselwirkungen (R = H) die Delokalisationseffekte {\"u}berwiegen, wodurch eine planare Struktur beg{\"u}nstigt wird. Die Bevorzugung einer offenschaligen Singulettkonfiguration anstelle eines Tripletts ist auf den großen Energieunterschied der beiden einfach besetzten Molek{\"u}lorbitale zur{\"u}ckzuf{\"u}hren. Durch die Berechnung der L{\"o}sungsmitteleffekte mithilfe von polarisierbaren Kontinuumsmodellen kann gefolgert werden, dass mit zunehmender statischer Dielektrizit{\"a}tskonstante eine planare und geschlossenschalige Struktur st{\"a}rker stabilisiert wird als eine verdrehte Anordnung. Ein weiteres Thema dieser Dissertation befasst sich mit der quantenchemischen Analyse eines makrozyklischen Perylenbisimid-Trimersystems, welches eingebettet in einer Polymethylmethacrylat-Matrix bei Temperaturen nahe dem absoluten Nullpunkt eine Lokalisierung der ersten drei angeregten Zust{\"a}nde zeigt. Quantenchemische Vakuumberechnungen ergeben, dass unabh{\"a}ngig von der gegenseitigen geometrischen Orientierung der drei Perylenbisimid-Chromophore der {\"U}bergang vom Grundzustand in den S1-Zustand verboten ist und dass die ersten drei angeregten Zust{\"a}nde delokalisiert vorliegen. Mithilfe von expliziten L{\"o}sungsmittelmodellen kann jedoch gezeigt werden, dass das Auftreten dieser Lokalisierungen auf eine inhomogene Polymethylmethacrylat-Umgebung zur{\"u}ckzuf{\"u}hren ist, die zu einem Symmetriebruch und somit zu einer Zunahme der Oszillatorst{\"a}rke f{\"u}r S1 und der Lokalisierungsgrade f{\"u}r S1, S2 und S3 f{\"u}hrt. Dar{\"u}ber hinaus wird der L{\"o}sungsmitteleinfluss auf die angeregten Zust{\"a}nde des Azulens mittels impliziter und expliziter L{\"o}sungsmittelmodelle berechnet. Bei einer Erh{\"o}hung der dynamischen Dielektrizit{\"a}tskonstante im impliziten Modell nehmen die Anregungsenergien der vertikalen Singulettzust{\"a}nde ab, wobei der Effekt mit steigender Oszillatorst{\"a}rke zunimmt. Die Auswirkung der statischen Dielektrizit{\"a}tskonstante auf die Anregungsenergien ist dagegen deutlich schw{\"a}cher ausgepr{\"a}gt. Im expliziten Modell bewirkt das L{\"o}sungsmittel ebenfalls eine Abnahme der Anregungsenergie des hellen Singulettzustands, wenn auch in geringerem Umfang als im impliziten Modell. Als letztes Thema wird der Inhibitionsmechanismus der Cysteinprotease Rhodesain durch zwei modifizierte 1,4-Naphthoquinone untersucht. W{\"a}hrend beide Naphthoquinone an der 2-Position eine Dipeptideinheit aufweisen, besitzen sie an der 3-Position entweder einen Nitril- oder Chloridsubstituenten. Zwar erfolgt bei beiden Derivaten die Inhibition {\"u}ber einen kovalent-reversiblen Mechanismus, jedoch verl{\"a}uft die Hemmung im Falle des Nitrilderivats erheblich effektiver. Die quantenchemischen Berechnungen eines vereinfachten Modells zeigen, dass die Cysteineinheit (HS-R) bevorzugt in einer exothermen und reversiblen Additionsreaktion an die elektronenarme C-C-Doppelbindung der Naphthoquinone anlagert. Dabei werden kleinere Reaktionsenergien f{\"u}r die Reaktion des Chlorderivats als f{\"u}r die Reaktion des Nitrilderivats erhalten. Durch die Ber{\"u}cksichtigung von Wasser in einem impliziten L{\"o}sungsmittelmodell kommt es bei fast allen Reaktionsprodukten zu einer Energiezunahme, die bei der Reaktion des Nitrilderivats st{\"a}rker ausf{\"a}llt als bei der Reaktion des Chlorderivats.}, subject = {Umgebungseinfluss}, language = {de} } @phdthesis{Kaiser2022, author = {Kaiser, Dustin}, title = {Non-standard computational approaches applied to molecular systems}, doi = {10.25972/OPUS-27664}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-276641}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2022}, abstract = {In this thesis, several contributions to the understanding and modeling of chemical phenomena using computational approaches are presented. These investigations are characterized by the usage of non-standard computational modeling techniques, which is necessitated by the complex nature of the electronic structure or atomic fluctuations of the target molecules. Multiple biradical-type molecules and their spectroscopic properties were modeled. In the course of the investigation, it is found that especially the impact of correct molecular geometries on the computationally predicted absorption properties may be critical. In order to find the correct minimum geometries, Multi-Reference methods may have to be invoked. The impact of geometry relaxation on the excitonic properties of Perylene Bisimide dimers were investigated. Oftentimes, these geometry factors are neglected in Organic Semiconductor modeling as an approximation. This present investigation suggests that this approximation is not always valid, as certain regimes are identified where geometrical parameters have critical impact on the localization and energetic properties of excitons. The mechanism of the Triazolinedione (TAD) tyrosine bioconjugation reaction is investigated using quantum-chemical methods. By comparison of different conceivable mechanisms and their energetic ordering, the TAD tyrosine bioconjugation is found to proceed by means of a base-mediated electrophilic aromatic substitution reaction. The kth nearest neighbor entropy estimation protocol is investigated. This estimator promises accurate entropy estimates even for flexible molecules with multiple structural minima. Our granular investigation of formal and practical properties of the estimator suggests that the uneven variance of a molecule's vibrational modes is the cause of the observed slow convergence of the estimator. A rescaling procedure to reestablish fast convergence is suggested and benchmarks are performed.}, subject = {Quantenchemie}, language = {en} } @phdthesis{Wirsing2023, author = {Wirsing, Sara}, title = {Computational Spectroscopic Studies with Focus on Organic Semiconductor Systems}, doi = {10.25972/OPUS-28655}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-286552}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2023}, abstract = {This work presents excited state investigations on several systems with respect to experimental spectroscopic work. The majority of projects covers the temporal evolution of excitations in thin films of organic semiconductor materials. In the first chapters, thinfilm and interface systems are build from diindeno[1,2,3-cd:1',2',3'-lm]perylene (DIP) and N,N'-bis-(2-ethylhexyl)-dicyanoperylene-3,4:9,10-bis(dicarboximide) (PDIR-CN2) layers, in the third chapter bulk systems consist of 4,4',4"-tris[(3-methylphenyl)phenylamino] triphenylamine (m-MTDATA), 4,7-diphenyl-1,10-phenanthroline (BPhen) and tris-(2,4,6-trimethyl-3-(pyridin-3-yl)phenyl)borane (3TPYMB). These were investigated by aggregate-based calculations. Careful selection of methods and incorporation of geometrical relaxation and environmental effects allows for a precise energetical assignment of excitations. The biggest issue was a proper description of charge-transfer excitations, which was resolved by the application of ionization potential tuning on aggregates. Subsequent characterization of excitations and their interplay condenses the picture. Therefore, we could assign important features of the experimental spectroscopic data and explain differences between systems. The last chapter in this work covers the analysis of single molecule spectroscopy on methylbismut. This poses different challenges for computations, such as multi-reference character of low-lying excitations and an intrinsic need for a relativistic description. We resolved this by combining complete active space self-consistent field based methods with scalarrelativistic density-functional theory. Thus we were able to confidently assign the spectroscopic features and explain underlying processes.}, subject = {Theoretische Chemie}, language = {en} } @phdthesis{Herok2024, author = {Herok, Christoph}, title = {Quantum Chemical Exploration of Potential Energy Surfaces: Reaction Cycles and Luminescence Phenomena}, doi = {10.25972/OPUS-35218}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-352185}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2024}, abstract = {This work aims at elucidating chemical processes involving homogeneous catalysis and photo-physical relaxation of excited molecules in the solid state. Furthermore, compounds with supposedly small singlet-triplet gaps and therefore biradicaloid character are investigated with respect to their electro-chemical behavior. The work on hydroboration catalysis via a reduced 9,10-diboraanthracene (DBA) was preformed in collaboration with the Wagner group in Frankfurt, more specifically Dr. Sven Prey, who performed all laboratory experiments. The investigation of delayed luminescence properties in arylboronic esters in their solid state was conducted in collaboration with the Marder group in W{\"u}rzburg. The author of this work took part in the synthesis of the investigated compounds while being supervised by Dr. Zhu Wu. The final project was a collaboration with the group of Anukul Jana from Hyderabad, India who provided the experimental data.}, subject = {Simulation}, language = {en} } @phdthesis{Weickert2023, author = {Weickert, Anastasia}, title = {Theoretische Untersuchungen zur Aufkl{\"a}rung der reversiblen Hemmung durch kovalente Inhibitoren der Enzyme Golgi-alpha-Mannosidase und Rhodesain}, doi = {10.25972/OPUS-32818}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-328181}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2023}, abstract = {Die Entwicklung maßgeschneiderter Proteinliganden ist ein integraler Bestandteil unterschiedlicher wissenschaftlicher Disziplinen, wie z.B. Wirkstoffentwicklung. Die vorliegende Arbeit befasst sich mit der reversiblen Inhibition in Form von kovalent gebundenen Enzym-Inhibitor-Komplexen der humanen Golgi-alpha-Mannosidase II (GM II) und der Cysteinprotease Rhodesain. Beide Enzyme sind erfolgversprechende Targets in der Bek{\"a}mpfung von zwei sehr unterschiedlichen Erkrankungen. Einerseits die Golgi-alpha-Mannosidase bei der Behandlung der Tumorprogression und andererseits die Cysteinprotease Rhodesain bei der Behandlung der Afrikanischen Schlafkrankheit. Die Arbeit an den zwei Enzymen unterteilt sich in zwei Teilprojekte. Die Entwicklung von maßgeschneiderten kovalent-reversiblen Inhibitoren f{\"u}r die genannten Enzyme wurde im Rahmen eines in-house entwickeltes Protokolls zwecks des rationalen Designs kovalenter Inhibitoren, durchgef{\"u}hrt. Dieses Protokoll basiert auf einer sich gegenseitig unterst{\"u}tzenden Zusammenarbeit zwischen Theorie und Experiment. Die vorliegende Arbeit befasst sich mit den theoretischen Untersuchungen mit Hilfe der quantenmechanischen (QM) als auch mit Hilfe der kombinierten quantenmechanisch/molek{\"u}lmechanischen (QM/MM) Methoden zu den genannten Enzymen. In einem ersten Schritt des Protokolls geht es um die Anwendung von Screeningverfahren. In einem Screening werden Leitstrukturen, zun{\"a}chst in L{\"o}sung (Schritt I), f{\"u}r eine weitere Untersuchung im Enzym (Schritt II) evaluiert. So k{\"o}nnen die Inhibitoren, f{\"u}r die experimentelle Mess- oder theoretische Dockingdaten vorhanden sind, als eine Leitstruktur betrachtet werden. Durch das Screening unter Verwendung der quantenmechanischen (QM-Modell) Methode kann eine Reihe von Inhibitoren nach einem sich konsistent ver{\"a}nderndem Muster erstellt werden und auf Bindungsparameter hin untersucht werden (Schritt I). Diese Parameter sind Reaktionsenergien und H{\"o}hen der Reaktionsbarriere einer Inhibitionsreaktion. Reaktionsenergien werden in dieser Betrachtung quantenmechanisch innerhalb der Born-Oppenheimer-(BO)-N{\"a}herung und im Rahmen des Konzeptes der Potentialhyperfl{\"a}chen (PES) als relative Energien zwischen den optimierten Geometrien der Produkte und der Edukte auf einer Potentialhyperfl{\"a}che f{\"u}r die Inhibitionsreaktion ermittelt. Die H{\"o}hen der Reaktionsbarrieren werden durch die relativen Energien zwischen den Geometrien der Edukte und der Zwischenstufen oder {\"U}bergangszust{\"a}nde abgesch{\"a}tzt. Unter Inhibitionsreaktion wird eine chemische Reaktion verstanden, bei der eine kovalente Bindung zwischen dem Inhibitormolek{\"u}l und den Aminos{\"a}uren in der aktiven Tasche eines Enzyms ausgebildet wird. F{\"u}r den Schritt I werden die Aminos{\"a}uren der aktiven Tasche durch kleine Molek{\"u}le, wie Essigs{\"a}ure und Methanthiol, angen{\"a}hert. Die kovalent-reversiblen Inhibitoren sollten in dieser Betrachtung nur leicht exotherme Reaktionen mit den relativen Energien im Bereich -5 bis -10 kcal/mol aufweisen. Der experimentelle Teil liefert w{\"a}hrenddessen die Synthese der neuen Inhibitoren und die Nachweise zur kovalenten Bindung mit Hilfe massenspektrometrischer Messungen (Schritt I). Die passenden Kandidaten aus dem ersten Schritt des Protokolls, d.h. Inhibitoren mit gew{\"u}nschten Bindungsparametern, werden durch die QM/MM-Berechnungen im Enzym (Schritt II) und durch die experimentellen Messungen an den Enzym-Inhibitor-Komplexen in Assays (Schritt II) analysiert. Die Untersuchungen f{\"u}r die Stufe II des Protokolls umfassen die Berechnungen der Reaktionsprofile und Minimumenergiereaktionspfade f{\"u}r die chemischen Reaktionen von Inhibitoren im Zielenzym. Ein Pfad minimaler potentieller Energie, der zwei Minima (hier Edukt und Produkt) verbindet, stellt ein Reaktionsprofil f{\"u}r eine chemische Reaktion dar. In der vorliegenden Arbeit wird dies auch als Minimumenergiepfad (MEP) bezeichnet. Der Letztere l{\"a}sst sich durch die Nudged Elastic Band (NEB)-Methode und mittels Potentialhyperfl{\"a}chen darstellen. Die Reversibilit{\"a}t der Inhibitoren wurde anhand der berechneten chemischen Reaktionen in Form von erstellten Reaktionsprofilen analysiert und diskutiert. Durch Protein-Ligand Docking (Schritt III) wird ein Screening von variierbaren Erkennungseinheiten der neuen Inhibitoren durchgef{\"u}hrt. Die Ergebnisse der Untersuchungen aus dem dritten Schritt liefern Hinweise zur Weiterentwicklung der ausgew{\"a}hlten Inhibitoren. Die letzte Stufe des in-house Protokolls besteht in der erneuten Untersuchung der optimierten Inhibitoren mit Hilfe von Theorie und Experiment (Schritt IV). Die theoretische Untersuchung anhand von QM/MM-Berechnungen {\"u}berpr{\"u}ft, ob die Inhibitionsreaktion der reaktiven Kopfgruppe nach der {\"A}nderung der Erkennungseinheit des Inhibitors weiterhin effektiv und nach dem gleichen Mechanismus mit der aktiven Seite des Enzyms ablaufen kann. Die experimentelle Untersuchung liefert, {\"a}hnlich wie im Schritt II, die messbaren Ergebnisse der Inhibition in Hinblick auf die Bindungseigenschaften und die Entstehung der Nebenprodukte. Die Untersuchungen am System Mannosidase GM II wurden in Zusammenarbeit mit den Arbeitskreisen von Prof.Dr. J. Seibel und Prof.Dr. T. Schirmeister durchgef{\"u}hrt. Die Leitstruktur zur Entwicklung des kovalent-reversiblen Inhibitors stellt der cyclische O,O-Acetal-Inhibitor (bestimmt anhand von Dockingexperimenten an der beta-L-Anhydrogulose durch Arbeitskreis Prof.Dr. J. Seibel) dar. Die Ergebnisse der theoretischen Studie liefern f{\"u}r den ersten Schritt im Rahmen des Protokolls den geeigneten Kandidaten aus einer Menge von insgesamt 22 modellierten Inhibitoren f{\"u}r die reversible Inhibition der Mannosidase GM II durch die Ausbildung einer kovalenten Bindung. Hierzu z{\"a}hlen zun{\"a}chst die thermodynamischen Modellberechnungen der Inhibitionsreaktion, welche die Reaktionsenergien f{\"u}r alle Kandidaten des Screenings liefern. Die Inhibitionsreaktion wird in diesem Schritt als Additionsreaktion von Essigs{\"a}ure an den Inhibitor-Kandidaten modelliert. F{\"u}r die Leitstruktur resultiert eine thermoneutrale Beschreibung der Reaktion mit Essigs{\"a}ure und dient im Weiteren als Referenz. Der Inhibitor Nr.7 (der cyclische N,O-Acetal-Inhibitor) zeigt mit -7,7 kcal/mol eine leicht exotherme Reaktion und somit eine bessere Triebkraft der untersuchten Reaktion im Vergleich zur Referenz. Die beiden Inhibitoren wurden dann f{\"u}r Stufe 2 des Protokolls untersucht, in der eine Analyse der Reaktionsprofile im Enzym mit Hilfe der QM/MM-Methodik durchgef{\"u}hrt wurde. Die Ergebnisse des zweiten Teils der Studie zeigen, dass der cyclische N,O-Acetal-Inhibitor eine deutlich bessere Affinit{\"a}t zur aktiven Seite der GM II im Vergleich zu seiner Leitstruktur aufweisen sollte. Dies zeigt sich auch in der deutlich h{\"o}heren Triebkraft der Inhibitionsreaktion von ca. -13 kcal/mol. Dieser Energiebeitrag ist klein genug, um eine Reversibilit{\"a}t der Inhibitionsreaktion gew{\"a}hrleisten zu k{\"o}nnen. Das bedeutet auch, dass der N,O-Acetal-Inhibitor im Vergleich zur Referenzstruktur eine deutlich st{\"a}rkere Inhibition bedingen sollte. Ber{\"u}cksichtigt man dann noch, dass die Reaktion laut Berechnungen nur leicht exotherm sein sollte, erh{\"a}lt man die M{\"o}glichkeit einer reversibel stattfinden kovalenten Hemmung. Zusammenfassend liefert dieser Teil der Arbeit, der mit Hilfe der QM- und QM/MM-Berechnungen durchgef{\"u}hrt wurde, ein reaktives molekulares Ger{\"u}st mit den gew{\"u}nschten Eigenschaften. Durch die theoretischen Untersuchungen (MD-Simulationen am Enzym-Inhibitor-Komplex) konnte außerdem eine zur Komplexbildung geeignete Konformation der Leitstruktur sowie des neuen Inhibitors gefunden werden. Die reversibel agierenden Acetal-Inhibitoren befinden sich in der aktiven Tasche in einer energetisch h{\"o}her liegenden Twist-Boot-Konformation und beg{\"u}nstigen mit zwei entstehenden Bindungen zum Zn2+-Ion die oktaedrische Koordination im Enzym. Als Teil dieser Arbeit wurden NEB-Berechnungen zur Bestimmung von Minimumenergiepfadendurchgef{\"u}hrt. Dies lieferte erweiterte Einblicke in der Berechnung von Reaktionsmechanismen jeweils auch in Kombination von 2- bzw. 3-dimensionalen Scans. Auch in der Beschreibung von Protonen{\"u}bertragungsreaktionen nach Grotthus, die einem Umklappen der kovalenten Bindungen entsprechen, erh{\"a}lt man hier Geometrien f{\"u}r die Teilschritte und somit eine detaillierte Beschreibung des Vorgangs. Der Mechanismus der Inhibition von GM II durch die Leitstruktur beinhaltet einen Wasser-katalysierten (oder auch Wasser-vermittelten) Ring{\"o}ffnungsschritt in der Tasche des Enzyms. Die Testrechnungen zum Protontransfer haben gezeigt, dass der Protontransfer {\"u}ber ein oder mehrere Wassermolek{\"u}le unter Verwendung von Standard-PES-Berechnungen nicht spontan stattfindet. Die Berechnung des MEP durch das Erstellen einer 3-dimensionalen Potentialhyperfl{\"a}che kann nur dann sinnvolle Ergebnisse liefern, wenn der Protontransfer vom Aspartat Asp341 zum Inhibitor {\"u}ber zwei Wassermolek{\"u}le explizit ber{\"u}cksichtigt wird. In diesem Fall ist die Berechnung der PES kein Standard und erfordert eine zus{\"a}tzliche Variation der Bindungsabst{\"a}nde O-H der beteiligten Molek{\"u}le des Protontransfers. Die Details f{\"u}r die zus{\"a}tzliche Variation der Bindungsabst{\"a}nde O-H bei der Berechnung der 3-dimensionalen PES haben die NEB-Berechnungen geliefert. Der NEB-Formalismus hat sich in der Beschreibung dieser komplexen Reaktionskoordinaten als besser geeignet erwiesen und wurde in dieser Arbeit aus diesem Grund haupts{\"a}chlich verwendet. Die Berechnung des Protonentransfers w{\"a}hrend einer Hemmungsreaktion durch zwei Wassermolek{\"u}le mit der NEB-Methode hat den MEP ermittelt, welcher zun{\"a}chst nicht auf der Grundlage eines 3-dimensionalen Scans ermittelt werden konnte. Solche QM/MM-Rechnungen wurden im Rahmen des in-house Protokolls zum ersten Mal durchgef{\"u}hrt. Dieser Protontransfer ist mit dem Grotthus-Mechanismus konform und kann plausibel anhand einer Klapp-Mechanismus-Betrachtung nachvollzogen werden. Mit Hilfe der NEB-Methode ist es m{\"o}glich MEPs effektiv und relativ schnell zu ermitteln. Es werden sowohl die Geometrien entlang des Pfades wie auch die einzelnen relativen Energien erhalten. Zur {\"U}berpr{\"u}fung der gefundenen {\"U}bergangszust{\"a}nde wurden die einzelnen Strukturen mit Hilfe der Normalmodenanalyse weiter untersucht und konnten verifiziert werden. Die MEP-Berechnungen f{\"u}r den Inhibitor Nr.1 erm{\"o}glichen die Etablierung eines Protokolls zur Berechnung eines Reaktionspfades {\"u}ber mehrere Molek{\"u}le, welches anschließend zur Berechnung des MEP f{\"u}r den Inhibitor Nr.7 angewendet wird. Das Protokoll beinhaltet in seiner einfachen Form die Ermittlung der Two-End-Komponenten einer chemischen Reaktion - Geometrien von Reaktant und Produkt. Betrifft dies eine Reaktion, die {\"u}ber mehrere Molek{\"u}le, z.B. Wassermolek{\"u}le oder deren Netzwerk, stattfindet, wird die Aufgabe komplexer. In diesem Fall ist eine Berechnung mit Hilfe des NEB-Moduls wesentlich produktiver als die Charakterisierung mit Hilfe der 3-dimensionalen PES. Der Vorteil liegt in der kollektiven Beschreibung der Reaktionskoordinaten, sodass die entscheidenden Reaktionskoordinaten und Variablen f{\"u}r die Durchf{\"u}hrung von Scans nicht einzeln bestimmt werden m{\"u}ssen. Dennoch kann es hier bei komplexen Reaktionskoordinaten auch zu Konvergenzproblemen bzw. zu langwierigen Optimierungszyklen kommen. Als weiteres Resultat liefern die durchgef{\"u}hrten MEP-Berechnungen Einblicke in die katalytischen Eigenschaften der Wassermolek{\"u}le f{\"u}r den Proton{\"u}bertragungsmechanismus nach Grotthus. Die Daten zeigen, dass die Barriere am niedrigsten wird, wenn zwei Wassermolek{\"u}le beim Protontransfer beteiligt sind. Wenn nur ein oder gar kein Wassermolek{\"u}l die Ring{\"o}ffnung katalysiert, steigt die Barriere auf 12 und 17 kcal/mol. Die Untersuchung in diesem Teil der Arbeit l{\"a}sst zudem Einblicke in die nukleophile Substitution der Vollacetale in der Enzym-Tasche der GM II erlangen. Die Rechnungen deuten darauf hin, dass die Vollacetal-Inhibitoren durch Wassermolek{\"u}le in der Tasche aktiviert werden. Die ausgebildeten Wasserstoffbr{\"u}ckenbindungen beg{\"u}nstigen die Geometrie des Enzym-Inhibitor-Komplexes. Dies bef{\"o}rdert die Ring{\"o}ffnungreaktion gleichzeitig mit dem nucleophilen Angriff des Aspartatrestes an dem C1-Atom des Inhibitors. Im Falle des gemischten Acetal-Inhibitors hingegen wird die Treibkraft bereits durch die Einf{\"u}hrung des Stickstoffatoms deutlich erh{\"o}ht. Durch die richtig angeordneten Grotthus-Wassermolek{\"u}le ist in diesem Fall die Barriere der Proton{\"u}bertragung durch das Aspartat-Aspartat-System der GM II (Asp341/Asp240) sekund{\"a}r. Betrachtet man die Schwingungsbewegung entlang der imagin{\"a}ren Moden der {\"U}bergangszust{\"a}nde, sind diese in beiden E-I-Komplexen {\"a}hnlich. Hierbei wird eine synergistische Bewegung der Bindungsabst{\"a}nde OD2-C1-O6 (Inhibitor Nr.1) bzw. OD2-C1-N (Inhibitor Nr.7) beobachtet. Die Entwicklung der kovalent-reversiblen Inhibitoren f{\"u}r das Enzym Rhodesain wurde in Zusammenarbeit mit dem Arbeitskreis Prof.Dr. T. Schirmeister durchgef{\"u}hrt. Als Leitstruktur zur Entwicklung des neuen kovalent-reversiblen Vinylsulfon-Inhibitors 4-Pyridyl-Phenylalanyl-Homophenylalanyl-alpha-Fluor-Phenylvinylsulfon dient in diesem Projekt der kovalent-irreversibel bindende Inhibitor K777, f{\"u}r den kristallographische Daten bekannt sind. Im Rahmen des Protokolls wurde eine Reihe von Inhibitoren untersucht, in denen ein alpha-H-Atom der Vinylsulfon-Einheit (im Weiteren VS f{\"u}r Vinylsulfon) durch verschiedene Gruppen X substituiert wurde. F{\"u}r den zun{\"a}chst vorgeschlagenen Cyano-Substituent (CN) in einem VS-Inhibitor ergab sich bei Berechnungen in einem polaren L{\"o}sungsmittel eine relativ niedrige Reaktionsenergie, d.h. es wurde eine reversible Reaktion vorhergesagt. Dies wurde experimentell best{\"a}tigt. Die theoretischen und experimentellen Ergebnisse von Schritt II widersprechen sich aber. W{\"a}hrend experimentell eine schwache reversible Hemmung gefunden wurde, sagten die Berechnungen keine Hemmung voraus. Tats{\"a}chlich zeigte sich im Nachhinein, dass die experimentell gefundene Hemmung nicht-kompetitiv ist, d.h. nicht in der aktiven Tasche stattfindet. Im Rahmen des Protokolls werden dagegen nur die kompetitiv interagierenden Inhibitoren ausgewertet. An dieser Stelle lassen sich die anhand theoretischen Methoden erhaltenen Daten {\"u}ber die Reversibilit{\"a}t der Hemmung (Reaktion im aktiven Zentrum) mit den Ergebnissen aus den experimentellen Untersuchungen (Reaktion außerhalb des aktiven Zentrums) nicht vergleichen. Durch den Wechsel von CN zu Halogenen wurde schließlich eine neue Reihe von Inhibitoren auf VS-Basis entwickelt. Die Berechnungen von Reaktionsenergien in L{\"o}sung und im Enzym haben f{\"u}r diese Inhibitoren eine reversible Hemmung vorhergesagt. Allerdings findet man eine einfache Additionsreaktion an der Doppelbindung der Vinylsulfon-Gruppe. F{\"u}r X=CN wurde von einer SN2-Reaktion ausgegangen. F{\"u}r X=Br fand man, dass sich nach der Addition ein HBr-Molek{\"u}l abspaltet, sodass die Hemmung insgesamt irreversibel ist. Da die Substitutionsreaktion ein irreversibler Prozess ist und die Freisetzung von Bromwasserstoff durch die experimentellen Untersuchungen best{\"a}tigt werden konnte, scheint Fluor ein geeigneter Substituent zu sein (X=F). Hier konnte man auch experimentell eine deutlich bessere Hemmung messen. Es wurden daher die Berechnungen im Enzym f{\"u}r Systeme mit den Inhibitoren K777-X mit X=F und X=H (K777-F- und K777-H-Inhibitor) durchgef{\"u}hrt und analysiert. In der vorliegenden Arbeit wurde versucht, die Reversibilit{\"a}t des K777-F-Inhibitors gegen die Irreversibilit{\"a}t des K777-H-Inhibitors durch die quantenmechanischen Berechnungen im Rahmen des Protokolls darzulegen. Die QM/MM-Berechnungen unterteilen sich in zwei Bereiche. Zun{\"a}chst wurde das Reaktionsprofil (auch Reaktionspfad) der Additionsreaktion des K777-X-Inhibitors an die aktive Tasche von Rhodesain ausgehend von der vorhandenen Kristallstruktur (PDB-Datei) berechnet. Im Folgenden wird dieses Teilergebnis als XP-Berechnung (im Weiteren XP f{\"u}r X-ray-Pfad) bezeichnet. Alle vier PES (X=H, F, Br und Cl) weisen prinzipiell die gleiche Form auf. Es ergeben sich aber Unterschiede in den berechneten Reaktionsenergien der Additionsreaktion (-20, -16, -10 und -11 kcal/mol). Die Reaktionsenergien der Substituenten Brom und Chlor entsprechen dem Bereich f{\"u}r reversible Reaktionen (ca. -10 kcal/mol), wobei Fluor mit -16 kcal/mol einen Grenzfall darstellt. Die Konturen der beiden PES (X=H vs. X=F) sind allerdings sehr {\"a}hnlich: In beiden F{\"a}llen findet sich f{\"u}r das anionische Intermediat kein Minimum. In der Potentialhyperfl{\"a}che f{\"u}r X=F steigt die Barriere der R{\"u}ckreaktion zwischen dem Intermediat und dem nicht-kovalenten Komplex auf etwa 5 kcal/mol an, die R{\"u}ckreaktion ist im Vergleich zu dem X=H mit ca. 1,5 kcal/mol leicht exotherm. Das ver{\"a}nderte Verh{\"a}ltnis zwischen der H{\"o}he der Reaktionsbarriere und dem Betrag der Reaktionsenergie (der {\"U}bergang von der endothermen zur exothermen Reaktion) auf diesem Abschnitt der PES k{\"o}nnte dazu beitragen, dass die Gesamtreaktion insgesamt reversibel ablaufen kann. Die Reversibilit{\"a}t des Inhibitors mit dem Substituenten Fluor l{\"a}sst sich auf diesem Schritt der Untersuchung durch die Absenkung der Reaktionsenergie der Additionsreaktion auf etwa 16 kcal/mol erkl{\"a}ren, denn die irreversible Reaktionen wurden bisher mit deutlich h{\"o}heren Reaktionsenergien assoziiert. Die erhaltenen nicht-kovalenten Enzym-Inhibitor-Komplexe der XP-Berechnungen wurden in einem zweiten Teilergebnis weiter verwendet, indem der Reaktionspfad der Additionsreaktion des K777-X-Inhibitors vom nicht-kovalenten Enzym-Inhibitor-Komplex zum kovalenten Enzym-Inhibitor-Komplex hin berechnet wurde. Im Folgenden wird dieses Teilergebnis als NP-Berechnung (NP f{\"u}r Nicht-kovalente-Pfad) bezeichnet. Die Berechnung der Reaktionsprofile der Additionsreaktion des VS-Inhibitors f{\"u}r X=H und X=F am alpha-Kohlenstoffatom der VS-Kopfgruppe lieferte konsistente Ergebnisse in Bezug auf die Reaktionsenergien. {\"A}hnlich den XP-Berechnungen, wurde ebenfalls die Tendenz der Absenkung der Reaktionsenergie von -7 kcal/mol (X=H) auf -4,3 (X=F) und -0,9 kcal/mol (X=Br) beobachtet. Die Thermodynamik der Additionsreaktion wurde durch einen Wechsel des Substituenten X von H nach F in der VS-Kopfgruppe des K777-X Inhibitors beeinflusst, indem die niedrigere Energiedifferenz zwischen den Edukten und den Produkten erzielt werden konnte. F{\"u}r beide Teile der Arbeit (XP- und NP-Berechnungen) implizierte dies einen Wechsel von einem irreversiblen zum einem reversiblen Verlauf in der Beschreibung der Reaktionsprofile. Die Ergebnisse des zweiten Teils der Arbeit (NP-Berechnungen) liefern nicht nur die konsistent geringeren Reaktionsenergien (Thermodynamik) sondern auch die h{\"o}heren Reaktionsbarrieren der Additionsreaktion im Vergleich zu den Ergebnissen der XP-Berechnungen. Die {\"A}nderungen der Reaktionsbarrieren im NP-Ansatz weisen zus{\"a}tzliche Diskrepanzen auf, wenn diese jeweils mittels der PES-Scan- und CI-NEB-Dimer-Methodik berechnet werden. W{\"a}hrend die Barriere des irreversiblen Inhibitors K11777 mit dem NEB-Ansatz ca. 11 kcal/mol betr{\"a}gt und durch den PES-Scan nur um 4 kcal/mol h{\"o}her liegt, ergibt sich eine umgekehrte Situation beim {\"U}bergang zu Fluor als Substituent: Durch die NEB-Berechnung liegt die Barriere bei ca. 18 kcal/mol und durch den PES-Scan ergibt sich eine um 4 kcal/mol niedrig liegende Barriere. Um die Ergebnisse der NP-Berechnungen zu {\"u}berpr{\"u}fen, wurden diese QM/MM-Rechnungen wiederholt durchgef{\"u}hrt. In den beiden neu durchgef{\"u}hrten Berechnungen f{\"u}r die Inhibitoren K777-X mit X=H und X=F wurden nur sehr kleine Abweichungen gefunden, die kleiner als die Fehler der Berechnung sind. Die Startstrukturen f{\"u}r die Berechnung des MEP stammten aus der erneut durchgef{\"u}hrten MD-Simulation an der Geometrie des nicht-kovalenten Enzym-Inhibitor-Komplexes, welche die XP-Berechnung resultierte. Die Gesamtdauer der MD-Simulation wurde zu einem Wert von 9 ns gew{\"a}hlt, welche insgesamt 900 Startstrukturen entlang der Simulation lieferte. Die Berechnung ergab die Reaktionsenergie von -8,4 kcal/mol (-7,0 kcal/mol als erstes Ergebnis) und die relative Energie des Int-Komplexes von 13,2 kcal/mol. Somit betr{\"a}gt die Barriere der R{\"u}ckreaktion zur Freisetzung des Inhibitors K11777 (X=H) in Form von einem nicht-kovalenten Enzym-Inhibitor-Komplex einen Wert von 21,6 kcal/mol. In analoger Vorgehensweise wurde die Evaluierung der NP-Berechnung f{\"u}r den Inhibitor K777-X mit X=F durchgef{\"u}hrt. Die Reaktionsenergien in den beiden Berechnungen unterscheiden sich in einem marginalen Abstieg zu den Werten von -2,9 kcal/mol und -3,3 kcal/mol (-4,3 kcal/mol als erstes Ergebnis). Beide Berechnungen liefern zudem die relativ kleinen Anstiege der Reaktionsbarriere zu den Werten von 19,8 kcal/mol und 20,9 kcal/mol. F{\"u}r die Inhibitoren K777-X mit X=H und X=F entsprechen die gefundenen Barrieren einer verzerrten Konformation des nicht-kovalenten Enzym-Inhibitor-Komplexes, die als eine bioaktive Konformation bezeichnet werden kann. Der anionische {\"U}bergangszustand Int*, der oft in der Literatur als ein anionisches Intermediat der Additionsreaktion bezeichnet wird, wurde nur f{\"u}r den Inhibitor mit dem Substituenten Brom (K777-X mit X=Br) identifiziert. Da der {\"U}bergangszustand (Int* mit der relativen Energie von 11,1 kcal/mol) nur 1,5 kcal/mol {\"u}ber der bioaktiven Konformation (Int mit der relativen Energie von 12,6 kcal/mol) liegt und die NEB-Reaktionspfade alleine die Barrieren {\"u}bersch{\"a}tzen, besitzen die anionischen {\"U}bergangszust{\"a}nde der Additionsreaktion der Inhibitoren mit X=H und X=F eine gesch{\"a}tzte relative Energie mit vergleichbaren Abweichungen von ca. 2 kcal/mol zu den identifizierten Int-Geometrien. Die durchgef{\"u}hrten Berechnungen zeigen, dass die Substituenten X=H und X=F im Vergleich zum Brom die anionischen Geometrien der nicht-kovalenten Enzym-Inhibitor-Komplexe jedoch mangelhaft bis ausreichend stabilisieren k{\"o}nnen. Zus{\"a}tzlich liegt die gesch{\"a}tzte Energiedifferenz zwischen den Geometrien Int* und Int unter der m{\"o}glichen Fehlergrenze der Berechnungen (ca. 3-4 kcal/mol). Aus diesem Grund misslang die Optimierung in Richtung der metastabilen anionischen Geometrien Int* mit Hilfe der CI-NEB-Dimer-Methodik im Fall der VS-Inhibitoren K777-X mit X=H und X=F. Der direkte Vergleich der geometrischen Parameter der nicht-kovalenten Enzym-Inhibitor-Komplexe f{\"u}r den Inhibitor K777-F aus den XP-Berechnungen mit solchen aus den NP-Berechnungen l{\"a}sst darauf schließen, dass die Geometrien der Enzym-Inhibitor-Komplexe der XP-Berechnung nur die lokalen Minima mit der verzerrten Geometrie des Inhibitors auf der PES darstellen und die Gesamtinformation {\"u}ber die Barrieren der Reaktion durch die Ergebnisse aus der NP-Berechnung erg{\"a}nzt werden sollten. Zusammenfassend sagen die Berechnungen f{\"u}r die reaktiven Kopfgruppen der Substanzklasse der halogenierten Vinylsulfone K777-X (X=Br, Cl und F) im Vergleich zur Leitstruktur des Vinylsulfon-Inhibitors K11777 deutlich geringere exotherme Additionsreaktionen im aktiven Zentrum von Rhodesain voraus. Dar{\"u}berhinaus konnte anhand der QM/MM-Berechnungen ein experimentell gemessenen verlangsamten Verlauf der reversiblen Inhibition im Falle von X=F (Inhibitor K777-X) durch die relativ erh{\"o}hte Reaktionsbarriere im Vergleich zur Leitstruktur erkl{\"a}rt werden. Dieser Inhibitor dient zun{\"a}chst als ein erfolgreich selektiertes reaktives Ger{\"u}st des neuen Inhibitors K777-X-S3 mit X=F und S3=4-Pyridyl (K777-F-Pyr), welcher mit Hilfe des Docking-Experiments (Schritt III durch die Arbeitsgruppe Prof.Dr. T. Schirmeister) deutlich verbessert werden konnte. Die Affinit{\"a}t des durch Docking verbesserten VS-Inhibitors mit Fluor als Substituent durch die eingef{\"u}hrte Seitenkette S3=4-Pyridyl (4-Pyridyl-Phenylalanyl-Homophenylalanyl-(Phenyl)-alpha-F-Vinylsulfon) stieg im Rhodesain von 190 nM zu 32 nM (Schritt IV, experimenteller Teil). Gleichzeitig konnte durch die QM/MM-Berechnungen in Schritt IV gezeigt werden, dass die Reaktion der reaktiven Kopfgruppe im neuen Inhibitor immer noch eine kovalent-reversible Hemmung von Rhodesain darstellt, auch wenn die Erkennungseinheit ge{\"a}ndert wurde. Hierf{\"u}r kann man die beiden Reaktionsprofile der NP-Berechnungen vergleichen. Die beiden fluorierten VS-Inhibitoren weisen eine {\"A}hnlichkeit bez{\"u}glich der Barrierenh{\"o}he und der Reaktionsenergie auf. Der fluorierte Vinylsulfon-Inhibitor K777-F wurde somit als ein neuer kovalent-reversibler Vinylsulfon-Inhibitor der Cysteinprotease Rhodesain erfolgreich eingef{\"u}gt.}, subject = {Computational Chemistry}, language = {de} } @phdthesis{Maier2021, author = {Maier, Jan Richard}, title = {Investigations of Metal-free Cannibalistic Hexadehydro-Diels-Alder and Pt-catalyzed Dimerization Reactions of Linked Aryl Bisdiynes}, doi = {10.25972/OPUS-24041}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-240411}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2021}, abstract = {The introductory chapter reviews the current state of mechanistic understanding of the hexadehydro-Diels-Alder (HDDA) reaction. With the rapid development of the HDDA reaction from its first discovery in 1997, the question of whether a concerted or stepwise mechanism better describes the thermally activated formation of ortho-benzyne from a diyne and a diynophile has been debated. Mechanistic and kinetic investigations were able to show that this is not a black or white situation, as minor changes can tip the balance. In chapter 2 of this thesis, the catalytic process leading from 1,11-bis(p-tolyl)undeca-1,3,8,10-tetrayne to fully-substituted naphthalene and azulene derivatives, by two different platinum-catalyzed dimerization pathways, was investigated. In chapter 3, the cannibalistic self-trapping reaction of an ortho-benzyne derivative generated from 1,11-bis(p-tolyl)undeca-1,3,8,10-tetrayne in an HDDA reaction was investigated. Without adding any specific trapping agent, the highly reactive benzyne is trapped by another bisdiyne molecule in at least three different modes. In chapter 4 direct UV/VIS spectroscopic evidence for the existence of an o-benzyne in solution is reported, and the dynamics of its formation in a photo-induced reaction are established. For this purpose, 1,11-bis(p-tolyl)undeca-1,3,8,10-tetrayne was investigated, using femtosecond transient absorption spectroscopy in the ultraviolet/visible region. In chapter 5, following the isolation and characterization of the reaction products discussed in chapter 3, further species resulting from reactions of the highly reactive ortho-benzyne derivative were identified.}, subject = {Diels-Alder-Reaktion}, language = {en} } @phdthesis{Sauer2021, author = {Sauer, Susanne}, title = {Implementation and Application of QM/MM Hybrid Methods}, doi = {10.25972/OPUS-24321}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-243213}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2021}, abstract = {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.}, subject = {Quantenmechanik}, language = {en} } @phdthesis{Hirsch2021, author = {Hirsch, Florian}, title = {IR/UV Double-Resonance Spectroscopy of Reactive Hydrocarbon Species and their Reaction Products in Cold Molecular Jets}, doi = {10.25972/OPUS-25175}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-251755}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2021}, abstract = {Reactive hydrocarbon species are important in a multitude of different scientific areas. In this thesis, the vibrational spectra of hydrocarbon radicals, biradicals and their reaction product have been studied in a gas-phase environment. The specific molecules investigated here, are of particular importance in the field of combustion and astrochemistry. They were produced from suitable precursors in a pyrolytically heated micro-reactor and subsequently seeded in an appropriate carrier gas. As methodology, IR/UV ion dip spectroscopy has been utilized, which delivers massselected gas-phase IR spectra of all ionizable species detectable in the molecular beam. These, with the help of DFT calculations, allow for determination of the fingerprint IR spectra, identification of mass carriers and formulation of potential reaction mechanisms. All studies have been conducted in collaboration with the group of Prof. Dr. Anouk M. Rjis and the necessary potent IR radiation has been provided by the free-electron laser FELIX. Thus, the IR/UV measurements have been executed at the FELIX Laboratory of the Radboud University in Nijmegen. The first study presented in this thesis is the investigation of ortho-benzyne in Chapter 3.1. This molecule is of particular interest due to its uncommon electronic structure and its role in high-temperature reactions. Although, the infrared spectrum of o-C6H4 was not accessible, a number of reaction products were identified via their fingerprint spectra. Masses in the range from 78 - 228 were assigned to their respective carrier. The identified species include typical PAHs like naphthalene, phenanthrene, up to triphenylene. The identified masses further suggest a PAH growth heavily influenced by diradical 1,4-cycloaddition followed by fragmentation, as well as by classical HACA- and PAC-like mechanisms. These results were augmented by threshold photoionization measurements from Engelbert Reusch, who identified lighter reaction products, which have insufficient IR absorption or unsuitable ionization characteristics to be identified in the IR/UV experiment. An interesting observation is the identification of m/z = 152. This carrier has been assigned differently by the IR and TPES experiments. Whereas the IR spectrum clearly identifies the species as 2-ethynylnaphthalene, the TPES evidently is in great agreement with biphenylene. This is a good example how different experimental methodologies can benefit from each other to gain a deeper insight into the actual science of a particular system. Probably, the prime example for an aromatically resonance stabilized radical is benzyl. This radical is of high importance for many combustion studies, as it represents the primary high-temperature decomposition product of toluene. The goal of the study was the identification of the benzyl self reaction products and the results are discussed in Section 3.2. The radical was pyrolytically produced by its respective nitrite precursor. The mass spectrum showed that the benzyl self reaction formed two products with C11 and three with C14 constitution. All mass peaks were evenly spaced by two mass units, respectively, which suggests a close relation in formation. Indeed, the C11 products were identified as diphenylmethane and fluorene, which are simply connected via cyclization. The heaviest product was identified as phenanthrene, which is formed via the cyclization of bibenzyl to 9,10-dihydrophenanthrene and subsequent elimination of hydrogen. This result was quiet interesting as the intermediate of this reaction was often assumed to be stilbene, which was not observed in the study. Hence, the reaction seems to undergo cyclization first before phenanthrene is finally formed via hydrogen elimination. Expanding the molecular frame of benzyl by an additional methyl group leads to the xylyl radicals and its decomposition product the xylylenes. Also important in combustion research, xylyl radicals represent the preferred decomposition products of xylene, a frequently used anti-knock agent in modern gasoline blends. After further hydrogen elimination the xylyl radicals can then form their respective xylylenes. The results of the xylyl experiments are discussed in Section 3.3. Here the gas-phase vibrational spectrum in the fingerprint region for all three isomers has been recorded for the first time in isolation. Although, all isomers have a very similar structure and symmetry, and consequently similar vibrational bands, the resolution of the experimental data was exceedingly sufficient for a clear assignment. Additionally, the dimerization products of meta- and para-xylyl could also be identified. A similar approach was taken to determine the fingerprint spectra for the xylylenes. Here, only para-xylylene could be unambiguously identified as the carrier of mass 104. For both ortho- and meta-xylylene precursors, only isomerization products were observed as the carriers of mass 104; benzocyclobutene and styrene, respectively. A possible explanation is elaborated upon in the troubleshooting Sec- tion 3.4.3.5. In the final experimental section a study on the decomposition of phthalide is presented. The objective of this experiment was mainly focused around the formation of C7 species, particularly the fulvenallenyl radical C7H5. In fact, the first experimental fingerprint spectrum of isolated C7H5 in the gas-phase was measured and is displayed in Fig. 3.45. Furthermore, the experiment demonstrates that the pyrolysis products of phthalide are excellent soot precursors, as many heavier reaction products have been identified. These include typical PAH species like naphthalene and phenanthrene as well as their methylated isomers. A large number of molecules with terminal ethynyl moieties indicate a strong influence of HACA growth in the experimental environment. However, many formation pathways of products have been discussed, which are formed involving experiment specific species, like C5H5 and C7H5, and often include expansion steps from 5- to 6-membered rings.}, subject = {Infrarotspektroskopie}, language = {en} } @phdthesis{Kudriashova2019, author = {Kudriashova, Liudmila}, title = {Photoluminescence Reveals Charge Carrier Recombination in Organic and Hybrid Semiconductors}, doi = {10.25972/OPUS-19343}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-193437}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2019}, abstract = {In this work, we elucidated recombination kinetics in organic and hybrid semiconductors by steady-state and time-resolved PL spectroscopy. Using these simple and very flexible experimental techniques, we probed the infrared emission from recombining free charge carriers in metal-halide perovskites, as well as the deep blue luminescence from intramolecular charge-transfer states in novel OLED emitters. We showed that similar state diagrams and kinetic models accurately describe the dynamics of excited species in these very different material systems. In Chapters 4 and 5, we focused on lead iodide perovskites (MAPI and FAPI), whose comparatively developed deposition techniques suited the systematic material research. In MAPI, we harnessed the anomalous dependence of transient PL on the laser repetition rate in order to investigate the role of interfaces with the commonly used charge-selective layers: PC60BM, spiro-MeOTAD, and P3HT. The film was deposited on a large precut substrate and separated into several parts, which were then covered with the charge-selective layers. Thereby, the same bulk perovskite structure was maintained for all samples. Consequently, we were able to isolate interface-affected and bulk carrier recombination. The first one dominated the fast component of PL decay up to 300 ns, whereas the last was assigned to the remaining slow component. The laser repetition rate significantly prolonged PL decay in MAPI with additional interfaces while shortening the charge carrier lifetime in the pristine film. We qualitatively explained this effect by a kinetic model that included radiative electron-hole recombination and nonradiative trap-assisted recombination. All in all, we showed that the apparent PL lifetime in MAPI is to large extend defined by the laser repetition rate and by the adjacent interfaces. Further, we studied photon recycling in MAPI and FAPI. We monitored how the microscopic PL transforms while propagating through the thin perovskite film. The emission was recorded within 5orders of magnitude in intensity up to 70μm away from the excitation spot. The Beer-Lambert law previously failed to describe the complex interplay of the intrinsic PL spectrum and the additional red-shifted peak. Therefore, we developed a general numerical model that accounts for self-absorption and diffusion of the secondary charge carriers. A simulation based on this model showed excellent agreement with the experimental spatially resolved PL maps. The proposed model can be applied to any perovskite film, because it uses easily measurable intrinsic PL spectrum and macroscopic absorption coefficient as seeding parameters. In Chapter 6, we conducted an extensive photophysical study of a novel compact deep blue OLED emitter, SBABz4, containing spiro-biacridine and benzonitrile units. We also considered its single-donor monomer counterpart, DMABz4, in order to highlight the structure-property relationships. Both compounds exhibited thermally activated delayed fluorescence (TADF), which was independently proven by oxygen quenching and temperature-dependent transient PL measurements. The spiro-linkage in the double-donor core of SBABz4 rendered its luminescence pure blue compared to the blue-green emission from the single-donor DMABz4. Thus, the core-donor provided desirable color tuning in the deep blue region, as opposed to the common TADF molecular design with core-acceptor. Using PL lifetimes and efficiencies, we predicted EQEmax = 7.1\% for SBABz4-based OLED, whereas a real test device showed EQEmax = 6.8\%. Transient PL was recorded from the solutions and solid films in the unprecedentedly broad dynamic range covering up to 6orders of magnitude in time and 8orders of magnitude in intensity. The stretched exponent was shown to fit the transient PL in the films very well, whereas PL decay in dilute solution was found purely exponential. When the emitter was embedded in the host matrix that prevented aggregation, its TADF properties were superior in comparison with the pure SBABz4 film. Finally, using temperature-dependent transient PL data, we calculated the TADF activation energy of 70 meV. To sum up, this Thesis contributes to the two fascinating topics of the last decade's material research: perovskite absorbers for photovoltaics and TADF emitters for OLEDs. We were lucky to work with the emerging systems and tailor for them new models out of the well-known physical concepts. This was both exciting and challenging. In the end, science of novel materials is always a mess. We hope that we brought there a bit of clarity and light.}, subject = {Time-resolved photoluminescence}, language = {en} } @phdthesis{Suess2021, author = {S{\"u}ß, Jasmin}, title = {Theoretische Untersuchungen an molekularen Aggregaten: 2D-Spektroskopie und Exzitonendynamik}, doi = {10.25972/OPUS-24713}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-247136}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2021}, abstract = {Diese Dissertation besch{\"a}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{\"u}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{\"u}nfter Ordnung in Abh{\"a}ngigkeit der Populationszeit T zu. Anhand des Dimer-Modells als kleinstm{\"o}gliches Aggregat lassen sich bereits gute Vorhersagen auch {\"u}ber das Verhalten gr{\"o}ßerer molekularer Aggregate treffen. Der Zerfall des oszillierenden Signals f{\"u}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{\"a}nden, der essentiell f{\"u}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{\"a}chst zueinander bewegen m{\"u}ssen, bevor die Startbedingung des Annihilierungsprozesses gegeben ist. Dieser zus{\"a}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{\"o}ßeren Modellsystems, bei dem sich der Effekt der zeitverz{\"o}gert eintretenden EEA deutlich in der Zerfallsdynamik manifestieren kann.}, subject = {Molekulardynamik}, language = {de} } @phdthesis{Endres2024, author = {Endres, Erik}, title = {Kovalente Inhibitoren: Modellierung und Design}, doi = {10.25972/OPUS-35933}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-359330}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2024}, abstract = {Kovalente Inhibition stellt einen effektiven Weg dar, die Verweildauer des Liganden innerhalb einer Bindetasche zu erh{\"o}hen. In dieser Arbeit wurden theoretische Methoden angewendet, um die Reaktivit{\"a}t und den nichtkovalenten Zustand vor der Reaktion zu modellieren. Im Rahmen einer Fallstudie zu Cathepsin K wurden nichtkovalente Modelle von kovalenten Inhibitoren generiert. F{\"u}r verschiedene Komplexe aus Cathepsin K und einem kovalent gebundenem Liganden wurde der Zustand vor der Reaktion modelliert und dessen Stabilit{\"a}t im Rahmen einer klassischen MD-Simulation {\"u}berpr{\"u}ft. Die Stabilit{\"a}t des Warheads in der Bindetasche hing haupts{\"a}chlich vom gew{\"a}hlten Protonierungszustand der katalytischen Aminos{\"a}uren ab. F{\"u}r eine Reihe von Inhibitoren der ChlaDUB1 wurde ein Protokoll aus quantenmechanischen Rechnungen genutzt, um die Reaktivit{\"a}t verschiedener Warheads abzusch{\"a}tzen. Die erhaltenen Aktivierungsenergien korrelierten mit experimentell bestimmten Raten zur Inaktivierung des Enzyms. Im Rahmen eines Wirkstoffdesign-Projektes zur Deubiquitinase USP28 wurden von unpublizierten Kristallstrukturen ausgehend erste Docking-Experimente durchgef{\"u}hrt. Es konnte gezeigt werden, dass ein literaturbekannter Inhibitor von USP28 mit einem Warhead so modifiziert werden kann, dass die reaktive Einheit in direkter Nachbarschaft zu einem Cystein positioniert wird. F{\"u}r diese Warheads wurden ebenfalls quantenmechanische Rechnungen zur Bestimmung der Aktivierungsenergie durchgef{\"u}hrt. Um besser nachvollziehen zu k{\"o}nnen, warum bei einem Photoswitch-Inhibitor der Butyrylcholin-Esterase der cis-Zustand des Molek{\"u}ls besser inhibiert als der trans-Zustand, wurde eine Docking-Studie des Zustandes vor der Reaktion durchgef{\"u}hrt. Es konnte ein qualitatives Modell aufgestellt werden, das zeigt, dass der trans-Zustand aufgrund seiner l{\"a}ngeren Form mit wichtigen Aminos{\"a}uren am Eingang der Bindungstasche kollidiert.}, subject = {Molekulardynamik}, language = {de} }