Einzelmolekülspektroskopische und quantenchemische Untersuchungen zum elektronischen Energietransfer

In der vorliegenden Arbeit wurden Untersuchungen zum Mechanismus, der Dynamik und der Kontrolle des elektronischen Energietransfers in multichromophoren Modellsystemen durchgeführt. Als Untersuchungsmethoden wurden hauptsächlich die konfokale Einzelmolekülspektroskopie und die Quantenchemie eingesetzt. Der Aufbau des Einzelmolekülmikroskops wurde bezüglich der Anregungs- und Detektionskomponenten variiert, um die unterschiedlichen Experimente durchzuführen. Die quantenchemischen Rechnungen wurden auf Dichtefunktional- und Coupled-Cluster-Niveau durchgeführt. Die aus den Rechnungen erhaltenen zusätzlichen Informationen über experimentell zum Teil schwer zugängliche Eigenschaften der Farbstoffe unterstützten die Interpretation der experimentellen Befunde. rnIn früheren Untersuchungen der AG Basché wurden die Energietransfer-Raten von Donor-Akzeptor-Systemen gemessen, die erhebliche Abweichungen von nach der Förster-Theorie berechneten Raten zeigten. Daher war ein Ziel der vorliegenden Arbeit, diese Abweichungen zu erklären. Zu diesem Zweck wurde die Geometrie der Diaden experimentell untersucht, sowie die elektronische Kopplung zwischen den Chromophoren quantenchemisch berechnet. Die relative Orientierung der Chromophore in den Diaden wurde in einem Einzelmolekül-Experiment mit rotierender Anregungspolarisation abgefragt. Die erhaltenen Winkelverteilungen konnten schließlich eindeutig auf die Flexibilität der die Chromophore verbrückenden Oligophenyl-Einheiten zurückgeführt werden. Die Unterschiede der gemessenen Energietransfer-Raten zu den nach der Förster-Theorie ermittelten Werten konnten jedoch nicht über die molekulare Flexibilität der Systeme erklärt werden. Aufklärung über die Diskrepanzen zur Förster-Theorie ergaben die quantenchemischen Rechnungen. In Rahmen dieser Arbeit wurde zum ersten Mal die Coupled-Cluster-Theorie zur Berechnung der elektronischen Kopplung eingesetzt. Die Betrachtung der isolierten Chromophore reichte aber nicht aus, um die gemessenen Abweichungen von der Förster-Theorie zu erklären. Erst über die Berücksichtigung der molekularen Brücke konnten die gefunden Abweichungen erklärt werden. Die deutliche Verstärkung der elektronischen Kopplung ist auf die Polarisierbarkeit der Brücke zurückzuführen.rnNach diesen Betrachtungen stand die Kontrolle des Energietransfers im Fokus der weiteren Untersuchungen. In den durchgeführten Einzelmolekülexperimenten wurden die Chromophore der Donor-Akzeptor-Systeme selektiv mit zwei Laserpulsen unterschiedlicher Wellenlänge angeregt. Beim gleichzeitigen Anregen beider Chromophore wurde Singulett-Singulett-Annihilation (SSA) induziert, ein Energietransferprozess, bei dem die Anregungsenergie vom vorigen Akzeptor zum vorigen Donor übertragen wird. Da über SSA Fluoreszenzphotonen gelöscht wurden, konnte über den Abstand der Laserpulse die Fluoreszenzintensität des einzelnen Moleküls moduliert werden. Konzeptionell verwandte Einzelmolekülexperimente wurden an einem weiteren molekularen System durchgeführt, das aus einem Kern und einer Peripherie bestand. Fluoreszenzauszeiten des Gesamtsystems bei selektiver Anregung des Kerns wurden auf die Population eines Triplett-Zustandes zurückgeführt, der die Fluoreszenz der Peripherie löschte. rnAbschließend wurde der SSA-Prozess zwischen zwei gleichartigen Chromophoren untersucht. Es wurde eine Methode entwickelt, die es zum ersten Mal erlaubte, die SSA-Zeitkonstante individueller Moleküle zu bestimmen. Hierfür wurden die Daten der gemessenen Photonen-Koinzidenzhistogramme mittels eines im Rahmen dieser Arbeit hergeleiteten analytischen Zusammenhangs ausgewertet, der über Monte-Carlo-Simulationen bestätigt wurde.

Scale bridging concepts in molecular simulation : coarse-graining and thermodynamic coupling

