Biaxial nematic phase in the Maier-Saupe model for a mixture of discs and cylinders

We analyze the global phase diagram of a Maier-Saupe lattice model with the inclusion of shape-disordered degrees of freedom to mimic a mixture of oblate and prolate molecules (discs and cylinders). In the neighborhood of a Landau multicritical point, solutions of the statistical problem can be written as a Landau-de Gennes expansion for the free energy. If the shape-disordered degrees of freedom are quenched, we confirm the existence of a biaxial nematic structure. If orientational and disorder degrees of freedom are allowed to thermalize, this biaxial solution becomes thermodynamically unstable. Also, we use a two-temperature formalism to mimic the presence of two distinct relaxation times, and show that a slight departure from complete thermalization is enough to stabilize a biaxial nematic phase.

The r-matrix of the Alday-Arutyunov-Frolov model

We investigate the classical integrability of the Alday-Arutyunov-Frolov model, and show that the Lax connection can be reduced to a simpler 2 x 2 representation. Based on this result, we calculate the algebra between the L-operators and find that it has a highly non-ultralocal form. We then employ and make a suitable generalization of the regularization technique proposed by Mail let for a simpler class of non-ultralocal models, and find the corresponding r- and s-matrices. We also make a connection between the operator-regularization method proposed earlier for the quantum case, and the Mail let's symmetric limit regularization prescription used for non-ultralocal algebras in the classical theory.

QUASI-TOPOLOGICAL QUANTUM FIELD THEORIES AND Z(2) LATTICE GAUGE THEORIES

We consider a two-parameter family of Z(2) gauge theories on a lattice discretization T(M) of a three-manifold M and its relation to topological field theories. Familiar models such as the spin-gauge model are curves on a parameter space Gamma. We show that there is a region Gamma(0) subset of Gamma where the partition function and the expectation value h < W-R(gamma)> i of the Wilson loop can be exactly computed. Depending on the point of Gamma(0), the model behaves as topological or quasi-topological. The partition function is, up to a scaling factor, a topological number of M. The Wilson loop on the other hand, does not depend on the topology of gamma. However, for a subset of Gamma(0), < W-R(gamma)> depends on the size of gamma and follows a discrete version of an area law. At the zero temperature limit, the spin-gauge model approaches the topological and the quasi-topological regions depending on the sign of the coupling constant.

On Equilibrium Distribution of a Reversible Growth Model

We study a probabilistic model of interacting spins indexed by elements of a finite subset of the d-dimensional integer lattice, da parts per thousand yen1. Conditions of time reversibility are examined. It is shown that the model equilibrium distribution converges to a limit distribution as the indexing set expands to the whole lattice. The occupied site percolation problem is solved for the limit distribution. Two models with similar dynamics are also discussed.

A stochastic spatially structured epidemic model with diffusive processes

We developed a stochastic lattice model to describe the vector-borne disease (like yellow fever or dengue). The model is spatially structured and its dynamical rules take into account the diffusion of vectors. We consider a bipartite lattice, forming a sub-lattice of human and another occupied by mosquitoes. At each site of lattice we associate a stochastic variable that describes the occupation and the health state of a single individual (mosquito or human). The process of disease transmission in the human population follows a similar dynamic of the Susceptible-Infected-Recovered model (SIR), while the disease transmission in the mosquito population has an analogous dynamic of the Susceptible-Infected-Susceptible model (SIS) with mosquitos diffusion. The occurrence of an epidemic is directly related to the conditional probability of occurrence of infected mosquitoes (human) in the presence of susceptible human (mosquitoes) on neighborhood. The probability of diffusion of mosquitoes can facilitate the formation of pairs Susceptible-Infected enabling an increase in the size of the epidemic. Using an asynchronous dynamic update, we study the disease transmission in a population initially formed by susceptible individuals due to the introduction of a single mosquito (human) infected. We find that this model exhibits a continuous phase transition related to the existence or non-existence of an epidemic. By means of mean field approximations and Monte Carlo simulations we investigate the epidemic threshold and the phase diagram in terms of the diffusion probability and the infection probability.

Irreversibility by competition on a Glauber-Ising model by means of the entropy production.

