Refine
Has Fulltext
- yes (12)
Is part of the Bibliography
- yes (12)
Document Type
- Doctoral Thesis (12)
Language
- English (12) (remove)
Keywords
- Femtosekundenspektroskopie (4)
- Ultrakurzzeitspektroskopie (4)
- Ultraschnelle Photochemie (3)
- ultrafast photochemistry (3)
- Pump-Probe-Technik (2)
- QM/MM (2)
- Verbindungen (2)
- theoretical chemistry (2)
- Benchmark (1)
- Biradical (1)
In this thesis, we apply the information-theoretic approach in the context of quantum dynamics and wave packet motion: Information-theoretic measures are calculated from position and momentum densities, which are obtained from time-dependent quantum wave functions. The aim of this thesis is to benchmark, analyze and interpret these quantities and relate their features to the wave packet dynamics. Firstly, this is done for the harmonic oscillator (HO) with and without static disorder. In the unperturbed HO, the analytical study of coherent and squeezed states reveals time-dependent entropy expressions related to the localization of the wave function. In the disordered HO, entropies from classical and quantum dynamics are compared for short and long times. In the quantum case, imprints of wave packet revivals are found in the entropy. Then, the energy dependence of the entropy for very long times is discussed. Secondly, this is donefor correlated electron-nuclear motion. Here, entropies derived from the total, electronic and nuclear density, respectively, are calculated in position and momentum space for weak and strong adiabatic electronic coupling. The correlation between electron and nucleus is investigated using different correlation measures, where some of these functions are sensitive to the nodal structure of the wave function. An analytic ansatz to interpret the information-theoretical quantities is applied as well.
This work aims at elucidating chemical processes involving homogeneous catalysis and photo–physical relaxation of excited molecules in the solid state. Furthermore, compounds with supposedly small singlet–triplet gaps and therefore biradicaloid character are investigated with respect to their electro–chemical behavior. The work on hydroboration catalysis via a reduced 9,10–diboraanthracene (DBA) was preformed in collaboration with the Wagner group in Frankfurt, more specifically Dr. Sven Prey, who performed all laboratory experiments. The investigation of delayed luminescence properties in arylboronic esters in their solid state was conducted in collaboration with the Marder group in Würzburg. The author of this work took part in the synthesis of the investigated compounds while being supervised by Dr. Zhu Wu. The final project was a collaboration with the group of Anukul Jana from Hyderabad, India who provided the experimental data.
Nowadays, computational-aided investigations become an essential part in the chemical, biochemical or pharmaceutical research. With increasing computing power, the calculation of larger biological systems becomes feasible. In this work molecular mechanical (MM) and quantum mechanical approaches (QM) and the combination of both (QM/MM) have been applied to study several questions which arose from different working groups. Thus, this work comprises eight different subjects which deals with chemical reactions or proton transfer in enzymes, conformational changes of ligands or proteins and verification of experimental data.
This work firstly deals with reaction mechanisms of aromatic inhibitors of cysteine proteases which can be found in many organisms. These enzymes are responsible for various cancer or diseases as for example Human African Trypanosomiasis (HAT) or the Chagas disease. Aromatic SNAr-type electrophiles might offer a new possibility to covalently modify these proteases. Quantum mechanical calculations have been performed to gain insights into the energetics and possible mechanisms.
The next chapter also deals with Trypanosomiasis but the focus was set on a different enzyme. The particularity of Trypanosomiasis is the thiol metabolism which can also be modified by covalent inhibitors. In this context, the wild type and point mutations of the enzyme tryparedoxin have been investigated via molecular dynamic (MD) simulations to examine the influence of specific amino acids in regard to the inhibitor. Experimental data showed that a dimerization of the enzyme occurs if the inhibitor is present. Simulations revealed that the stability of the dimer decreases in absence of the inhibitor and thus confirms these experiments.
Further investigations concerning cysteine proteases such as cruzain and rhodesain have been conducted with respect to experimental kinetic data of covalent vinylsulfone inhibitors. Several approaches such as QM or QM/MM calculations and docking, MD or MMPBSA/MMGBSA simulations have been applied to reproduce these data. The utilization of force field approaches resulted in a qualitatively accurate prediction.
The kinase AKT is involved in a range of diseases and plays an important role in the formation of cancer. Novel covalent-allosteric inhibitors have been developed and crystallized in complex with AKT. It was shown that depending on the inhibitor a different cysteine residue is modified. To investigate these differences in covalent modification computational simulations have been applied.
Enoyl-(acyl carrier) (ENR) proteins are essential in the last step of the fatty acid biosynthesis II (FAS) and represent a good target for inhibition. The diphenylether inhibitor SKTS1 which was originally designed to target the ENR’s of Staphylococcus aureus was also crystallized in InhA, the ENR of Mycobacterium tuberculosis (TB). Crystal structures indicate a change of the inhibitor's tautomeric form. This subject was investigated via MD simulations. Results of these simulations confirmed the tautomerization of the inhibitor.
This work also deals with the development of a covalent inhibitor originating from a non-covalent ligand. The target FadA5 is an essential enzyme for the degradation of steroids in TB and is responsible for chronic tuberculosis. This enzyme was crystallized in complex with a non-covalent ligand which served as starting point for this study. Computations on QM or QM/MM level and docking and MD simulations have been applied to evaluate potential candidates.
The next chapter focuses on the modification of the product spectrum of Bacillus megaterium levansucrase, a polymerase which catalyzes the biosynthesis of fructans. The covalent modification of the wild type or mutants of the enzyme lead to an accumulation of oligosaccharides but also to polymers with higher polymerization degree. To understand these changes in product spectra MD simulations have been performed.
Finally, the proton transfer in catalytic cysteine histidine dyads was investigated. The focus was set on the influence of the relaxation of the protein environment to the reaction. Calculations of the enzymes FadA5 and rhodesain revealed that the preferred protonation state of the dyade depends on the protein environment and has an impact on the reaction barrier. Furthermore, the adaptation of the environment to a fixed protonation state was analyzed via MD simulations.
In the context of quantum mechanical calculations, the properties of non-adiabatic coupling in a small system, the Shin-Metiu model, is investigated.
The transition from adiabatic to non-adiabatic dynamics is elucidated in modifying the electron-nuclear interaction. This allows the comparison of weakly correlated electron-nuclear motion with the case where the strong correlations determine the dynamics.
The studies of the model are extended to include spectroscopical transitions being present in two-dimensional and degenerate four-wave mixing spectroscopy.
Furthermore, the quantum and classical time-evolution of the coupled motion in the complete electron-nuclear phase space is compared for the two coupling cases.
Additionally, the numerically exact electron flux within the weak coupling case is compared to the Born-Oppenheimer treatment.
In the last part of the thesis, the model is extended to two dimensions.
The system then possesses potential energy surfaces which exhibit a typical 'Mexican hat'-like structure and a conical intersection in the adiabatic representation.
Thus, it is possible to map properties of the system onto a vibronic coupling (Jahn-Teller) hamiltonian. Exact wave-packet propagations as well as nuclear wave-packet dynamics in the adiabatic and diabatic representation are performed.
The invention of laser pulse shapers allowed for various quantum control experiments, where a chemical reaction is guided by specifically tailored laser pulses. However, despite of the prominent role of the liquid phase in chemistry, no successful attempt for controlling the selectivity of a bond-fission reaction has yet been reported in this state of matter. Promising candidates for such an experiment are C$_{\infty\mathrm{v}}$-symmetric trihalide anions with two different chemical bonds like $\ce{I2Cl-}$, because these molecules notionally offer the most simplest selectivity-control scenario of breaking either the one or the other bond and they are expected to dissociate under ultraviolet (UV) irradiation like it is known for the most-studied trihalide $\ce{I3-}$.
In order to investigate in this thesis the possibility that the dissociation reaction of such trihalides branches into two different photofragments, the ultrafast photodissociation dynamics of $\ce{I3-}$, $\ce{Br3-}$, $\ce{IBr2-}$ and $\ce{ICl2-}$ (point group D$_{\infty\mathrm{h}}$) as well as of $\ce{I2Br-}$ and $\ce{I2Cl-}$ (point group C$_{\infty\mathrm{v}}$) in dichloromethane solution were measured with broadband transient absorption spectroscopy in magic-angle configuration. The identification of the reaction pathway(s) relies on vibrational wavepacket oscillations, which survive the dissociation process and therefore carry not only informations about the reactant trihalides but also about the fragment dihalides.
These characteristic vibrational wavenumbers were extracted from the measured transient absorption spectra by globally fitting the population dynamics together with the wavepacket dynamics. Until recently, such a combined model function was not available in the well-established fitting tool Glotaran. This made it inevitable to develop a custom implementation of the underlying variable-projection fitting algorithm, for which the computer-algebra software Mathematica was chosen. Mathematica's sophisticated built-in functions allow not only for a high flexibility in constructing arbitrary model functions, but also offer the possibility to automatically calculate the derivative(s) of a model function. This allows the fitting procedure to use the exact Jacobian matrix instead of approximating it with the finite difference method.
Against the expectation, only one of the two thinkable photodissociation channels was found for each of the investigated C$_{\infty\mathrm{v}}$ trihalides. Since the photofragments recombine, their absorption signal as well as the reactant ground state bleach recover. This happens in a biexponential manner, which in the case of $\ce{I3-}$ was interpreted by Ruhman and coworkers with the direct formation of a neutral dihalogen fragment $\ce{I2}$ beside the negatively charged dihalide fragment $\ce{I2-}$. In this thesis, such a direct reaction channel was not found and instead the fast component of the biexponential decay is explained with vibrational excess energy mediating the recombination-preceding electron transfer process $\ce{I2- + I -> I2 + I-}$, while the slow component is attributed to cooled-down fragments.
In addition to the trihalide experiments, the possibility of a magic-angle configuration for polarization-shaping control experiments was theoretically investigated in this thesis by deriving magic-angle conditions for the third-order electric-dipole response signal of arbitrarily polarized laser pulses. Furthermore, the subtleties of anisotropy signals violating the well-known range of \numrange{-0.2}{0.4} were studied.
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.
Time-resolved spectroscopy allows for analyzing light-induced energy conversion and
chromophore–chromophore interactions in molecular systems, which is a prerequisite in
the design of new materials and for improving the efficiency of opto-electronic devices.
To elucidate photo-induced dynamics of complex molecular systems, transient absorption
(TA) and coherent two-dimensional (2D) spectroscopy were employed and combined
with additional experimental techniques, theoretical approaches, and simulation models
in this work.
A systematic series of merocyanines, synthetically varied in the number of chromophores
and subsitution pattern, attached to a benzene unit was investigated in cooperation with
the group of Prof. Dr. Frank Würthner at the University of Würzburg. The global analysis
of several TA experiments, and additional coherent 2D spectroscopy experiments, provided
the basis to elaborate a relaxation scheme which was applicable for all merocyanine
systems under investigation. This relaxation scheme is based on a double minimum on the
excited-state potential energy surface. One of these minima is assigned to an intramolecular
charge-transfer state which is stabilized in the bis- and tris-chromophoric dyes by
chromphore–chromophore interactions, resulting in an increase in excited-state lifetime.
Electro-optical absorption and density functional theory (DFT) calculations revealed a
preferential chromophore orientation which compensates most of the dipole moment of
the individual chromophores. Based on this structural assignment the conformationdependent
exciton energy splitting was calculated. The linear absorption spectra of the
multi-chromophoric merocyanines could be described by a combination of monomeric and
excitonic spectra.
Subsequently, a structurally complex polymeric squaraine dye was studied in collaboration
with the research groups of Prof. Dr. Christoph Lambert and Prof. Dr. Roland Mitric
at the University of Würzburg. This polymer consists of a superposition of zigzag and
helix structures depending on the solvent. High-level DFT calculations confirmed the previous
assignment that zigzag and helix structures can be treated as J- and H-aggregates,
respectively. TA experiments revealed that in dependence on the solvent as well as the
excitation energy, ultrafast energy transfer within the squaraine polymer proceeds from
initially excited helix segments to zigzag segments or vice versa. Additionally, 2D spectroscopy
confirmed the observed sub-picosecond dynamics. In contrast to other conjugated
polymers such as MEH-PPV, which is investigated in the last chapter, ultrafast
energy transfer in squaraine polymers is based on the matching of the density of states
between donor and acceptor segments due to the small reorganization energy in cyanine-like
chromophores.
Finally, the photo-induced dynamics of the aggregated phase of the conjugated polymer
MEH-PPV was investigated in cooperation with the group of Prof. Dr. Anna Köhler at the University of Bayreuth. Our collaborators had previously described the aggregation of MEH-PPV upon cooling by the formation of so-called HJ-aggregates based on exciton
theory. By TA measurements and by making use of an affiliated band analysis distinct
relaxation processes in the excited state and to the ground state were discriminated. By
employing 2D spectroscopy the energy transfer between different conjugated segments
within the aggregated polymer was resolved. The initial exciton relaxation within the
aggregated phase indicates a low exciton mobility, in contrast to the subsequent energy
transfer between different chromophores within several picoseconds.
This work contributes by its systematic study of structure-dependent relaxation dynamics
to the basic understanding of the structure-function relationship within complex
molecular systems. The investigated molecular classes display a high potential to increase
efficiencies of opto-electronic devices, e.g., organic solar cells, by the selective choice of
the molecular morphology.
Within the framework of this thesis, photolysis reactions in the liquid phase were investigated by means of ultrafast optical spectroscopy. Apart from molecular studies dealing with the highly spin-dependent reactivity of diphenylcarbene (DPC) in binary solvent
mixtures and ligand dissociation reactions of so-called CO-releasing molecules (CORMs),
special emphasis was put on the implementation and characterization of methods improving
and extending the signal detection in conventional pump–probe transient absorption setups.
The assumption of DPC being an archetypal triplet-ground-state arylcarbene was recently questioned by matrix-isolation studies at low temperatures. DPC embedded in argon matrices revealed a hitherto unknown reactivity when the carbene environment was modified by small amounts of methanol dopant molecules. To complement these findings with liquid-phase experiments at room temperature, femtosecond pump–probe transient absorption spectroscopy with probing in the visible and ultraviolet regime was employed to unravel primary reaction processes of DPC in solvent mixtures. Supported by quantum chemical simulations conducted by our collaborators, it was shown that a competition between the reaction pathways occurs that not only depends on the solvent molecule near-by but also on its interaction with other solvent molecules. In-depth analysis of the solvation dynamics and the amount of nascent intermediates corroborates the importance of a hydrogen-bonded complex with a protic solvent molecule, in striking analogy to complexes found at cryogenic temperatures.
Probing the transient absorption of molecules in the mid-infrared spectral range benefits from the high chemical specificity of molecules’ vibrational signatures. The technique of chirped-pulse upconversion (CPU) constitutes a promising alternative to standard direct multichannel MCT detection when accessing this spectral detection window. Hence, one chapter of this thesis is dedicated to a direct comparison between both detection methods. By conducting an exemplary pump–probe transient absorption experiment, it became evident, that the additional nonlinear interaction step is responsible for increased noise levels when using CPU. However, a correction procedure capable of removing these additional noise contributions—stemming from the fundamental laser radiation used for upconversion—was successfully tested. Perhaps most importantly for various spectroscopic applications, CPU scored with a significantly extended detection bandwidth owing to the high pixel numbers of modern CCD cameras.
Transition-metal complexes capable of releasing small molecular messengers upon photoactivation are promising sources of gasotransmitters such as carbon monoxide (CO) or nitric oxide (NO) in biological applications. However, only little is known about the characteristic time scales of ligand dissociation in this class of molecules. For this purpose, two complexes were investigated with femtosecond time resolution: [Mn(CO)3(tpm)]Cl with tpm=tris(2-pyrazolyl)methane, a manganese tricarbonyl complex which has proven to be selective and cytotoxic to cancer cells, and [Mo(CO)2(NO)(iPr3tacn)]PF6 with iPr3tacn=1,4,7-triisopropyl-1,4,7-triazacyclononane, a molybdenum complex containing both carbonyl and nitrosyl ligands. By conducting pump–probe transient absorption measurements in different spectral probing windows supported by quantum chemical calculations and linear absorption spectroscopy, it was shown that both complexes are able to release one CO ligand within the first few picoseconds after UV excitation. The results complement existing studies which focused on the molecules’ ligand-releasing properties upon long-term exposure. The additional information gained on an ultrafast time scale provides a comprehensive understanding of individual reaction steps connected with ligand release in this class of molecules. Hence, the studies might create new incentives to develop modified molecules for specific applications.
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.