In dieser Arbeit wurden Simulation von Flüssigkeiten auf molekularer Ebene durchgeführt, wobei unterschiedliche Multi-Skalen Techniken verwendet wurden. Diese erlauben eine effektive Beschreibung der Flüssigkeit, die weniger Rechenzeit im Computer benötigt und somit Phänomene auf längeren Zeit- und Längenskalen beschreiben kann.rnrnEin wesentlicher Aspekt ist dabei ein vereinfachtes (“coarse-grained”) Modell, welches in einem systematischen Verfahren aus Simulationen des detaillierten Modells gewonnen wird. Dabei werden ausgewählte Eigenschaften des detaillierten Modells (z.B. Paar-Korrelationsfunktion, Druck, etc) reproduziert.rnrnEs wurden Algorithmen untersucht, die eine gleichzeitige Kopplung von detaillierten und vereinfachten Modell erlauben (“Adaptive Resolution Scheme”, AdResS). Dabei wird das detaillierte Modell in einem vordefinierten Teilvolumen der Flüssigkeit (z.B. nahe einer Oberfläche) verwendet, während der Rest mithilfe des vereinfachten Modells beschrieben wird.rnrnHierzu wurde eine Methode (“Thermodynamische Kraft”) entwickelt um die Kopplung auch dann zu ermöglichen, wenn die Modelle in verschiedenen thermodynamischen Zuständen befinden. Zudem wurde ein neuartiger Algorithmus der Kopplung beschrieben (H-AdResS) der die Kopplung mittels einer Hamilton-Funktion beschreibt. In diesem Algorithmus ist eine zur Thermodynamischen Kraft analoge Korrektur mit weniger Rechenaufwand möglich.rnrnAls Anwendung dieser grundlegenden Techniken wurden Pfadintegral Molekulardynamik (MD) Simulationen von Wasser untersucht. Mithilfe dieser Methode ist es möglich, quantenmechanische Effekte der Kerne (Delokalisation, Nullpunktsenergie) in die Simulation einzubeziehen. Hierbei wurde zuerst eine Multi-Skalen Technik (“Force-matching”) verwendet um eine effektive Wechselwirkung aus einer detaillierten Simulation auf Basis der Dichtefunktionaltheorie zu extrahieren. Die Pfadintegral MD Simulation verbessert die Beschreibung der intra-molekularen Struktur im Vergleich mit experimentellen Daten. Das Modell eignet sich auch zur gleichzeitigen Kopplung in einer Simulation, wobei ein Wassermolekül (beschrieben durch 48 Punktteilchen im Pfadintegral-MD Modell) mit einem vereinfachten Modell (ein Punktteilchen) gekoppelt wird. Auf diese Weise konnte eine Wasser-Vakuum Grenzfläche simuliert werden, wobei nur die Oberfläche im Pfadintegral Modell und der Rest im vereinfachten Modell beschrieben wird.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Efficient and Accurate Characterization of the Bergman Cyclization for an Expanded Substructure of Esperamicin A1

Incorporation of enediynes into anticancer drugs remains an intriguing yet elusive strategy for the design of therapeutically active agents. Density functional theory was used to locate reactants, products, and transition states along the Bergman cyclization pathways connecting enediynes to reactive para-biradicals. Sum method correction to low-level calculations confirmed B3LYP/6-31G(d,p) as the method of choice in investigating enediynes. Herein described as MI:Sum, calculated reaction enthalpies differed from experiment by an average of 2.1 kcal·mol−1 (mean unsigned error). A combination of strain energy released across the reaction coordinate and the critical intramolecular distance between reacting diynes explains reactivity differences. Where experimental and calculated barrier heights are in disagreement, higher level multireference treatment of the enediynes confirms lower level estimates. Previous work concerning the chemically reactive fragment of esperamcin, MTC, is expanded to our model system MTC2.

Au@Hg nanoalloy formation through direct amalgamation: Structural, spectroscopic, and computational evidence for slow nanoscale diffusion

Dynamic core-shell nanoparticles have received increasing attention in recent years. This paper presents a detailed study of Au-Hg nanoalloys, whose composing elements show a large difference in cohesive energy. A simple method to prepare Au@Hg particles with precise control over the composition up to 15 atom% mercury is introduced, based on reacting a citrate stabilized gold sol with elemental mercury. Transmission electron microscopy shows an increase of particle size with increasing mercury content and, together with X-ray powder diffraction, points towards the presence of a core-shell structure with a gold core surrounded by an Au-Hg solid solution layer. The amalgamation process is described by pseudo-zero-order reaction kinetics, which indicates slow dissolution of mercury in water as the rate determining step, followed by fast scavenging by nanoparticles in solution. Once adsorbed at the surface, slow diffusion of Hg into the particle lattice occurs, to a depth of ca. 3 nm, independent of Hg concentration. Discrete dipole approximation calculations relate the UV-vis spectra to the microscopic details of the nanoalloy structure. Segregation energies and metal distribution in the nanoalloys were modeled by density functional theory calculations. The results indicate slow metal interdiffusion at the nanoscale, which has important implications for synthetic methods aimed at core-shell particles.

