Refine
Has Fulltext
- yes (23)
Is part of the Bibliography
- yes (23)
Year of publication
Document Type
- Doctoral Thesis (23)
Keywords
- Theoretische Chemie (5)
- Exziton (4)
- QM/MM (4)
- Computational Chemistry (3)
- Quantenchemie (3)
- organische Halbleiter (3)
- Benchmark (2)
- Berechnung (2)
- Biradikal (2)
- CAAC (2)
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 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.
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. ...
Im Rahmen dieser Arbeit konnten nasschemische Synthesen fü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üssig zur Synthese von B2I4 gelang die Darstellung aller vier Dibortetrahalogenide mittels einfacher Schlenktechnik basierend auf der Synthese von B2Br4 durch Nö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ührte zunächst zur Ausbildung des Bisadduktes 53, allerdings konnte in Lösung eine Umlagerung beobachtet werden, deren Verlauf 11B-NMR-spektroskopisch verfolgt wurde, jedoch nicht final aufgeklärt werden konnte. Durch die Zugabe eines Unterschusses der Lewis-Basen IDipp bzw. PCy3 sollten zunächst Monoaddukte von B2Cl4 dargestellt werden, deren Umsetzung mit einer weiteren Lewis-Base die Synthese asymmetrischer Lewis-Basen-Addukte von B2Cl4 ermöglichen sollte. Die sp2-sp3-Diborane 56 und 57 konnten bei tiefen Temperaturen 11B-NMR-spektroskopisch nachgewiesen werden, allerdings führte eine Abfangreaktion mit diversen Lewis-Basen nicht zu den gewünschten asymmetrischen Addukten. Bei Raumtemperatur konnte eine Folgereaktion von 56 zur Chlorid-verbrü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ät von B2Cl4 gegenüber ungesä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ächst ü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ü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ür eine weitere Reduktion zur B-B-Doppelbindung nötig wä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ührte zu den chloridverbrü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ür die CAAC-stabilisierten Verbindungen auch röntgenkristallografisch nachgewiesen werden konnten. Theoretische Untersuchungen bestätigten die Relevanz des verbrückenden Chloratoms zur Stabilisierung dieser Verbindungen. Für die Stammverbindung der Diborylalkene (59 (R = H)) konnte bei der Umsetzung mit MeCAAC eine unlösliche Verbindung erhalten werden, deren Struktur als Bisaddukt 82 mittels NMR-spektroskopischen Untersuchungen im Festkörper und durch Verbrennungsanalyse bestätigt werden konnte.
Die Reduktion der CAAC-stabilisierten Diborylalkene 79 und 80 in Gegenwart von MeCAAC führte zu den captodativ-stabilisierten Diborylradikalen 83 und 84, deren Strukturanalyse eine orthogonale Anordnung der C2-Brücke zur B(CAAC)-Einheit offenlegt. Ausführliche EPR-spektroskopische Untersuchungen bei variabler Temperatur und theoretische Berechnungen bestä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ülorbitale (SOMO) bestätigen eine Delokalisierung der Spindichte über die NCAAC-CCAAC-B-Einheit. Der Austausch der verbrückenden Einheit und die somit einhergehende Verringerung der Sterik führt zu einer Planarisierung des Moleküls im Festkörper (87). Theoretische Untersuchungen und die Auswertung der strukturellen Parameter ergeben eine Delokalisierung der Elektronendichte ü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ä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ührt, die in der Ausbildung der Diborylene 99−104 resultierten. Die Darstellung einer para-Phenylen-verbrückten Donor-Akzeptor-Verbindung (D: Borylen, A: BMes2) gelang nicht.
Im Rahmen dieser Arbeit wurden die elektronischen Eigenschaften von Graphen auf Metalloberflächen mittels Rastertunnelmikroskopie und Quasiteilcheninterferenz (englisch quasiparticle interference, QPI)-Messungen untersucht. Durch das Verwenden schwerer Substrate sollte die Spin-Bahn-Wechselwirkung des Graphen verstärkt werden und damit eine Bandlücke am K-Punkt der Bandstruktur mittels QPI beobachtet werden. Um das Messen von QPI auf Graphen zu testen, wurde auf der Oberfläche eines SiC(0001)-Kristalls durch Erhitzen Graphen erzeugt und mit dem Rastertunnelmikroskop untersucht. Dieses System wurde schon ausfü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ä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ü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ätzlich eine Temperaturabhängigkeit nachgewiesen werden konnte. Beim Versuch, Gadolinium zu interkalieren, wurden zwei verschieden Oberflächenstrukturen beobachtet, die auf eine unterschiedlich Anordnung bzw. Menge des interkalierten Gadoliniums zurückzuführen sein kö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ängige Austauschaufspaltung reproduziert werden. Darüber hinaus konnten sechs verschieden magnetische Domänen beobachtet werden. Zusätzlich sind auf der Oberfläche magnetische Streifen auszumachen, die möglicherweise auf einer Spinspirale basieren. Als Grundlage für die mögliche zukünftige Erzeugung Graphen-artiger Molekülgitter wurde das Wachstum von H-TBTQ und Me-TBTQ auf Ag(111) untersucht. Die Moleküle richten sich dabei nach der Oberflächenstruktur des Silber aus und bilden längliche Inseln, deren Kanten in drei Vorzugsrichtungen verlaufen. Auf H-TBTQ wurde zudem eine zweite, Windmühlen-artige Ausrichtung der Moleküle auf der Oberfläche beobachtet. Auf den mit den Molekülen bedeckten Stellen der Oberfläche wurde eine Verschiebung des Ag-Oberflächenzustands beobachtet, die mit einem Ladungstransfer vom Ag(111)-Substrat auf die TBTQ-Moleküle zu erklären sein könnte.
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.
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.
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.
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.
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.