The Dynamics of a Small Model Glass Former as Viewed from Its Potential Energy Landscape

Despite intensive research during the last decades, thetheoreticalunderstanding of supercooled liquids and the glasstransition is stillfar from being complete. Besides analytical investigations,theso-called energy-landscape approach has turned out to beveryfruitful. In the literature, many numerical studies havedemonstratedthat, at sufficiently low temperatures, all thermodynamicquantities can be predicted with the help of the propertiesof localminima in the potential-energy-landscape (PEL). The main purpose of this thesis is to strive for anunderstanding ofdynamics in terms of the potential energy landscape. Incontrast to the study of static quantities, this requirestheknowledge of barriers separating the minima.Up to now, it has been the general viewpoint that thermallyactivatedprocesses ('hopping') determine the dynamics only belowTc(the critical temperature of mode-coupling theory), in thesense that relaxation rates follow from local energybarriers.As we show here, this viewpoint should be revisedsince the temperature dependence of dynamics is governed byhoppingprocesses already below 1.5Tc.At the example of a binary mixture of Lennard-Jonesparticles (BMLJ),we establish a quantitative link from the diffusioncoefficient,D(T), to the PEL topology. This is achieved in three steps:First, we show that it is essential to consider wholesuperstructuresof many PEL minima, called metabasins, rather than singleminima. Thisis a consequence of strong correlations within groups of PELminima.Second, we show that D(T) is inversely proportional to theaverageresidence time in these metabasins. Third, the temperaturedependenceof the residence times is related to the depths of themetabasins, asgiven by the surrounding energy barriers. We further discuss that the study of small (but not toosmall) systemsis essential, in that one deals with a less complex energylandscapethan in large systems. In a detailed analysis of differentsystemsizes, we show that the small BMLJ system consideredthroughout thethesis is free of major finite-size-related artifacts.

Picard-Fuchs equations of dimensionally regulated Feynman integrals

Zusammenfassung In der vorliegenden Arbeit besch¨aftige ich mich mit Differentialgleichungen von Feynman– Integralen. Ein Feynman–Integral h¨angt von einem Dimensionsparameter D ab und kann f¨ur ganzzahlige Dimension als projektives Integral dargestellt werden. Dies ist die sogenannte Feynman–Parameter Darstellung. In Abh¨angigkeit der Dimension kann ein solches Integral divergieren. Als Funktion in D erh¨alt man eine meromorphe Funktion auf ganz C. Ein divergentes Integral kann also durch eine Laurent–Reihe ersetzt werden und dessen Koeffizienten r¨ucken in das Zentrum des Interesses. Diese Vorgehensweise wird als dimensionale Regularisierung bezeichnet. Alle Terme einer solchen Laurent–Reihe eines Feynman–Integrals sind Perioden im Sinne von Kontsevich und Zagier. Ich beschreibe eine neue Methode zur Berechnung von Differentialgleichungen von Feynman– Integralen. ¨ Ublicherweise verwendet man hierzu die sogenannten ”integration by parts” (IBP)– Identit¨aten. Die neue Methode verwendet die Theorie der Picard–Fuchs–Differentialgleichungen. Im Falle projektiver oder quasi–projektiver Variet¨aten basiert die Berechnung einer solchen Differentialgleichung auf der sogenannten Griffiths–Dwork–Reduktion. Zun¨achst beschreibe ich die Methode f¨ur feste, ganzzahlige Dimension. Nach geeigneter Verschiebung der Dimension erh¨alt man direkt eine Periode und somit eine Picard–Fuchs–Differentialgleichung. Diese ist inhomogen, da das Integrationsgebiet einen Rand besitzt und daher nur einen relativen Zykel darstellt. Mit Hilfe von dimensionalen Rekurrenzrelationen, die auf Tarasov zur¨uckgehen, kann in einem zweiten Schritt die L¨osung in der urspr¨unglichen Dimension bestimmt werden. Ich beschreibe außerdem eine Methode, die auf der Griffiths–Dwork–Reduktion basiert, um die Differentialgleichung direkt f¨ur beliebige Dimension zu berechnen. Diese Methode ist allgemein g¨ultig und erspart Dimensionswechsel. Ein Erfolg der Methode h¨angt von der M¨oglichkeit ab, große Systeme von linearen Gleichungen zu l¨osen. Ich gebe Beispiele von Integralen von Graphen mit zwei und drei Schleifen. Tarasov gibt eine Basis von Integralen an, die Graphen mit zwei Schleifen und zwei externen Kanten bestimmen. Ich bestimme Differentialgleichungen der Integrale dieser Basis. Als wichtigstes Beispiel berechne ich die Differentialgleichung des sogenannten Sunrise–Graphen mit zwei Schleifen im allgemeinen Fall beliebiger Massen. Diese ist f¨ur spezielle Werte von D eine inhomogene Picard–Fuchs–Gleichung einer Familie elliptischer Kurven. Der Sunrise–Graph ist besonders interessant, weil eine analytische L¨osung erst mit dieser Methode gefunden werden konnte, und weil dies der einfachste Graph ist, dessen Master–Integrale nicht durch Polylogarithmen gegeben sind. Ich gebe außerdem ein Beispiel eines Graphen mit drei Schleifen. Hier taucht die Picard–Fuchs–Gleichung einer Familie von K3–Fl¨achen auf.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Computer simulations of charged systems in partially periodic geometries

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

