Materials properties from electron density distributions: Molecular magnetism and linear optical properties

The general goal of this thesis is correlating observable properties of organic and metal-organic materials with their ground-state electron density distribution. In a long-term view, we expect to develop empirical or semi-empirical approaches to predict materials properties from the electron density of their building blocks, thus allowing to rationally engineering molecular materials from their constituent subunits, such as their functional groups. In particular, we have focused on linear optical properties of naturally occurring amino acids and their organic and metal-organic derivatives, and on magnetic properties of metal-organic frameworks. For analysing the optical properties and the magnetic behaviour of the molecular or sub-molecular building blocks in materials, we mostly used the more traditional QTAIM partitioning scheme of the molecular or crystalline electron densities, however, we have also investigated a new approach, namely, X-ray Constrained Extremely Localized Molecular Orbitals (XC-ELMO), that can be used in future to extracted the electron densities of crystal subunits. With the purpose of rationally engineering linear optical materials, we have calculated atomic and functional group polarizabilities of amino acid molecules, their hydrogen-bonded aggregates and their metal-organic frameworks. This has enabled the identification of the most efficient functional groups, able to build-up larger electric susceptibilities in crystals, as well as the quantification of the role played by intermolecular interactions and coordinative bonds on modifying the polarizability of the isolated building blocks. Furthermore, we analysed the dependence of the polarizabilities on the one-electron basis set and the many-electron Hamiltonian. This is useful for selecting the most efficient level of theory to estimate susceptibilities of molecular-based materials. With the purpose of rationally design molecular magnetic materials, we have investigated the electron density distributions and the magnetism of two copper(II) pyrazine nitrate metal-organic polymers. High-resolution X-ray diffraction and DFT calculations were used to characterize the magnetic exchange pathways and to establish relationships between the electron densities and the exchange-coupling constants. Moreover, molecular orbital and spin-density analyses were employed to understand the role of different magnetic exchange mechanisms in determining the bulk magnetic behaviour of these materials. As anticipated, we have finally investigated a modified version of the X-ray constrained wavefunction technique, XC-ELMOs, that is not only a useful tool for determination and analysis of experimental electron densities, but also enables one to derive transferable molecular orbitals strictly localized on atoms, bonds or functional groups. In future, we expect to use XC-ELMOs to predict materials properties of large systems, currently challenging to calculate from first-principles, such as macromolecules or polymers. Here, we point out advantages, needs and pitfalls of the technique. This work fulfils, at least partially, the prerequisites to understand materials properties of organic and metal-organic materials from the perspective of the electron density distribution of their building blocks. Empirical or semi-empirical evaluation of optical or magnetic properties from a preconceived assembling of building blocks could be extremely important for rationally design new materials, a field where accurate but expensive first-principles calculations are generally not used. This research could impact the community in the fields of crystal engineering, supramolecular chemistry and, of course, electron density analysis.

Treatment planning aspects and Monte Carlo methods in proton therapy

Over the last years, the interest in proton radiotherapy is rapidly increasing. Protons provide superior physical properties compared with conventional radiotherapy using photons. These properties result in depth dose curves with a large dose peak at the end of the proton track and the finite proton range allows sparing the distally located healthy tissue. These properties offer an increased flexibility in proton radiotherapy, but also increase the demand in accurate dose estimations. To carry out accurate dose calculations, first an accurate and detailed characterization of the physical proton beam exiting the treatment head is necessary for both currently available delivery techniques: scattered and scanned proton beams. Since Monte Carlo (MC) methods follow the particle track simulating the interactions from first principles, this technique is perfectly suited to accurately model the treatment head. Nevertheless, careful validation of these MC models is necessary. While for the dose estimation pencil beam algorithms provide the advantage of fast computations, they are limited in accuracy. In contrast, MC dose calculation algorithms overcome these limitations and due to recent improvements in efficiency, these algorithms are expected to improve the accuracy of the calculated dose distributions and to be introduced in clinical routine in the near future.

From Complex to Stochastic Potential: Heavy Quarkonia in the Quark-Gluon Plasma

