Development of a multi-species method by GC-coupling with inductively coupled plasma isotope dilution mass spectrometry for the simultaneous determination of alkylated lead, mercury and tin compounds

An accurate and sensitive species-specific GC-ICP-IDMS (gas chromatography inductively coupled plasma isotope dilution mass spectrometry) method for the determination of trimethyllead and a multi-species-specific GC-ICP-IDMS method for the simultaneous determination of trimethyllead, methylmercury, and butyltins in biological and environmental samples were developed. They allow the determination of corresponding elemental species down to the low ng g-1 range. The developed synthesis scheme for the formation of isotopically labeled Me3206Pb+ can be used for future production of this spike. The novel extraction technique, stir bar sorptive extraction (SBSE), was applied for the first time in connection with species-specific isotope dilution GC-ICP-MS for the determination of trimethyllead, methylmercury and butyltins. The results were compared with liquid-liquid extraction. The developed methods were validated by the analysis of certified reference materials. The liquid-liquid extraction GC-ICP-IDMS method was applied to seafood samples purchased from a supermarket. The methylated lead fraction in these samples, correlated to total lead, varied in a broad range of 0.01-7.6 %. On the contrary, the fraction of methylmercury is much higher, normally in the range of 80-98 %. The highest methylmercury content of up to 12 µg g-1 has been determined in shark samples, an animal which is at the end of the marine food chain, whereas in other seafood samples a MeHg+ content of less than 0.2 µg g-1 was found. Butyltin species could only be determined in samples, where anthropogenic contaminations must be assumed. This explains the observed broad variation of the butylated tin fraction in the range of <0.3-49 % in different seafood samples. Because all isotope-labelled spike compounds, except trimethyllead, are commercially available, the developed multi-species-specific GC-ICP-IDMS method has a high potential in future for routine analysis.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

Numerical coupling of thermal-electric network models and energy-transport equations including optoelectronic semiconductor devices

In this work the numerical coupling of thermal and electric network models with model equations for optoelectronic semiconductor devices is presented. Modified nodal analysis (MNA) is applied to model electric networks. Thermal effects are modeled by an accompanying thermal network. Semiconductor devices are modeled by the energy-transport model, that allows for thermal effects. The energy-transport model is expandend to a model for optoelectronic semiconductor devices. The temperature of the crystal lattice of the semiconductor devices is modeled by the heat flow eqaution. The corresponding heat source term is derived under thermodynamical and phenomenological considerations of energy fluxes. The energy-transport model is coupled directly into the network equations and the heat flow equation for the lattice temperature is coupled directly into the accompanying thermal network. The coupled thermal-electric network-device model results in a system of partial differential-algebraic equations (PDAE). Numerical examples are presented for the coupling of network- and one-dimensional semiconductor equations. Hybridized mixed finite elements are applied for the space discretization of the semiconductor equations. Backward difference formluas are applied for time discretization. Thus, positivity of charge carrier densities and continuity of the current density is guaranteed even for the coupled model.

Computational Modeling of Cardiac Excitation-Contraction Coupling in Physiological and Pathological Conditions

The cardiomyocyte is a complex biological system where many mechanisms interact non-linearly to regulate the coupling between electrical excitation and mechanical contraction. For this reason, the development of mathematical models is fundamental in the field of cardiac electrophysiology, where the use of computational tools has become complementary to the classical experimentation. My doctoral research has been focusing on the development of such models for investigating the regulation of ventricular excitation-contraction coupling at the single cell level. In particular, the following researches are presented in this thesis: 1) Study of the unexpected deleterious effect of a Na channel blocker on a long QT syndrome type 3 patient. Experimental results were used to tune a Na current model that recapitulates the effect of the mutation and the treatment, in order to investigate how these influence the human action potential. Our research suggested that the analysis of the clinical phenotype is not sufficient for recommending drugs to patients carrying mutations with undefined electrophysiological properties. 2) Development of a model of L-type Ca channel inactivation in rabbit myocytes to faithfully reproduce the relative roles of voltage- and Ca-dependent inactivation. The model was applied to the analysis of Ca current inactivation kinetics during normal and abnormal repolarization, and predicts arrhythmogenic activity when inhibiting Ca-dependent inactivation, which is the predominant mechanism in physiological conditions. 3) Analysis of the arrhythmogenic consequences of the crosstalk between β-adrenergic and Ca-calmodulin dependent protein kinase signaling pathways. The descriptions of the two regulatory mechanisms, both enhanced in heart failure, were integrated into a novel murine action potential model to investigate how they concur to the development of cardiac arrhythmias. These studies show how mathematical modeling is suitable to provide new insights into the mechanisms underlying cardiac excitation-contraction coupling and arrhythmogenesis.

