Computer simulations of charged systems in partially periodic geometries

This work presents algorithms for the calculation of the electrostatic interaction in partially periodic systems. The framework for these algorithms is provided by the simulation package ESPResSo, of which the author was one of the main developers. The prominent features of the program are listed and the internal structure is described. In the following, algorithms for the calculation of the Coulomb sum in three dimensionally periodic systems are described. These methods are the foundations for the algorithms for partially periodic systems presented in this work. Starting from the MMM2D method for systems with one non-periodic coordinate, the ELC method for these systems is developed. This method consists of a correction term which allows to use methods for three dimensional periodicity also for the case of two periodic coordinates. The computation time of this correction term is neglible for large numbers of particles. The performance of MMM2D and ELC are demonstrated by results from the implementations contained in ESPResSo. It is also discussed, how different dielectric constants inside and outside of the simulation box can be realized. For systems with one periodic coordinate, the MMM1D method is derived from the MMM2D method. This method is applied to the problem of the attraction of like-charged rods in the presence of counterions, and results of the strong coupling theory for the equilibrium distance of the rods at infinite counterion-coupling are checked against results from computer simulations. The degree of agreement between the simulations at finite coupling and the theory can be characterized by a single parameter gamma_RB. In the special case of T=0, one finds under certain circumstances flat configurations, in which all charges are located in the rod-rod plane. The energetically optimal configuration and its stability are determined analytically, which depends on only one parameter gamma_z, similar to gamma_RB. These findings are in good agreement with results from computer simulations.

Nanoscale effects and applications of self-organized nanostructured thin films

A nanostructured thin film is a thin material layer, usually supported by a (solid) substrate, which possesses subdomains with characteristic nanoscale dimensions (10 ~ 100 nm) that are differentiated by their material properties. Such films have captured vast research interest because the dimensions and the morphology of the nanostructure introduce new possibilities to manipulating chemical and physical properties not found in bulk materials. Block copolymer (BCP) self-assembly, and anodization to form nanoporous anodic aluminium oxide (AAO), are two different methods for generating nanostructures by self-organization. Using poly(styrene-block-methyl methacrylate) (PS-b-PMMA) nanopatterned thin films, it is demonstrated that these polymer nanopatterns can be used to study the influence of nanoscale features on protein-surface interactions. Moreover, a method for the directed assembly of adsorbed protein nanoarrays, based on the nanoscale juxtaposition of the BCP surface domains, is also demonstrated. Studies on protein-nanopattern interactions may inform the design of biomaterials, biosensors, and relevant cell-surface experiments that make use of nanoscale structures. In addition, PS-b-PMMA and AAO thin films are also demonstrated for use as optical waveguides at visible wavelengths. Due to the sub-wavelength nature of the nanostructures, scattering losses are minimized, and the optical response is amenable to analysis with effective medium theory (EMT). Optical waveguide measurements and EMT analysis of the films’ optical anisotropy enabled the in situ characterization of the PS-b-PMMA nanostructure, and a variety of surface processes within the nanoporous AAO involving (bio)macromolecules at high sensitivity.

Curvature induced nanoscale rafts in lipid membranes

We present a coarse grained model for computer simulations of lipid mixtures, which we use to study generic mechanisms for the formation of nanoscale membrane structures (lipid rafts). We observe that even a two component system can separate into rafts of finite size, and we study these rafts and other membrane structures in detail. We look at the characteristics of our model that enable these phenomena and how they may relate to lipid-cholesterol or lipid-lipid mixtures. We propose an explanation for our findings using elastic theory to describe a possible mechanism of raft stabilization via curvature differences between coexisting lipid phases and we investigate whether this theory can be used to explain the results of our computer simulations.

Periodic Higgs-de Rham flows and representations of algebraic fundamental groups

Let k := bar{F}_p for p > 2, W_n(k) := W(k)/p^n and X_n be a projective smooth W_n(k)-scheme which is W_{n+1}(k)-liftable. For all n > 1, we construct explicitly a functor, which we call the inverse Cartier functor, from a subcategory of Higgs bundles over X_n to a subcategory of flat Bundles over X_n. Then we introduce the notion of periodic Higgs-de Rham flows and show that a periodic Higgs-de Rham flow is equivalent to a Fontaine-Faltings module. Together with a p-adic analogue of Riemann-Hilbert correspondence established by Faltings, we obtain a coarse p-adic Simpson correspondence.