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.

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.

The S-1/S-2 exciton interaction in 2-pyridone center dot 6-methyl-2-pyridone: Davydov splitting, vibronic coupling, and vibronic quenching

The excitonic splitting between the S-1 and S-2 electronic states of the doubly hydrogen-bonded dimer 2-pyridone center dot 6-methyl-2-pyridone (2PY center dot 6M2PY) is studied in a supersonic jet, applying two-color resonant two-photon ionization (2C-R2PI), UV-UV depletion, and dispersed fluorescence spectroscopies. In contrast to the C-2h symmetric (2-pyridone) 2 homodimer, in which the S-1 <- S-0 transition is symmetry-forbidden but the S-2 <- S-0 transition is allowed, the symmetry-breaking by the additional methyl group in 2PY center dot 6M2PY leads to the appearance of both the S-1 and S-2 origins, which are separated by Delta(exp) = 154 cm(-1). When combined with the separation of the S-1 <- S-0 excitations of 6M2PY and 2PY, which is delta = 102 cm(-1), one obtains an S-1/S-2 exciton coupling matrix element of V-AB, el = 57 cm(-1) in a Frenkel-Davydov exciton model. The vibronic couplings in the S-1/S-2 <- S-0 spectrum of 2PY center dot 6M2PY are treated by the Fulton-Gouterman single-mode model. We consider independent couplings to the intramolecular 6a' vibration and to the intermolecular sigma' stretch, and obtain a semi-quantitative fit to the observed spectrum. The dimensionless excitonic couplings are C(6a') = 0.15 and C(sigma') = 0.05, which places this dimer in the weak-coupling limit. However, the S-1/S-2 state exciton splittings Delta(calc) calculated by the configuration interaction singles method (CIS), time-dependent Hartree-Fock (TD-HF), and approximate second-order coupled-cluster method (CC2) are between 1100 and 1450 cm(-1), or seven to nine times larger than observed. These huge errors result from the neglect of the coupling to the optically active intra-and intermolecular vibrations of the dimer, which lead to vibronic quenching of the purely electronic excitonic splitting. For 2PY center dot 6M2PY the electronic splitting is quenched by a factor of similar to 30 (i.e., the vibronic quenching factor is Gamma(exp) = 0.035), which brings the calculated splittings into close agreement with the experimentally observed value. The 2C-R2PI and fluorescence spectra of the tautomeric species 2-hydroxypyridine center dot 6-methyl-2-pyridone (2HP center dot 6M2PY) are also observed and assigned. (C) 2011 American Institute of Physics.

A novel isolation method of Brucella species and molecular tracking of Brucella suis biovar 2 in domestic and wild animals

Brucella suis biovar 2 is the most common aetiological agent of porcine brucellosis in Europe. B. suis biovar 2 is considered to have low zoonotic potential, but is a causative agent of reproductive losses in pigs, and it is thus economically important. The multilocus variable-number of tandem repeats genotyping analysis of 16 loci (MLVA-16) has proven to be highly discriminatory and is the most suitable assay for simultaneously identifying B. suis and tracking infections. The aim of this study was to investigate the relatedness between isolates of B. suis biovar 2 obtained during a brucellosis outbreak in domestic pigs and isolates from wild boars and hares collected from proximal or remote geographical areas by MLVA-16. A cluster analysis of the MLVA-16 data revealed that most of the isolates obtained from Switzerland clustered together, with the exception of one isolate. The outbreak isolates constituted a unique subcluster (with a genetic similarity >93.8%) distinct from that of the isolates obtained from wild animals, suggesting that direct transmission of the bacterium from wild boars to domestic pigs did not occur in this outbreak. To obtain a representative number of isolates for MLVA-16, alternative methods of Brucella spp. isolation from tissue samples were compared with conventional direct cultivation on a Brucella-selective agar. We observed an enhanced sensitivity when mechanical homogenisation was followed by host cell lysis prior to cultivation on the Brucella-selective agar. This work demonstrates that MLVA-16 is an excellent tool for both monitoring brucellosis and investigating outbreaks. Additionally, we present efficient alternatives for the isolation of Brucella spp.

A novel computer simulation method for simulating the multiscale transduction dynamics of signal proteins

Signal proteins are able to adapt their response to a change in the environment, governing in this way a broad variety of important cellular processes in living systems. While conventional molecular-dynamics (MD) techniques can be used to explore the early signaling pathway of these protein systems at atomistic resolution, the high computational costs limit their usefulness for the elucidation of the multiscale transduction dynamics of most signaling processes, occurring on experimental timescales. To cope with the problem, we present in this paper a novel multiscale-modeling method, based on a combination of the kinetic Monte-Carlo- and MD-technique, and demonstrate its suitability for investigating the signaling behavior of the photoswitch light-oxygen-voltage-2-Jα domain from Avena Sativa (AsLOV2-Jα) and an AsLOV2-Jα-regulated photoactivable Rac1-GTPase (PA-Rac1), recently employed to control the motility of cancer cells through light stimulus. More specifically, we show that their signaling pathways begin with a residual re-arrangement and subsequent H-bond formation of amino acids near to the flavin-mononucleotide chromophore, causing a coupling between β-strands and subsequent detachment of a peripheral α-helix from the AsLOV2-domain. In the case of the PA-Rac1 system we find that this latter process induces the release of the AsLOV2-inhibitor from the switchII-activation site of the GTPase, enabling signal activation through effector-protein binding. These applications demonstrate that our approach reliably reproduces the signaling pathways of complex signal proteins, ranging from nanoseconds up to seconds at affordable computational costs.