Coarse-graining and quantum-classical adaptive coupling in soft matter

Die vorliegende Arbeit untersucht den Zusammenhang zwischen Skalen in Systemen weicher Materie, der für Multiskalen-Simulationen eine wichtige Rolle spielt. Zu diesem Zweck wurde eine Methode entwickelt, die die Approximation der Separierbarkeit von Variablen für die Molekulardynamik und ähnliche Anwendungen bewertet. Der zweite und größere Teil dieser Arbeit beschäftigt sich mit der konzeptionellen und technischen Erweiterung des Adaptive Resolution Scheme'' (AdResS), einer Methode zur gleichzeitigen Simulation von Systemen mit mehreren Auflösungsebenen. Diese Methode wurde auf Systeme erweitert, in denen klassische und quantenmechanische Effekte eine Rolle spielen.rnrnDie oben genannte erste Methode benötigt nur die analytische Form der Potentiale, wie sie die meisten Molekulardynamik-Programme zur Verfügung stellen. Die Anwendung der Methode auf ein spezielles Problem gibt bei erfolgreichem Ausgang einen numerischen Hinweis auf die Gültigkeit der Variablenseparation. Bei nicht erfolgreichem Ausgang garantiert sie, dass keine Separation der Variablen möglich ist. Die Methode wird exemplarisch auf ein zweiatomiges Molekül auf einer Oberfläche und für die zweidimensionale Version des Rotational Isomer State (RIS) Modells einer Polymerkette angewandt.rnrnDer zweite Teil der Arbeit behandelt die Entwicklung eines Algorithmus zur adaptiven Simulation von Systemen, in denen Quanteneffekte berücksichtigt werden. Die Quantennatur von Atomen wird dabei in der Pfadintegral-Methode durch einen klassischen Polymerring repräsentiert. Die adaptive Pfadintegral-Methode wird zunächst für einatomige Flüssigkeiten und tetraedrische Moleküle unter normalen thermodynamischen Bedingungen getestet. Schließlich wird die Stabilität der Methode durch ihre Anwendung auf flüssigen para-Wasserstoff bei niedrigen Temperaturen geprüft.

Kinematic models of interseismic deformation from inversion of GPS and InSAR measurements to estimate fault parameters and coupling degree

We have used kinematic models in two Italian regions to reproduce surface interseismic velocities obtained from InSAR and GPS measurements. We have considered a Block modeling, BM, approach to evaluate which fault system is actively accommodating the occurring deformation in both considered areas. We have performed a study for the Umbria-Marche Apennines, obtaining that the tectonic extension observed by GPS measurements is explained by the active contribution of at least two fault systems, one of which is the Alto Tiberina fault, ATF. We have estimated also the interseismic coupling distribution for the ATF using a 3D surface and the result shows an interesting correlation between the microseismicity and the uncoupled fault portions. The second area analyzed concerns the Gargano promontory for which we have used jointly the available InSAR and GPS velocities. Firstly we have attached the two datasets to the same terrestrial reference frame and then using a simple dislocation approach, we have estimated the best fault parameters reproducing the available data, providing a solution corresponding to the Mattinata fault. Subsequently we have considered within a BM analysis both GPS and InSAR datasets in order to evaluate if the Mattinata fault may accommodate the deformation occurring in the central Adriatic due to the relative motion between the North-Adriatic and South-Adriatic plates. We obtain that the deformation occurring in that region should be accommodated by more that one fault system, that is however difficult to detect since the poor coverage of geodetic measurement offshore of the Gargano promontory. Finally we have performed also the estimate of the interseismic coupling distribution for the Mattinata fault, obtaining a shallow coupling pattern. Both of coupling distributions found using the BM approach have been tested by means of resolution checkerboard tests and they demonstrate that the coupling patterns depend on the geodetic data positions.