An out of equilibrium Ising model subjected to an irreversible dynamics is analyzed by means of a stochastic dynamics, on a effort that aims to understand the observed critical behavior as consequence of the intrinsic microscopic characteristics. The study focus on the kinetic phase transitions that take place by assuming a lattice model with inversion symmetry and under the influence of two competing Glauber dynamics, intended to describe the stationary states using the entropy production, which characterize the system behavior and clarifies its reversibility conditions. Thus, it is considered a square lattice formed by two sublattices interconnected, each one of which is in contact with a heat bath at different temperature from the other. Analytical and numerical treatments are faced, using mean-field approximations and Monte Carlo simulations. For the one dimensional model exact results for the entropy production were obtained, though in this case the phase transition that takes place in the two dimensional counterpart is not observed, fact which is in accordance with the behavior shared by lattice models presenting inversion symmetry. Results found for the stationary state show a critical behavior of the same class as the equilibrium Ising model with a phase transition of the second order, which is evidenced by a divergence with an exponent µ ¼ 0:003 of the entropy production derivative.

Statistical models for experimental model membranes.

Biological membranes are constituted from lipid bilayers and proteins. Investigation of protein-membrane interaction, essential for biological function of cells, must rest upon solid knowledge of lipid bilayer behavior. Thus, extensive studies of an experimental model for membranes, lipid bilayers in water solution, have been undertaken in the last decades. These systems present structural, thermal and electrical properties which depend on temperature, ionic strength or concentration. In this talk, we shall discuss statistical models for lipid bilayers, as well as the relation between their properties and results for properties of lipid dispersions investigated by the laboratories supervised by Teresa Lamy (IF-USP) and Amando Ito (FFCL-USP).

Comparison of model potentials for molecular dynamics simulation of crystalline silica

Simulationen von SiO2 mit dem von van Beest, Kramer und vanSanten (BKS) entwickelten Paarpotenzial erzeugen vielezufriedenstellende Ergebnisse, aber auch charakteristischeSchwachstellen. In dieser Arbeit wird das BKS-Potenzial mitzwei kürzlich vorgeschlagenen Potenzialen verglichen, dieeffektiv Mehrteilchen-Wechselwirkungen beinhalten. Der ersteAnsatz erlaubt dazu fluktuierende Ladungen, der zweiteinduzierbare Polarisierungen auf den Sauerstoffatomen. Die untersuchten Schwachstellen des BKS Potenzialsbeinhalten das Verhältnis der zwei Gitterkonstanten a und cim Quarzübergang, das von BKS falsch beschrieben wird.Cristobalit und Tridymit erscheinen instabil mit BKS.Weiterhin zeigt die BKS-Zustandsdichte charakteristischeAbweichungen von der wahren Zustandsdichte. DerÜbergangsdruck für den Stishovit I-II Übergang wird deutlichüberschätzt. Das Fluktuierende-Ladungs-Modell verbesserteinige der genannten Punkte, reproduziert aber viele andereEigenschaften schlechter als BKS. DasFluktierende-Dipol-Modell dagegen behebt alle genanntenArtefakte. Zusätzlich wird der druckinduzierte Phasenübergang imalpha-Quarz untersucht. Alle Potentiale finden die selbeStruktur für Quarz II. Bei anschliessender Dekompressionerzeugt BKS eine weitere Phase, während die beiden anderenPotentiale wieder zum alpha-Quarz zurückkehren. Weiterhinwerden zwei Methoden entwickelt, um die piezoelektrischenKonstanten bei konstantem Druck zu bestimmen. Die Ergebnissegeben Hinweise auf eine möglicherweisenicht-elektrostatische Natur der Polarisierungen imFluktuierende-Dipole-Modell. Mit dieser Interpretation scheint das Fluktuierende-DipolPotential alle verfügbaren experimentellen Daten am bestenvon allen drei untersuchten Ansätzen zu reproduzieren.

Simulations of physics and chemistry of polar stratospheric clouds with a general circulation model

A polar stratospheric cloud submodel has been developed and incorporated in a general circulation model including atmospheric chemistry (ECHAM5/MESSy). The formation and sedimentation of polar stratospheric cloud (PSC) particles can thus be simulated as well as heterogeneous chemical reactions that take place on the PSC particles. For solid PSC particle sedimentation, the need for a tailor-made algorithm has been elucidated. A sedimentation scheme based on first order approximations of vertical mixing ratio profiles has been developed. It produces relatively little numerical diffusion and can deal well with divergent or convergent sedimentation velocity fields. For the determination of solid PSC particle sizes, an efficient algorithm has been adapted. It assumes a monodisperse radii distribution and thermodynamic equilibrium between the gas phase and the solid particle phase. This scheme, though relatively simple, is shown to produce particle number densities and radii within the observed range. The combined effects of the representations of sedimentation and solid PSC particles on vertical H2O and HNO3 redistribution are investigated in a series of tests. The formation of solid PSC particles, especially of those consisting of nitric acid trihydrate, has been discussed extensively in recent years. Three particle formation schemes in accordance with the most widely used approaches have been identified and implemented. For the evaluation of PSC occurrence a new data set with unprecedented spatial and temporal coverage was available. A quantitative method for the comparison of simulation results and observations is developed and applied. It reveals that the relative PSC sighting frequency can be reproduced well with the PSC submodel whereas the detailed modelling of PSC events is beyond the scope of coarse global scale models. In addition to the development and evaluation of new PSC submodel components, parts of existing simulation programs have been improved, e.g. a method for the assimilation of meteorological analysis data in the general circulation model, the liquid PSC particle composition scheme, and the calculation of heterogeneous reaction rate coefficients. The interplay of these model components is demonstrated in a simulation of stratospheric chemistry with the coupled general circulation model. Tests against recent satellite data show that the model successfully reproduces the Antarctic ozone hole.