Diversity of Mycoplasma hyopneumoniae in pig farms revealed by direct molecular typing of clinical material

Mycoplasma hyopneumoniae is the etiological agent of enzootic pneumonia in swine. Various reports indicate that different strains are circulating in the swine population. We investigated the variety of M. hyopneumoniae strains by a newly developed genetic typing method based on the polyserine repeat motif of the LppS homolog P146. PCR amplification using M. hyopneumoniae specific, conserved primers flanking the region encoding the repeat motif, followed by sequencing and cluster analysis was carried out. The study included strains isolated from different geographic regions as well as lysates from lung swabs from a series of pig farms in Switzerland. High diversity of M. hyopneumoniae was observed but farms being in close geographic or operative contact generally seemed to be affected by the same strains. Moreover, analysis of multiple samples from single pig farms indicated that these harbored the same, farm-specific strain. The results indicate that multiple strains of M. hyopneumoniae are found in the swine population but that specific strains or clones are responsible for local outbreaks. The method presented is a highly reproducible epidemiologic tool allowing direct typing of M. hyopneumoniae from clinical material without prior isolation and cultivation of strains.

An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120--800 ka

An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120–800 ka. In this framework, new measurements of δ18Oatm over Marine Isotope Stage (MIS) 11–12 on EDC and a complete δ18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from δ18Oatm, δO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one.

Calculation of crystal optical properties from molecular electron density

We have recently developed a method to obtain distributed atomic polarizabilities adopting a partitioning of the molecular electron density (for example, the Quantum Theory of Atoms in Molecules, [1]), calculated with or without an applied electric field. The procedure [2] allows to obtained atomic polarizability tensors, which are perfectly exportable, because quite representative of an atom in a given functional group. Among the many applications of this idea, the calculation of crystal susceptibility is easily available, either from a rough estimation (the polarizability of the isolated molecule is used) or from a more precise estimation (the polarizability of a molecule embedded in a cluster representing the first coordination sphere is used). Lorentz factor is applied to include the long range effect of packing, which is enhancing the molecular polarizability. Simple properties like linear refractive index or the gyration tensor can be calculated at relatively low costs and with good precision. This approach is particularly useful within the field of crystal engineering of organic/organometallic materials, because it would allow a relatively easy prediction of a property as a function of the packing, thus allowing "reverse crystal engineering". Examples of some amino acid crystals and salts of amino acids [3] will be illustrated, together with other crystallographic or non-crystallographic applications. For example, the induction and dispersion energies of intermolecular interactions could be calculated with superior precision (allowing anisotropic van der Waals interactions). This could allow revision of some commonly misunderstood intermolecular interactions, like the halogen bonding (see for example the recent remarks by Stone or Gilli [4]). Moreover, the chemical reactivity of coordination complexes could be reinvestigated, by coupling the conventional analysis of the electrostatic potential (useful only in the circumstances of hard nucleophilic/electrophilic interaction) with the distributed atomic polarizability. The enhanced reactivity of coordinated organic ligands would be better appreciated. [1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory. Oxford Univ. Press, 1990. [2] A. Krawczuk-Pantula, D. Pérez, K. Stadnicka, P. Macchi, Trans. Amer. Cryst. Ass. 2011, 1-25 [3] A. S. Chimpri1, M. Gryl, L. H.R. Dos Santos1, A. Krawczuk, P. Macchi Crystal Growth & Design, in the press. [4] a) A. J. Stone, J. Am. Chem. Soc. 2013, 135, 7005−7009; b) V. Bertolasi, P. Gilli, G. Gilli Crystal Growth & Design, 2013, 12, 4758-4770.

Method for the determination of specific molecular markers of biomass burning in lake sediments

Fire has an influence on regional to global atmospheric chemistry and climate. Molecular markers of biomass burning archived in lake sediments are becoming increasingly important in paleoenvironmental reconstruction and may help determine the interaction between climate and fire activity. Here, we present a high performance anion exchange chromatography–mass spectrometry method to allow separation and analysis of levoglucosan, mannosan and galactosan in lake sediments, with implications for reconstructing past biomass burning events. Determining mannosan and galactosan in Lake Kirkpatrick, New Zealand (45.03°S, 168.57°E) sediment cores and comparing these isomers with the more abundant biomass burning markers levoglucosan and charcoal represents a significant advancement in our ability to analyze past fire activity. Levoglucosan, mannosan and galactosan concentrations correlated significantly with macroscopic charcoal concentration. Levoglucosan/mannosan and levoglucosan/(mannosan + galactosan) ratios may help determine not only when fires occurred, but also if changes in the primary burned vegetation occurred.

Development of a method for fast and automatic radiocarbon measurement of aerosol samples by online coupling of an elemental analyzer with a MICADAS AMS

A fast and automatic method for radiocarbon analysis of aerosol samples is presented. This type of analysis requires high number of sample measurements of low carbon masses, but accepts precisions lower than for carbon dating analysis. The method is based on online Trapping CO2 and coupling an elemental analyzer with a MICADAS AMS by means of a gas interface. It gives similar results to a previously validated reference method for the same set of samples. This method is fast and automatic and typically provides uncertainties of 1.5–5% for representative aerosol samples. It proves to be robust and reliable and allows for overnight and unattended measurements. A constant and cross contamination correction is included, which indicates a constant contamination of 1.4 ± 0.2 μg C with 70 ± 7 pMC and a cross contamination of (0.2 ± 0.1)% from the previous sample. A Real-time online coupling version of the method was also investigated. It shows promising results for standard materials with slightly higher uncertainties than the Trapping online approach.