The grain-scale distribution and behaviour of melt and fluid in crystalline analogue systems

The production, segregation and migration of melt and aqueous fluids (henceforth called liquid) plays an important role for the transport of mass and energy within the mantle and the crust of the Earth. Many properties of large-scale liquid migration processes such as the permeability of a rock matrix or the initial segregation of newly formed liquid from the host-rock depends on the grain-scale distribution and behaviour of liquid. Although the general mechanisms of liquid distribution at the grain-scale are well understood, the influence of possibly important modifying processes such as static recrystallization, deformation, and chemical disequilibrium on the liquid distribution is not well constrained. For this thesis analogue experiments were used that allowed to investigate the interplay of these different mechanisms in-situ. In high-temperature environments where melts are produced, the grain-scale distribution in “equilibrium” is fully determined by the liquid fraction and the ratio between the solid-solid and the solid-liquid surface energy. The latter is commonly expressed as the dihedral or wetting angle between two grains and the liquid phase (Chapter 2). The interplay of this “equilibrium” liquid distribution with ongoing surface energy driven recrystallization is investigated in Chapter 4 and 5 with experiments using norcamphor plus ethanol liquid. Ethanol in contact with norcamphor forms a wetting angle of about 25°, which is similar to reported angles of rock-forming minerals in contact with silicate melt. The experiments in Chapter 4 show that previously reported disequilibrium features such as trapped liquid lenses, fully-wetted grain boundaries, and large liquid pockets can be explained by the interplay of the liquid with ongoing recrystallization. Closer inspection of dihedral angles in Chapter 5 reveals that the wetting angles are themselves modified by grain coarsening. Ongoing recrystallization constantly moves liquid-filled triple junctions, thereby altering the wetting angles dynamically as a function of the triple junction velocity. A polycrystalline aggregate will therefore always display a range of equilibrium and dynamic wetting angles at raised temperature, rather than a single wetting angle as previously thought. For the deformation experiments partially molten KNO3–LiNO3 experiments were used in addition to norcamphor–ethanol experiments (Chapter 6). Three deformation regimes were observed. At a high bulk liquid fraction >10 vol.% the aggregate deformed by compaction and granular flow. At a “moderate” liquid fraction, the aggregate deformed mainly by grain boundary sliding (GBS) that was localized into conjugate shear zones. At a low liquid fraction, the grains of the aggregate formed a supporting framework that deformed internally by crystal plastic deformation or diffusion creep. Liquid segregation was most efficient during framework deformation, while GBS lead to slow liquid segregation or even liquid dispersion in the deforming areas.

Rheology and Fourier-transform rheology on water-based systems