The in-medium physics of heavy quarkonium is an ideal proving ground for our ability to connect knowledge about the fundamental laws of physics to phenomenological predictions. One possible route to take is to attempt a description of heavy quark bound states at finite temperature through a Schrödinger equation with an instantaneous potential. Here we review recent progress in devising a comprehensive approach to define such a potential from first principles QCD and extract its, in general complex, values from non-perturbative lattice QCD simulations. Based on the theory of open quantum systems we will show how to interpret the role of the imaginary part in terms of spatial decoherence by introducing the concept of a stochastic potential. Shortcomings as well as possible paths for improvement are discussed.

Redox control of thermopower and figure of merit in phase-coherent molecular wires

We demonstrate how redox control of intra-molecular quantum interference in phase-coherent molecular wires can be used to enhance the thermopower (Seebeck coefficient) S and thermoelectric figure of merit ZT of single molecules attached to nanogap electrodes. Using first principles theory, we study the thermoelectric properties of a family of nine molecules, which consist of dithiol-terminated oligo (phenylene-ethynylenes) (OPEs) containing various central units. Uniquely, one molecule of this family possesses a conjugated acene-based central backbone attached via triple bonds to terminal sulfur atoms bound to gold electrodes and incorporates a fully conjugated hydroquinonecentral unit. We demonstrate that both S and the electronic contribution Z el T to the figure of merit ZT can be dramatically enhanced by oxidizing the hydroquinone to yield a second molecule, which possesses a cross-conjugated anthraquinone central unit. This enhancement originates from the conversion of the pi-conjugation in the former to cross-conjugation in the latter, which promotes the appearance of a sharp anti-resonance at the Fermi energy. Comparison with thermoelectric properties of the remaining seven conjugated molecules demonstrates that such large values of S and Z el T are unprecedented. We also evaluate the phonon contribution to the thermal conductance, which allows us to compute the full figure of merit ZT = Z el T/(1 + κ p/κ el), where κ p is the phonon contribution to the thermal conductance and κ el is the electronic contribution. For unstructured gold electrodes, κ p/κ el Gt⃒ 1 and therefore strategies to reduce κ p are needed to realize the highest possible figure of merit.

A first look at maximally twisted mass lattice QCD calculations at the physical point

In this contribution, a first look at simulations using maximally twisted mass Wilson fermions at the physical point is presented. A lattice action including clover and twisted mass terms is presented and the Monte Carlo histories of one run with two mass-degenerate flavours at a single lattice spacing are shown. Measurements from the light and heavy-light pseudoscalar sectors are compared to previous Nf = 2 results and their phenomenological values. Finally, the strategy for extending simulations to Nf = 2+1+1 is outlined.

From Structure to Function: Characterization of Cu(I) Adducts in Leveler Additives by DFT Calculations

We present the first molecular model of the coordination complex formed by Cu(I) and imidazole-epichlorohydrin polymers. Our calculations show that the Cu(I) ion has linear coordination and the whole complex has neutral charge. Our model suggests salt couple pairing as the driving force for the formation of the surface-confined precipitation, which is crucial to obtain flat surfaces in industrial copper deposition processes, required for mass fabrication of state-of-the-art electronic and memory devices.

Experimental dissolution of molybdenum-sulphides at low oxygen concentrations: A first-order approximation of late Archean atmospheric conditions

The abundance of atmospheric oxygen and its evolution through Earth's history is a highly debated topic. The earliest change of the Mo concentration and isotope composition of marine sediments are interpreted to be linked to the onset of the accumulation of free O2 in Earth's atmosphere. The O2 concentration needed to dissolve significant amounts of Mo in water is not yet quantified, however. We present laboratory experiments on pulverized and surface-cleaned molybdenite (MoS2) and a hydrothermal breccia enriched in Mo-bearing sulphides using a glove box setup. Duration of an experiment was 14 days, and first signs of oxidation and subsequent dissolution of Mo compounds start to occur above an atmospheric oxygen concentration of 72 ± 20 ppmv (i.e., 2.6 to 4.6 × 10−4 present atmospheric level (PAL)). This experimentally determined value coincides with published model calculations supporting atmospheric O2 concentrations between 1 × 10−5 to 3 × 10−4 PAL prior to the Great Oxidation Event and sets an upper limit to the molecular oxygen needed to trigger Mo accumulation and Mo isotope variations recorded in sediments. In combination with the published Mo isotope composition of the rock record, this result implies an atmospheric oxygen concentration prior to 2.76 Ga of below 72 ± 20 ppmv.