Photoemissions-Elektronenmikroskopie an Adsorbatsystemen

Die Photoemissions-Elektronenmikroskopie ist eine hervorragend geeignete Methode zur Untersuchung dynamischer Vorgänge auf realen polykristallinen Oberflächen im sub-μm Bereich. Bei der Anwendung auf Adsorbatsysteme lassen sich geringe Bedeckungsunterschiede, sowie Adsorbatstrukturen und -phasen unterscheiden. Die Methode erlaubt dabei ein breites Anwendungsspektrum über weite Temperaturbereiche und Systeme unterschiedlichster Bindungsenergie. Bei der Chemisorption von Sauerstoff auf polykristallinen Metallen wird unterschiedliches Aufwachsverhalten in den Helligkeitswerten im Mikroskopbild widergespiegelt. Bei Kupferproben zeigen Oberflächen mit unterschiedlicher kristalliner Richtung aufgrund der Symmetrie des fcc-Gitters ein ähnliches Verhalten. Das hexagonale Gitter des Titans zeigt dagegen große Unterschiede im Adsorptionsverhalten in Abhängigkeit der kristallinen Richtung. Diese Unterschiede konnten auf verschiedene Haftkoeffizienten und Oxidationsstufen der Metalle zurückgeführt werden. In einem Modell zur Photostromanalyse konnte beim Kupfer der Übergang von verschiedenen Überstrukturen bei wachsender Bedeckung gezeigt und die Übergänge ermittelt werden.. Auf den Titanoberflächen wurde so das Wachstum der Oxide TiO, TiO₂ und Ti₂O₃ unterschieden und die Übergänge des unterschiedlichen Wachstums ermittelt. Bei der thermischen Desorption der Schichten konnten unterschiedliche Haftkoeffizient auf einzelnen Kristalliten qualitativ gezeigt werden. Diese erstmalig eingesetzte Analysemethode weist Ähnlichkeiten zur Thermo-Desorptions-Spektroskopie (TDS) auf, zeigt jedoch ortsaufgelöst lokale Unterschiede auf polykristallinen Oberflächen. Bei thermisch gestützten Oberflächenreaktionen ließen sich die Reaktionskeime deutlich identifizieren und mit einer Grauwertanalyse konnte die Oxidation der karbidischen Lagen zu Kohlenmonoxid und die Metalloxidation unterschieden werden. Dabei konnte gezeigt werden, daß die Reaktionskeime nur an Plattengrenzen auftreten, nicht jedoch auf der Oberfläche. Durch die Aufrauhung der Plattengrenzen mit zunehmender Reaktionsdauer nimmt die Zahl der Reaktionskeime kontinuierlich zu, die laterale Ausdehnung der Einzelreaktionen bleibt aber konstant. Bei der Physisorption von Xenon auf Graphit wurde erstmals für die Photoemissionsmikroskopie die resonanten Anregung ausgenutzt. Die verschiedenen Phasen des Adsorbats können dabei deutlich unterschieden werden; bei niedrigen Temperaturen (40K) findet ein gleichmäßiges Wachstum auf der gesamten Oberfläche statt, bei höheren Temperaturen von 60-65K ist dagegen ein Inselwachstum in verschiedenen Phasen zu beobachten. Die zeitliche Entwicklung des Wachstums, die örtliche Lage der Phasen und die Phasenübergänge (gas, fest inkommensurabel, fest kommensurabel) konnten bestimmt werden. Bei der Desorption der Schichten konnten die einzelnen Phasen ebenfalls getrennt werden und das unterschiedliche Desorptionsverhalten sowie die Phasenübergänge selber verifiziert werden.

Lattice preferred orientation in omphacite: Examples from the Tauern Window, Austria and the Western Gneiss Region, Norway