The influence of shear fields on water-based systems was investigated within this thesis. The non-linear rheological behaviour of spherical and rod-like particles was examined with Fourier-Transform rheology under LAOS conditions. As a model system for spherical particles two different kinds of polystyrene dispersions, with a solid content higher than 0.3 each, were synthesised within this work. Due to the differences in polydispersity and Debye-length, differences were also found in the rheology. In the FT-rheology both kinds of dispersions showed a similar rise in the intensities of the magnitudes of the odd higher harmonics, which were predicted by a model. The in some cases additionally appearing second harmonics were not predicted. A novel method to analyse the time domain signal was developed, that splits the time domain signal up in four characteristic functions. Those characteristic functions correspond to rheological phenomena. In some cases the intensities of the Fourier components can interfere negatively. FD-virus particles were used as a rod-like model system, which already shows a highly non-linear behaviour at concentrations below 1. % wt. Predictions for the dependence of the higher harmonics from the strain amplitude described the non-linear behaviour well at large, but no so good at small strain amplitudes. Additionally the trends of the rheological behaviour could be described by a theory for rod-like particles. An existing rheo-optical set-up was enhanced by reducing the background birefringence by a factor of 20 and by increasing the time resolution by a factor of 24. Additionally a combination of FT-rheology and rheo-optics was achieved. The influence of a constant shear field on the crystallisation process of zinc oxide in the presence of a polymer was examined. The crystallites showed a reduction in length by a factor of 2. The directed addition of polymers in combination with a defined shear field can be an easy way for a defined change of the form of crystallites.

Influences of deep convective cloud systems on tropospheric trace gases and photochemistry over the tropical West Pacific

A numerical model for studying the influences of deep convective cloud systems on photochemistry was developed based on a non-hydrostatic meteorological model and chemistry from a global chemistry transport model. The transport of trace gases, the scavenging of soluble trace gases, and the influences of lightning produced nitrogen oxides (NOx=NO+NO2) on the local ozone-related photochemistry were investigated in a multi-day case study for an oceanic region located in the tropical western Pacific. Model runs considering influences of large scale flows, previously neglected in multi-day cloud resolving and single column model studies of tracer transport, yielded that the influence of the mesoscale subsidence (between clouds) on trace gas transport was considerably overestimated in these studies. The simulated vertical transport and scavenging of highly soluble tracers were found to depend on the initial profiles, reconciling contrasting results from two previous studies. Influences of the modeled uptake of trace gases by hydrometeors in the liquid and the ice phase were studied in some detail for a small number of atmospheric trace gases and novel aspects concerning the role of the retention coefficient (i.e. the fraction of a dissolved trace gas that is retained in the ice phase upon freezing) on the vertical transport of highly soluble gases were illuminated. Including lightning NOx production inside a 500 km 2-D model domain was found to be important for the NOx budget and caused small to moderate changes in the domain averaged ozone concentrations. A number of sensitivity studies yielded that the fraction of lightning associated NOx which was lost through photochemical reactions in the vicinity of the lightning source was considerable, but strongly depended on assumptions about the magnitude and the altitude of the lightning NOx source. In contrast to a suggestion from an earlier study, it was argued that the near zero upper tropospheric ozone mixing ratios which were observed close to the study region were most probably not caused by the formation of NO associated with lightning. Instead, it was argued in agreement with suggestions from other studies that the deep convective transport of ozone-poor air masses from the relatively unpolluted marine boundary layer, which have most likely been advected horizontally over relatively large distances (both before and after encountering deep convection) probably played a role. In particular, it was suggested that the ozone profiles observed during CEPEX (Central Equatorial Pacific Experiment) were strongly influenced by the deep convection and the larger scale flow which are associated with the intra-seasonal oscillation.

Synthetical engineering of supramolecular properties of large polycyclic aromatic hydrocarbons