A Computational Analysis of the Collision-Induced Dissociation Mechanisms of the Dipeptide Glycine-Serine Under Mass Spectrometry

Collision-induced dissociation (CID) of peptides using tandem mass spectrometry (MS) has been used to determine the identity of peptides and other large biological molecules. Mass spectrometry (MS) is a useful tool for determining the identity of molecules based on their interaction with electromagnetic fields. If coupled with another method like infrared (IR) vibrational spectroscopy, MS can provide structural information, but in its own right, MS can only provide the mass-to-charge (m/z) ratio of the fragments produced, which may not be enough information to determine the mechanism of the collision-induced dissociation (CID) of the molecule. In this case, theoretical calculations provide a useful companion for MS data and yield clues about the energetics of the dissociation. In this study, negative ion electrospray tandem MS was used to study the CID of the deprotonated dipeptide glycine-serine (Gly-Ser). Though negative ion MS is not as popular a choice as positive ion MS, studies by Bowie et al. show that it yields unique clues about molecular structure which complement positive ion spectroscopy, such as characteristic fragmentations like the loss of formaldehyde from the serine residue.2 The increase in the collision energy in the mass spectrometer alters the flexibility of the dipeptide backbone, enabling isomerizations (reactions not resulting in a fragment loss) and dissociations to take place. The mechanism of the CID of Gly-Ser was studied using two computational methods, B3LYP/6-311+G* and M06-2X/6-311++G**. The main pathway for molecular dissociation was analyzed in 5 conformers in an attempt to verify the initial mechanism proposed by Dr. James Swan after examination of the MS data. The results suggest that the loss of formaldehyde from serine, which Bowie et al. indicates is a characteristic of the presence of serine in a protein residue, is an endothermic reaction that is made possible by the conversion of the translational energy of the ion into internal energy as the ion collides with the inert collision gas. It has also been determined that the M06-2X functional¿s improved description of medium and long-range correlation makes it more effective than the B3LYP functional at finding elusive transition states. M06-2X also more accurately predicts the energy of those transition states than does B3LYP. A second CID mechanism, which passes through intermediates with the same m/z ratio as the main pathway for molecular dissociation, but different structures, including a diketopiperazine intermediate, was also studied. This pathway for molecular dissociation was analyzed with 3 conformers and the M06-2X functional, due to its previously determined effectiveness. The results suggest that the latter pathway, which meets the same intermediate masses as the first mechanism, is lower in overall energy and therefore a more likely pathway of dissociation than the first mechanism.

Spin-strain coupling in a 3-D transition metal oxides

ab-initio Hartree Fock (HF), density functional theory (DFT) and hybrid potentials were employed to compute the optimized lattice parameters and elastic properties of perovskite 3-d transition metal oxides. The optimized lattice parameters and elastic properties are interdependent in these materials. An interaction is observed between the electronic charge, spin and lattice degrees of freedom in 3-d transition metal oxides. The coupling between the electronic charge, spin and lattice structures originates due to localization of d-atomic orbitals. The coupling between the electronic charge, spin and crystalline lattice also contributes in the ferroelectric and ferromagnetic properties in perovskites. The cubic and tetragonal crystalline structures of perovskite transition metal oxides of ABO3 are studied. The electronic structure and the physics of 3-d perovskite materials is complex and less well considered. Moreover, the novelty of the electronic structure and properties of these perovskites transition metal oxides exceeds the challenge offered by their complex crystalline structures. To achieve the objective of understanding the structure and property relationship of these materials the first-principle computational method is employed. CRYSTAL09 code is employed for computing crystalline structure, elastic, ferromagnetic and other electronic properties. Second-order elastic constants (SOEC) and bulk moduli (B) are computed in an automated process by employing ELASTCON (elastic constants) and EOS (equation of state) programs in CRYSTAL09 code. ELASTCON, EOS and other computational algorithms are utilized to determine the elastic properties of tetragonal BaTiO3, rutile TiO2, cubic and tetragonal BaFeO3 and the ferromagentic properties of 3-d transition metal oxides. Multiple methods are employed to crosscheck the consistency of our computational results. Computational results have motivated us to explore the ferromagnetic properties of 3-d transition metal oxides. Billyscript and CRYSTAL09 code are employed to compute the optimized geometry of the cubic and tetragonal crystalline structure of transition metal oxides of Sc to Cu. Cubic crystalline structure is initially chosen to determine the effect of lattice strains on ferromagnetism due to the spin angular momentum of an electron. The 3-d transition metals and their oxides are challenging as the basis functions and potentials are not fully developed to address the complex physics of the transition metals. Moreover, perovskite crystalline structures are extremely challenging with respect to the quality of computations as the latter requires the well established methods. Ferroelectric and ferromagnetic properties of bulk, surfaces and interfaces are explored by employing CRYSTAL09 code. In our computations done on cubic TMOs of Sc-Fe it is observed that there is a coupling between the crystalline structure and FM/AFM spin polarization. Strained crystalline structures of 3-d transition metal oxides are subjected to changes in the electromagnetic and electronic properties. The electronic structure and properties of bulk, composites, surfaces of 3-d transition metal oxides are computed successfully.

