Partial reconstruction of the trajectories of a discretely observed branching diffusion with immigration and an application to inference

In this treatise we consider finite systems of branching particles where the particles move independently of each other according to d-dimensional diffusions. Particles are killed at a position dependent rate, leaving at their death position a random number of descendants according to a position dependent reproduction law. In addition particles immigrate at constant rate (one immigrant per immigration time). A process with above properties is called a branching diffusion withimmigration (BDI). In the first part we present the model in detail and discuss the properties of the BDI under our basic assumptions. In the second part we consider the problem of reconstruction of the trajectory of a BDI from discrete observations. We observe positions of the particles at discrete times; in particular we assume that we have no information about the pedigree of the particles. A natural question arises if we want to apply statistical procedures on the discrete observations: How can we find couples of particle positions which belong to the same particle? We give an easy to implement 'reconstruction scheme' which allows us to redraw or 'reconstruct' parts of the trajectory of the BDI with high accuracy. Moreover asymptotically the whole path can be reconstructed. Further we present simulations which show that our partial reconstruction rule is tractable in practice. In the third part we study how the partial reconstruction rule fits into statistical applications. As an extensive example we present a nonparametric estimator for the diffusion coefficient of a BDI where the particles move according to one-dimensional diffusions. This estimator is based on the Nadaraya-Watson estimator for the diffusion coefficient of one-dimensional diffusions and it uses the partial reconstruction rule developed in the second part above. We are able to prove a rate of convergence of this estimator and finally we present simulations which show that the estimator works well even if we leave our set of assumptions.

Pathomechanische Analysen zum Untergang retinaler Ganglienzellen im experimentellen autoimmunen Glaukom-Modell

In den vergangenen Jahren konnten zahlreiche Studien die Veränderung des natürlichen Autoantikörperrepertoirs bei Glaukompatienten aufzeigen. Zu den Antigenen zählen verschiedenen Hitzeschockproteinen, aber ebenso neuronal assoziierte Strukturproteine wie das Myelin basische Protein (MBP) oder das sauren Gliafaserprotein und einige neuropyhsiologische Proteine aus der Retina und dem Sehnerven. Da bei den Glaukompatienten nicht einzelne Antikörperreaktionen verändert sind, sondern vielmehr komplexe Autoantikörpermuster vorliegen, bestand das primäre Ziel der Dissertation zu zeigen, ob eine systemische Immunisierung mit MBP, Homogenaten opticus-assoziierter Antigene (ONA) und Antigenen der retinalen Ganglienzellschicht (RGA) den Verlust von retinalen Ganglienzellen (RGZ) in einem Experimentellen Autoimmunen Glaukom (EAG) Tiermodell auslösen können. Die systemische Injektion von MBP, ONA oder RGA induzierten ophthalmopathologische Veränderungen in der Retina, gekennzeichnet durch retinalen Ganglienzellverlust mitsamt Zerstörung der Axone im Sehnerv. Unter der Annahme, dass die Neurodegeneration durch Autoantiköper vermittelt ist, wurde ebenfalls untersucht, ob sich die Antikörperreaktivität gegen okulare Strukturen oder den Sehnerv im Verlauf der Studie verändern. Getestet wurde die Antikörperreaktivität gegen Gewebsschnitte gesunder Tiere mit dem Ergebnis einer signifikanten und zeitabhängigen Zunahme der Immunreaktivität. Darüber hinaus war es erstmals möglich die Ablagerung von IgG Autoantikörpern in der Retina und dem Sehnerv nachzuweisen sowie die Caspase mediierte Apoptose zu untersuchen. Ebenfalls konnte die Verteilung von aktivierten Mikroglia im optischen System evaluiert werden, wobei diese mehrmals in Kolokalisation mit den IgG-Autoantikörpern auftraten. Diese Beobachtungen lassen den Schluss zu, dass die Immunreaktionen von Autoantikörpern alleine und im Zusammenspiel mit der Mikroglia im Zusammenhang mit der Neurodegeneration der retinalen Ganglienzelle im EAG Modell stehen könnten.

Atomic-scale characterization of diamond surfaces and fullerene self-assembly

