Refine
Has Fulltext
- yes (16)
Is part of the Bibliography
- yes (16)
Year of publication
Document Type
- Doctoral Thesis (14)
- Journal article (1)
- Preprint (1)
Language
- English (16) (remove)
Keywords
- Exziton (16) (remove)
Institute
EU-Project number / Contract (GA) number
- 614623 (2)
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.
We present a theoretical study on exciton–exciton annihilation (EEA) in a molecular dimer. This process is monitored using a fifth-order coherent two-dimensional (2D) spectroscopy as was recently proposed by Dostál et al. [Nat. Commun. 9, 2466 (2018)]. Using an electronic three-level system for each monomer, we analyze the different paths which contribute to the 2D spectrum. The spectrum is determined by two entangled relaxation processes, namely, the EEA and the direct relaxation of higher lying excited states. It is shown that the change of the spectrum as a function of a pulse delay can be linked directly to the presence of the EEA process.
We present a theoretical study on exciton–exciton annihilation (EEA) in a molecular dimer. This process is monitored using a fifth-order coherent two-dimensional (2D) spectroscopy as was recently proposed by Dostál et al. [Nat. Commun. 9, 2466 (2018)]. Using an electronic three-level system for each monomer, we analyze the different paths which contribute to the 2D spectrum. The spectrum is determined by two entangled relaxation processes, namely, the EEA and the direct relaxation of higher lying excited states. It is shown that the change of the spectrum as a function of a pulse delay can be linked directly to the presence of the EEA process.
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.
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.
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
The chirality of the interlocked bay-arylated perylene motif is investigated upon its material prospect and the enhancement of its chiroptical response to the NIR spectral region. A considerable molecular library of inherently chiral perylene bisimides (PBIs) was utilized as acceptors in organic solar cells to provide decent device performances and insights into the structure-property relationship of PBI materials within a polymer blend. For the first time in the family of core-twisted PBIs, the effects of enantiopurity on the device performance was thoroughly investigated. The extraordinary structural sensitivity of CD spectroscopy served as crucial analytical tool to bridge the highly challenging gap between molecular properties and device analytics by proving the excitonic chirality of a helical PBI dimer. The chirality of this perylene motif could be further enhanced on a molecular level by both the expansion and the enhanced twisting of the π-scaffold to achieve a desirable strong chiroptical NIR response introducing a new family of twisted QBI-based nanoribbons. These achievements could be substantially further developed by expanding this molecular concept to a supramolecular level. The geometrically demanding supramolecular arrangement necessary for the efficient excitonic coupling was carefully encoded into the molecular design. Accordingly, the QBIs could form the first J-type aggregate constituting a fourfold-stranded superhelix of a rylene bisimide with strong excitonic chirality. Therefore, this thesis has highlighted the mutual corroboration of experimental and theoretical data from the molecular to the supramolecular level. It has demonstrated that for rylene bisimide dyes, the excitonic contribution to the overall chiroptical response can be designed and rationalized. This can help to pave the way for new organic functional materials to be used for
chiral sensing or chiral organic light-emitting devices.
The goal of this thesis was the development and application of higher-order spectroscopic techniques. In contrast to ordinary pump–probe (PP) and two-dimensional (2D) spectroscopy, higher-order coherently detected spectroscopic methods measure a polarization that has an order of nonlinearity higher than three. The key idea of the techniques in this thesis is to isolate the higher-order signals from the lower-order signals either by their excitation frequency or by their excitation intensity dependence. Due to the increased number of interactions in higher-order spectroscopy, highly excited states can be probed. For excitonic systems such as aggregates and polymers, the fifth-order signal allows one to directly measure exciton–exciton annihilation (EEA). In polymers and aggregates, the exciton transport is not connected to a change of the absorption and can therefore not be investigated with conventional third-order techniques. In contrast, EEA can be used as a probe to study exciton diffusion in these isonergetic systems. As a part of this thesis, anisotropy in fifth-order 2D spectroscopy was investigated and was used to study geometric properties in polymers.
In 2D spectroscopy, the multi-quantum signals are separated from each other by their spectral position along the excitation axis. This concept can be extended systematically to higher signals. Another approach to isolate multi-quantum signals in PP spectroscopy utilizes the excitation intensity. The PP signal is measured at specific excitation intensities and linear combinations of these measurements result in different signal contributions. However, these signals do not correspond to clean nonlinear signals because the higher-order signals contaminate the lower-order multi-quantum signals. In this thesis, a correction protocol was derived that uses the isolated multiquantum signals, both from 2D spectroscopy and from PP spectroscopy, to remove the contamination of higher-order signals resulting in clean nonlinear signals. Using the correction on the third-order signal allows one to obtain annihilation-free signals at high excitation intensities, i.e., with high signal-to-noise ratio. Isolation and correction in PP and 2D spectroscopy were directly compared by measuring the clean third-order signals of squaraine oligomers at high excitation intensities. Furthermore, higher-order PP spectroscopy was used to isolate up to the 13th nonlinear order of squaraine polymers.
The demonstrated spectroscopic techniques represent general procedures to isolate clean signals in terms of perturbation theory. The technique of higher-order PP spectroscopy needs only small modifications of ordinary PP setups which opens the field of higher-order spectroscopy to the broad scientific community. The technique to obtain clean nonlinear signals allows one to systematically increase the number of interacting (quasi)particles in a system and to characterize their interaction energies and dynamics.
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 present thesis demonstrates the importance of the solid state packing of dipolar merocyanine dyes with regard to charge transport and exciton coupling.
Due to the charge transport theory for disordered materials, it is expected that high ground state dipole moments in amorphous thin films lead to low mobility values due to a broadening of the density of states. However, due to their inherent dipolarity, merocyanine dyes usually align in antiparallel dimers in an ordered fashion. The examination of twenty different molecules with ground state dipole moments up to 15.0 D shows that by a high dipolarity and well-defined sterics, the molecules pack in a highly regular two-dimensional brickwork-type structure, which is beneficial for hole transport. Utilization of these molecules for organic thin-film transistors (OTFTs) leads to hole mobility values up to 0.21 cm²/Vs. By fabrication of single crystal field-effect transistors (SCFETs) for the derivative showing the highest mobility values in OTFTs, even hole mobilities up to 2.34 cm²/Vs are achieved. Hence, merocyanine based transistors show hole mobility values comparable to those of conventional p-type organic semiconductors and therefore high ground state dipole moments are not necessarily disadvantageous regarding high mobility applications.
By examination of a different series of ten merocyanine dyes with the same chromophore backbone but different donor substituents, it is demonstrated that the size of the donor has a significant influence on the optical properties of thin films. For small and rigid donor substituents, a hypsochromic shift of the absorption compared to the monomer absorption in solution is observed due to the card stack like packing of the molecules in the solid state. By utilization of sterical demanding or flexible donor substituents, a zig-zag type packing is observed, leading to a bathochromical shift of the absorption. These packing motifs and spectral shifts with an offset of 0.93 eV of the H- and J-bands comply with the archetype examples of H- and J-aggregates from Kasha’s exciton theory.