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.

Helix-coil-transition in solid polypeptides

Abstract Poly(L-glutamic acid) (PLGA) was synthesized by living anionic ring-opening polymerization of the NCA monomer, which was obtained by reacting diphosgene with an amino acid derivative. The chemical structures and thermal properties were characterized by 1H-NMR, 13C-NMR, TGA and DSC. XRD powder patterns found to be amorphous for all polymers obtained. The molecular weights could be determined under severe limitations due to low solubility and high aggregation tendency. The secondary structure of the PLGA films was analyzed in the solid state by IR spectroscopy; the order was determined mainly by XRD. Uniform bulk films (1-5 µm) were produced by drop-casting of PLGA solutions in TFA on silica. The XRD film analysis indicated the absence of a long range order or an orientation even if a helical microstructure was confirmed by IR spectroscopy. The coil solvent TFA delivered constantly a helical or a β-sheet structure in the solid state depending on the water content of the solvent which was observed for the first time to exhibit a high influence on the crystallization process for PLGA. Temperature dependent in-situ IR measurements were examined to analyze if a helix-coil transition occurs, but there could be no solvent system determined, which resulted in a disordered coil structure in the solid state. General parameters like solvent systems, evaporation conditions, concentration, substrates etc. were analyzed. New crystallizations were obtained on silica prepared by drop-casting of solutions of PLGA in DMF, DMA, TMU, NMP, and pyridine/water mixtures, respectively. PSCBC in DMF, CDCl3/TFA-d, and PSBC in CDCl3/TFA-d exhibited the same crystalline diffraction patterns like PLGA. The long range order in the X-ray diffraction pattern is proven by extremely sharp crystalline signals, which are not changing the shape or the position of the peak by increasing the temperature up to 160°C. The substrate seems to play a decisive role because the crystalline structures were not obtainable on glass. The crystal structure consists probably of two different layered structures based on the intensity ratios of the two series of crystalline signals in the X-ray diffraction patterns. The source of the layered structure remains unclear and needs further studies to investigate the spatial arrangement of the chains in more detail. The secondary structure was still not changing upon heating even if a highly crystalline diffraction pattern occurs. Concluding that even the newly investigated crystallization did not show a helix-coil transition in the solid state by annealing, the phenomenon known in solution has to be claimed as unachievable in the solid state based on the results of this work. A remaining open question represents the observation that the same crystalline pattern can be reproducibly prepared with exhibiting two different ordered secondary structures (helix and β-sheet). After the investigation that the evaporation time cannot be decisive for the crystal growth, the choice of a strong hydrogen bonding interrupting solvent is most probably the key to support and induce the crystallization process.

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.