In this thesis, elemental research towards the implantation of a diamond-based molecular quantum computer is presented. The approach followed requires linear alignment of endohedral fullerenes on the diamond C(100) surface in the vicinity of subsurface NV-centers. From this, four fundamental experimental challenges arise: 1) The well-controlled deposition of endohedral fullerenes on a diamond surface. 2) The creation of NV-centers in diamond close to the surface. 3) Preparation and characterization of atomically-flat diamondsurfaces. 4) Assembly of linear chains of endohedral fullerenes. First steps to overcome all these challenges were taken in the framework of this thesis. Therefore, a so-called “pulse injection” technique was implemented and tested in a UHV chamber that was custom-designed for this and further tasks. Pulse injection in principle allows for the deposition of molecules from solution onto a substrate and can therefore be used to deposit molecular species that are not stable to sublimation under UHV conditions, such as the endohedral fullerenes needed for a quantum register. Regarding the targeted creation of NV-centers, FIB experiments were carried out in cooperation with the group of Prof. Schmidt-Kaler (AG Quantum, Physics Department, Johannes Gutenberg-Universität Mainz). As an entry into this challenging task, argon cations were implanted into (111) surface-oriented CaF2 crystals. The resulting implantation spots on the surface were imaged and characterized using AFM. In this context, general relations between the impact of the ions on the surface and their valency or kinetic energy, respectively, could be established. The main part of this thesis, however, is constituted by NCAFM studies on both, bare and hydrogen-terminated diamond C(100) surfaces. In cooperation with the group of Prof. Dujardin (Molecular Nanoscience Group, ISMO, Université de Paris XI), clean and atomically-flat diamond surfaces were prepared by exposure of the substrate to a microwave hydrogen plasma. Subsequently, both surface modifications were imaged in high resolution with NC-AFM. In the process, both hydrogen atoms in the unit cell of the hydrogenated surface were resolved individually, which was not achieved in previous STM studies of this surface. The NC-AFM images also reveal, for the first time, atomic-resolution contrast on the clean, insulating diamond surface and provide real-space experimental evidence for a (2×1) surface reconstruction. With regard to the quantum computing concept, high-resolution NC-AFM imaging was also used to study the adsorption and self-assembly potential of two different kinds of fullerenes (C60 and C60F48) on aforementioned diamond surfaces. In case of the hydrogenated surface, particular attention was paid to the influence of charge transfer doping on the fullerene-substrate interaction and the morphology emerging from self-assembly. Finally, self-assembled C60 islands on the hydrogen-terminated diamond surface were subject to active manipulation by an NC-AFM tip. Two different kinds of tip-induced island growth modes have been induced and were presented. In conclusion, the results obtained provide fundamental informations mandatory for the realization of a molecular quantum computer. In the process it was shown that NC-AFM is, under proper circumstances, a very capable tool for imaging diamond surfaces with highest resolution, surpassing even what has been achieved with STM up to now. Particular attention was paid to the influence of transfer doping on the morphology of fullerenes on the hydrogenated diamond surface, revealing new possibilities for tailoring the self-assembly of molecules that have a high electron affinity.

Chiral properties of two-flavour QCD at zero and non-zero temperature

Lattice Quantum Chromodynamics (LQCD) is the preferred tool for obtaining non-perturbative results from QCD in the low-energy regime. It has by nowrnentered the era in which high precision calculations for a number of phenomenologically relevant observables at the physical point, with dynamical quark degrees of freedom and controlled systematics, become feasible. Despite these successes there are still quantities where control of systematic effects is insufficient. The subject of this thesis is the exploration of the potential of todays state-of-the-art simulation algorithms for non-perturbativelyrn$\mathcal{O}(a)$-improved Wilson fermions to produce reliable results in thernchiral regime and at the physical point both for zero and non-zero temperature. Important in this context is the control over the chiral extrapolation. Thisrnthesis is concerned with two particular topics, namely the computation of hadronic form factors at zero temperature, and the properties of the phaserntransition in the chiral limit of two-flavour QCD.rnrnThe electromagnetic iso-vector form factor of the pion provides a platform to study systematic effects and the chiral extrapolation for observables connected to the structure of mesons (and baryons). Mesonic form factors are computationally simpler than their baryonic counterparts but share most of the systematic effects. This thesis contains a comprehensive study of the form factor in the regime of low momentum transfer $q^2$, where the form factor is connected to the charge radius of the pion. A particular emphasis is on the region very close to $q^2=0$ which has not been explored so far, neither in experiment nor in LQCD. The results for the form factor close the gap between the smallest spacelike $q^2$-value available so far and $q^2=0$, and reach an unprecedented accuracy at full control over the main systematic effects. This enables the model-independent extraction of the pion charge radius. The results for the form factor and the charge radius are used to test chiral perturbation theory ($\chi$PT) and are thereby extrapolated to the physical point and the continuum. The final result in units of the hadronic radius $r_0$ is rn$$ \left\langle r_\pi^2 \right\rangle^{\rm phys}/r_0^2 = 1.87 \: \left(^{+12}_{-10}\right)\left(^{+\:4}_{-15}\right) \quad \textnormal{or} \quad \left\langle r_\pi^2 \right\rangle^{\rm phys} = 0.473 \: \left(^{+30}_{-26}\right)\left(^{+10}_{-38}\right)(10) \: \textnormal{fm} \;, $$rn which agrees well with the results from other measurements in LQCD and experiment. Note, that this is the first continuum extrapolated result for the charge radius from LQCD which has been extracted from measurements of the form factor in the region of small $q^2$.rnrnThe order of the phase transition in the chiral limit of two-flavour QCD and the associated transition temperature are the last unkown features of the phase diagram at zero chemical potential. The two possible scenarios are a second order transition in the $O(4)$-universality class or a first order transition. Since direct simulations in the chiral limit are not possible the transition can only be investigated by simulating at non-zero quark mass with a subsequent chiral extrapolation, guided by the universal scaling in the vicinity of the critical point. The thesis presents the setup and first results from a study on this topic. The study provides the ideal platform to test the potential and limits of todays simulation algorithms at finite temperature. The results from a first scan at a constant zero-temperature pion mass of about 290~MeV are promising, and it appears that simulations down to physical quark masses are feasible. Of particular relevance for the order of the chiral transition is the strength of the anomalous breaking of the $U_A(1)$ symmetry at the transition point. It can be studied by looking at the degeneracies of the correlation functions in scalar and pseudoscalar channels. For the temperature scan reported in this thesis the breaking is still pronounced in the transition region and the symmetry becomes effectively restored only above $1.16\:T_C$. The thesis also provides an extensive outline of research perspectives and includes a generalisation of the standard multi-histogram method to explicitly $\beta$-dependent fermion actions.

