Use of computer algebra in the calculation of Feynman diagrams

The increasing precision of current and future experiments in high-energy physics requires a likewise increase in the accuracy of the calculation of theoretical predictions, in order to find evidence for possible deviations of the generally accepted Standard Model of elementary particles and interactions. Calculating the experimentally measurable cross sections of scattering and decay processes to a higher accuracy directly translates into including higher order radiative corrections in the calculation. The large number of particles and interactions in the full Standard Model results in an exponentially growing number of Feynman diagrams contributing to any given process in higher orders. Additionally, the appearance of multiple independent mass scales makes even the calculation of single diagrams non-trivial. For over two decades now, the only way to cope with these issues has been to rely on the assistance of computers. The aim of the xloops project is to provide the necessary tools to automate the calculation procedures as far as possible, including the generation of the contributing diagrams and the evaluation of the resulting Feynman integrals. The latter is based on the techniques developed in Mainz for solving one- and two-loop diagrams in a general and systematic way using parallel/orthogonal space methods. These techniques involve a considerable amount of symbolic computations. During the development of xloops it was found that conventional computer algebra systems were not a suitable implementation environment. For this reason, a new system called GiNaC has been created, which allows the development of large-scale symbolic applications in an object-oriented fashion within the C++ programming language. This system, which is now also in use for other projects besides xloops, is the main focus of this thesis. The implementation of GiNaC as a C++ library sets it apart from other algebraic systems. Our results prove that a highly efficient symbolic manipulator can be designed in an object-oriented way, and that having a very fine granularity of objects is also feasible. The xloops-related parts of this work consist of a new implementation, based on GiNaC, of functions for calculating one-loop Feynman integrals that already existed in the original xloops program, as well as the addition of supplementary modules belonging to the interface between the library of integral functions and the diagram generator.

Polymer controlled mineralization of zinc phosphate hydrates and applications in corrosion protection, catalysis and biomedicine