It is lively debated how eclogites find their way from deep to mid-crustal levels during exhumation. Different exhumation models for high-pressure and ultrahigh-pressure rocks were suggested in previous studies, based mainly on field observations and less on microstructural studies on the exhumed rocks. The development and improvement of electron microscopy techniques allows it, to focus interest on direct investigations of microstructures and crystallographic properties in eclogites. In this case, it is of importance to study the applicability of crystallographic measurements on eclogites for exhumation processes and to unravel which processes affect eclogite textures. Previous studies suggested a strong relationship between deformation and lattice preferred orientation (LPO) in omphacite but it is still unclear if the deformation is related to the exhumation of eclogites. This study is focused on the questions which processes affect omphacite LPO and if textural investigations of omphacite are applicable for studying eclogite exhumation. Therefore, eclogites from two examples in the Alps and in the Caledonides were collected systematically and investigated with respect to omphacite LPO by using the electron backscattered diffraction (EBSD) technique. Omphacite textures of the Tauern Window (Austria) and the Western Gneiss Region (Norway) were studied to compare lattice preferred orientation with field observations and suggested exhumation models from previous studies. The interpretation of omphacite textures, regarding the deformation regime is mainly based on numerical simulations in previous studies. Omphacite LPO patterns of the Eclogite Zone are clearly independent from any kind of exhumation process. The textures were generated during omphacite growth on the prograde path of eclogite development until metamorphic peak conditions. Field observations in the Eclogite Zone show that kinematics in garnet mica schist, surrounding the eclogites, strongly indicate an extrusion wedge geometry. Stretching lineations show top-N thrusting at the base and a top-S normal faulting with a sinistral shear component at the top of the Eclogite Zone. The different shear sense on both sides of the unit does not affect the omphacite textures in any way. The omphacite lattice preferred orientation patterns of the Western Gneiss Region can not be connected with any exhumation model. The textures were probably generated during the metamorphic peak and reflect the change from subduction to exhumation. Eclogite Zone and Western Gneiss Region differ significantly in size and especially in metamorphic conditions. While the Eclogite Zone is characterized by constant P-T conditions (600-650°C, 20-25 kbar), the Western Gneiss Region contains a wide P-T range from high- to ultrahigh pressure conditions (400-800°C, 20-35 kbar). In contrast to this, the omphacite textures of both units are very similar. This means that omphacite LPO is independent from P-T conditions and therefore from burial depth. Further, in both units, omphacite LPO is independent from grain and subgrain size as well as from any shape preferred orientation (SPO) on grain and subgrain scale. Overall, omphacite lattice preferred orientation are generated on the prograde part of omphacite development. Therefore, textural investigations on omphacite LPO are not applicable to study eclogite exhumation.

Single polymers at surfaces: adsorption and detachment

This thesis is concerned with the adsorption and detachment of polymers at planar, rigid surfaces. We have carried out a systematic investigation of adsorption of polymers using analytical techniques as well as Monte Carlo simulations with a coarse grained off-lattice bead spring model. The investigation was carried out in three stages. In the first stage the adsorption of a single multiblock AB copolymer on a solid surface was investigated by means of simulations and scaling analysis. It was shown that the problem could be mapped onto an effective homopolymer problem. Our main result was the phase diagram of regular multiblock copolymers which shows an increase in the critical adsorption potential of the substrate with decreasing size of blocks. We also considered the adsorption of random copolymers which was found to be well described within the annealed disorder approximation. In the next phase, we studied the adsorption kinetics of a single polymer on a flat, structureless surface in the regime of strong physisorption. The idea of a ’stem-flower’ polymer conformation and the mechanism of ’zipping’ during the adsorption process were used to derive a Fokker-Planck equation with reflecting boundary conditions for the time dependent probability distribution function (PDF) of the number of adsorbed monomers. The numerical solution of the time-dependent PDF obtained from a discrete set of coupled differential equations were shown to be in perfect agreement with Monte Carlo simulation results. Finally we studied force induced desorption of a polymer chain adsorbed on an attractive surface. We approached the problem within the framework of two different statistical ensembles; (i) by keeping the pulling force fixed while measuring the position of the polymer chain end, and (ii) by measuring the force necessary to keep the chain end at fixed distance above the adsorbing plane. In the first case we treated the problem within the framework of the Grand Canonical Ensemble approach and derived analytic expressions for the various conformational building blocks, characterizing the structure of an adsorbed linear polymer chain, subject to pulling force of fixed strength. The main result was the phase diagram of a polymer chain under pulling. We demonstrated a novel first order phase transformation which is dichotomic i.e. phase coexistence is not possible. In the second case, we carried out our study in the “fixed height” statistical ensemble where one measures the fluctuating force, exerted by the chain on the last monomer when a chain end is kept fixed at height h over the solid plane at different adsorption strength ε. The phase diagram in the h − ε plane was calculated both analytically and by Monte Carlo simulations. We demonstrated that in the vicinity of the polymer desorption transition a number of properties like fluctuations and probability distribution of various quantities behave differently, if h rather than the force, f, is used as an independent control parameter.