Mode-coupling theory: generalizations, high dimensions and microscopic dynamics

This work contains several applications of the mode-coupling theory (MCT) and is separated into three parts. In the first part we investigate the liquid-glass transition of hard spheres for dimensions d→∞ analytically and numerically up to d=800 in the framework of MCT. We find that the critical packing fraction ϕc(d) scales as d²2^(-d), which is larger than the Kauzmann packing fraction ϕK(d) found by a small-cage expansion by Parisi and Zamponi [J. Stat. Mech.: Theory Exp. 2006, P03017 (2006)]. The scaling of the critical packing fraction is different from the relation ϕc(d)∼d2^(-d) found earlier by Kirkpatrick and Wolynes [Phys. Rev. A 35, 3072 (1987)]. This is due to the fact that the k dependence of the critical collective and self nonergodicity parameters fc(k;d) and fcs(k;d) was assumed to be Gaussian in the previous theories. We show that in MCT this is not the case. Instead fc(k;d) and fcs(k;d), which become identical in the limit d→∞, converge to a non-Gaussian master function on the scale k∼d^(3/2). We find that the numerically determined value for the exponent parameter λ and therefore also the critical exponents a and b depend on the dimension d, even at the largest evaluated dimension d=800. In the second part we compare the results of a molecular-dynamics simulation of liquid Lennard-Jones argon far away from the glass transition [D. Levesque, L. Verlet, and J. Kurkijärvi, Phys. Rev. A 7, 1690 (1973)] with MCT. We show that the agreement between theory and computer simulation can be improved by taking binary collisions into account [L. Sjögren, Phys. Rev. A 22, 2866 (1980)]. We find that an empiric prefactor of the memory function of the original MCT equations leads to similar results. In the third part we derive the equations for a mode-coupling theory for the spherical components of the stress tensor. Unfortunately it turns out that they are too complex to be solved numerically.

Uptake mechanism, intracellular trafficking and endo-lysosomal pH monitoring of polystyrene nanoparticles

