Computer simulations of two-dimensional colloidal crystals in confinement

Monte Carlo simulations are used to study the effect of confinement on a crystal of point particles interacting with an inverse power law potential in d=2 dimensions. This system can describe colloidal particles at the air-water interface, a model system for experimental study of two-dimensional melting. It is shown that the state of the system (a strip of width D) depends very sensitively on the precise boundary conditions at the two ``walls'' providing the confinement. If one uses a corrugated boundary commensurate with the order of the bulk triangular crystalline structure, both orientational order and positional order is enhanced, and such surface-induced order persists near the boundaries also at temperatures where the system in the bulk is in its fluid state. However, using smooth repulsive boundaries as walls providing the confinement, only the orientational order is enhanced, but positional (quasi-) long range order is destroyed: The mean-square displacement of two particles n lattice parameters apart in the y-direction along the walls then crosses over from the logarithmic increase (characteristic for $d=2$) to a linear increase (characteristic for d=1). The strip then exhibits a vanishing shear modulus. These results are interpreted in terms of a phenomenological harmonic theory. Also the effect of incommensurability of the strip width D with the triangular lattice structure is discussed, and a comparison with surface effects on phase transitions in simple Ising- and XY-models is made

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.

Mineral chemistry and structural relationships of inclusions in diamond samples

Diamant ist das härteste Mineral – und dazu ein Edelstein -, das unter höchstem Druck und hohen Temperaturen in tiefen kontinentalen Regionen der Erde kristallisiert. Die Mineraleinschlüsse in Diamanten werden durch die physikalische Stabilität und chemische Beständigkeit der umgebenden – eigentlich metastabilen -Diamant-Phase geschützt. Aufgrund der koexistierenden Phasenkombination ermöglichen sie, die Mineral-Entwicklung zu studieren, während deren der Einschlüssen und die Diamanten kristallisierten. rnDie Phasenkombinationen von Diamant und Chrom-Pyrop, Chrom-Diopsid, Chromit, Olivin, Graphit und Enstatit nebeneinander (teilweise in Berührungsexistenz) mit Chrom-Pyrop Einschlüssen wurden von neunundzwanzig Diamant-Proben von sechs Standorten in Südafrika (Premier, Koffiefontein, De Beers Pool, Finsch, Venetia und Koingnaas Minen) und Udachnaya (Sibirien/Russland) identifiziert und charakterisiert. Die Mineraleinschlüsse weisen z.T. kubo-oktaedrische Form auf, die unabhängig von ihren eigenen Kristallsystemen ausgebildet werden können. Das bedeutet, dass sie syngenetische Einschlüsse sind, die durch die sehr hohe Formenergie des umgebenden Diamanten morphologisch unter Zwang stehen. Aus zweidiemnsionalen Messungen der ersten Ordnung von charakteristischen Raman-Banden lassen sich relative Restdrucke in Diamanten zwischen Diamant und Einschlussmineral gewinnen; sie haben charakteristische Werte von ca. 0,4 bis 0,9 GPa um Chrom-Pyrop-Einschlüsse, 0,6 bis 2,0 GPa um Chrom-Diopsid-Einschlüsse, 0,3 bis 1,2 GPa um Olivin-Einschlüsse, 0,2 bis 1,0 GPa um Chromit-Einschlüsse, beziehungsweise 0,5 GPa um Graphit Einschlüsse.rnDie kristallstrukturellen Beziehung von Diamanten und ihren monomineralischen Einschlüssen wurden mit Hilfe der Quantifizierung der Winkelkorrelationen zwischen der [111] Richtung von Diamanten und spezifisch ausgewählten Richtungen ihrer mineralischen Einschlüsse untersucht. Die Winkelkorrelationen zwischen Diamant [111] und Chrom-Pyrop [111] oder Chromit [111] zeigen die kleinsten Verzerrungen von 2,2 bis zu 3,4. Die Chrom-Diopsid- und Olivin-Einschlüsse zeigen die Missorientierungswerte mit Diamant [111] bis zu 10,2 und 12,9 von Chrom-Diopsid [010] beziehungsweise Olivin [100].rnDie chemische Zusammensetzung von neun herausgearbeiteten (orientiertes Anschleifen) Einschlüssen (drei Chrom-Pyrop-Einschlüsse von Koffiefontein-, Finsch- und Venetia-Mine (zwei von drei koexistieren nebeneinander mit Enstatit), ein Chromit von Udachnaya (Sibirien/Russland), drei Chrom-Diopside von Koffiefontein, Koingnaas und Udachnaya (Sibirien/Russland) und zwei Olivin Einschlüsse von De Beers Pool und Koingnaas) wurden mit Hilfe EPMA und LA-ICP-MS analysiert. Auf der Grundlage der chemischen Zusammensetzung können die Mineraleinschlüsse in Diamanten in dieser Arbeit der peridotitischen Suite zugeordnet werden.rnDie Geothermobarometrie-Untersuchungen waren aufgrund der berührenden Koexistenz von Chrom-Pyrop- und Enstatit in einzelnen Diamanten möglich. Durchschnittliche Temperaturen und Drücke der Bildung sind mit ca. 1087 (± 15) C, 5,2 (± 0,1) GPa für Diamant DHK6.2 von der Koffiefontein Mine beziehungsweise ca. 1041 (± 5) C, 5,0 (± 0,1) GPa für Diamant DHF10.2 von der Finsch Mine zu interpretieren.rn