Molecular genetic investigation of the Neolithic population history in the western Carpathian Basin

Der Fokus dieser Dissertation ist die populationsgenetische Analyse der neolithischen Bevölkerungswechsel in den 6.-5. Jahrtausende vor Christus, die im westlichen Karpatenbecken stattfanden. Die Zielsetzung der Studie war, mittels der Analyse von mitochondrialer und Y-chromosomaler aDNA, den Genpool der sechs neolithischen und kupferzeitlichen Populationen zu untersuchen und die daraus resultierenden Ergebnisse mit anderen prähistorischen und modernen genetischen Daten zu vergleichen.rnInsgesamt wurden 323 Individuen aus 32 ungarischen, kroatischen und slowakischen Fundplätzen beprobt und bearbeitet in den archäogenetischen Laboren der Johannes Gutenberg-Universität in Mainz. Die DNA Ergebnisse wurden mit verschiedenen populationsgenetischen Methoden ausgewertet. Vergleichsdaten von prähistorischen und modernen eurasiatischen Populationen wurden dazu gesammelt.rnDie HVS-I Region der mitochondrialen DNA konnten bei 256 Individuen reproduziert und authentifiziert werden (mit einer Erfolgsrate von 85.9%). Die Typisierung der HVS-II Region war in 80 Fällen erfolgreich. Testend alle gut erhaltene Proben, die Y-chromosomale Haplogruppe konnte in 33 männlichen Individuen typisiert werden.rnDie neolithischen, mitochondrialen Haplogruppen deuten auf eine hohe Variabilität des maternalen Genpools hin. Sowohl die mitochondrialen als auch die Y-chromosomalen Daten lassen Rückschlüsse auf eine nah-östliche bzw. südwestasiatische Herkunft der frühen Bauern zu. Die Starčevo- und linearbandkermaischen-Populationen in westlichem Karpatenbecken (letztere abgekürzt als LBKT) und die linearbandkermaischen-Population in Mitteleuropa (LBK) haben so starke genetische Ähnlichkeit, dass die Verbreitung der LBK nach Mitteleuropa mit vorangegangenen Wanderungsereignissen zu erklären ist. Die Transdanubische aDNA Daten zeigen hohe Affinität zu den publizierten prähistorischen aDNA Datensätzen von Mitteleuropa aus den 6.-4. Jahrtausende vor Chr. Die maternal-genetische Variabilität der Starčevo-Population konnte auch innerhalb der nachfolgenden Populationen Transdanubiens festgestellt werden. Nur kleinere Infiltrationen und Immigrationsereignissen konnten während der Vinča-, LBKT-, Sopot- und Balaton-Lasinja-Kultur in Transdanubien identifiziert werden. Zwischen den transdanubischen Regionen konnten mögliche genetische Unterschiede nur in der LBKT und in der Lengyel-Periode beobachtet werden, als sich die nördlichen Gruppen von den südlichen Populationen trennten. rnDie festgestellte Heterogenität der mtDNA in Zusammenhang mit der Y-chromosomalen Homogenität in den Starčevo- und LBK-Populationen, weisen auf patrilokale Residenzregeln und patrilineare Abstammungsregeln in den ersten Bauergemeinschaften hin. rnObwohl die hier präsentierten Daten einen großen Fortschritt in der Forschung von aDNA und Neolithikum des Karpatenbeckens und Mitteleuropas bedeuten, werfen sie auch mehrere Fragen auf, deren Beantwortung durch zukünftige Genomforschungen erbracht werden könnte.