@phdthesis{Graefe2005,
  author    = {Gr{\"a}fe, Stefanie},
  title     = {Laser-control of molecular dynamics},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-13388},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2005},
  abstract  = {In this work a new algorithm to determine quantum control fields from the instantaneous response of systems has been developed. The derived fields allow to establish a direct connection between the applied perturbation and the molecular dynamics. The principle is most easily illustrated in regarding a classical forced oscillator. A particle moving inside the respective potential is accelerated if an external field is applied acting in the same direction as its momentum (heating). In contrary, a deceleration is achieved by a field acting in the opposite direction as the momentum (cooling). Furthermore, when the particle reaches a classical turning point and then changes its direction, the sign of the field has to be changed to further drive the system in the desired way. The frequency of the field therefore is in resonance with the oscillator. This intuitively clear picture of a driven classical oscillator can be used for directing (or controlling) quantum mechanical wave packet motion. The efficiency of the instantaneous dynamics algorithm was demonstrated in treating various model problems, the population transfer in double well potentials, excitation and dissociation of selective modes, and the population transfer between electronic states. Although it was not tried to optimize the fields to gain higher yields, the control was found to be very efficient. Driving population transfer in a double well potential could be shown to take place with nearly 100\% efficiency. It was shown that selective dissociation within the electronic ground state of HOD can be performed by either maximizing a selected coordinate's differential momentum change or the energy absorption. Concerning the population transfer into excited electronic states, a direct comparison with common control algorithms as optimal control theory and genetic algorithms was accomplished using a one-dimensional representation of methyl iodide. The fields derived from the various control theories were effective in transferring population into the chosen target state but the underlying physical background of the derived optimal fields was not obvious to explain. The instantaneous dynamics algorithm allowed to establish a direct relation between the derived fields and the underlying molecular dynamics. Bound-to-bound transitions could be handled more effectively. This was demonstrated on the sodium dimer in a representation of 3 electronic states being initially in its vibronic ground state. The objective was to transfer population into a predefined excited state. Choosing the first or the second state as a target, the control fields exhibited quite different features. The pulse-structure is related to the excited state wave packet, moving in, and out of the Franck-Condon region. Changing the control objective, the derived control field performed pure electronic transitions on a fast time-scale via a two-step transition. Futhermore, orientational effects have been investigated. The overall-efficiency of the population transfer for differently oriented molecules was about 70 \% or more if applying a control field derived for a 45° orientation. Spectroscopic methods to gain information about the outcome of the control process have been investigated. It was shown that pump/probe femtosecond ionization spectroscopy is suited to monitor time-dependent molecular probability distributions. In particular, time-dependent photoelectron spectra are able to monitor the population in the various electronic states. In the last chapter a different possibility of controlling molecules was regarded by investigating molecular iodine with a setup similar to the STIRAP ("Stimulated Raman Adiabatic passage") scenario. The possibility to extend this technique to a fs-time scale was examined in theory as well as in experiments, the latter being performed by Dr. Torsten Siebert in the Kiefer group, University of W{\"u}rzburg. It was shown that off-resonant excitation with implementation of the pulses with a higher intensity of the Stokes pulse as compared to the pump pulse - describing a so-called f-STIRAP like configuration - was shown to effectively transfer population into excited ground-state vibrational levels. This was theoretically underlined by comparing the numerically exact coupling case with the adiabatic picture. The process was described to run in the vicinity of adibaticity. A new model explaining the process by the system's vector rotating around the dressed state vector will be adopted in future calculations. Altogether, a new promising algorithm to control dynamical processes based on the instantaneous response has been developed. Because the derived control fields have been shown to be very efficient in selectively influencing molecules, it is to be expected that farther reaching applications can be realized in future investigations.},
  subject      = {Laserstrahlung},
  language  = {en}
}
@phdthesis{Papastathopoulos2005,
  author    = {Papastathopoulos, Evangelos},
  title     = {Adaptive control of electronic excitation utilizing ultrafast laser pulses},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-12533},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2005},
  abstract  = {The subject of this work has been the investigation of dynamical processes that occur during and after the interaction of matter with pulses of femtosecond laser radiation. The experiments presented here were performed in the gas phase and involve one atomic and several model molecular systems. Absorption of femtosecond laser radiation by these systems induces an electronic excitation, and subsequently their ionization, photofragmentation or isomerization. The specific adjustment of the excitation laser field properties offers the possibility to manipulate the induced electronic excitation and to influence the formation of the associated photoproducts. From the perspective of the employed spectroscopic methods, the development of photoelectron spectroscopy and its implementation in laser control experiments has been of particular interest in this thesis. This technique allows for a most direct and intuitive observation of electronic excitation dynamics in atomic as well as in complex polyatomic molecular systems. The propagation of an intermediate electronic transient state, associated to the formation of a particular photoproduct, can be interrogated by means of its correlation to a specific state of the atomic or molecular continuum. Such correlations involve the autoionization of the transient state, or by means of a second probe laser field, a structural correlation, as summarized by the Koopman's theorem (section 2.4.1). The technique of adaptive femtosecond quantum control has been the subject of development in our group for many years. The basic method, by which the temporal profile of near-infrared laser pulses at a central wavelength of 800 nm, can be adjusted, is a programmable femtosecond pulse-shaper that comprises of a zero dispersion compressor and a commercial liquid crystal modulator (LCD). This experimental arrangement was realized prior to this thesis and served as a starting point to extend the pulse-shaping technique to the ultraviolet spectral region. This technological development was realized for the purposes of the experiments presented in Chapter 5. It involves a combination of the LCD-pulse-shaper with frequency up-conversion techniques on the basis of producing specifically modulated laser pulses of central wavelength 266 nm. Furthermore, the optical method X-FROG had to be developed in order to characterize the often complex structure of generated ultraviolet pulses. In the adaptive control experiments presented in this work, the generated femtosecond laser pulses could be automatically adjusted by means of specifically addressing the 128 independent voltage parameters of the programmable liquid-crystal modulator. Additionally a machine learning algorithm was employed for the cause of defining laser pulse-shapes that delivered the desired (optimal) outcome in the investigated laser interaction processes. In Chapter 4, the technique of feedback-controlled femtosecond pulse shaping was combined with time-of-flight mass spectroscopy as well as photoelectron spectroscopy in order to investigate the multiphoton double ionization of atomic calcium. A pronounced absolute enhancement of the double ionization yield was obtained with optimized femtosecond laser pulses. On the basis of the measured photoelectron spectra and of the electron optimization experiments, a non-sequential process was found, which plays an important role in the formation of doubly charged Calcium ions. Then in Chapter 5, the dynamics following the pp* excitation of ethylene-like molecules were investigated. In this context, the model molecule stilbene was studied by means of femtosecond photoelectron spectroscopy. Due to the simplicity of its chemical structure, stilebene is one of the most famous models used in experimental as well as theoretical studies of isomerization dynamics. From the time-resolved experiments described in that chapter, new spectroscopic data involving the second excited electronic state S2 of the molecule were acquired. The second ethylenic product was the molecule tetrakis (dimethylamino) ethylene (TDMAE). Due to the presence of numerous lone pair electrons on the four dimethylamino groups, TDMAE exhibits a much more complex structure than stilbene. Nevertheless, previously reported studies on the dynamics of TDMAE provided vital information for planning and conducting a successful optimisation control experiment of the wavepacket propagation upon the (pp*) S1 excited potential surface of the molecule. Finally, in Chapter 6 the possibility of employing femtosecond laser pulses as an alternative method for activating a metallocene molecular catalyst was addressed. By means of an adaptive laser control scheme, an optimization experiment was realized. There, the target was the selective cleavage of one methyl-ligand of the model catalyst (Cp)^2Zr(CH3)^2, which induces a catalytic coordination position on the molecule. The spectroscopic studies presented in that chapter were performed in collaboration to the company BASF A.G. and constitute a proof-of principle attempt for a commercial application of the adaptive femtosecond quantum control technique.},
  subject      = {Ultrakurzer Lichtimpuls},
  language  = {en}
}
@phdthesis{Merget2015,
  author    = {Merget, Benjamin},
  title     = {Computational methods for assessing drug-target residence times in bacterial enoyl-ACP reductases and predicting small-molecule permeability for the \(Mycobacterium\) \(tuberculosis\) cell wall},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-127386},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2015},
  abstract  = {\textbf{Molecular Determinants of Drug-Target Residence Times of Bacterial Enoyl-ACP Reductases.} Whereas optimization processes of early drug discovery campaigns are often affinity-driven, the drug-target residence time \$t_R\$ should also be considered due to an often strong correlation with \textit{in vivo} efficacy of compounds. However, rational optimization of \$t_R\$ is not straightforward and generally hampered by the lack of structural information about the transition states of ligand association and dissociation. The enoyl-ACP reductase FabI of the fatty acid synthesis (FAS) type II is an important drug-target in antibiotic research. InhA is the FabI enzyme of \textit{Mycobacterium tuberculosis}, which is known to be inhibited by various compound classes. Slow-onset inhibition of InhA is assumed to be associated with the ordering of the most flexible protein region, the substrate binding loop (SBL). Diphenylethers are one class of InhA inhibitors that can promote such SBL ordering, resulting in long drug-target residence times. Although these inhibitors are energetically and kinetically well characterized, it is still unclear how the structural features of a ligand affect \$t_R\$. Using classical molecular dynamics (MD) simulations, recurring conformational families of InhA protein-ligand complexes were detected and structural determinants of drug-target residence time of diphenyl\-ethers with different kinetic profiles were described. This information was used to deduce guidelines for efficacy improvement of InhA inhibitors, including 5'-substitution on the diphenylether B-ring. The validity of this suggestion was then analyzed by means of MD simulations. Moreover, Steered MD (SMD) simulations were employed to analyze ligand dissociation of diphenylethers from the FabI enzyme of \textit{Staphylococcus aureus}. This approach resulted in a very accurate and quantitative linear regression model of the experimental \$ln(t_R)\$ of these inhibitors as a function of the calculated maximum free energy change of induced ligand extraction. This model can be used to predict the residence times of new potential inhibitors from crystal structures or valid docking poses. Since correct structural characterization of the intermediate enzyme-inhibitor state (EI) and the final state (EI*) of two-step slow-onset inhibition is crucial for rational residence time optimization, the current view of the EI and EI* states of InhA was revisited by means of crystal structure analysis, MD and SMD simulations. Overall, the analyses affirmed that the EI* state is a conformation resembling the 2X23 crystal structure (with slow-onset inhibitor \textbf{PT70}), whereas a twist of residues Ile202 and Val203 with a further opened helix \$\alpha 6\$ corresponds to the EI state. Furthermore, MD simulations emphasized the influence of close contacts to symmetry mates in the SBL region on SBL stability, underlined by the observation that an MD simulation of \textbf{PT155} chain A with chain B' of a symmetry mate in close proximity of the SBL region showed significantly more stable loops, than a simulation of the tetrameric assembly. Closing Part I, SMD simulations were employed which allow the delimitation of slow-onset InhA inhibitors from rapid reversible ligands. \textbf{Prediction of \textit{Mycobacterium tuberculosis} Cell Wall Permeability.} The cell wall of \textit{M. tuberculosis} hampers antimycobacterial drug design due to its unique composition, providing intrinsic antibiotic resistance against lipophilic and hydrophilic compounds. To assess the druggability space of this pathogen, a large-scale data mining endeavor was conducted, based on multivariate statistical analysis of differences in the physico-chemical composition of a normally distributed drug-like chemical space and a database of antimycobacterial--and thus very likely permeable--compounds. The approach resulted in the logistic regression model MycPermCheck, which is able to predict the permeability probability of small organic molecules based on their physico-chemical properties. Evaluation of MycPermCheck suggests a high predictive power. The model was implemented as a freely accessible online service and as a local stand-alone command-line version. Methodologies and findings from both parts of this thesis were combined to conduct a virtual screening for antimycobacterial substances. MycPermCheck was employed to screen the chemical permeability space of \textit{M. tuberculosis} from the entire ZINC12 drug-like database. After subsequent filtering steps regarding ADMET properties, InhA was chosen as an exemplary target. Docking to InhA led to a principal hit compound, which was further optimized. The quality of the interaction of selected derivatives with InhA was subsequently evaluated using MD and SMD simulations in terms of protein and ligand stability, as well as maximum free energy change of induced ligand egress. The results of the presented computational experiments suggest that compounds with an indole-3-acethydrazide scaffold might constitute a novel class of InhA inhibitors, worthwhile of further investigation.},
  subject      = {Computational chemistry},
  language  = {en}
}
@phdthesis{Becker2015,
  author    = {Becker, Johannes},
  title     = {Development and implementation of new simulation possibilities in the CAST program package},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-132032},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2015},
  abstract  = {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.},
  subject      = {Molekulardynamik},
  language  = {en}
}
@phdthesis{Lindner2019,
  author    = {Lindner, Joachim Oliver},
  title     = {Multistate Metadynamics with Electronic Collective Variables},
  doi       = {10.25972/OPUS-19163},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-191638},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2019},
  abstract  = {The aim of this thesis was to develop new automatic enhanced sampling methods by extending the idea of Parrinello's metadynamics to multistate problems and by introducing new quantum-mechanical electronic collective variables. These methods open up a rich perspective for applications to the photophysical processes in complex molecular systems, which play a major role in many natural processes such as vision and photosynthesis, but also in the development of new materials for organic electronics, whose function depends on specific electronic properties such as biradicalicity.},
  subject      = {Theoretische Chemie},
  language  = {en}
}