Particle Monte Carlo simulation of string-like colloidal assembly in 2 dimensions

We simulate structural phase behavior of polymer-grafted colloidal particles by molecular Monte Carlo technique. Interparticle potential, which has a finite repulsive square-step outside a rigid core of the colloid, was previously confirmed via numerical self-consistent field calculation. This model potential is purely repulsive. We simulate these model colloids in the canonical ensemble in 2 dimensions and find that these particles containing no interparticle attraction self-assemble and align in a string-like assembly, at low temperature and high density. This string-like colloidal assembly is related to percolation phenomena. Analyzing the cluster size distribution and the average string length, we build phase diagrams and discover that the average string length diverges around the region where the melting transition line and the percolation transition line cross. This result is similar to Ising spin systems, in which the percolation transition line and the order-disorder line meet at a critical point.Comment: 11 pages, 14 figure

Messung der Freien Energie von spontan strukturbildenden Systemen mittels Computersimulation

Die Fusion von Membranen ist Bestandteil einer Vielzahl von biologischen Prozessen, wie z.B. der Endo- und Exozytose, Virusinfektion oder synaptischen Ausschüttung. Der detaillierte Mechanismus ist jedoch nicht vollständig verstanden, weil die mesoskopischen Zeit- und Längenskalen von Mikrosekunden und einigen Nanometern eine direkte experimentelle Beobachtung erschweren. Die Untersuchung dieser kollektiven Prozesse, an denen viele Lipidmoleküle teilnehmen, ist auch für atomistische Simulationen schwierig. Während die Details des Fusionspfades immer noch strittig sind, beginnen die meisten Fusionsszenarien mit der Bildung einer sanduhrförmigen Verbindung ( Stalk ) zwischen den beiden gegenüberliegenden Membranen, welche fusionieren. Das Verständnis der Eigenschaften dieser ersten Übergangsstruktur ist ein Schlüssel zur Kontrolle der Fusion von Membranen. Wir verwenden ein vergröbertes Modell für Doppelschichtmembranen. Die Lipidmoleküle werden durch ein einfaches Kugel-Feder-Modell beschrieben und die Freiheitsgrade des Lösungsmittels sind herausintegriert. Die effektiven, nichtgebundenenWechselwirkungen zwischen den Segmenten der Lipidmoleküle haben die Form einer Virialentwicklung. Im Rahmen der Molekularfeldnäherung sind die Entwicklungskoeffizienten mit der Dichte und Kompressibilität des hydrophoben Membraninnerens und der Abstossung zwischen hydrophilien und hydrophoben Segmenten verknüpft. Um ein solches Dichtefunktional der ungebundenen Wechselwirkungen in Teilchensimulationen zu verwenden, werden die lokalen Dichten aus den expliziten Teilchenkoordinaten auf ein Gitter abgebildet. Dieses lösungsmittelfreie vergröberte Modell für Membranen mit weichen Wechselwirkungen erlaubt die effiziente Simulation von Membraneigenschaften. Dieses vergröberte Modell wurde verwendet, um die freie Energie eines Stalks zwischen zwei gegenüberliegenden Membranen als Funktion der molekularen Asymmetrie der Lipide und der Membranspannung zu bestimmen. Zu diesem Zweck haben wir ein allgemeines Verfahren zur Berechnung von freien Energien in selbstordnenden Systemen entwickelt. Die Methode beruht auf einem thermodynamischen Pfad, welcher das System von zwei unverbundenen Membranen mit der Stalk Struktur reversibel verbindet. Dieser Pfad verläuft in einem erweiterten Zustandsraum mit einem inhomogenen externen Feld. Dieses ist entworfen, um das System ohne ungebundene Wechselwirkungen zur Bildung der zwei gegenüberliegenden Membranen bzw. der Stalk Struktur zu veranlassen. Mittels expandedensemble - Simulationen wird gezeigt, dass der Pfad tatsächlich reversibel ist und die Helmholtz sche freie Energie wird mit großer Genauigkeit bestimmt. Durch zusätzliche großkanonische Simulationen haben wir die freie Energie des Stalks als Funktion der Membranspannung gemessen. Um die Abhängigkeit der freien Energie von der molekularen Architektur zu berechnen, verwenden wir ein semi-großkanonisches Ensemble, in dem Monte- Carlo-Schritte Lipide von einer Sorte in eine andere Lipidsorte mit verschiedener Architektur mutieren und umgekehrt. In diesem Ensemble wird die Zusammensetzung der Membran durch die Differenz der chemischen Potentiale der Spezies bestimmt und die Änderung der freien Energie bei dem Austausch einer Lipidsorte durch die andere kann leicht ermittelt werden. Mit diesen Rechentechniken haben wir die Stabilität von Stalks systematisch untersucht. Die Simulationen zeigen, dass die freie Energie eines Stalks in der Größenordnung von etwa 10 kBT liegt, wobei kBT die thermische Energieeinheit ist. Stalks bestehen aus 10-100 Lipidmolekülen und die freie Energie eines Stalks nimmt mit der Membranspannung zu. Die Stabilität hängt sehr stark von der Molekülarchitektur ab, Amphiphile mit grossen Kopfgruppen führen zu höchst metastabilen Stalks , während sehr asymmetrische Lipide die freie Energie von Stalks sogar zu negativenWerten verschieben, welche thermodynamisch stabilen Strukturen entsprechen

Conducting transition analysis of thin films composed of long flexible macromolecules: Percolation study

By simulating percolation and critical phenomena of labelled species inside films composed of single-component linear homogeneous macromolecules using the molecular Monte Carlo method in 3 dimensions, we study the dependence of these conducting transition and critical phenomena upon both thermal movements, i.e. spontaneous mobility, and extra-molecular topological constraints of the molecules. Systems containing topological constraints and/or composed of immobile particles, e.g. lattice models and chemical gelation, were studied in conventional works on percolation. Coordinates of the randomly distributed particles in the conventional lattice models are limited to discrete lattice points. Moreover, each particle is spatially fixed at the distributed position, which results in a temporally unchanged network structure. Although each polymer in the chemical gels can spontaneously move in the continuous space, the network structure is fixed when cross-linking reaction ends. By contrast to these conventional systems, all the molecules in the present system freely move and spontaneously diffuse in the continuous space. The network structure of the present molecules continues changing dynamically. The percolation and critical phenomena of such dynamic network structures are examined here. We reveal that these phenomena also occur in the present system, and that both the universality class and percolation threshold are independent of the extra-molecular topological constraints