Ein inverses Rückstreuproblem der elektrischen Impedanztomographie

Die vorliegende Arbeit untersucht das inverse Hindernisproblem der zweidimensionalen elektrischen Impedanztomographie (EIT) mit Rückstreudaten. Wir präsentieren und analysieren das mathematische Modell für Rückstreudaten, diskutieren das inverse Problem für einen einzelnen isolierenden oder perfekt leitenden Einschluss und stellen zwei Rekonstruktionsverfahren für das inverse Hindernisproblem mit Rückstreudaten vor. Ziel des inversen Hindernisproblems der EIT ist es, Inhomogenitäten (sogenannte Einschlüsse) der elektrischen Leitfähigkeit eines Körpers aus Strom-Spannungs-Messungen an der Körperoberfläche zu identifizieren. Für die Messung von Rückstreudaten ist dafür nur ein Paar aus an der Körperoberfläche nahe zueinander angebrachten Elektroden nötig, das zur Datenerfassung auf der Oberfläche entlang bewegt wird. Wir stellen ein mathematisches Modell für Rückstreudaten vor und zeigen, dass Rückstreudaten die Randwerte einer außerhalb der Einschlüsse holomorphen Funktion sind. Auf dieser Grundlage entwickeln wir das Konzept des konvexen Rückstreuträgers: Der konvexe Rückstreuträger ist eine Teilmenge der konvexen Hülle der Einschlüsse und kann daher zu deren Auffindung dienen. Wir stellen einen Algorithmus zur Berechnung des konvexen Rückstreuträgers vor und demonstrieren ihn an numerischen Beispielen. Ferner zeigen wir, dass ein einzelner isolierender Einschluss anhand seiner Rückstreudaten eindeutig identifizierbar ist. Der Beweis dazu beruht auf dem Riemann'schen Abbildungssatz für zweifach zusammenhängende Gebiete und dient als Grundlage für einen Rekonstruktionsalgorithmus, dessen Leistungsfähigkeit wir an verschiedenen Beispielen demonstrieren. Ein perfekt leitender Einschluss ist hingegen nicht immer aus seinen Rückstreudaten rekonstruierbar. Wir diskutieren, in welchen Fällen die eindeutige Identifizierung fehlschlägt und zeigen Beispiele für unterschiedliche perfekt leitende Einschlüsse mit gleichen Rückstreudaten.