540 Chemie und zugeordnete Wissenschaften
Refine
Has Fulltext
- yes (7)
Is part of the Bibliography
- yes (7)
Document Type
- Doctoral Thesis (7)
Language
- English (7)
Keywords
- Theoretische Chemie (7) (remove)
Institute
- Institut für Physikalische und Theoretische Chemie (7) (remove)
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.
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.
Theoretical Investigations on the Interactions of Small Compounds with their Molecular Environments
(2015)
In the first part of this work, a combination of theoretical methods for the rational design of covalent inhibitor is presented. Starting from the crystal structure of the covalent complex of a lead compound, quantum mechanical and QM/MM calculations were used to derive the exact geometry of the preceeding non-covalent enzyme inhibitor complex. The geometry of the latter mainly determines the reactivity of the inhibitor against its target enzyme concerning the formation of the covalent bond towards an active site residue. Therefore, this geometry was used as starting point for the optimization of the substitution pattern of the inhibitor such as to increase its binding affinity without loosing its ability to covalently bind to the target protein. The optimization of the chemical structure was supported by using docking procedures, which are best suited to estimate binding affinities that arise from the introduced changes. A screening of the novel substitution patterns resulted in a first generation of model compounds which were further tested for their reactivity against the target. Dynamic simulations on the novel compounds revealed that the orientation that compounds adopt within the active site are such that a covalent interaction with the enzyme is no longer possible. Hence, the chemical structure was further modified, including not only changes in the substituents but also within the core of the molecule. Docking experiments have been conducted to assure sufficiently high binding affinities and to obtain the most favored binding poses. Those have then again been used for dynamic simulations which resulted in structures, for which the bond formation process appeared feasible. A final series of QM/MM calculations considering various protonation states was computed to estimate the reaction energies for the covalent attachment of the inhibitor to the enzyme. The theoretical results indicate a reasonable high inhibition potency of the novel compounds.
The second part concentrates on the environmental influences on the electron density of an inhibitor molecule. Therefore, a vinylsulfone-based model compound was selected for which an experimental crystal structure for the pure compound as well as a theoretically determined enzyme-inhibitor complex have been available. To provide reference data for the larger systems, the conformational space of the isolated molecule was screened for favorable geometries which were later compared to those within the crystal and protein surrounding. The geometry of the crystal structure could readily be taken from the experimental data whereas calculations on the protein complex revealed four potential non-covalent complexes exhibiting different arrangements of the molecule within the active site of the protein as well as two possible protonation states of the catalytic dyad. Hence, all four protein complexes have been compared to the crystal structure of the molecule as well as against the more favorable geometries of the isolated molecule being determined within vacuum or aqueous surrounding. Whereas the molecule itself was found to adopt comparable geometries within all investigated environments, the interactions pattern between the crystal surrounding and the protein differed largely from each other. The favorable formation of dimers within the crystal has a strong stabilizing effect and explains the extraordinarily good quality of the crystal. Within the protein however, repulsive forces have been found between the protein and the inhibitor. The origin of the repulsion could be traced back to effect of on of the substituents to the vinyl scaffold. The difference in the chemical structure in comparison to a well known inhibitor might also explain the experimentally found loss of activity for the model compound in comparison to K11777.
The present work presents investigations on energy and charge transport properties in organic crystals. Chapter 4 treats exciton transport in anthracene, which is an example for weakly coupled π-systems. The electronic coupling parameter is evaluated by the monomer transition density approach. With these and the reorganization energy hopping rates are calculated in the framework of the Marcus theory. Together with the knowledge of the crystal structure, these allow us to calculate the experimental accessible exciton diffusion lengths, whose isotropic part fits nicely within the scattering of experimental values found in the literature. Furthermore, the anisotropy of the exciton diffusion lengths is reproduced qualitatively and quantitatively correct. This chapter also contains studies about electron and hole transport in both polymorphs (α and β) of perylene. Reorganization energies as well as diffusion coefficients for both crystal structures and types of charge transport were calculated. The best transport is hole transport in β-perylene, but it is strongly isotropic. The preferred transport direction is along the b-axis of the unit cell with couplings of greater than 100 meV. However, there is no transport along the c-axis. The diffusion constant in b-direction is bigger by two orders of magnitude than in c-direction (62.7•10-6 m2/s vs. 0.4•10-6 m2/s). Charge transport is calculated to be strongly anisotropic for holes as well as electrons in both modifications. To verify these results experimental electron mobilities have been compared to the simulations. Good agreement was found with errors of less than 27%. As it was shown above, the calculation and measurement of transport properties between weakly coupled systems is possible. However, it is difficult to exactly determine the quality of the electronic coupling. For this reason a collaboration about strongly interacting π-systems was started between us and the research group of Prof. Ingo Fischer. There, [2.2]paracyclophanes and its derivates were investigated to show how hydroxyl substitution influences absorption properties. Overall, a combination of SCS-MP2 and SCS-CC2 performs best to address the description of geometric and electronic structures for both ground and excited states of these model systems as well as their parent compounds benzene and phenol. Only [2.2]paracyclophane shows a double minimum potential regarding a twist and shift motion between the benzene/phenol subunits towards each other. All other systems are less flexible due to their substitution pattern. Almost all [2.2]paracyclophanes display minor changes in their geometric structure upon excitation to the S1 state: The inter-ring distance shortens, but qualitatively they keep their shift and twist characteristics, although the extent of these deformations diminishes. The exception is p-DHPC, which turns from a shifted ground state structure into a twisted excited state structure. Consequently, the intensity of the 0-0 transition cannot be observed experimentally due to small Franck-Condon factors and impurities of o-DHPC. In the present thesis, the structures and their changes due to excitation are explained by electrostatic potentials as well as antibonding (bonding) HOMO (LUMO) orbitals. Adiabatic excitation energies have been corrected by ZPEs and result in accuracies with errors smaller than 0.1 eV. Note that corrections on the B3LYP level worsen the results and one has to apply SCS-CC2 to achieve this accuracy. These calculations allow an interpretation of the experimental [1+1]REMPI spectra. Band progressions of the twist, shift and breathing of the [2.2]paracyclophane skeleton vibrations have been identified and show good agreement to the experiment. This work shows that the substitution pattern in [2.2]paracyclophanes can have a significant impact on spectroscopic properties. Because these properties are directly linked to the transport properties of these materials, the hereby gained insight can be used to design materials with customized transport properties. It was shown that the SCS-CC2 method is very appropriate to predict the interaction between the π-systems
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.
Theory and simulation of ultrafast autodetachment dynamics and nonradiative relaxation in molecules
(2024)
In this thesis, theoretical approaches for the simulation of electron detachment processes in molecules following vibrational or electronic excitation are developed and applied. These approaches are based on the quantum-classical surface-hopping methodology, in which nuclear motion is treated classically as an ensemble of trajectories in the potential of quantum-mechanically described electronic degrees of freedom.
The visualization of energy functions is based on the possibility of separating different degrees of freedom. The most important one is the Born-Oppenheimer-approximation, which separates nucleus and electron movements. This allows the illustration of the potential energy as a function of the nuclei coordinates. Minima of the surface correspond to stable points like isomers or conformers. They are important for predicting the stability or thermodynamical of a system. Stationary points of first order correspond to transition points. They describe phase transitions, chemical reaction, or conformational changes. Furthermore, the partition function connects the potential hypersurface to the free energy of the system. The aim of the present work is the development and application of new approaches for the efficient exploration of multidimensional hypersurfaces. Initially, the Conformational Analysis and Search Tool (CAST) program was developed to create a basis for the new methods and algorithms. The development of CAST in object oriented C++ included, among other things, the implementation of a force field, different interfaces to external programs, analysis tools, and optimization libraries. Descriptions of an energy landscape require knowledge about the most stable minima. The Gradient Only Tabu Search (GOTS) has been shown to be very efficient in the optimization of mathematical test functions. Therefore, GOTS was taken as a starting point. Tabu-Search is based on the steepest descent - modest ascent strategy. The steepest descent is used for finding local minima, while the modest ascent is taken for leaving a minimum quickly. Furthermore, Tabu-Search is combined with an adaptive memory design to avoid cycling or returning. The highly accurate exploration of the phase space by Tabu-Search is often too expensive for complex optimization problems. Therefore, an algorithm for diversification of the search is required. After exploration of the proximity of the search space, the algorithm would guide the search to new and hopefully promising parts of the phase space. First application of GOTS to conformational search revealed weaknesses in the diversification search and the modest ascent part. On the one hand, the original methodology for diversification is insufficiently diverse. The algorithm is considerably improved by combining the more local GOTS with the wider searching Basin Hopping (BH) approach. The second weak point is a too inaccurate and inefficient modest ascent strategy. Analysis of common transition state search algorithms lead to the adaption of the Dimer-method to the Tabu-Search approach. The Dimer-method only requires the first derivatives for locating the closest transition state. For conformational search, dihedral angles are usually the most flexible degrees of freedom. Therefore, only those are used in the Dimer-method for leaving a local minimum. Furthermore, the exact localization of the reaction pathway and the transition state is not necessary as the local minimum position should only be departed as fast as possible. This allows for larger step sizes during the Dimer-search. In the following optimization step, all coordinates are relaxed to remove possible strains in the system. The new Tabu-Search method with Dimer-search delivers more and improved minima. Furthermore, the approach is faster for larger systems. For a system with approximately 1200 atoms, an acceleration of 40 was measured. The new approach was compared to Molecular Dynamics with optimization (MD), Simulated Annealing (SA), and BH with the help of conformational search problems of bio-organic systems. In all cases, a better performance was found. A comparison to the Monte Carlo Multiple Minima/Low Mode Sampling (MCMM/LM) method proved the outstanding performance of the new Tabu-Search approach. The solvation of the chignolin protein further revealed the possibility of uncovering discrepancies between the employed theoretical model and the experimental starting structure. Ligand optimization for improvement of x-ray structures was one further new application field. Besides the global optimization, the search for transition states and reaction pathways is also of paramount importance. These points describe different transitions of stable states. Therefore, a new approach for the exploration of such cases was developed. The new approach is based on a global minimization of a hyperplane being perpendicular to the reaction coordinate. Minima of this reduced phase space belong to traces of transition states between reactant and product states on the unchanged hypersurface. Optimization to the closest transition state using the Dimer-method delivers paths lying between the initial and the final state. An iterative approach finally yields complex reaction pathways with many intermediate local minima. The PathOpt algorithm was tested by means of rearrangements of argon clusters showing very promising results.