An important property for devices is the charge-carrier mobility values for discotic organic materials like hexa-peri-hexabenzocoronenes. A close relation exists between the degree of their columnar self-arrangement of the molecules and their mobilities. Within this first step an induction of a higher order via hydrogen-bonding was considered, which mainly pointed towards the improvement of the intracolumnar stacking of the materials. For the analytics a broad range of methods was used including differential scanning calorimetry (DSC), wide-angle X-ray diffractometry (WAXS), solid-state NMR spectroscopy and scanning tunneling microscopy (STM). Indeed, a specific influence of the hydrogen-bonds could be identified, although in several cases by the cost of a severe reduction of solubility and processability. This effect was dampened by the addition of a long alkyl chain next to the hydrogen-bond exerting functional group, which resulted in an improved columnar arrangement by retention of processability. In contrast to the before mentioned example of inducing a higher intracolumnar order by hydrogen-bonding, the focus was also be set upon larger aromatic systems. The charge-carrier mobility is also in close relation to the size of the aromatic core and larger π-areas are expected to lead to improved mobilities. For photovoltaic applications a high extinction coefficient over a broad range of the spectrum is favorable, which can also be achieved by enlarging the aromatic core component. In addition the stronger π-interactions between the aromatic core components should yield an improved columnar stability and order. However the strengthening of the π-interactions between the aromatic core components led to a reduction of the solubility and the processability due to the stronger aggregation of the molecules. This required the introduction of efficiently solubilizing features in terms of long alkyl chains in the corona of the aromatic entity, in combination of a distortion of the aromatic core moiety by bulky tert-butyl groups. By this approach not only the processing and cleaning of the materials with standard laboratory techniques became possible, but moreover the first structure-rich UV/vis and a resolved 1H-NMR spectra for an aromatic system two times larger than hexa-peri-hexabenzocoronene were recorded. The bulk properties in an extruded fiber as well as on the surface showed a columnar self-assembly including a phase in which a homeotropic alignment on a substrate was observed, which turns the material into an interesting candidate for future applications in electronic devices.

Existence of solutions and asymtptotic limits of the Euler-Poisson and the quantum drift-diffusion systems

My work concerns two different systems of equations used in the mathematical modeling of semiconductors and plasmas: the Euler-Poisson system and the quantum drift-diffusion system. The first is given by the Euler equations for the conservation of mass and momentum, with a Poisson equation for the electrostatic potential. The second one takes into account the physical effects due to the smallness of the devices (quantum effects). It is a simple extension of the classical drift-diffusion model which consists of two continuity equations for the charge densities, with a Poisson equation for the electrostatic potential. Using an asymptotic expansion method, we study (in the steady-state case for a potential flow) the limit to zero of the three physical parameters which arise in the Euler-Poisson system: the electron mass, the relaxation time and the Debye length. For each limit, we prove the existence and uniqueness of profiles to the asymptotic expansion and some error estimates. For a vanishing electron mass or a vanishing relaxation time, this method gives us a new approach in the convergence of the Euler-Poisson system to the incompressible Euler equations. For a vanishing Debye length (also called quasineutral limit), we obtain a new approach in the existence of solutions when boundary layers can appear (i.e. when no compatibility condition is assumed). Moreover, using an iterative method, and a finite volume scheme or a penalized mixed finite volume scheme, we numerically show the smallness condition on the electron mass needed in the existence of solutions to the system, condition which has already been shown in the literature. In the quantum drift-diffusion model for the transient bipolar case in one-space dimension, we show, by using a time discretization and energy estimates, the existence of solutions (for a general doping profile). We also prove rigorously the quasineutral limit (for a vanishing doping profile). Finally, using a new time discretization and an algorithmic construction of entropies, we prove some regularity properties for the solutions of the equation obtained in the quasineutral limit (for a vanishing pressure). This new regularity permits us to prove the positivity of solutions to this equation for at least times large enough.

From pyrene to large polycyclic aromatic hydrocarbons