ZUSAMMENFASSUNG Die Tauglichkeit von Hybridmaterialien auf der Basis von Zinkphosphathydrat-Zementen zum Einsatz als korrosionshemmende anorganische Pigmente oder zur prothetischen und konservierenden Knochen- und Zahntherapie wird weltweit empirisch seit den neunziger Jahren intensiv erforscht. In der vorliegenden Arbeit wurden zuerst Referenzproben, d.h. alpha-und beta-Hopeite (Abk. a-,b-ZPT) dank eines hydrothermalen Kristallisationsverfahrens in wässerigem Milieu bei 20°C und 90°C hergestellt. Die Kristallstruktur beider Polymorphe des Zinkphosphattetrahydrats Zn3(PO4)2  4 H2O wurde komplett bestimmt. Einkristall-strukturanalyse zeigt, daß der Hauptunterschied zwischen der alpha-und beta-Form des Zinkphosphattetrahydrats in zwei verschiedenen Anordnungen der Wasserstoffbrücken liegt. Die entsprechenden drei- und zweidimensionalen Anordnungen der Wasserstoffbrücken der a-und b-ZPT induzieren jeweils unterschiedliches thermisches Verhalten beim Aufwärmen. Während die alpha-Form ihr Kristallwasser in zwei definierten Stufen verliert, erzeugt die beta-Form instabile Dehydratationsprodukt. Dieses entspricht zwei unabhängigen, aber nebeneinander ablaufenden Dehydratationsmechanismen: (i) bei niedrigen Heizraten einen zweidimensionalen Johnson-Mehl-Avrami (JMA) Mechanismus auf der (011) Ebene, der einerseits bevorzugt an Kristallkanten stattfindet und anderseits von existierenden Kristalldefekten auf Oberflächen gesteuert wird; (ii) bei hohen Heizraten einem zweidimensionalen Diffusionsmechanismus (D2), der zuerst auf der (101) Ebene und dann auf der (110) Ebene erfolgt. Durch die Betrachtung der ZPT Dehydratation als irreversibele heterogene Festkörperstufenreaktion wurde dank eines „ähnlichen Endprodukt“-Protokolls das Dehydratationsphasendiagramm aufgestellt. Es beschreibt die möglichen Zusammenhänge zwischen den verschiedenen Hydratationszuständen und weist auf die Existenz eines Übergangszustandes um 170°C (d.h. Reaktion b-ZPT  a-ZPT) hin. Daneben wurde auch ein gezieltes chemisches Ätzverfahren mit verdünnten H3PO4- und NH3 Lösungen angewendet, um die ersten Stufe des Herauslösens von Zinkphosphat genau zu untersuchen. Allerdings zeigen alpha- und beta-Hopeite charakteristische hexagonale und kubische Ätzgruben, die sich unter kristallographischer Kontrolle verbreitern. Eine zuverlässige Beschreibung der Oberfächenchemie und Topologie konnte nur durch AFM und FFM Experimente erfolgen. Gleichzeitig konnte in dieser Weise die Oberflächendefektdichte und-verteilung und die Volumenauflösungsrate von a-ZPT und b-ZPT bestimmt werden. Auf einem zweiten Weg wurde eine innovative Strategie zur Herstellung von basischen Zinkphosphatpigmenten erster und zweiter Generation (d.h. NaZnPO4  1H2O und Na2ZnPO4(OH)  2H2O) mit dem Einsatz von einerseits oberflächenmodifizierten Polystyrolatices (z.B. produziert durch ein Miniemulsionspolymerisationsverfahren) und anderseits von Dendrimeren auf der Basis von Polyamidoamid (PAMAM) beschritten. Die erhaltene Zeolithstruktur (ZPO) hat in Abhängigkeit von steigendem Natrium und Wassergehalt unterschiedliche kontrollierte Morphologie: hexagonal, würfelförmig, herzförmig, sechsarmige Sterne, lanzettenförmige Dendrite, usw. Zur quantitativen Evaluierung des Polymereinbaus in der Kristallstruktur wurden carboxylierte fluoreszenzmarkierte Latices eingesetzt. Es zeigt sich, daß Polymeradditive nicht nur das Wachstum bis zu 8 µm.min-1 reduzierten. Trotzdem scheint es auch als starker Nukleationsbeschleuniger zu wirken. Dank der Koordinationschemie (d.h. Bildung eines sechszentrigen Komplexes L-COO-Zn-PO4*H2O mit Ligandenaustausch) konnten zwei einfache Mechanismen zur Wirkung von Latexpartikeln bei der ZPO Kristallisation aufgezeigt werden: (i) ein Intrakorona- und (ii) ein Extrakorona-Keimbildungsmechanismus. Weiterhin wurde die Effizienz eines Kurzzeit- und Langzeitkorrosionschutzes durch maßgeschneiderte ZPO/ZPT Pigmente und kontrollierte Freisetzung von Phosphationen in zwei Näherungen des Auslösungsgleichgewichts abgeschätzt: (i) durch eine Auswaschungs-methode (thermodynamischer Prozess) und (ii) durch eine pH-Impulsmethode (kinetischer Prozess. Besonders deutlich wird der Ausflösungs-Fällungsmechanismus (d.h. der Metamorphismus). Die wesentliche Rolle den Natriumionen bei der Korrosionshemmung wird durch ein passendes zusammensetzungsabhängiges Auflösungsmodell (ZAAM) beschrieben, das mit dem Befund des Salzsprühteste und der Feuchtigkeitskammertests konsistent ist. Schließlich zeigt diese Arbeit das herausragende Potential funktionalisierter Latices (Polymer) bei der kontrollierten Mineralisation zur Herstellung maßgeschneiderter Zinkphosphat Materialien. Solche Hybridmaterialien werden dringend in der Entwicklung umweltfreundlicher Korrosionsschutzpigmente sowie in der Dentalmedizin benötigt.

Regularisierung und Renormierung in der Quantenfeldtheorie : Resultate aus dem Vergleich konsistenter und praktikabler Methoden

Die Arbeit beginnt mit dem Vergleich spezieller Regularisierungsmethoden in der Quantenfeldtheorie mit dem Verfahren zur störungstheoretischen Konstruktion der S-Matrix nach Epstein und Glaser. Da das Epstein-Glaser-Verfahren selbst als Regularisierungsverfahren verwandt werden kann und darüberhinaus ausschließlich auf physikalisch motivierten Postulaten basiert, liefert dieser Vergleich ein Kriterium für die Zulässigkeit anderer Regularisierungsmethoden. Zusätzlich zur Herausstellung dieser Zulässigkeit resultiert aus dieser Gegenüberstellung als weiteres wesentliches Resultat ein neues, in der Anwendung praktikables sowie konsistentes Regularisierungsverfahren, das modifizierte BPHZ-Verfahren. Dieses wird anhand von Ein-Schleifen-Diagrammen aus der QED (Elektronselbstenergie, Vakuumpolarisation und Vertexkorrektur) demonstriert. Im Gegensatz zur vielverwandten Dimensionalen Regularisierung ist dieses Verfahren uneingeschränkt auch für chirale Theorien anwendbar. Als Beispiel hierfür dient die Berechnung der im Rahmen einer axialen Erweiterung der QED-Lagrangedichte auftretenden U(1)-Anomalie. Auf der Stufe von Mehr-Schleifen-Diagrammen zeigt der Vergleich der Epstein-Glaser-Konstruktion mit dem bekannten BPHZ-Verfahren an mehreren Beispielen aus der Phi^4-Theorie, darunter das sog. Sunrise-Diagramm, daß zu deren Berechnung die nach der Waldformel des BPHZ-Verfahrens zur Regularisierung beitragenden Unterdiagramme auf eine kleinere Klasse eingeschränkt werden können. Dieses Resultat ist gleichfalls für die Praxis der Regularisierung bedeutsam, da es bereits auf der Stufe der zu berücksichtigenden Unterdiagramme zu einer Vereinfachung führt.

Systematische Untersuchung von Instabilitäten an lasergekühlten 40 Ca + - Ionen in einer linearen Paulfalle

In linearen Paulfallen gespeicherte und lasergekühlte Ionen stellen in weiten Bereichen der Physik ideale Objekte hinsichtlich störungsfreier und präziser Messungen atomarer Übergangsfrequenzen und der gezielten Manipulation von Quantenzuständen dar. Eine Einschränkung dieser optimalen Bedingungen erfolgt durch Heizmechanismen, die aus Abweichungen des Speicherpotentials von der idealen Quadrupolform resultieren. Höhere Potentialordnungen führen zu einer Kopplung der radialen Bewegungsmoden und bei bestimmten Speicherparametern zu nichtlinearen Resonanzen. Hierbei werden die Ionenbahnen durch eine Energieaufnahme aus dem Speicherfeld destabilisiert. Dieses kann zu Linienverbreiterungen, einer Limitierung der Kohärenzzeiten und unter Umständen zu einem Ionenverlust führen. Die systematische Untersuchung dieser Instabilitäten in einer linearen Paulfalle erfolgt durch Spektroskopie an einer kleinen Anzahl lasergekühlter ^40Ca^+ - Ionen. Der experimentell zugängliche Speicherbereich wird mit hoher Auflösung abgetastet. Durch eine eingehende Quantifizierung der Falleneigenschaften werden die nichtlinearen Resonanzen eindeutig den erzeugenden Potentialtermen zugeordnet. Die Resonanzlinien zeigen eine charakteristische Aufspaltung, deren Größe vom angelegten Axialpotential bestimmt wird. Diese zusätzliche Kopplung der Radialbewegung an die Axialbewegung führt zu einer modifizierten Resonanzbedingung. Nichtlineare Resonanzen treten massenspezifisch auf. Da eine präzise Kontrolle der Axialpotentiale sehr einfach ist, könnten die beobachteten radial-axial koppelnden Resonanzen eine Anwendung in der Massenspektrometrie finden.

Statistical mechanics of protein complexed and condensed DNA

In this thesis I treat various biophysical questions arising in the context of complexed / ”protein-packed” DNA and DNA in confined geometries (like in viruses or toroidal DNA condensates). Using diverse theoretical methods I consider the statistical mechanics as well as the dynamics of DNA under these conditions. In the first part of the thesis (chapter 2) I derive for the first time the single molecule ”equation of state”, i.e. the force-extension relation of a looped DNA (Eq. 2.94) by using the path integral formalism. Generalizing these results I show that the presence of elastic substructures like loops or deflections caused by anchoring boundary conditions (e.g. at the AFM tip or the mica substrate) gives rise to a significant renormalization of the apparent persistence length as extracted from single molecule experiments (Eqs. 2.39 and 2.98). As I show the experimentally observed apparent persistence length reduction by a factor of 10 or more is naturally explained by this theory. In chapter 3 I theoretically consider the thermal motion of nucleosomes along a DNA template. After an extensive analysis of available experimental data and theoretical modelling of two possible mechanisms I conclude that the ”corkscrew-motion” mechanism most consistently explains this biologically important process. In chapter 4 I demonstrate that DNA-spools (architectures in which DNA circumferentially winds on a cylindrical surface, or onto itself) show a remarkable ”kinetic inertness” that protects them from tension-induced disruption on experimentally and biologically relevant timescales (cf. Fig. 4.1 and Eq. 4.18). I show that the underlying model establishes a connection between the seemingly unrelated and previously unexplained force peaks in single molecule nucleosome and DNA-toroid stretching experiments. Finally in chapter 5 I show that toroidally confined DNA (found in viruses, DNAcondensates or sperm chromatin) undergoes a transition to a twisted, highly entangled state provided that the aspect ratio of the underlying torus crosses a certain critical value (cf. Eq. 5.6 and the phase diagram in Fig. 5.4). The presented mechanism could rationalize several experimental mysteries, ranging from entangled and supercoiled toroids released from virus capsids to the unexpectedly short cholesteric pitch in the (toroidaly wound) sperm chromatin. I propose that the ”topological encapsulation” resulting from our model may have some practical implications for the gene-therapeutic DNA delivery process.

Silicon-containing block copolymers with well-defined, branched architectures

In the past decade, block copolymers (BCPs) have attracted increasing scientific and technological interest because of their inherent capability to spontaneously self-assemble into ordered arrays of nanostructures. The importance of nanostructures in a number of applications has fostered the need for well-defined, complex macromolecular architectures. In this thesis, the influence of macromolecular architecture on the bulk morphologies of novel linear-hyperbranched and linear brush-like diblock copolymer structure is investigated. An innovative, generally applicable strategy for the preparation of these defined diblock copolymers, consisting of linear polystyrene and branched polycarbosilane blocks, is demonstrated. Furthermore, complete characterization and solid-state morphological studies are provided. Finally, the concept is extended to linear-hyperbrached and linear brush-like polyalkoxysilanes. A shift of the classical phase boundaries to higher PS weight fractions as well as the appearance of new morphologies confirms the dramatic effect that polymer topology has on the morphology of BCPs.

Microscopic Modeling on complex networks

The field of complex systems is a growing body of knowledge, It can be applied to countless different topics, from physics to computer science, biology, information theory and sociology. The main focus of this work is the use of microscopic models to study the behavior of urban mobility, which characteristics make it a paradigmatic example of complexity. In particular, simulations are used to investigate phase changes in a finite size open Manhattan-like urban road network under different traffic conditions, in search for the parameters to identify phase transitions, equilibrium and non-equilibrium conditions . It is shown how the flow-density macroscopic fundamental diagram of the simulation shows,like real traffic, hysteresis behavior in the transition from the congested phase to the free flow phase, and how the different regimes can be identified studying the statistics of road occupancy.

Higher-order calculations in manifestly Lorentz-invariant baryon chiral perturbation theory

This thesis is concerned with calculations in manifestly Lorentz-invariant baryon chiral perturbation theory beyond order D=4. We investigate two different methods. The first approach consists of the inclusion of additional particles besides pions and nucleons as explicit degrees of freedom. This results in the resummation of an infinite number of higher-order terms which contribute to higher-order low-energy constants in the standard formulation. In this thesis the nucleon axial, induced pseudoscalar, and pion-nucleon form factors are investigated. They are first calculated in the standard approach up to order D=4. Next, the inclusion of the axial-vector meson a_1(1260) is considered. We find three diagrams with an axial-vector meson which are relevant to the form factors. Due to the applied renormalization scheme, however, the contributions of the two loop diagrams vanish and only a tree diagram contributes explicitly. The appearing coupling constant is fitted to experimental data of the axial form factor. The inclusion of the axial-vector meson results in an improved description of the axial form factor for higher values of momentum transfer. The contributions to the induced pseudoscalar form factor, however, are negligible for the considered momentum transfer, and the axial-vector meson does not contribute to the pion-nucleon form factor. The second method consists in the explicit calculation of higher-order diagrams. This thesis describes the applied renormalization scheme and shows that all symmetries and the power counting are preserved. As an application we determine the nucleon mass up to order D=6 which includes the evaluation of two-loop diagrams. This is the first complete calculation in manifestly Lorentz-invariant baryon chiral perturbation theory at the two-loop level. The numerical contributions of the terms of order D=5 and D=6 are estimated, and we investigate their pion-mass dependence. Furthermore, the higher-order terms of the nucleon sigma term are determined with the help of the Feynman-Hellmann theorem.

Thin nanostructured crystalline TiO 2 films and their applications in solar cells

Research on thin nanostructured crystalline TiO2 films has attracted considerable interests because of their intriguing physical properties and potential applications in photovoltaics. Nanostructured TiO2 film plays an important role in the TiO2 based dye-sensitized solar cells because they act as a substrate for the adsorption of dye molecules and a matrix for the transportation of electrons as well. Thus they can influence the solar cell performance significantly. Consequently, the control of the morphology including the shape, size and size distribution of the TiO2 nanostructures is critical to tune and optimize the performance of the solar cells. To control the TiO2 morphology, a strategy using amphiphilic block copolymer as templating agent coupled with sol-gel chemistry has been applied. Especially, a good-poor solvent pair induced phase separation process has been developed to guide the microphase separation behavior of the block copolymers. The amphiphilic block copolymers used include polystyrene-block-poly (ethylene oxide) (PS-b-PEO), poly (methyl methacrylate)-block-poly (ethylene oxide) (PMMA-b-PEO), and poly (ethylene oxide)-block-polystyrene-block-poly (ethylene oxide) (PEO-b-PS-b-PEO). The block copolymer undergoes a good-poor-solvent pair induced phase separation in a mixed solution of 1, 4-dioxane or N, N’-dimethyl formamide (DMF), concentrated hydrochloric acid (HCl) and Titanium tetraisopropoxide (TTIP). Specifically, in the system of PS-b-PEO, a morphology phase diagram of the inorganic-copolymer composite films was mapped by adjusting the weight fractions among 1, 4-dioxane, HCl, and TTIP in solution. The amorphous TiO2 within the titania-block copolymer composite films was crystallized by calcination at temperatures above 400C, where the organic block copolymer was simultaneously burned away. This strategy is further extended to other amphiphilic block copolymers of PMMA-b-PEO and PEO-b-PS-b-PEO, where the morphology of TiO2 films can also be controlled. The local and long range structures of the titania films were investigated by the combination of imaging techniques (AFM, SEM) and x-ray scattering techniques (x-ray reflectivity and grazing incidence small-angle x-ray scattering). Based on the knowledge of the morphology control, the crystalline TiO2 nanostructured films with different morphologies were introduced into solid state dye-sensitized solar cells. It has been found that all of the morphologies help to improve the performance of the solar cells. Especially, clustered nanoparticles, worm-like structures, foam-like structures, large collapsed nanovesicles show more pronounced performance improvement than other morphologies such as nanowires, flakes, and nanogranulars.

Virtual Compton scattering in baryon chiral perturbation theory

This thesis is concerned with the calculation of virtual Compton scattering (VCS) in manifestly Lorentz-invariant baryon chiral perturbation theory to fourth order in the momentum and quark-mass expansion. In the one-photon-exchange approximation, the VCS process is experimentally accessible in photon electro-production and has been measured at the MAMI facility in Mainz, at MIT-Bates, and at Jefferson Lab. Through VCS one gains new information on the nucleon structure beyond its static properties, such as charge, magnetic moments, or form factors. The nucleon response to an incident electromagnetic field is parameterized in terms of 2 spin-independent (scalar) and 4 spin-dependent (vector) generalized polarizabilities (GP). In analogy to classical electrodynamics the two scalar GPs represent the induced electric and magnetic dipole polarizability of a medium. For the vector GPs, a classical interpretation is less straightforward. They are derived from a multipole expansion of the VCS amplitude. This thesis describes the first calculation of all GPs within the framework of manifestly Lorentz-invariant baryon chiral perturbation theory. Because of the comparatively large number of diagrams - 100 one-loop diagrams need to be calculated - several computer programs were developed dealing with different aspects of Feynman diagram calculations. One can distinguish between two areas of development, the first concerning the algebraic manipulations of large expressions, and the second dealing with numerical instabilities in the calculation of one-loop integrals. In this thesis we describe our approach using Mathematica and FORM for algebraic tasks, and C for the numerical evaluations. We use our results for real Compton scattering to fix the two unknown low-energy constants emerging at fourth order. Furthermore, we present the results for the differential cross sections and the generalized polarizabilities of VCS off the proton.

Mixtures of colloidal rods and spheres in bulk and in confinement

When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. Monte Carlo simulations were performed to study the phase diagram of such rod-polymer mixtures. The colloidal rods were modelled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while overlap of polymers and rods is forbidden. In this thesis the emphasis was on the depletion effects caused by the addition of spheres on the isotropic phase of rod-like particles. Although most of the present experimental studies consider systems close to or beyond the isotropic-nematic transition, the isotropic phase with depletion interactions turns out to be a not less interesting topic. First, the percolation problem was studied in canonical simulations of a system of hard rods and soft spheres, where the amount of depletant was kept low to prevent phase separation of the mixture. The lowering of the percolation threshold seen in experiment is confirmed to be due to the depletion interactions. The local changes in the structure of the fluid of rods, which were measured in the simulations, indicated that the depletion forces enhance local alignment and aggregation of the rods. Then, the phase diagram of isotropic-isotropic demixing of short spherocylinders was calculated using grand canonical ensemble simulations with successive umbrella sampling. Finite size scaling analysis allowed to estimate the location of the critical point. Also, estimates for the interfacial tension between the coexisting isotropic phases and analyses of its power-law behaviour on approach of the critical point are presented. The obtained phase diagram was compared to the predictions of the free volume theory. After an analysis of the bulk, the phase behaviour in confinement was studied. The critical point of gas-liquid demixing is shifted to higher concentrations of rods and smaller concentrations of spheres due to the formation of an orientationally ordered surface film. If the separation between the walls becomes very small, the critical point is shifted back to smaller concentrations of rods because the surface film breaks up. A method to calculate the contact angle of the liquid-gas interface with the wall is introduced and the wetting behaviour on the approach to the critical point is analysed.

Multiple stellar populations in globular clusters with photometry and low resolution spectroscopy

Our view of Globular Clusters has deeply changed in the last decade. Modern spectroscopic and photometric data have conclusively established that globulars are neither coeval nor monometallic, reopening the issue of the formation of such systems. Their formation is now schematized as a two-step process, during which the polluted matter from the more massive stars of a first generation gives birth, in the cluster innermost regions, to a second generation of stars with the characteristic signature of fully CNO-processed matter. To date, star-to-star variations in abundances of the light elements (C, N, O, Na) have been observed in stars of all evolutionary phases in all properly studied Galactic globular clusters. Multiple or broad evolutionary sequences have also been observed in nearly all the clusters that have been observed with good signal-to-noise in the appropriate photometric bands. The body of evidence suggests that spreads in light-element abundances can be fairly well traced by photometric indices including near ultraviolet passbands, as CNO abundance variations affect mainly wavelengths shorter than ~400 nm owing to the rise of some NH and CN molecular absorption bands. Here, we exploit this property of near ultraviolet photometry to trace internal chemical variations and combined it with low resolution spectroscopy aimed to derive carbon and nitrogen abundances in order to maximize the information on the multiple populations. This approach has been proven to be very effective in (i) detecting multiple population, (ii) characterizing their global properties (i.e., relative fraction of stars, location in the color-magnitude diagram, spatial distribution, and trends with cluster parameters) and (iii) precisely tagging their chemical properties (i.e., extension of the C-N anticorrelation, bimodalities in the N content).

Proton-conducting copolymers, blends and composites with phosphonic acid as protogenic group

ABSTRACT: In this work, proton conducting copolymers, polymer blends and composites containing phosphonic acid groups have been prepared. Proton conduction mechanisms in these materials are discussed respectively in both, the anhydrous and humidified state. Atom transfer radical copolymerization (ATRCP) of diisopropyl-p-vinylbenzyl phosphonate (DIPVBP) and 4-vinyl pyridine (4VP) is studied for the first time in this work. The kinetic parameters are obtained by using the 1H-NMR online technique. Proton conduction in poly(vinylbenzyl phosphonic acid) (PVBPA) homopolymer and its statistical copolymers with 4-vinyl pyridine (poly(VBPA-stat-4VP)s) are comprehensively studied in both, the “dry” and “wet” state. Effects of temperature, water content and polymer composition on proton conductivities are studied and proton transport mechanisms under various conditions are discussed. The proton conductivity of the polymers is in the range of 10-6-10-8 S/cm in nominally dry state at 150 oC. However, proton conductivity of the polymers increases rapidly with water content in the polymers which can reach 10-2 S/cm at the water uptake of 25% in the polymers. The highest proton conductivity obtained from the polymers can even reach 0.3 S/cm which was measured at 85oC with 80% relative humidity in the measuring atmosphere. Poly(4-vinyl pyridine) was grafted from the surface of SiO2 nanoparticles using ATRP in this work for the first time. Following this approach, silica nanoparticles with a shell of polymeric layer are used as basic particles in a polymeric acidic matrix. The proton conductivities of the composites are studied under both, humidified and dry conditions. In dry state, the conductivity of the composites is in the range of 10-10~10-4 S/cm at 150 oC. While in humid state, the composites show much higher proton conductivity. The highest proton conductivity obtained with the composites is 0.5 S/cm measured at 85oC with 80% relative humidity in the measuring atmosphere. The miscibility of poly (vinyl phosphonic acid) and PEO is studied for the first time in this work and a phase diagram is plotted based on a DSC study and optical microscopy. With this knowledge, homogeneous PVPA/PEO mixtures are prepared as proton-conducting polymer blends. The mobility of phosphonic acid groups and PEO in the blends is determined by 1H-MAS-NMR in temperature dependent measurements. The effect of composition and the role of PEO on proton conduction are discussed.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Crystal phases and glassy dynamics in monodisperse hard ellipsoids

We have performed Monte Carlo and molecular dynamics simulations of suspensions of monodisperse, hard ellipsoids of revolution. Hard-particle models play a key role in statistical mechanics. They are conceptually and computationally simple, and they offer insight into systems in which particle shape is important, including atomic, molecular, colloidal, and granular systems. In the high density phase diagram of prolate hard ellipsoids we have found a new crystal, which is more stable than the stretched FCC structure proposed previously . The new phase, SM2, has a simple monoclinic unit cell containing a basis of two ellipsoids with unequal orientations. The angle of inclination is very soft for length-to-width (aspect) ratio l/w=3, while the other angles are not. A symmetric state of the unit cell exists, related to the densest-known packings of ellipsoids; it is not always the stable one. Our results remove the stretched FCC structure for aspect ratio l/w=3 from the phase diagram of hard, uni-axial ellipsoids. We provide evidence that this holds between aspect ratios 3 and 6, and possibly beyond. Finally, ellipsoids in SM2 at l/w=1.55 exhibit end-over-end flipping, warranting studies of the cross-over to where this dynamics is not possible. Secondly, we studied the dynamics of nearly spherical ellipsoids. In equilibrium, they show a first-order transition from an isotropic phase to a rotator phase, where positions are crystalline but orientations are free. When over-compressing the isotropic phase into the rotator regime, we observed super-Arrhenius slowing down of diffusion and relaxation, and signatures of the cage effect. These features of glassy dynamics are sufficiently strong that asymptotic scaling laws of the Mode-Coupling Theory of the glass transition (MCT) could be tested, and were found to apply. We found strong coupling of positional and orientational degrees of freedom, leading to a common value for the MCT glass-transition volume fraction. Flipping modes were not slowed down significantly. We demonstrated that the results are independent of simulation method, as predicted by MCT. Further, we determined that even intra-cage motion is cooperative. We confirmed the presence of dynamical heterogeneities associated with the cage effect. The transit between cages was seen to occur on short time scales, compared to the time spent in cages; but the transit was shown not to involve displacements distinguishable in character from intra-cage motion. The presence of glassy dynamics was predicted by molecular MCT (MMCT). However, as MMCT disregards crystallization, a test by simulation was required. Glassy dynamics is unusual in monodisperse systems. Crystallization typically intervenes unless polydispersity, network-forming bonds or other asymmetries are introduced. We argue that particle anisometry acts as a sufficient source of disorder to prevent crystallization. This sheds new light on the question of which ingredients are required for glass formation.