What is the intracellular fate of nanoparticles (NPs) taken up by the cells? This question has been investigated for polystyrene NPs of different sizes with a set of molecular biological and biophysical techniques.rnTwo sets of fluorescent NPs, cationic and non-ionic, were synthesized with three different polymerization techniques. Non-ionic particles (132 – 846 nm) were synthesized with dispersion polymerization in an ethanol/water solution. Cationic NPs with 120 nm were synthesized by miniemulsion polymerization Particles with 208, 267 and 603 nm were produced by seeding the 120 nm particle obtained by miniemulsion polymerization with drop-wise added monomer and polymerization of such. The colloidal characterization of all particles showed a comparable amount of the surface groups. In addition, particles were characterized with regard to their size, morphology, solid content, amount of incorporated fluorescent dye and zeta potential. The fluorescent intensities of all particles were measured by fluorescence spectroscopy for calibration in further cellular experiments. rnThe uptake of the NPs to HeLa cells after 1 – 24 h revealed a much higher uptake of cationic NPs in comparison to non-ionic NPs. If the same amount of NPs with different sizes is introduced to the cell, a different amount of particles is present in the cell medium, which complicates a comparison of the uptake. The same conclusion is valid for the particles’ overall surface area. Therefore, HeLa cells were incubated with the same concentration, amount and surface area of NPs. It was found that with the same concentration always the same polymer amount is taking up by cells. However, the amount of particles taken up decreases for the biggest. A correlation to the surface area could not be found. We conclude that particles are endocytosed by an excavator-shovel like mechanism, which does not distinguish between different sizes, but is only dependent on the volume that is taken up. For the decreased amount of large particles, an overload of this mechanism was assumed, which leads to a decrease in the uptake. rnThe participation of specific endocytotic processes has been determined by the use of pharmacological inhibitors, immunocytological staining and immunofluorescence. The uptake of NPs into the endo-lysosomal machinery is dominated by a caveolin-mediated endocytosis. Other pathways, which include macropinocytosis and a dynamin-dependent mechanism but exclude clathrin mediated endocytosis, also occur as competing processes. All particles can be found to some extent in early endosomes, but only bigger particles were proven to localize in late endosomes. No particles were found in lysosomes; at least not in lysosomes that are labeled with Lamp1 and cathepsin D. However, based on the character of the performed experiment, a localization of particles in lysosomes cannot be excluded.rnDuring their ripening process, vesicles undergo a gradual acidification from early over late endosomes to lysosomes. It is hypothesized that NPs in endo-lysosomal compartments experience the same change in pH value. To probe the environmental pH of NPs after endocytosis, the pH-sensitive dye SNARF-4F was grafted onto amino functionalized polystyrene NPs. The pH value is a ratio function of the two emission wavelengths of the protonated and deprotonated form of the dye and is hence independent of concentration changes. The particles were synthesized by the aforementioned miniemulsion polymerization with the addition of the amino functionalized copolymer AEMH. The immobilization of SNARF-4F was performed by an EDC-coupling reaction. The amount of physically adsorbed dye in comparison to covalently bonded dye was 15% as determined by precipitation of the NPs in methanol, which is a very good solvent for SNARF-4F. To determine influences of cellular proteins on the fluorescence properties, a intracellular calibration fit was established with platereader measurements and cLSM imaging by the cell-penetrable SNARF-4F AM ester. Ionophores equilibrated the extracellular and intracellular pH.rnSNARF-4F NPs were taken up well by HeLa cells and showed no toxic effects. The pH environment of SNARF-4F NPs has been qualitatively imaged as a movie over a time period up to 1 h in pseudo-colors by a self-written automated batch program. Quantification revealed an acidification process until pH value of 4.5 over 24 h, which is much slower than the transport of nutrients to lysosomes. NPs are present in early endosomes after min. 1 h, in late endosomes at approx. 8 h and end up in vesicles with a pH value typical for lysosomes after > 24 h. We therefore assume that NPs bear a unique endocytotic mechanism, at least with regards to the kinetic involvedrn

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.

Synthesis, characterization and catalytic performances of ruthenium-based catalysts for the acceptorless dehydrogenative coupling of butanol

A growing interest towards new sources of energy has led in recent years to the development of a new generation of catalysts for alcohol dehydrogenative coupling (ADC). This green, atom-efficient reaction is capable of turning alcohol derivatives into higher value and chemically more attractive ester molecules, and it finds interesting applications in the transformation of the large variety of products deriving from biomass. In the present work, a new series of ruthenium-PNP pincer complexes are investigated for the transformation of 1-butanol, one of the most challenging substrates for this type of reactions, into butyl butyrate, a short-chain symmetrical ester widely used in flavor industries. Since the reaction kinetics depends on hydrogen diffusion, the study aimed at identifying proper reactor type and right catalyst concentration to avoid mass transfer interferences and to get dependable data. A comparison between catalytic activities and productivities has been made to establish the role of the different ligands bonded both to the PNP binder and to the ruthenium metal center, and hence to find the best catalyst for this type of reaction.

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

The bulk transition of many-flavour QCD and the search for a UVFP at strong coupling

We explore the nature of the bulk transition observed at strong coupling in the SU(3) gauge theory with Nf=12 fermions in the fundamental representation. The transition separates a weak coupling chirally symmetric phase from a strong coupling chirally broken phase and is compatible with the scenario where conformality is restored by increasing the flavour content of a non abelian gauge theory. We explore the intriguing possibility that the observed bulk transition is associated with the occurrence of an ultraviolet fixed point (UVFP) at strong coupling, where a new theory emerges in the continuum.