Investigation of a superthermal ultracold neutron source based on a solid deuterium converter for the TRIGA Mainz reactor

Research in fundamental physics with the free neutron is one of the key tools for testing the Standard Model at low energies. Most prominent goals in this field are the search for a neutron electric dipole moment (EDM) and the measurement of the neutron lifetime. Significant improvements of the experimental performance using ultracold neutrons (UCN) require reduction of both systematic and statistical errors.rnThe development and construction of new UCN sources based on the superthermal concept is therefore an important step for the success of future fundamental physics with ultracold neutrons. rnSignificant enhancement of today available UCN densities strongly correlates with an efficient use of an UCN converter material. The UCN converter here is to be understood as a medium which reduces the velocity of cold neutrons (CN, velocity of about 600 m/s) to the velocity of UCN (velocity of about 6 m/s).rnSeveral big research centers around the world are presently planning or constructing new superthermal UCN sources, which are mainly based on the use of either solid deuterium or superfluid helium as UCN converter.rnThanks to the idea of Yu.Pokotilovsky, there exists the opportunity to build competitive UCN sources also at small research reactors of the TRIGA type. Of course these smaller facilities don't promise high UCN densities of several 1000 UCN/cm³, but they are able to provide densities around 100 UCN/cm³ for experiments.rnIn the context of this thesis, it was possible to demonstrate succesfully the feasibility of a superthermal UCN source at the tangential beamport C of the research reactor TRIGA Mainz. Based on a prototype for the future UCN source at the Forschungs-Neutronenquelle Heinz Maier-Leibnitz (FRMII) in munich, which was planned and built in collaboration with the Technical University of Munich, further investigations and improvements were done and are presented in this thesis. rnIn parallel, a second UCN source for the radial beamport D was designed and built. The comissioning of this new source is foreseen in spring 2010.rnAt beamport D with its higher thermal neutron flux, it should be possible to increase the available UCN densities of 4 UCN/cm³ by minimum one order of magnitude.

Solid-state NMR on columnar supramolecular structures

The present work deals with the characterisation of three columnar self-assembled systems, that is, benzene-1,3,5-tricarboxamides, a peripherally thioalkyl-substituted phthalocyanine, and several oligo-(p-phenylenevinylene)s. In order to probe the supramolecular organisation solid-state NMR has been used as the main technique, supported by X-ray measurements, theoretical methods, and thermal analysis. rnrnBenzene-1,3,5-tricarboxamides (BTAs) turned out to be well suited model compounds to study various fundamental supramolecular interactions, such as π-π-interactions, hydrogen bonding, as well as dynamic and steric effects of attached side chains. Six BTAs have been investigated in total, five with a CO-centred amide group bearing different side chains and one with an inverted N-centred amide group. The physical properties of these BTAs have been investigated as a function of temperature. The results indicated that in case of the CO-centred BTAs the stability of the columnar mesophase depends strongly on the nature of the side chains. Further experiments revealed a coplanar orientation of adjacent BTA molecules in the columnar assembly of CO-centred BTAs, whereas the N-centred BTA, showed a deviating not fully coplanar arrangement. These differences were ascribed to distinct hydrogen bonding schemes, involving a parallel alignment of hydrogen bonds in case of CO-centred BTAs and an antiparallel alignment in case of the N-centred counterpart.rnrn The fundamental insights of the supramolecular organisation of BTAs could be partially adapted to an octa-substituted phthalocyanine with thiododecyl moieties. Solid-state NMR in combination with chemical shift calculations determined a tilted herringbone arrangement of phthalocyanine rings in the crystalline phase as well as in the mesophase. Moreover, 1H NMR measurements in the mesophase of this compound suggested an axial rotation of molecules, which is inhibited in the crystalline phase.rnrnAs a third task, the supramolecular assembly of oligo-(p-phenylenevinylene)s of varying length and with different polar head groups have been investigated by a combined X-ray and solid-state NMR study. The results revealed a columnar structure formation of these compounds, being promoted by phase separation of alkyl side chains and aromatic rigid rods. In this system solid-state NMR yielded meaningful insight into the isotropisation process of butoxy and 2-S-methylbutoxy substituted oligo-(p-phenylenevinylene) rods.rn

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.