Studies of polycyclic aromatic hydrocarbons have shown that the overall size, periphery, and functionalization of PAHs are crucial parameters which significantly alter their electronic structure and chemical reactivity. Therefore, the major direction of this thesis is the synthesis and characterization of extended PAHs: (i) with different functional groups improving their processability, (ii) with different periphery changing their chemical reactivity, (iii) with inclusions of different metal ions, which influence their physical properties. • The cyclodehydrogenation reaction has been proposed for to synthesise polyphenylene ribbons with preplanarized (dibenzo[e,l]pyrene) moieties in the aromatic core with up to 10 nm linear size. The synthetic strategy employed is discussed in Chapter 2 and is based on stoichiometrically controlled DIELS-ALDER cycloaddition. All molecules possessed very good solubility in common organic solvents allowing their characterization by standard analytical techniques. • A new concept was developed to extend PAH’s core. Here the introduction of “zigzag” sites, discussed in Chpater 3 was shown to lower the HOMO-LUMO gap and to result in higher chemical reactivities. This allowed, in Chapters 3, 4 and 5, further functionalization of PAH and enlargement of their aromatic cores up to 224 atoms. Despite the size of these novel molecules, extraordinary solubilities in common organic solvents were obtained due to distortions from planarity of the aromatic cores by bulky tert-butyl groups, which hampered the usually very pronounced aggregation tendency of extended π-systems. All extended PAHs posses the small HOMO-LUMO gap together with good electron affinities making them potential candidates for application in organic FETs. • Another alternative synthetic route has been proposed to obtain extended the metal-PAH complexes. Using the quinoxaline methodology in Chapter 5 three new phenanthroline ligands (up to 60 skeletal atoms) have been synthesized and characterized. Four different (Ru(II), Cu(II) and Pt(II)) complexes were synthesized, allowing to construct a range of large metal complexes by varying the metal as well as the number and nature of ligands.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Effect of drug physicochemical properties on the release from liposomal systems in vitro and in vivo

Liposomes were discovered about 40 years ago by A. Bangham and since then they became very versatile tools in biology, biochemistry and medicine. Liposomes are the smallest artificial vesicles of spherical shape that can be produced from natural untoxic phospholipids and cholesterol. Liposome vesicles can be used as drug carriers and become loaded with a great variety of molecules, such as small drug molecules, proteins, nucleotides and even plasmids. Due to the variability of liposomal compositions they can be used for a large number of applications. In this thesis the β-adrenoceptor antagonists propranolol, metoprolol, atenolol and pindolol, glucose, 18F-Fluorodeoxyglucose (FDG) and Er-DTPA were used for encapsulation in liposomes, characterization and in vitro release studies. Multilamellar vesicles (MLV), large unilamellar vesicles (LUV) and smaller unilamellar vesicles (SUV) were prepared using one of the following lipids: 1,2-Dimyristoyl-sn-Glycero-3-Phosphocholine (DMPC), 1,2-Distearoyl-sn-Glycero-3-Phosphocholine (DSPC), Phospholipone 90H (Ph90H) or a mixture of DSPC and DMPC (1:1). The freeze thawing method was used for preparation of liposomes because it has three advantages (1) avoiding the use of chloroform, which is used in other methods and causes toxicity (2) it is a simple method and (3) it gives high entrapping efficiency. The percentage of entrapping efficiencies (EE) was different depending on the type and phase transition temperature (Tc) of the lipid used. The average particle size and particle size distribution of the prepared liposomes were determined using both dynamic light scattering (DLS) and laser diffraction analyzer (LDA). The average particle size of the prepared liposomes differs according to both liposomal type and lipid type. Dispersion and dialysis techniques were used for the study of the in vitro release of β-adrenoceptor antagonists. The in vitro release rate of β-adrenoceptor antagonists was increased from MLV to LUV to SUV. Regarding the lipid type, β-adrenoceptor antagonists exhibited different in vitro release pattern from one lipid to another. Two different concentrations (50 and 100mg/ml) of Ph90H were used for studying the effect of lipid concentration on the in vitro release of β-adrenoceptor antagonists. It was found that liposomes made from 50 mg/ml Ph90H exhibited higher release rates than liposomes made at 100 mg/ml Ph90H. Also glucose was encapsulated in MLV, LUV and SUV using 1,2-Dimyristoyl-sn-Glycero-3-Phosphocholine (DMPC), 1,2-Distearoyl-sn-Glycero-3-Phosphocholine (DSPC), Phospholipone 90H (Ph90H), soybean lipid (Syb) or a mixture of DSPC and DMPC (1:1). The average particle size and size distribution were determined using laser diffraction analysis. It was found that both EE and average particle size differ depending on both lipid and liposomal types. The in vitro release of glucose from different types of liposomes was performed using a dispersion method. It was found that the in vitro release of glucose from different liposomes is dependent on the lipid type. 18F-FDG was encapsulated in MLV 1,2-Dimyristoyl-sn-Glycero-3-Phosphocholine (DMPC), 1,2-Distearoyl-sn-Glycero-3-Phosphocholine (DSPC), Phospholipone 90H (Ph90H), soybean lipid (Syb) or a mixture of DSPC and DMPC (1:1). FDG-containing LUV and SUV were prepared using Ph90H lipid. The in vitro release of FDG from the different types of lipids was accomplished using a dispersion method. Results similar to that of glucose release were obtained. In vivo imaging of FDG in both uncapsulated FDG and FDG-containing MLV was performed in the brain and the whole body of rats using PET scanner. It was found that the release of FDG from FDG-containing MLV was sustained. In vitro-In vivo correlation was studied using the in vitro release data of FDG from liposomes and in vivo absorption data of FDG from injected liposomes using microPET. Erbium, which is a lanthanide metal, was used as a chelate with DTPA for encapsulation in SUV liposomes for the indirect radiation therapy of cancer. The liposomes were prepared using three different concentrations of soybean lipid (30, 50 and 70 mg/ml). The stability of Er-DTPA SUV liposomes was carried out by storage of the prepared liposomes at three different temperatures (4, 25 and 37 °C). It was found that the release of Er-DTPA complex is temperature dependent, the higher the temperature, the higher the release. There was an inverse relationship between the release of the Er-DTPA complex and the concentration of lipid.

