Monte-Carlo simulation of hard spherocylinders under confinement

When a liquid crystal is confined to a cavity its director field becomes subject to competing forces: on the one hand, the surface of the cavity orients the director field (``surface anchoring''), on the other hand deformations of the director field cost elastic energy. Hence the equilibrium director field is determined by a compromise between surface anchoring and elasticity. One example of a confined liquid crystal that has attracted particular interest from physicists is the nematic droplet. In this thesis a system of hard rods is considered as the simplest model for nematic liquid crystals consisting of elongated molecules. First, systems of hard spherocylinders in a spherical geometry are investigated by means of canonical Monte Carlo simulations. In contrast to previous simulation work on this problem, a continuum model is used. In particular, the effects of ordering near hard curved walls are studied for the low-density regime. With increasing density, first a uniaxial surface film forms and then a biaxial surface film, which eventually fills the entire cavity. We study how the surface order, the adsorption and the shape of the director field depend on the curvature of the wall. We find that orientational ordering at a curved wall in a cavity is stronger than at a flat wall, while adsorption is weaker. For densities above the isotropic-nematic transition, we always find bipolar configurations. As a next step, an extension of the Asakura-Oosawa-Vrij model for colloid-polymer mixtures to anisotropic colloids is considered. By means of computer simulations we study how droplets of hard, rod-like particles optimize their shape and structure under the influence of the osmotic compression caused by the presence of spherical particles that act as depletion agents. At sufficiently high osmotic pressures the rods that make up the drops spontaneously align to turn them into uniaxial nematic liquid crystalline droplets. The nematic droplets or ``tactoids'' that so form are not spherical but elongated, resulting from the competition between the anisotropic surface tension and the elastic deformation of the director field. In agreement with recent theoretical predictions we find that sufficiently small tactoids have a uniform director field, whilst large ones are characterized by a bipolar director field. From the shape and director-field transformation of the droplets we estimate the surface anchoring strength.

Quantenkorrekturen an magnetischen Domänenwänden auf Gitterstrukturen

Wir betrachten die eindimensionale Heisenberg-Spinkette aus einem neuen und aktuelleren Blickwinkel. Experimentelle Techniken der Herstellung und selbstverständlich auch experimentelle Meßmethoden erlauben nicht nur die Herstellung von Nanopartikeln und Nanodrähten, sondern gestatten es auch, Domänenwände in diesen Strukturen auszumessen. Die meisten heute verwendeten Theorien und Simulationsmethoden haben ihre Grundlage im mikromagnetischen Kontinuumsmodell, daß schon über Jahrzehnte hinweg erforscht und erprobt ist. Wir stellen uns jedoch die Frage, ob die innere diskrete Struktur der Substrate und die quantenmechanischen Effekte bei der Genauigkeit heutiger Messungen in Betracht gezogen werden müssen. Dazu wählen wir einen anderen Ansatz. Wir werden zunächst den wohlbekannten klassischen Fall erweitern, indem wir die diskrete Struktur der Materie in unseren Berechnungen berücksichtigen. Man findet in diesem Formalismus einen strukturellen Phasenübergang zwischen einer Ising-artigen und einer ausgedehnten Wand. Das führt zu bestimmten Korrekturen im Vergleich zum Kontinuumsfall. Der Hauptteil dieser Arbeit wird sich dann mit dem quantenmechanischen Fall beschäftigen. Wir rotieren das System zunächst mit einer Reihe lokaler Transformationen derart, daß alle Spins in die z-Richtung ausgerichtet sind. Im Rahmen einer 1/S-Entwicklung läßt sich der erhaltene neue Hamilton-Operator diagonalisieren. Setzt man hier die klassische Lösung ein, so erhält man Anregungsmoden in diesem Grenzfall. Unsere Resultate erweitern und bestätigen frühere Berechnungen. Mit Hilfe der Numerik wird schließlich der Erwartungswert der Energie minimiert und somit die Form der Domänenwand im quantenmechanischen Fall berechnet. Hieraus ergeben sich auch bestimmte Korrekturen zum kritischen Verhalten des Systems. Diese Ergebnisse sind vollkommen neu.

Exact quantum Monte Carlo methods for correlated electron systems

Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.

Vortex dynamics in superconductors in the presence of anisotropic pinning

Vortex dynamics in two different classes of superconductors with anisotropic unidirected pinning sites was experimentally investigated by magnetoresistivity measurements: YBCO−films with unidirected twins and Nb-films deposited on faceted $mathrm Al_2O_3$ substrate surfaces. For the interpretation of the experimental results a theoretical model based on the Fokker-Planck equation was used. It was proved by X-ray measurements that YBCO films prepared on (001) $mathrm NdGaO_3$ substrates exhibit only one twin orientation in contrast to YBCO films grown on (100) $mathrm SrTiO_$3 substrates. The magnetoresistivity measurements of the YBCO films with unidirected twin boundaries revealed the existence of two new magnetoresistivity components, which is a characteristic feature of a guided vortex motion: an odd longitudinal component with respect to the magnetic field sign reversal and an even transversal component. However, due to the small coherence length in YBCO and the higher density of point-like defects comparing to high-quality YBCO single crystals, the strength of the isotropic point pinning was comparable with the strength of the pinning produced by twins. This smeared out all effects caused by the pinning anisotropy. The behaviour of the odd longitudinal component was found to be independent of the transport current direction with respect to the twin planes. The magnetoresistivity measurements of faceted Nb films demonstrated the appearance of an odd longitudinal and even transversal component of the magnetoresistivity. The temperature and magnetic field dependences of all relevant magnetoresistivity components were measured. The angles between the average vortex velocity vector and the transport current direction calculated from the experimental data for the different transport current orientations with respect to the facet ridges showed that the vortices moved indeed along the facet ridges. An anomalous Hall effect, i.e. a sign change of the odd transversal magnetoresistivity, has been found in the temperature and magnetic field dependences of the Hall resisitivity of the samples. The theory developed by V.~A.~Shklovskij was used for the explanation of the experimental data. It shows very good agreement with the experiment. The temperature dependence of the even longitudinal magnetoresistivity component of the samples could be very well fitted within the theoretical approach, using for the isotropic and anisotropic pinning potential simple potential with a symmetric triangular potential wells whose depths were estimated from the experimental data.