Distributed Atomic Polarizabilities of Amino Acids and their Hydrogen-Bonded Aggregates

With the purpose of rational design of optical materials, distributed atomic polarizabilities of amino acid molecules and their hydrogen-bonded aggregates are calculated in order to identify the most efficient functional groups, able to buildup larger electric susceptibilities in crystals. Moreover, we carefully analyze how the atomic polarizabilities depend on the one-electron basis set or the many-electron Hamiltonian, including both wave function and density functional theory methods. This is useful for selecting the level of theory that best combines high accuracy and low computational costs, very important in particular when using the cluster method to estimate susceptibilities of molecular-based materials.

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

Experimental and Theoretical Electron Density Analysis of Copper Pyrazine Nitrate Quasi-Low-Dimensional Quantum Magnets

The accurate electron density distribution and magnetic properties of two metal-organic polymeric magnets, the quasi-one-dimensional (1D) Cu(pyz)(NO3)2 and the quasi-two-dimensional (2D) [Cu(pyz)2(NO3)]NO3·H2O, have been investigated by high-resolution single-crystal X-ray diffraction and density functional theory calculations on the whole periodic systems and on selected fragments. Topological analyses, based on quantum theory of atoms in molecules, enabled the characterization of possible magnetic exchange pathways and the establishment of relationships between the electron (charge and spin) densities and the exchange-coupling constants. In both compounds, the experimentally observed antiferromagnetic coupling can be quantitatively explained by the Cu-Cu superexchange pathway mediated by the pyrazine bridging ligands, via a σ-type interaction. From topological analyses of experimental charge-density data, we show for the first time that the pyrazine tilt angle does not play a role in determining the strength of the magnetic interaction. Taken in combination with molecular orbital analysis and spin density calculations, we find a synergistic relationship between spin delocalization and spin polarization mechanisms and that both determine the bulk magnetic behavior of these Cu(II)-pyz coordination polymers.

Oxygen Reduction Reaction on Pt Overlayers Deposited onto a Gold Film: Ligand, Strain, and Ensemble Effect

We study the oxygen reduction reaction (ORR), the catalytic process occurring at the cathode in fuel cells, on Pt layers prepared by electrodeposition onto an Au substrate. Using a nominal Pt layer by layer deposition method previously proposed, imperfect layers of Pt on Au are obtained. The ORR on deposited Pt layers decreases with increasing Pt thickness. In the submonolayer region, however, the ORR activity is superior to that of bulk Pt. Using density functional theory (DFT) calculations, we correlate the observed activity trend to strain, ligand, and ensemble effects. At submonolayer coverage certain atom configurations weaken the binding energies of reaction intermediates due to a ligand and ensemble effect, thus effectively increasing the ORR activity. At higher Pt coverage the activity is governed by a strain effect, which lowers the activity by decreasing the oxidation potential of water. This study is a nice example of how the influence of strain, ligand, and ensemble effects on the ORR can be deconvoluted.

Theoretical study of band alignment in nano-porous ZnO interacting with substituted Phthalocyanines

The aim of this work is the theoretical study of the band alignment between the two components of a hybrid organic-inorganic solar-cell. The working organic molecules are metal tetra-sulphonated phthalocyanines (M-Pc) and the inorganic material is nano-porous ZnO growth in the 001 direction. The theoretical calculations are being made using the density functional theory (DFT) using a GGA functional with the SIESTA code, which projects electron wave functions and density onto a real space grid and uses as basis set a linear combination of numerical, finite-range localized atomic orbitals. We also used the DFT+U method included in the code that allows a semi-empirical inclusion of electronic correlations in the description of electronic spectra for systems such as zinc oxide.