Structure property relations in complex copolymer systems

A thorough investigation was made of the structure-property relation of well-defined statistical, gradient and block copolymers of various compositions. Among the copolymers studied were those which were synthesized using isobornyl acrylate (IBA) and n-butyl acrylate (nBA) monomer units. The copolymers exhibited several unique properties that make them suitable materials for a range of applications. The thermomechanical properties of these new materials were compared to acrylate homopolymers. By the proper choice of the IBA/nBA monomer ratio, it was possible to tune the glass transition temperature of the statistical P(IBA-co-nBA) copolymers. The measured Tg’s of the copolymers with different IBA/nBA monomer ratios followed a trend that fitted well with the Fox equation prediction. While statistical copolymers showed a single glass transition (Tg between -50 and 90 ºC depending on composition), DSC block copolymers showed two Tg’s and the gradient copolymer showed a single, but very broad, glass transition. PMBL-PBA-PMBL triblock copolymers of different composition ratios were also studied and revealed a microphase separated morphology of mostly cylindrical PMBL domains hexagonally arranged in the PBA matrix. DMA studies confirmed the phase separated morphology of the copolymers. Tensile studies showed the linear PMBL-PBA-PMBL triblock copolymers having a relatively low elongation at break that was increased by replacing the PMBL hard blocks with the less brittle random PMBL-r-PMMA blocks. The 10- and 20-arm PBA-PMBL copolymers which were studied revealed even more unique properties. SAXS results showed a mixture of cylindrical PMBL domains hexagonally arranged in the PBA matrix, as well as lamellar. Despite PMBL’s brittleness, the triblock and multi-arm PBA-PMBL copolymers could become suitable materials for high temperature applications due to PMBL’s high glass transition temperature and high thermal stability. The structure-property relation of multi-arm star PBA-PMMA block copolymers was also investigated. Small-angle X-ray scattering revealed a phase separated morphology of cylindrical PMMA domains hexagonally arranged in the PBA matrix. DMA studies found that these materials possess typical elastomeric behavior in a broad range of service temperatures up to at least 250°C. The ultimate tensile strength and the elastic modulus of the 10- and 20-arm star PBA-PMMA block copolymers are significantly higher than those of their 3-arm or linear ABA type counterparts with similar composition, indicating a strong effect of the number of arms on the tensile properties. Siloxane-based copolymers were also studied and one of the main objectives here was to examine the possibility to synthesize trifluoropropyl-containing siloxane copolymers of gradient distribution of trifluoropropyl groups along the chain. DMA results of the PDMS-PMTFPS siloxane copolymers synthesized via simultaneous copolymerization showed that due to the large difference in reactivity rates of 2,4,6-tris(3,3,3-trifluoropropyl)-2,4,6-trimethylcyclotrisiloxane (F) and hexamethylcyclotrisiloxane (D), a copolymer of almost block structure containing only a narrow intermediate fragment with gradient distribution of the component units was obtained. A more dispersed distribution of the trifluoropropyl groups was obtained by the semi-batch copolymerization process, as the DMA results revealed more ‘‘pure gradient type’’ features for the siloxane copolymers which were synthesized by adding F at a controlled rate to the polymerization of the less reactive D. As with trifluoropropyl-containing siloxane copolymers, vinyl-containing polysiloxanes may be converted to a variety of useful polysiloxane materials by chemical modification. But much like the trifluoropropyl-containing siloxane copolymers, as a result of so much difference in the reactivities between the component units 2,4,6-trivinyl-2,4,6-trimethylcyclotrisiloxane (V) and hexamethylcyclotrisiloxane (D), thermal and mechanical properties of the PDMS-PMVS copolymers obtained by simultaneous copolymerization was similar to those of block copolymers. Only the copolymers obtained by semi-batch method showed properties typical for gradient copolymers.

