Refine
Has Fulltext
- yes (7)
Is part of the Bibliography
- yes (7) (remove)
Document Type
- Doctoral Thesis (7)
Keywords
- Quantenchemie (7) (remove)
Institute
- Institut für Physikalische und Theoretische Chemie (7) (remove)
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.
Quantenchemische Untersuchungen von Umgebungseinflüssen bei offen- und geschlossenschaligen Systemen
(2022)
In dieser Dissertation werden die Umgebungseinflü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ättigte C2R2-Brücke (R = Et, Me) besitzen und durch cyclische (Alkyl)(amino)carbene (CAACs) stabilisiert werden. Quantenchemische Berechnungen zeigen, dass hauptsächlich sterische Effekte für die Ausbildung einer verdrehten Molekülanordnung verantwortlich sind, während bei geringen sterischen Wechselwirkungen (R = H) die Delokalisationseffekte überwiegen, wodurch eine planare Struktur begünstigt wird. Die Bevorzugung einer offenschaligen Singulettkonfiguration anstelle eines Tripletts ist auf den großen Energieunterschied der beiden einfach besetzten Molekülorbitale zurückzuführen. Durch die Berechnung der Lösungsmitteleffekte mithilfe von polarisierbaren Kontinuumsmodellen kann gefolgert werden, dass mit zunehmender statischer Dielektrizitätskonstante eine planare und geschlossenschalige Struktur stä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ände zeigt. Quantenchemische Vakuumberechnungen ergeben, dass unabhängig von der gegenseitigen geometrischen Orientierung der drei Perylenbisimid-Chromophore der Übergang vom Grundzustand in den S1-Zustand verboten ist und dass die ersten drei angeregten Zustände delokalisiert vorliegen. Mithilfe von expliziten Lösungsmittelmodellen kann jedoch gezeigt werden, dass das Auftreten dieser Lokalisierungen auf eine inhomogene Polymethylmethacrylat-Umgebung zurückzuführen ist, die zu einem Symmetriebruch und somit zu einer Zunahme der Oszillatorstärke für S1 und der Lokalisierungsgrade für S1, S2 und S3 führt.
Darüber hinaus wird der Lösungsmitteleinfluss auf die angeregten Zustände des Azulens mittels impliziter und expliziter Lösungsmittelmodelle berechnet. Bei einer Erhöhung der dynamischen Dielektrizitätskonstante im impliziten Modell nehmen die Anregungsenergien der vertikalen Singulettzustände ab, wobei der Effekt mit steigender Oszillatorstärke zunimmt. Die Auswirkung der statischen Dielektrizitätskonstante auf die Anregungsenergien ist dagegen deutlich schwächer ausgeprägt. Im expliziten Modell bewirkt das Lö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ä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 über einen kovalent-reversiblen Mechanismus, jedoch verlä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ür die Reaktion des Chlorderivats als für die Reaktion des Nitrilderivats erhalten. Durch die Berücksichtigung von Wasser in einem impliziten Lösungsmittelmodell kommt es bei fast allen Reaktionsprodukten zu einer Energiezunahme, die bei der Reaktion des Nitrilderivats stärker ausfällt als bei der Reaktion des Chlorderivats.
Diese Arbeit befasst sich mit verschiedenen Aspekten der Dynamik von Kernen, Elektronen und gekoppelten Kern-Elektron-Systemen, wobei je nach System unterschiedliche Herangehensweisen gewählt wurden.
Zentrale Punkte sind bei allen drei Kapiteln einerseits die Lokalisierung von Teilchen und Energie und andererseits eine hohe Sensitivität in Bezug auf die Wahl der Anfangsbedingungen.
Im ersten Teil wurden von der Carrier-Envelope-Phase (CEP) abhängende, laser-induzierte Lokalisierungen betrachtet. Das zentrale Element ist dabei das entwickelte Doppelpulsschema, mit welchem eine CEP-Abhängigkeit in beobachtbaren Größen erzeugt wird. Als Beispielsysteme wurden die Fragmentation im D₂⁺-Modellsystem und eine Isomerisierung im Doppelminimumpotential (DMP) untersucht. Als Observable wird die Asymmetrie betrachtet
Im DMP kann die Asymmetrie mit dem Entantiomeren/Isomerenüberschuss gleich gesetzt werden kann und im D₂⁺-Modellsystem mit der Lokalisierung des Elektrons auf einem der beiden dissoziierenden Kerne.
Eine Phasenabhängigkeit der Asymmetrien besteht nur für die CEP des zweiten Pulses φ₂, für welchen keine Begrenzungen für die Anzahl an Laserzyklen auftreten.
Im DMP wurde die CEP-Abhängigkeit der Asymmetrien auch bei unterschiedlichen Startkonfigurationen untersucht. Für alle untersuchten Startkonfigurationen konnte ein Laserparametersatz gefunden werden, der für zumindest eine der beiden Asymmetrien eine CEP-Abhängigkeit liefert. Aufgrund der aufgehobenen energetischen Entartung der Paare gerader und ungerader Symmetrie ist die resultierende Lokalisierung zeitabhängig. Zur Messung der vorhergesagten Dynamiken ist z.B. die Aufnahme eines Photoelektronen-Spektrums denkbar.
In nächsten Kapitel wurden unterschiedliche Dynamiken innerhalb eines 4d Kern-Elektron-Modells in der Nähe einer konischen Durchschneidung (CI) zweier Potentiale betrachtet. Hierbei ist hervorzuheben, dass eine solche gleichzeitige Untersuchung von Kern- und Elektron-Dynamik in Systemen mit CIs in der Literatur, nach Wissen des Autors, bisher nicht veröffentlicht ist. Das 4d-Potential wurde mit Hilfe des sogenannten Potfit-Algorithmus gefittet. Dieser Fit wurde anschließend verwendet, um die Dynamik des gekoppelten Systems mit Hilfe der ”Multi-Configuration Time-Dependent Hartree”(MCTDH)-Methode zu berechnen. Aus der Analyse der gekoppelten Kern-Elektron-Wellenfunktion ergaben sich zwei grundlegend unterschiedliche Klassen von Dynamiken:
• Diabatisch: Kern- und Elektrondynamik sind nahezu entkoppelt. Der Kern bewegt sich und das Elektron bleibt statisch.
• Adiabatisch: Kern- und Elektrondynamik sind stark gekoppelt. Die Kerndynamik findet auf Kreisbahnen statt. Mit der Rotation der Kerndichte um den Winkel φ geht eine Rotation der Elektron-Dichte einher.
Die diabatische Bewegung entspricht der Dynamik durch die konische Durchschneidung und die adiabatische Bewegung der Dynamik auf der unteren Potentialfläche. Welche der beiden Dynamiken stattfindet, wird durch die Wahl der Anfangsbedingung bestimmt. Der wesentliche Unterschied zwischen den beiden Startzuständen ist dabei die Lage des Knotens im elektronischen Anteil der Wellenfunktion. In den diabatischen Bewegungen bleibt z.B. der pₓ -artige Charakter der elektronischen Wellenfunktion konstant, wohingegen sich bei der adiabatischen Dynamik der Charakter mit der Kernbewegung ändert. Die Zeitersparnis durch die Verwendung des MCTDH-Ansatzes im Vergleich zur Split-Operator-Methode liegt etwa bei einem Faktor 5.
Das letzte Kapitel widmet sich der mikroskopischen Beschreibung von Exziton-Exziton-
Annihilierung (EEA). Dabei werden numerische Lösungen der aus einem mikro-
skopischen Modell hergeleiteten Ratengleichungen mit Messungen ( transienter Absorption) verglichen. Es wurden zwei Systeme untersucht: ein Squarain-basiertes Heteropolymer (SQA-SQB)ₙ und ein [2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenvinylen]-Polymer, auch bekannt als MEH-PPV. In beiden Fällen gelang die systematische Parameterbestimmung mit Hilfe einer Aufteilung in lokalisierte Subsysteme. Diese Subsysteme werden einzeln gewichtet und anschließend aufsummiert, wobei die Gewichte optimiert werden können. Aus den so erhaltenen Parametern ergibt sich für beide Systeme ein ähnliches Bild:
• Durch ultraschnelle Lokalisierung der Anregung im fs-Bereich auf kleinere Aggregateinheiten bilden sich voneinander getrennte Subsysteme.
• Die in den Subsystemen lokalisierten Exzitonen können sich nur innerhalb dieser Bereiche frei bewegen. Es ist ausreichend, direkt benachbarte Mono-, Bi-, Tri- und Tetra-Exzitonen in bis zu zwei Dimensionen zu berücksichtigen.
• Auf einer fs-Zeitskala annihilieren direkt benachbarte Exzitonen.
• Im MEH-PPV ergibt sich der Signalzerfall im fs-Bereich als Mittelwert aus einer schnellen (zwischen Ketten) und einer langsamen (innerhalb von Ketten) Annihilierung.
• Im ps- bis ns-Bereich wird sowohl durch Diffusion vermittelte Annihilierung, also auch der Zerfall der ersten angeregten Zustände bedeutsam.
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.
In the context of this dissertation very long ranged exciton diffusion lengths (LD) were simulated for perylene-based materials under ideal conditions. This leads to the conclusion that the short LD values in existing materials result from an extrinsic and intrinsic immobilization. The latter, which is a specific material property, is based on a relaxation of the exciton into self-trapping states. An in-depth understanding of the atomistic processes defining self-trapping is essential to developing materials with long LD in the future, in which intrinsic immobilization is prevented. For the development of such a mechanistic understanding it is crucial that a clear relationship between molecular structure and LD is available. This is given by single crystals of diindeno perylene (DIP) and α-perylene tetracarboxylic anhydride (α-PTCDA). An extraordinary large LD of 90 nm was measured for the first one, while the latter possesses only 22 nm. Part of this thesis was to deliver reasons for this discrepancy. Only self-trapping comes into question to explain the different LD values. One reason for the different self-trapping in DIP and α-PTCDA could lie in the electronic structure. However, it was possible to demonstrate that a wide range of perylene-based materials possess no significant differences in their electronic structures. Consequently, such differences can be neglected for the explanation of immobilization mechanisms for the exciton. A further possible explanation could be polarization effects in the crystal, which influences the electronic structure of perylene based materials differently. Especially their influence on charge transfer (CT) states, which are located above the optically bright Frenkel state, was in question because such states could be stabilized by a polarizable surrounding. A significant influence of polarization effects on all considered states were excluded by using a polarizable continuum model. Hence, the small LD values in α-PTCDA are an evidence for self-trapping, which produces a crystal structure built up by π-stacks, while the one of DIP is of herringbone type. Since polarization effects can be neglected, is the dimer only via steric restrictions influenced by the crystal. Hence, a method describing self-trapping has to consider such effects, so that a mechanical embedding QM/MM approach is sufficient. Now, potential energy surfaces were calculated, on which wave packet dynamics were subsequently performed. In this way, atomistic mechanisms for the immobilization of excitons were described for the first time in organic materials. Self-trapping was studied in crystals of α-PTCDA by potential energy surfaces, which map an intermolecular shift motion of the dimer in the crystal. An immobilization of excitons occurs within 500 fs, which results from an irreversible energy loss together with a local deformation of the crystal lattice. This prevents a further transport of the exciton. In the case of DIP, this immobilization does not proceed due to high barriers. These barriers result from the herringbone type packing motif in the DIP crystal. This discrepancy in the dynamics explains the different LD values in DIP and α-PTCDA. In a further example, an exciton immobilization was found in helical π-aggregates of perylene tetracarboxylic bisimide (PBI) molecules. Self-trapping is caused by a relaxation mechanism, in which the exciton is transferred by asymmetric vibrations of the aggregate from the bright to a dark Frenkel state within 200 fs, whereby the transition is mediated by a CT state. However, the CT state is almost non-populated during the whole mechanism so that its participation could not yet be proven experimentally. This entire procedure is solely possible in helical aggregates, because only for such structures is there a CT state located next to the bright Frenkel state. At the final Frenkel state a torsional motion around the π-stacking axis is possible so that the loss in energy and the local rearrangement of the aggregate structure occurs, which means a self-trapping of the exciton. This mechanism is in perfect agreement with all available experimental data. These insights allow the conclusion that in future materials for organic solar cells an irreversible and ultrafast deformation of aggregates after photo-absorption must be avoided. Only in this way long LD values can be achieved and exciton self-trapping can be prevented. However, small LD values are always predicted in helical aggregates of perylene-based materials, because exciton immobilization occurs already due to small molecular motions. For this reason such aggregates are inappropriate for the use in organic solar cells. Long LD values are expected for aggregate structures with long intermolecular shifts or molecules with bulky substituents.
This work focuses on theoretical approaches for predicting the valence and core excited states of aggregate systems. For the valence excitations, TD-HF and TD-DFT with different functionals have been tested at the Perylene bisimide (PBI) system. A simple character analysis method based on the calculated transition dipole moments is proposed. However, this method does not work for excited states without any transition dipole moment. Thus, we proposed a more general and more valid method based on a calculated CIS type wavefunction for the character analysis. Furthermore, a model Hamiltonian method is derived from a localized picture. The energies of the diabatic states and the corresponding coupling parameters were also determined on the basis of ab initio calculations. For the core excitation, three different methods were validated for C 1s-excited and ionized states if several small molecules. Also we tested the basis sets dependence of these core excited states. Based on those results, we chose the frozen core approximation method to evaluate the core excited states of NTCDA molecules. In order to explain the findings in the experiments, we developed an algorithm to evaluate the exciton coupling parameter where non-orthogonal MOs are used.
The spectroscopic properties of molecular aggregates have been investigated by means of quantum dynamical calculations. Thereby both linear and nonlinear spectroscopic techniques have been taken into account. For the simulation of absorption and CD-spectra, coupling effects were regarded as well as the relative orientation of the monomer units in order to determine the parameters by reproducing measured spectra. For a more detailled description, results from quantum chemical calculations have also been included. Furthermore, investigations on nonlinear spectroscopy of molecular dimers have been performed.