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.

Experimental investigation into the effect of stress on dissolution and growth of very soluble brittle salts in aqueous solution

Das Studium der Auflösungs- und Wachstumsprozesse an Feststoff-Flüssigkeits-Grenzflächen unter nicht-hydrostatischen Beanspruchungen ist wesentlich für das Verständnis von Defor-mationsprozessen, die in der Erde ablaufen. Unter diesen genannten Prozessen gehört die Drucklösung zu den wichtigsten duktilen Deformationsprozessen, von der Diagenese bishin zur niedrig- bis mittelgradigen metamorphen Bedingungen. Bisher ist allerdings wenig darüber bekannt, welche mechanischen, physikalischen oder chemischen Potentialenergie-Gradienten die Drucklösung steuern. I.a. wird angenommen, daß die Drucklösung durch Un-terschiede kristallplastischer Verformungsenergien oder aber durch Unterschiede der Normal-beanspruchung an Korngrenzen gesteuert wird. Unterschiede der elastischen Verformungs-energien werden dabei allerdings als zu gering erachtet, um einen signifikanten Beitrag zu leisten. Aus diesem Grund werden sie als mögliche treibende Kräfte für die Drucklösung vernachlässigt. Andererseits haben neue experimentelle und theoretische Untersuchungen gezeigt, daß die elastische Verformung in der Tat einen starken Einfluß auf Lösungs- und Wachstumsmechanismen von Kristallen in einer Lösung haben kann. Da die in der Erdkruste vorherrschenden Deformationsmechanismen überwiegend im elastischen Verformungsbereich der Gesteine ablaufen, ist es sehr wichtig, das Verständnis für die Effekte, die die elastische Verformung verursacht, zu erweitern, und ihre Rolle während der Deformation durch Drucklösung zu definieren. Die vorliegende Arbeit beschäftigt sich mit Experimenten, bei denen der Effekt der mechanisch kompressiven Beanspruchung auf Lösungs- und Wachstumsprozesse von Einzelkristallen unterschiedlicher, sehr gut löslicher, elastisch/spröder Salze untersucht wurde. Diese Salze wurden als Analoga gesteinsbildender Minerale wie Quarz und Calcit ausgewählt. Der Einfluß von Stress auf die Ausbildung der Oberflächenmikrostrukturen in einer untersättigten Lösung wurde an Kaliumalaun untersucht.Lösungsrillen (20 – 40 µm breit, 10 – 40 µm tief und 20 – 80 µm Abstand) entwickelten sich in den Bereichen, in denen die Beanspruchung im Kristall am größten war. Sie verschwanden wieder, sobald der Kristall entlastet wurde. Diese Rillen entwickelten sich parallel zu niedrig indizierten kristallographischen Richtungen und sub-perpendikular zu den Trajektorien, die der maximalen, lokalen kompressiven Beanspruchung entsprachen. Die Größe der Lösungsrillen hing von der lokalen Oberflächenbeanspruchung, der Oberflächenenergie und dem Untersättigungsgrad der wässrigen Lösung ab. Die mikrostrukturelle Entwicklung der Kristalloberflächen stimmte gut mit den theoretischen Vorhersagen überein, die auf den Modellen von Heidug & Leroy (1994) und Leroy & Heidug (1994) basieren. Der Einfluß der Beanspruchung auf die Auflösungsrate wurde an Natriumchlorat-Einzelkristallen untersucht. Dabei wurde herausgefunden, daß sich gestresste Kristalle schneller lösen als Kristalle, auf die keine Beanspruchung einwirkt. Der experimentell beobachtete Anstieg der Auflösungsrate der gestressten Kristalle war ein bis zwei Größenordnungen höher als theoretisch erwartet. Die Auflösungsrate stieg linear mit dem Stress an, und der Anstieg war um so größer, je stärker die Lösung untersättigt war. Außerdem wurde der Effekt der Bean-spruchung auf das Kristallwachstum an Kaliumalaun- und Kaliumdihydrogenphosphat-Ein-zelkristallen untersucht. Die Wachstumsrate der Flächen {100} und {110} von Kalium-alaun war bei Beanspruchung stark reduziert. Für all diese Ergebnisse spielte die Oberflächenrauhigkeit der Kristalle eine Schlüsselrolle, indem sie eine nicht-homogene Stressverteilung auf der Kristalloberfläche verursachte. Die Resultate zeigen, daß die elastische Verformung eine signifikante Rolle während der Drucklösung spielen kann, und eine signifikante Deformation in der oberen Kruste verursachen kann, bei Beanspruchungen, die geringer sind, als gemeinhin angenommen wird. Somit folgt, daß die elastische Bean-spruchung berücksichtigt werden muß, wenn mikrophysikalische Deformationsmodelle entwickelt werden sollen.

Simulation studies of correlation functions and relaxation in polymeric systems

We investigate the statics and dynamics of a glassy,non-entangled, short bead-spring polymer melt with moleculardynamics simulations. Temperature ranges from slightlyabove the mode-coupling critical temperature to the liquidregime where features of a glassy liquid are absent. Ouraim is to work out the polymer specific effects on therelaxation and particle correlation. We find the intra-chain static structure unaffected bytemperature, it depends only on the distance of monomersalong the backbone. In contrast, the distinct inter-chainstructure shows pronounced site-dependence effects at thelength-scales of the chain and the nearest neighbordistance. There, we also find the strongest temperaturedependence which drives the glass transition. Both the siteaveraged coupling of the monomer and center of mass (CM) andthe CM-CM coupling are weak and presumably not responsiblefor a peak in the coherent relaxation time at the chain'slength scale. Chains rather emerge as soft, easilyinterpenetrating objects. Three particle correlations arewell reproduced by the convolution approximation with theexception of model dependent deviations. In the spatially heterogeneous dynamics of our system weidentify highly mobile monomers which tend to follow eachother in one-dimensional paths forming ``strings''. Thesestrings have an exponential length distribution and aregenerally short compared to the chain length. Thus, arelaxation mechanism in which neighboring mobile monomersmove along the backbone of the chain seems unlikely.However, the correlation of bonded neighbors is enhanced. When liquids are confined between two surfaces in relativesliding motion kinetic friction is observed. We study ageneric model setup by molecular dynamics simulations for awide range of sliding speeds, temperatures, loads, andlubricant coverings for simple and molecular fluids. Instabilities in the particle trajectories are identified asthe origin of kinetic friction. They lead to high particlevelocities of fluid atoms which are gradually dissipatedresulting in a friction force. In commensurate systemsfluid atoms follow continuous trajectories for sub-monolayercoverings and consequently, friction vanishes at low slidingspeeds. For incommensurate systems the velocity probabilitydistribution exhibits approximately exponential tails. Weconnect this velocity distribution to the kinetic frictionforce which reaches a constant value at low sliding speeds. This approach agrees well with the friction obtaineddirectly from simulations and explains Amontons' law on themicroscopic level. Molecular bonds in commensurate systemslead to incommensurate behavior, but do not change thequalitative behavior of incommensurate systems. However,crossed chains form stable load bearing asperities whichstrongly increase friction.