Molecular dynamics simulations of sheared brush-like systems

Mit Hilfe von Molekulardynamik-Simulationen untersuchen wir bürstenartige Systeme unter guten Lösungsmittelbedingungen. Diese Systeme sind, dank ihren vielfältigen Beschaffenheiten, die von Molekularparametern und äußeren Bedingungen abhängig sind, wichtig für viele industrielle Anwendungen. Man vermutet, dass die Polymerbürsten eine entscheidende Rolle in der Natur wegen ihrer einzigartigen Gleiteigenschaften spielen. Ein vergröbertes Modell wird verwendet, um die strukturellen und dynamischen Eigenschaften zweier hochkomprimierter Polymerbürsten, die eine niedrige Reibung aufweisen, zu untersuchen. Allerdings sind die Lubrikationseigenschaften dieser Systeme, die in vielen biologischen Systemen vorhanden sind, beeinflußt. Wir untersuchen so-genannte "weiche Kolloide", die zwischen den beiden Polymerbürsten eingebettet sind, und wie diese Makroobjekte auf die Polymerbürsten wirken.rnrnNicht-Gleichgewichts-Molekulardynamik-Simulationen werden durchgeführt, in denen die hydrodynamischen Wechselwirkungen durch die Anwendung des DPD-Thermostaten mit expliziten Lösungsmittelmolekülen berücksichtigt werden. Wir zeigen, dass die Kenntnis der Gleichgewichtseigenschaften des Systems erlaubt, dynamische Nichtgleichgewichtsigenschaften der Doppelschicht vorherzusagen.rnrnWir untersuchen, wie die effektive Wechselwirkung zwischen kolloidalen Einschlüßen durch die Anwesenheit der Bürsten (in Abhängigkeit der Weichheit der Kolloide und der Pfropfdichte der Bürsten) beeinflußt wird. Als nächsten Schritt untersuchen wir die rheologische Antwort von solchen komplexen Doppelschichten auf Scherung. Wir entwickeln eine Skalen-Theorie, die die Abhängigkeit der makroskopischen Transporteigenschaften und der lateralen Ausdehnung der verankerten Ketten von der Weissenberg Zahl oberhalb des Bereichs, in dem die lineare Antwort-Theorie gilt, voraussagt. Die Vorhersagen der Theorie stimmen gut mit unseren und früheren numerischen Ergebnissen und neuen Experimenten überein. Unsere Theorie bietet die Möglichkeit, die Relaxationszeit der Doppelschicht zu berechnen. Wenn diese Zeit mit einer charakteristischen Längenskala kombiniert wird, kann auch das ''transiente'' (nicht-stationäre) Verhalten beschrieben werden.rnrnrnWir untersuchen die Antwort des Drucktensors und die Deformation der Bürsten während der Scherinvertierung für grosse Weissenberg Zahlen. Wir entwickeln eine Vorhersage für die charakteristische Zeit, nach der das System wieder den stationären Zustand erreicht.rnrnrnElektrostatik spielt eine bedeutende Rolle in vielen biologischen Prozessen. Die Lubrikationseigenschaften der Polymerbürsten werden durch die Anwesenheit langreichweitiger Wechselwirkungen stark beeinflusst. Für unterschiedliche Stärken der elektrostatischen Wechselwirkungen untersuchen wir rheologische Eigenschaften der Doppelschicht und vergleichen mit neutralen Systemen. Wir studieren den kontinuierlichen Übergang der Systemeigenschaften von neutralen zu stark geladenen Bürsten durch Variation der Bjerrumlänge und der Ladungsdichte.

