Computersimulation von Positronium-Annihilation in Polymeren

In dieser Arbeit wird eine neue Methode zur Berechnung derLebensdauern von Ortho-Positronium in Polymerenbeschrieben. Zur Auswertung vonPositronium-Annihilations-Experimenten wurde bislang eineinfaches Modell verwendet, dessen Gültigkeit im Falle vonPolymeren fragwürdig ist. Durch ein verbessertes Modell wirdeine realistischere Beschreibung der Lebensdauernmöglich. Mit verschiedenen quantenmechanischen Methoden wirdeine Näherung für das Wechselwirkungspotential desPositroniums mit der Polymermatrix bestimmt. DieElektronendichte im Polymer wird mit quantenchemischenRechnungen ermittelt, während die Polymerstruktur selbst mitMolekulardynamik-Simulationen auf atomistis cher Ebeneberechnet wird. Die Pfad-integral-Monte-Carlo-Methodegestattet es dann, die Aufenthaltswahrscheinlichkeit unddaraus die Lebensdauer von Ortho-Positronium in den zumfreien Volumen gehörenden Löchern im Polymer zuberechnen. Diese werden mit Ergebnissen ausPositroniu-Annihilations-Experimenten verglichen. Wenn sieübereinstimmen, läßt sich schließen, daß diezugrundeliegenden simulierten Polymerkonfigurationen inBezug auf das freie Volumen ein realistisches Abbild desrealen Polymers sind. Da die Positionen der Atome in der Simulation bekannt sind,wird das freie Volumen mit der Gittermethode bestimmt undmit einem Clusteralgorithmus zu Löchern zusammengefaßt, diedann geometrisch analysiert werden. So können Informationenüber Gestalt, Größe und Verteilung der Löcher gewonnenwerden, die weit über das bislang verwendete Standardmodellhinausgehen. Es zeigt sich, daß einige der Annahmen desStandardmodells nicht erfüllt sind. Die Möglichkeiten des neuen Verfahrens werden an Polystyrolund Bisphenol-A-Polycarbonat demonstriert. ExperimentelleLebensdauerspektren können gut reproduziert werden. Eswurden umfangreiche Analysen des freien Volumensdurchgeführt. Im Falle des Polystyrols wurde auch seineTemperaturabhängigkeit untersucht. Des weiteren werden Simulationen zur Diffusion vonKohlendioxid in Polystyrol bei verschiedenen Temperaturenund Drücken präsentiert. Die Berechnung vonDiffusionskoeffizienten in NpT-Simulationen wird diskutiert;es wird eine Methode dargestellt, die die Berechnung derVerschiebung der Teilchen in ungefalteten Koordinatenerlaubt, auch wenn die Größe der Simulationszellefluktuiert, weil eine automatische Druckkontrolle verwendetwird. Damit werden Diffusionskoeffizientenermittelt. Außerdem wird der Frage nachgegangen, welcheLöcher die Kohlendioxidmoleküle besetzen. Dies ist von derexperimentellen Seite interessant, weil häufig gasbeladenePolymere durch Positronium untersucht werden, um dieAnzahldichte der Löcher zu bestimmen; diese wird wiederumbenötigt, um aus dem mittleren Volumen eines Loches dasgesamte freie Volumen zu berechnen.

Computer simulations of two-dimensional colloidal crystals under confinement and shear

In this thesis we are presenting a broadly based computer simulation study of two-dimensional colloidal crystals under different external conditions. In order to fully understand the phenomena which occur when the system is being compressed or when the walls are being sheared, it proved necessary to study also the basic motion of the particles and the diffusion processes which occur in the case without these external forces. In the first part of this thesis we investigate the structural transition in the number of rows which occurs when the crystal is being compressed by placing the structured walls closer together. Previous attempts to locate this transition were impeded by huge hysteresis effects. We were able to determine the transition point with higher precision by applying both the Schmid-Schilling thermodynamic integration method and the phase switch Monte Carlo method in order to determine the free energies. These simulations showed not only that the phase switch method can successfully be applied to systems with a few thousand particles and a soft crystalline structure with a superimposed pattern of defects, but also that this method is way more efficient than a thermodynamic integration when free energy differences are to be calculated. Additionally, the phase switch method enabled us to distinguish between several energetically very similar structures and to determine which one of them was actually stable. Another aspect considered in the first result chapter of this thesis is the ensemble inequivalence which can be observed when the structural transition is studied in the NpT and in the NVT ensemble. The second part of this work deals with the basic motion occurring in colloidal crystals confined by structured walls. Several cases are compared where the walls are placed in different positions, thereby introducing an incommensurability into the crystalline structure. Also the movement of the solitons, which are created in the course of the structural transition, is investigated. Furthermore, we will present results showing that not only the well-known mechanism of vacancies and interstitial particles leads to diffusion in our model system, but that also cooperative ring rotation phenomena occur. In this part and the following we applied Langevin dynamics simulations. In the last chapter of this work we will present results on the effect of shear on the colloidal crystal. The shear was implemented by moving the walls with constant velocity. We have observed shear banding and, depending on the shear velocity, that the inner part of the crystal breaks into several domains with different orientations. At very high shear velocities holes are created in the structure, which originate close to the walls, but also diffuse into the inner part of the crystal.