Mechanical alterations in airway smooth muscle after interaction with indoor pollutant - A mathematical model

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully w1derstood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholinc (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle fibers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, ir1corporates the phenomenology of both MCh and FA and reproduces experirnental results observed with ir1 vitro exposure of smooth muscle to .FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells ir1 a tissue level model. The model can also be used in different biological scales.

Electronic Structure Calculations in a 2D SixGe1-x Alloy Under an Applied Electric Field.

The recent advances and promises in nanoscience and nanotechnology have been focused on hexagonal materials, mainly on carbon-based nanostructures. Recently, new candidates have been raised, where the greatest efforts are devoted to a new hexagonal and buckled material made of silicon, named Silicene. This new material presents an energy gap due to spin-orbit interaction of approximately 1.5 meV, where the measurement of quantum spin Hall effect(QSHE) can be made experimentally. Some investigations also show that the QSHE in 2D low-buckled hexagonal structures of germanium is present. Since the similarities, and at the same time the differences, between Si and Ge, over the years, have motivated a lot of investigations in these materials. In this work we performed systematic investigations on the electronic structure and band topology in both ordered and disordered SixGe1-x alloys monolayer with 2D honeycomb geometry by first-principles calculations. We show that an applied electric field can tune the gap size for both alloys. However, as a function of electric field, the disordered alloy presents a W-shaped behavior, similarly to the pure Si or Ge, whereas for the ordered alloy a V-shaped behavior is observed.

Probing the interaction of brain fatty acid binding protein (B-FABP) with model membranes

Brain fatty acid-binding protein (B-FABP) interacts with biological membranes and delivers polyunsaturated fatty acids (FAs) via a collisional mechanism. The binding of FAs in the protein and the interaction with membranes involve a motif called "portal region", formed by two small α-helices, A1 and A2, connected by a loop. We used a combination of site-directed mutagenesis and electron spin resonance to probe the changes in the protein and in the membrane model induced by their interaction. Spin labeled B-FABP mutants and lipidic spin probes incorporated into a membrane model confirmed that BFABP interacts with micelles through the portal region and led to structural changes in the protein as well in the micelles. These changes were greater in the presence of LPG when compared to the LPC models. ESR spectra of B-FABP labeled mutants showed the presence of two groups of residues that responded to the presence of micelles in opposite ways. In the presence of lysophospholipids, group I of residues, whose side chains point outwards from the contact region between the helices, had their mobility decreased in an environment of lower polarity when compared to the same residues in solution. The second group, composed by residues with side chains situated at the interface between the α-helices, experienced an increase in mobility in the presence of the model membranes. These modifications in the ESR spectra of B-FABP mutants are compatible with a less ordered structure of the portal region inner residues (group II) that is likely to facilitate the delivery of FAs to target membranes. On the other hand, residues in group I and micelle components have their mobilities decreased probably as a result of the formation of a collisional complex. Our results bring new insights for the understanding of the gating and delivery mechanisms of FABPs.

Structural and Interaction Studies of Bacterial Polysaccharides by NMR Spectroscopy

An introduction to bacterial polysaccharides and the methods for structural determination are described in the first two parts of the thesis. In a structural elucidation of bacterial polysaccharides NMR experiments are important as is component analysis. A short description of immunochemical methods such as enzyme immunoassays is included. Two NMR techniques used for interaction studies, trNOE and STD NMR, are also discussed. The third part of the thesis discusses and summarizes the results from the included papers. The structures of the exopolysaccharides produced by two lactic acid bacteria are determined by one- and two dimensional NMR experiments. One is a heteropolysaccharide produced by Streptococcus thermophilus and the other a homopolysaccharide produced by Propionibacterium freudenreichii. The structure of an acidic polysaccharide from a marine bacterium with two serine residues in the repeating unit is also investigated. The structural and immunological relationship between two O-antigenic polysaccharides from Escherichia coli strain 180/C3 and O5 is discussed and investigated. Finally, interaction studies of an octasaccharide derived from the Salmonella enteritidis O-antigen and a bacteriophage are described which were performed with NMR experiments.

Local Structure of Hydrogen-Bonded Liquids

Ordinary yet unique, water is the substance on which life is based. Water seems, at first sight, to be a very simple molecule, consisting of two hydrogen atoms attached to one oxygen. Its small size belies the complexity of its action and its numerous anomalies, central to a broad class of important phenomena, ranging from global current circulation, terrestrial water and CO2 cycles to corrosion and wetting. The explanation of this complex behavior comes from water's unique ability to form extensive three-dimensional networks of hydrogen-bonds, whose nature and structures, in spite of a great deal of efforts involving a plethora of experimental and theoretical techniques, still lacks a complete scientific understanding. This thesis is devoted to the study of the local structure of hydrogen-bonded liquids, with a particular emphasis on water, taking advantage of a combination of core-level spectroscopies and density functional theory spectra calculations. X-ray absorption, in particular, is found to be sensitive to the local hydrogen-bond environment, thus offering a very promising tool for spectroscopic identification of specific structural configurations in water, alcohols and aqueous solutions. More specifically, the characteristic spectroscopic signature of the broken hydrogen-bond at the hydrogen side is used to analyze the structure of bulk water, leading to the finding that most molecules are arranged in two hydrogen-bond configurations, in contrast to the picture provided by molecular dynamics simulations. At the liquid-vapor interface, an interplay of surface sensitive measurements and theoretical calculations enables us to distinguish a new interfacial species in equilibrium with the gas. In a similar approach the cluster form of the excess proton in highly concentrated acid solutions and the different coordination of methanol at the vacuum interface and in the bulk can also be clearly identified. Finally the ability of core-level spectroscopies, aided by sophisticated density functional theory calculations, to directly probe the valence electronic structure of a system is used to observe the nature of the interaction between water molecules and solvated ions in solution. Water around transition metal ions is found to interact with the solute via orbital mixing with the metal d-orbitals. The hydrogen-bond between water molecules is explained in terms of electrostatic interactions enhanced by charge rehybridization in which charge transfer between connecting molecules is shown to be fundamental.

On soil-structurae interaction in large non-slender partially buried structures.

[EN]This paper addresses the seismic analysis of a deeply embedded non-slender structure hosting the pumping unit of a reservoir. The dynamic response in this type of problems is usually studied under the assumption of a perfectly rigid structure using a sub-structuring procedure (three-step solution) proposed specifically for this hypothesis.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

Structure learning of graphical models for task-oriented robot grasping

In the collective imaginaries a robot is a human like machine as any androids in science fiction. However the type of robots that you will encounter most frequently are machinery that do work that is too dangerous, boring or onerous. Most of the robots in the world are of this type. They can be found in auto, medical, manufacturing and space industries. Therefore a robot is a system that contains sensors, control systems, manipulators, power supplies and software all working together to perform a task. The development and use of such a system is an active area of research and one of the main problems is the development of interaction skills with the surrounding environment, which include the ability to grasp objects. To perform this task the robot needs to sense the environment and acquire the object informations, physical attributes that may influence a grasp. Humans can solve this grasping problem easily due to their past experiences, that is why many researchers are approaching it from a machine learning perspective finding grasp of an object using information of already known objects. But humans can select the best grasp amongst a vast repertoire not only considering the physical attributes of the object to grasp but even to obtain a certain effect. This is why in our case the study in the area of robot manipulation is focused on grasping and integrating symbolic tasks with data gained through sensors. The learning model is based on Bayesian Network to encode the statistical dependencies between the data collected by the sensors and the symbolic task. This data representation has several advantages. It allows to take into account the uncertainty of the real world, allowing to deal with sensor noise, encodes notion of causality and provides an unified network for learning. Since the network is actually implemented and based on the human expert knowledge, it is very interesting to implement an automated method to learn the structure as in the future more tasks and object features can be introduced and a complex network design based only on human expert knowledge can become unreliable. Since structure learning algorithms presents some weaknesses, the goal of this thesis is to analyze real data used in the network modeled by the human expert, implement a feasible structure learning approach and compare the results with the network designed by the expert in order to possibly enhance it.

Systematic studies of the interaction between amyloid beta-protein and lipids

The amyloid peptide (Aß), a normal constituent of neuronal and non-neuronal cells, has been shown to be a major component of the extracellular plaque of Alzheimer’s disease (AD). The interaction of Aß peptides with the lipid matrix of neuronal cell membranes plays an important role in the pathogenesis of AD. In this study, we have developed peptide-tethered artificial lipid membranes by the Langmuir-Blodgett and Langmuir-Schaefer methods. Anti-Aß40-mAb labeled with a fluorophore was used to probe the Aß40 binding to the model membrane system. Systematic studies on the antibody or Aß-membrane interactions were carried out in our model systems by Surface Plasmon Field-Enhanced Fluorescence Spectroscopy (SPFS). Aß adsorption is critically determined by the lipid composition of the membranes. Aß specifically binds with membranes of sphingomyelin, and this preferential adsorption was markedly amplified by the addition of sterols (cholesterol or 25-OH-Chol). Fluorescence microscopy indicated that 25-OH-Chol could also form micro-domains with sphingomyelin as cholesterol does at the conditions used for the built-up of the model membranes. Our findings suggest that micro-domains composed of sphingomyelin and the sterols could be the binding sites of Aß and the role of sphingomyelin in AD should receive much more attention. The artificial membranes provide a novel platform for the study on AD, and SPFS is a potential tool for detecting Aß-membrane interaction. Numerous investigations indicate that the ability of Aß to form fibrils is considerably dependent upon the levels of ß-sheet structure adopted by Aß. Membrane-mediated conformational transition of Aß has been demonstrated. In this study, we focus on the interaction of Aß and the membranes composed of POPC/SM/25-OH-Chol (2:1:1). The artificial membrane system was established by the methods as described above. Immunoassy based on a pair of monoclonal antibodies (mAbs) against different epitopes was employed to detect the orientation of the Aß at the model membranes. Kinetics of antibody-Aß binding was determined by surface plasmon field-enhanced fluorescence spectroscopy (SPFS). The attempt has also been made to probe the change in the conformation of Aß using SPFS combined with immunoassay. Melatonin was employed to induce the conformational change of Aß. The orientation and the conformational change of Aß are evaluated by analysing kinetic/affinity parameters. This work provides novel insight into the investigation on the structure of Aß at the membrane surface.

Structure of the Southern Tyrrhenian subduction system: insights from seismological analysis of anisotropy and attenuation

The Southern Tyrrhenian subduction system shows a complex interaction among asthenospheric flow, subducting slab and overriding plate. To shed light on the deformations and mechanical properties of the slab and surrounding mantle, I investigated seismic anisotropy and attenuation properties through the subduction region. I used both teleseisms and slab earthquakes, analyzing shear-wave splitting on SKS and S phases, respectively. The fast polarization directions φ, and the delay time, δt, were retrieved using the method of Silver and Chan [1991. SKS and S φ reveal a complex anisotropy pattern across the subduction zone. SKS-rays sample primarily the sub-slab region showing rotation of fast directions following the curved shape of the slab and very strong anisotropy. S-rays sample mainly the slab, showing variable φ and a smaller δt. SKS and S splitting reveals a well developed toroidal flow at SW edge of the slab, while at its NE edge the pattern is not very clear. This suggests that the anisotropy is controlled by the slab rollback, responsible for about 100 km slab parallel φ in the sub-slab mantle. The slab is weakly anisotropic, suggesting the asthenosphere as main source of anisotropy. To investigate the physical properties of the slab and surrounding regions, I analyzed the seismic P and S wave attenuation. By inverting high-quality S-waves t* from slab earthquakes, 3D attenuation models down to 300 km were obtained. Attenuation results image the slab as low-attenuation body, but with heterogeneous QS and QP structure showing spot of high attenuation , between 100-200 km depth, which could be due dehydration associated to the slab metamorphism. A low QS anomaly is present in the mantle wedge beneath the Aeolian volcanic arc and could indicate mantle melting and slab dehydration.

Fiber beam-columns models with flexure-shear interaction for nonlinear analysis of reinforced concrete structures

The aim of this study was to develop a model capable to capture the different contributions which characterize the nonlinear behaviour of reinforced concrete structures. In particular, especially for non slender structures, the contribution to the nonlinear deformation due to bending may be not sufficient to determine the structural response. Two different models characterized by a fibre beam-column element are here proposed. These models can reproduce the flexure-shear interaction in the nonlinear range, with the purpose to improve the analysis in shear-critical structures. The first element discussed is based on flexibility formulation which is associated with the Modified Compression Field Theory as material constitutive law. The other model described in this thesis is based on a three-field variational formulation which is associated with a 3D generalized plastic-damage model as constitutive relationship. The first model proposed in this thesis was developed trying to combine a fibre beamcolumn element based on the flexibility formulation with the MCFT theory as constitutive relationship. The flexibility formulation, in fact, seems to be particularly effective for analysis in the nonlinear field. Just the coupling between the fibre element to model the structure and the shear panel to model the individual fibres allows to describe the nonlinear response associated to flexure and shear, and especially their interaction in the nonlinear field. The model was implemented in an original matlab® computer code, for describing the response of generic structures. The simulations carried out allowed to verify the field of working of the model. Comparisons with available experimental results related to reinforced concrete shears wall were performed in order to validate the model. These results are characterized by the peculiarity of distinguishing the different contributions due to flexure and shear separately. The presented simulations were carried out, in particular, for monotonic loading. The model was tested also through numerical comparisons with other computer programs. Finally it was applied for performing a numerical study on the influence of the nonlinear shear response for non slender reinforced concrete (RC) members. Another approach to the problem has been studied during a period of research at the University of California Berkeley. The beam formulation follows the assumptions of the Timoshenko shear beam theory for the displacement field, and uses a three-field variational formulation in the derivation of the element response. A generalized plasticity model is implemented for structural steel and a 3D plastic-damage model is used for the simulation of concrete. The transverse normal stress is used to satisfy the transverse equilibrium equations of at each control section, this criterion is also used for the condensation of degrees of freedom from the 3D constitutive material to a beam element. In this thesis is presented the beam formulation and the constitutive relationships, different analysis and comparisons are still carrying out between the two model presented.

Stress – structure correlations in grafted polymer films

In der vorliegenden Arbeit wurden experimentelle Untersuchungen zu gepfropften Polymerfilmen durchgeführt. Dabei wurden endgepfropfte poly-methyl-methacrylate (PMMA) Bürsten hergestellt durch „grafting from“ Methoden und polystyrol (PS)/ poly-vinyl-methyl-ether (PVME) Polymerfilme gepfropft auf UV sensitiven Oberflächen untersucht. Zur Strukturuntersuchung wurden die hergestellten Systeme wurden mit Rasterkraftmikroskopie (engl.: Surface Probe Microscopy, SPM), Röntgen - und Neutronenreflektivitätsmessungen, sowie mit Röntgenstreuung unter streifenden Einfall (engl.: Grazing Incidence Small Angle X-Ray Scattering, GISAXS) untersucht. rnEs wurde gezeigt, dass ein aus der Transmissionsstreuung bekanntes Model auch für auch für die GISAXS Analyse polydisperser Polymerdomänen und Kolloidsysteme verwendet werden kann. Der maximale Fehler durch die gemachten Näherungen wurde auf < 20% abgeschätzt.rnErgebnisse aus der Strukturanalyse wurden mit mechanischen Filmeigenschaften verknüpft. Dazu wurden mechanische Spannungsexperimente durchgeführt. Hierzu wurden die zu untersuchenden Filme selektiv auf einzelne Mikro-Federbalken-Sensoren (engl.: Micro Cantilever Sensor, MCS) der MCS Arrays aufgebracht. Dies wurde durch Maskierungstechniken und Mikro-Kontaktdrucken bewerkstelligt. rnPhasenübergansexperimente der gepfropften PS/PVME Filme haben gezeigt, dass die Möglichkeit einer Polymer/Polymer Phasenseparation stark von Propfpunktdichte der gebundenen Polymerketten mit der Oberfläche abhängt. PS/PVME Filmsysteme mit hohen Pfropfpunktdichten zeigten keinen Phasenübergang. Bei niedrig gepfropften Filmsystemen waren hingegen Polymer/Polymer Phasenseparationen zu beobachten. Es wurde geschlussfolgert, dass die gepfropften Polymersysteme einen hinreichenden Grad an entropischen Freiheitsgraden benötigen um eine Phasenseparation zu zeigen. Mechanische Spannungsexperimente haben dabei das Verstehen der Phasenseparationsmechanismen möglich gemacht.rnAus Quellexperimenten dichtgepfropfter PMMA Bürsten, wurden Lösungsmittel-Polymer Wechselwirkungsparameter (-Parameter) bestimmt. Dabei wurde festgestellt, dass sich die erhaltenen Parameter aufgrund von Filmbenetzung und entropischen Effekten maßgeblich von den errechneten Bulkwerten unterscheiden. Weiterhin wurden nicht reversible Kettenverschlaufungseffekt beobachtet.

Nuclear structure studies in the xenon and radon region and the discovery of a new radon isotope by Penning trap mass spectrometry

Heutzutage gewähren hochpräzise Massenmessungen mit Penning-Fallen tiefe Einblicke in die fundamentalen Eigenschaften der Kernmaterie. Zu diesem Zweck wird die freie Zyklotronfrequenz eines Ions bestimmt, das in einem starken, homogenen Magnetfeld gespeichert ist. Am ISOLTRAP-Massenspektrometer an ISOLDE / CERN können die Massen von kurzlebigen, radioaktiven Nukliden mit Halbwertszeiten bis zu einigen zehn ms mit einer Unsicherheit in der Größenordnung von 10^-8 bestimmt werden. ISOLTRAP besteht aus einem Radiofrequenz-Quadrupol zum akkumulieren der von ISOLDE gelieferten Ionen, sowie zwei Penning-Fallen zum säubern und zur Massenbestimmung der Ionen. Innerhalb dieser Arbeit wurden die Massen von neutronenreichen Xenon- und Radonisotopen (138-146Xe und 223-229Rn) gemessen. Für elf davon wurde zum ersten Mal die Masse direkt bestimmt; 229Rn wurde im Zuge dieses Experimentes sogar erstmalig beobachtet und seine Halbwertszeit konnte zu ungefähr 12 s bestimmt werden. Da die Masse eines Nuklids alle Wechselwirkungen innerhalb des Kerns widerspiegelt, ist sie einzigartig für jedes Nuklid. Eine dieser Wechselwirkungen, die Wechselwirkung zwischen Protonen und Neutronen, führt zum Beispiel zu Deformationen. Das Ziel dieser Arbeit ist eine Verbindung zwischen kollektiven Effekten, wie Deformationen und Doppeldifferenzen von Bindungsenergien, sogenannten deltaVpn-Werten zu finden. Insbesondere in den hier untersuchten Regionen zeigen deltaVpn-Werte ein sehr ungewöhnliches Verhalten, das sich nicht mit einfachen Argumenten deuten lässt. Eine Erklärung könnte das Auftreten von Oktupoldeformationen in diesen Gebieten sein. Nichtsdestotrotz ist eine quantitative Beschreibung von deltaVpn-Werten, die den Effekt von solchen Deformationen berücksichtigt mit modernen Theorien noch nicht möglich.

Effect of mixing, particle interaction, and spatial dimension on the glass transition

In this thesis, the influence of composition changes on the glass transition behavior of binary liquids in two and three spatial dimensions (2D/3D) is studied in the framework of mode-coupling theory (MCT).The well-established MCT equations are generalized to isotropic and homogeneous multicomponent liquids in arbitrary spatial dimensions. Furthermore, a new method is introduced which allows a fast and precise determination of special properties of glass transition lines. The new equations are then applied to the following model systems: binary mixtures of hard disks/spheres in 2D/3D, binary mixtures of dipolar point particles in 2D, and binary mixtures of dipolar hard disks in 2D. Some general features of the glass transition lines are also discussed. The direct comparison of the binary hard disk/sphere models in 2D/3D shows similar qualitative behavior. Particularly, for binary mixtures of hard disks in 2D the same four so-called mixing effects are identified as have been found before by Götze and Voigtmann for binary hard spheres in 3D [Phys. Rev. E 67, 021502 (2003)]. For instance, depending on the size disparity, adding a second component to a one-component liquid may lead to a stabilization of either the liquid or the glassy state. The MCT results for the 2D system are on a qualitative level in agreement with available computer simulation data. Furthermore, the glass transition diagram found for binary hard disks in 2D strongly resembles the corresponding random close packing diagram. Concerning dipolar systems, it is demonstrated that the experimental system of König et al. [Eur. Phys. J. E 18, 287 (2005)] is well described by binary point dipoles in 2D through a comparison between the experimental partial structure factors and those from computer simulations. For such mixtures of point particles it is demonstrated that MCT predicts always a plasticization effect, i.e. a stabilization of the liquid state due to mixing, in contrast to binary hard disks in 2D or binary hard spheres in 3D. It is demonstrated that the predicted plasticization effect is in qualitative agreement with experimental results. Finally, a glass transition diagram for binary mixtures of dipolar hard disks in 2D is calculated. These results demonstrate that at higher packing fractions there is a competition between the mixing effects occurring for binary hard disks in 2D and those for binary point dipoles in 2D.

A soft-tetramer model for diblock copolymer melts: structure formation and geometric characterization

The ability of block copolymers to spontaneously self-assemble into a variety of ordered nano-structures not only makes them a scientifically interesting system for the investigation of order-disorder phase transitions, but also offers a wide range of nano-technological applications. The architecture of a diblock is the most simple among the block copolymer systems, hence it is often used as a model system in both experiment and theory. We introduce a new soft-tetramer model for efficient computer simulations of diblock copolymer melts. The instantaneous non-spherical shape of polymer chains in molten state is incorporated by modeling each of the two blocks as two soft spheres. The interactions between the spheres are modeled in a way that the diblock melt tends to microphase separate with decreasing temperature. Using Monte Carlo simulations, we determine the equilibrium structures at variable values of the two relevant control parameters, the diblock composition and the incompatibility of unlike components. The simplicity of the model allows us to scan the control parameter space in a completeness that has not been reached in previous molecular simulations.The resulting phase diagram shows clear similarities with the phase diagram found in experiments. Moreover, we show that structural details of block copolymer chains can be reproduced by our simple model.We develop a novel method for the identification of the observed diblock copolymer mesophases that formalizes the usual approach of direct visual observation,using the characteristic geometry of the structures. A cluster analysis algorithm is used to determine clusters of each component of the diblock, and the number and shape of the clusters can be used to determine the mesophase.We also employ methods from integral geometry for the identification of mesophases and compare their usefulness to the cluster analysis approach.To probe the properties of our model in confinement, we perform molecular dynamics simulations of atomistic polyethylene melts confined between graphite surfaces. The results from these simulations are used as an input for an iterative coarse-graining procedure that yields a surface interaction potential for the soft-tetramer model. Using the interaction potential derived in that way, we perform an initial study on the behavior of the soft-tetramer model in confinement. Comparing with experimental studies, we find that our model can reflect basic features of confined diblock copolymer melts.