Theoretical study of correlated systems using hybrid functionals

Diese Arbeit unterstreicht das Potential von Hybridfunktionalen (B3LYP) für die Untersuchung einer großen Bandbreite von Systemen. Durch die Einbeziehung der exakten Hartree-Fock Austauschenergie kann B3LYP für molekulare und kristalline Systeme eingesetzt werden. Zum Beispiel können stark korrelierte Systeme mit B3LYP erfolgreich erforscht werden. Die elektronische Struktur von PAHs wurde mit B3LYP Hybriddichtefunktionalen untersucht. Mit der ∆SCF-Methode wurden Elektronenbindungsenergien bestimmt, welche die mit UPS gewonnenen experimentellen Resultate bestätigen und ergänzen. Symmetrieeigenschaften der molekularen Orbitale wurden analysiert, um eine Zuordnung und Einschätzung der zugehörigen Signalstärke zu ermöglichen. Während σ-artige Orbitale nur schwer durch UPS-Messungen an dünnen Filmen detektiert werden können, bieten Rechnungen eine detaillierte Einsicht in die verborgenen Teile der Spektren.rnWeiterhin wurden π−π-Komplexe untersucht, welche von verschiedenen Donor- und Akzeptor-Molekülen gebildet werden. Die Moleküle basieren auf polyzyklischen, aromatischen Kohlenwasserstoffen. Für Ladungstransferkomplexe finden DFT Rechnungen ein Minimum in der Oberfläche der potentiellen Energie. Diese attraktive Wechselwirkung wird durch Coulombanziehung verursacht. Allerdings ist die Coulombanziehung nicht die stärkste Wechselwirkung in Ladungstransferkomplexen. Die Einbeziehung von van der Waals-Korrekturen verbessert den intermolekularen Abstand und die Bindungsenergie.rnEine Verkleinerung der intermolekularen Abstände führt zu einer großen Verschiebung der HOMO- und LUMO-Energie.rnAus der Klasse der kristallinen korrelierten Systeme wurden Rb4O6 und FeSe untersucht. Im Falle von Rb4O6 führen Ladungsordnung und Korrelationen zu einem isolierenden Grundzustand. Das hypothetische druckabhängige Phasendiagramm wurde untersucht. Eine Erhöhung des Drucks führt zu einer vergrößerten Bandlücke. Bei etwa 75 GPa wird die Bandbreite W größer als der Bandabstand U und das System nimmt einen homogen gemischt valenten Zustand mit teilweise besetzten π−π-Orbitalen an. Für Drücke ab 160 GPa wird W sehr viel größer als U und das System wird metallisch.rnIm Fall von FeSe finden wir eine korrelierte und isolierende Phase bei hohen Drücken, während das System bei niedrigen Drücken supraleitendes Verhalten zeigt. Die Berechnungen der Elektronenstruktur mit dem Hybridfunktional B3LYP führt zum korrekten halbleitenden Grundzustand in der NiAs- und MnP-Struktur von FeSe. Die Rolle der Korrelationen, der Stöchiometrie und der Nähe zum Magnetismus wird besprochen. Im Speziellen wird gezeigt, dass die Phase mit NiAs-Struktur starke lokale Korrelationen aufweist, was zu einem halbleitenden Zustand in einem weiten Druckbereich führt.