Computer simulation of ordering and dynamics in liquid crystals in the bulk and close to the surface

The aim of this PhD thesis is to investigate the orientational and dynamical properties of liquid crystalline systems, at molecular level and using atomistic computer simulations, to reach a better understanding of material behavior from a microscopic point view. In perspective this should allow to clarify the relation between the micro and macroscopic properties with the objective of predicting or confirming experimental results on these systems. In this context, we developed four different lines of work in the thesis. The first one concerns the orientational order and alignment mechanism of rigid solutes of small dimensions dissolved in a nematic phase formed by the 4-pentyl,4 cyanobiphenyl (5CB) nematic liquid crystal. The orientational distribution of solutes have been obtained with Molecular Dynamics Simulation (MD) and have been compared with experimental data reported in literature. we have also verified the agreement between order parameters and dipolar coupling values measured in NMR experiments. The MD determined effective orientational potentials have been compared with the predictions of Maier­Saupe and Surface tensor models. The second line concerns the development of a correct parametrization able to reproduce the phase transition properties of a prototype of the oligothiophene semiconductor family: sexithiophene (T6). T6 forms two crystalline polymorphs largely studied, and possesses liquid crystalline phases still not well characterized, From simulations we detected a phase transition from crystal to liquid crystal at about 580 K, in agreement with available experiments, and in particular we found two LC phases, smectic and nematic. The crystal­smectic transition is associated to a relevant density variation and to strong conformational changes of T6, namely the molecules in the liquid crystal phase easily assume a bent shape, deviating from the planar structure typical of the crystal. The third line explores a new approach for calculating the viscosity in a nematic through a virtual exper- iment resembling the classical falling sphere experiment. The falling sphere is replaced by an hydrogenated silicon nanoparticle of spherical shape suspended in 5CB, and gravity effects are replaced by a constant force applied to the nanoparticle in a selected direction. Once the nanoparticle reaches a constant velocity, the viscosity of the medium can be evaluated using Stokes' law. With this method we successfully reproduced experimental viscosities and viscosity anisotropy for the solvent 5CB. The last line deals with the study of order induction on nematic molecules by an hydrogenated silicon surface. Gaining predicting power for the anchoring behavior of liquid crystals at surfaces will be a very desirable capability, as many properties related to devices depend on molecular organization close to surfaces. Here we studied, by means of atomistic MD simulations, the flat interface between an hydrogenated (001) silicon surface in contact with a sample of 5CB molecules. We found a planar anchoring of the first layers of 5CB where surface interactions are dominating with respect to the mesogen intermolecular interactions. We also analyzed the interface 5CB­vacuum, finding a homeotropic orientation of the nematic at this interface.

Modeling the gas and aqueous phase chemistry of the marine boundary layer

Ein eindimensionales numerisches Modell der maritimenGrenzschicht (MBL) wurde erweitert, um chemische Reaktionenin der Gasphase, von Aerosolpartikeln und Wolkentropfen zu beschreiben. Ein Schwerpunkt war dabei die Betrachtung derReaktionszyklen von Halogenen. Soweit Ergebnisse vonMesskampagnen zur Verfuegung standen, wurden diese zurValidierung des Modells benutzt. Die Ergebnisse von frueheren Boxmodellstudien konntenbestaetigt werden. Diese zeigten die saeurekatalysierteAktivierung von Brom aus Seesalzaerosolen, die Bedeutung vonHalogenradikalen fuer die Zerstoerung von O3, diepotentielle Rolle von BrO bei der Oxidation von DMS und dievon HOBr und HOCl in der Oxidation von S(IV). Es wurde gezeigt, dass die Beruecksichtigung derVertikalprofile von meteorologischen und chemischen Groessenvon grosser Bedeutung ist. Dies spiegelt sich darin wider,dass Maxima des Saeuregehaltes von Seesalzaerosolen und vonreaktiven Halogenen am Oberrand der MBL gefunden wurden.Darueber hinaus wurde die Bedeutung von Sulfataerosolen beidem aktiven Recyceln von weniger aktiven zu photolysierbarenBromspezies gezeigt. Wolken haben grosse Auswirkungen auf die Evolution und denTagesgang der Halogene. Dies ist nicht auf Wolkenschichtenbeschraenkt. Der Tagesgang der meisten Halogene ist aufgrundeiner erhoehten Aufnahme der chemischen Substanzen in die Fluessigphase veraendert. Diese Ergebnisse betonen dieWichtigkeit der genauen Dokumentation der meteorologischenBedingungen bei Messkampagnen (besonders Wolkenbedeckungsgrad und Fluessigwassergehalt), um dieErgebnisse richtig interpretieren und mit Modellresultatenvergleichen zu koennen. Dieses eindimensionale Modell wurde zusammen mit einemBoxmodell der MBL verwendet, um die Auswirkungen vonSchiffemissionen auf die MBL abzuschaetzen, wobei dieVerduennung der Abgasfahne parameterisiert wurde. DieAuswirkungen der Emissionen sind am staerksten, wenn sie insauberen Gebieten stattfinden, die Hoehe der MBL gering istund das Einmischen von Hintergrundluft schwach ist.Chemische Reaktionen auf Hintergrundaerosolen spielen nureine geringe Rolle. In Ozeangebieten mit schwachemSchiffsverkehr sind die Auswirkungen auf die Chemie der MBL beschraenkt. In staerker befahrenen Gebieten ueberlappensich die Abgasfahnen mehrerer Schiffe und sorgen fuerdeutliche Auswirkungen. Diese Abschaetzung wurde mitSimulationen verglichen, bei denen die Emissionen alskontinuierliche Quellen behandelt wurden, wie das inglobalen Chemiemodellen der Fall ist. Wenn die Entwicklungder Abgasfahne beruecksichtigt wird, sind die Auswirkungendeutlich geringer da die Lebenszeit der Abgase in der erstenPhase nach Emission deutlich reduziert ist.

Electrode Materials for Ionic Liquid Based-Supercapacitors

The development of safe, high energy and power electrochemical energy-conversion systems can be a response to the worldwide demand for a clean and low-fuel-consuming transport. This thesis work, starting from a basic studies on the ionic liquid (IL) electrolytes and carbon electrodes and concluding with tests on large-size IL-based supercapacitor prototypes demonstrated that the IL-based asymmetric configuration (AEDLCs) is a powerful strategy to develop safe, high-energy supercapacitors that might compete with lithium-ion batteries in power assist-hybrid electric vehicles (HEVs). The increase of specific energy in EDLCs was achieved following three routes: i) the use of hydrophobic ionic liquids (ILs) as electrolytes; ii) the design and preparation of carbon electrode materials of tailored morphology and surface chemistry to feature high capacitance response in IL and iii) the asymmetric double-layer carbon supercapacitor configuration (AEDLC) which consists of assembling the supercapacitor with different carbon loadings at the two electrodes in order to exploit the wide electrochemical stability window (ESW) of IL and to reach high maximum cell voltage (Vmax). Among the various ILs investigated the N-methoxyethyl-N-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide (PYR1(2O1)TFSI) was selected because of its hydrophobicity and high thermal stability up to 350 °C together with good conductivity and wide ESW, exploitable in a wide temperature range, below 0°C. For such exceptional properties PYR1(2O1)TFSI was used for the whole study to develop large size IL-based carbon supercapacitor prototype. This work also highlights that the use of ILs determines different chemical-physical properties at the interface electrode/electrolyte with respect to that formed by conventional electrolytes. Indeed, the absence of solvent in ILs makes the properties of the interface not mediated by the solvent and, thus, the dielectric constant and double-layer thickness strictly depend on the chemistry of the IL ions. The study of carbon electrode materials evidences several factors that have to be taken into account for designing performing carbon electrodes in IL. The heat-treatment in inert atmosphere of the activated carbon AC which gave ACT carbon featuring ca. 100 F/g in IL demonstrated the importance of surface chemistry in the capacitive response of the carbons in hydrophobic ILs. The tailored mesoporosity of the xerogel carbons is a key parameter to achieve high capacitance response. The CO2-treated xerogel carbon X3a featured a high specific capacitance of 120 F/g in PYR14TFSI, however, exhibiting high pore volume, an excess of IL is required to fill the pores with respect to that necessary for the charge-discharge process. Further advances were achieved with electrodes based on the disordered template carbon DTC7 with pore size distribution centred at 2.7 nm which featured a notably high specific capacitance of 140 F/g in PYR14TFSI and a moderate pore volume, V>1.5 nm of 0.70 cm3/g. This thesis work demonstrated that by means of the asymmetric configuration (AEDLC) it was possible to reach high cell voltage up to 3.9 V. Indeed, IL-based AEDLCs with the X3a or ACT carbon electrodes exhibited specific energy and power of ca. 30 Wh/kg and 10 kW/kg, respectively. The DTC7 carbon electrodes, featuring a capacitance response higher of 20%-40% than those of X3a and ACT, respectively, enabled the development of a PYR14TFSI-based AEDLC with specific energy and power of 47 Wh/kg and 13 kW/kg at 60°C with Vmax of 3.9 V. Given the availability of the ACT carbon (obtained from a commercial material), the PYR1(2O1)TFSI-based AEDLCs assembled with ACT carbon electrodes were selected within the EU ILHYPOS project for the development of large-size prototypes. This study demonstrated that PYR1(2O1)TFSI-based AEDLC can operate between -30°C and +60°C and its cycling stability was proved at 60°C up to 27,000 cycles with high Vmax up to 3.8 V. Such AEDLC was further investigated following USABC and DOE FreedomCAR reference protocols for HEV to evaluate its dynamic pulse-power and energy features. It was demonstrated that with Vmax of 3.7 V at T> 30 °C the challenging energy and power targets stated by DOE for power-assist HEVs, and at T> 0 °C the standards for the 12V-TSS and 42V-FSS and TPA 2s-pulse applications are satisfied, if the ratio wmodule/wSC = 2 is accomplished, which, however, is a very demanding condition. Finally, suggestions for further advances in IL-based AEDLC performance were found. Particularly, given that the main contribution to the ESR is the electrode charging resistance, which in turn is affected by the ionic resistance in the pores that is also modulated by pore length, the pore geometry is a key parameter in carbon design not only because it defines the carbon surface but also because it can differentially “amplify” the effect of IL conductivity on the electrode charging-discharging process and, thus, supercapacitor time constant.

The bursting behavior of gas slugs: laboratory and analytical insights into Strombolian volcanic eruptions.

The aim of this thesis is to study how explosive behavior and geophysical signals in a volcanic conduit are related to the development of overpressure in slug-driven eruptions. A first suite of laboratory experiments of gas slugs ascending in analogue conduits was performed. Slugs ascended into a range of analogue liquids and conduit diameters to allow proper scaling to the natural volcanoes. The geometrical variation of the slug in response to the explored variables was parameterised. Volume of gas slug and rheology of the liquid phase revealed the key parameters in controlling slug overpressure at bursting. Founded on these results, a theoretical model to calculate burst overpressure for slug-driven eruptions was developed. The dimensionless approach adopted allowed to apply the model to predict bursting pressure of slugs at Stromboli. Comparison of predicted values with measured data from Stromboli volcano showed that the model can explain the entire spectrum of observed eruptive styles at Stromboli – from low-energy puffing, through normal Strombolian eruptions, up to paroxysmal explosions – as manifestations of a single underlying physical process. Finally, another suite of laboratory experiments was performed to observe oscillatory pressure and forces variations generated during the expansion and bursting of gas slugs ascending in a conduit. Two end-member boundary conditions were imposed at the base of the pipe, simulating slug ascent in closed base (zero magma flux) and open base (constant flux) conduit. At the top of the pipe, a range of boundary conditions that are relevant at a volcanic vent were imposed, going from open to plugged vent. The results obtained illustrate that a change in boundary conditions in the conduit concur to affect the dynamic of slug expansion and burst: an upward flux at the base of the conduit attenuates the magnitude of the pressure transients, while a rheological stiffening in the top-most region of conduit changes dramatically the magnitude of the observed pressure transients, favoring a sudden, and more energetic pressure release into the overlying atmosphere. Finally, a discussion on the implication of changing boundary on the oscillatory processes generated at the volcanic scale is also given.

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

C 3-symmetric discotic liquid crystalline materials for molecular electronics: versatile synthesis and self-organization

This thesis presents the versatile synthesis and self-organization of C3-symmetric discotic nanographene molecules as well as their potential applications as materials in molecular electronics. The details can be described as follows: 1) A novel synthetic strategy towards properly designed C3 symmetric 1,3,5-tris-2’arylbenzene precursors has been developed. After the final planarization by treatment with FeCl3 under mild conditions, for the first time, it became possible to access a variety of new C3-symmetric hexa-peri-hexabenzocoronenes (HBCs) and a series of triangle-shaped nanographenes. D3 symmetric HBC with three alkyl substituents and C2 symmetric HBC with two alkyl substituents were synthesized and found to show the surprising decrease of isotropic points., the self-assembly at the liquid-solid interface displayed a unique zigzag and flower patterns. 2) Triangle-shaped discotics revealed a unique self-assembly behavior in solution, solid state as well as at the solution-substrate interface. A mesophase stability over the broad temperature range with helical supramoelcular arrangement were observed in the bulk state. The honeycomb pattern as the result of novel self-assembly was presented. Triangle-shaped discotics with swallow alkyl tails were fabricated into photovoltaic devices, the supramolecular arrangement upon thermal treatment was found to play a key role in the improvement of solar efficiency. 3) A novel class of C3 symmetric HBCs with alternating polar/apolar substituents was synthesized. Their peculiar self-assembly in solution, in the bulk and on the surface were investigated by NMR techniques, X-ray diffraction as well as different electron microscope techniques. 4) A novel concept for manipulating the intracolumnar stacking of discotics and thus for controlling the helical pitch was presented. A unique staggered stacking in the column was achieved for the first time. Theoretical simulations confirmed this self-organization and predicted that this packing should show the highest charge carrier mobility for all discotics.

Molecular organization of n-cyanobiphenyl liquid crystals on a molybdenite surface

The alignement and anchoring of liquid crystals on solid surfaces is a key problem for modern device technology that until now has been treated empirically, but that can now be tackled by atomistic computer simulations. Molecular dynamics (MD) simulations were used in this thesis work to study two films of 7 and 8 n-alkyl-4’cyanobiphenyl (7CB and 8CB) liquid crystals , with a thickness of 15 nm, confined between two (001) surfaces of MoS2 (molybdenite). The isotropic and nematic phases of both liquid crystals were simulated, and the resulting structures characterized structurally. A new force field was designed to model the interactions between the liquid crystal (LC) molecules and the surface of molybdenite, while an accurate force field developed previously was used to model the 7CB and 8CB molecules. The results show that the (001) molybdenite surface induces a planar orientation in both the liquid crystals. For the nematic phase of 8CB, one of the two solid/LC interfaces is composed of a first layer of molecules aligned parallel to the surface, followed by a second layer of molecules aligned perpendicular to the surface (also called, homeotropic). The effect of the surface appears to be local in nature as it is confined to the first 15 Angström of the LC film. Conversely, for the nematic phase of 7CB, a planar ordering is established into the LC film. The LC molecules at the interface with the molybdenite appear to align preferentially their alkyl chains toward the solid substrate. The resulting tilt angle of molecules was found to be in good agreement with experimental measurements available in literature. Despite the fact that the MD simulations spanned a time range of more than 100 ns, the nematic phases of both 7CB and 8CB were found not to be completely formed. In order to confirm the findings presented in this thesis, we propose to extend the current study.

Entwicklung einer Thermodesorptionseinheit für die GC/MS zur Bestimmung hochreaktiver, biogener Kohlenwasserstoffe und deren Anwendung im Rahmen von Labor- und Feldstudien

Flüchtige organische Verbindungen (volatile organic compounds, VOCs), besonders rnTerpene, gelten als Vorläufer des sekundären organischen Aerosols (secondary rnorganic aerosols, SOA). Terpene werden von Pflanzen zur Abwehr oder zur rnAttraktion von Bestäubern emittiert. Ungesättigten Verbindungen, wie Sesquiterpene, rnsind sehr ozonolyseempfindlich und weisen nur geringe Konzentrationen in der rnAtmosphäre auf. Zudem lassen sie sich mit handelsüblichen rnThermodesorptionseinheiten meist nicht ohne Artefakte nachweisen, da sie eine hohe rnReaktivität mit vielen Oberflächen zeigen. rnDiese Arbeit präsentiert Forschungsergebnisse der qualitativen und rn(semi)quantitativen Auswertungen flüchtiger organischer Verbindungen in rnLaborstudien und biogenen Emissionsproben (Feldmessungen) mittels rnThermodesorption-Gaschromatographie-Massenspektrometrie (TD-GC/MS). Speziell rnzur Analyse hochreaktiver Verbindungen wurde eine Thermodesorptionseinheit für die rnGC/MS entwickelt. Diese besteht aus einer Kryofokussierung, einem Desorptionsofen rnund einer Heizung. Die Steuerung erfolgt über eine eigens dafür geschaffene rnBedienoberf läche von Labview® an einem PC über eine nachgeschaltete SPS rn(speicherprogrammierbare Steuerung). Das komplette Desorptionsprogramm und der rnStart der GC/MS-Analyse wurden automatisiert. rnDie Quantifizierung alle Proben wurde über Diffusionsquellen und rnVergleichsmessungen durch auf Adsorptionsröhrchen aufgebrachte rnVerdünnungsreihen realisiert. Um Informationen über die mutmaßlichen Vorläufer des rnsekundären organischen Aerosols zu erhalten, wurde zudem eine Ozon-Scrubber-rnMethode basierend auf Propen entwickelt. Diese wurde anhand von Standards in einer rnReaktionskammer getestet und in Feldmessungen erfolgreich angewendet. rnQuantitative Analysen zeigen, dass die meisten Terpene so vollständig vor der rnOzonolyse bewahrt werden können. Für hochreaktive Analyte wie α-Humulen oder rnβ-Caryophyllen wurden Wiederfindungsraten von über 80 % erreicht. So konnte die rnTemperaturabhängigkeit der Terpen-Emissionen der Fichte (Picea abies) in rnFeldmessungen nachgewiesen werden. rnEine weitere Anwendungsmöglichkeit wurde mit der Unterscheidung verschiedener rnArten der gleichen Gattung anhand der Emissionsmuster und der möglichen rnAbgrenzung verschiedener Bestäubertypen am Beispiel der Gattung Salvia untersucht. rnDie Emissionsanalysen zeigen, dass eine Zuordnung der Verwandtschaftsverhältnisse rnzusätzlich zu anderen Vergleichen möglich ist. Das gleiche gilt für die Differenzierung rnvon Bestäubertypen. Die Ergebnisse der Feldmessungen wurden durch rnMethodenvergleich zwischen biogenen Emissionsmessungen mit anschließender rnTD-GC/MS-Analyse und Extraktionen der jeweiligen Blüten/Blätter mit rnanschließender GC/MS-Messung bestätigt.

Kinetic modeling and experiments on gas uptake and chemical transformation of organic aerosol in the atmosphere

Aerosolpartikel beeinflussen das Klima durch Streuung und Absorption von Strahlung sowie als Nukleations-Kerne für Wolkentröpfchen und Eiskristalle. Darüber hinaus haben Aerosole einen starken Einfluss auf die Luftverschmutzung und die öffentliche Gesundheit. Gas-Partikel-Wechselwirkunge sind wichtige Prozesse, weil sie die physikalischen und chemischen Eigenschaften von Aerosolen wie Toxizität, Reaktivität, Hygroskopizität und optische Eigenschaften beeinflussen. Durch einen Mangel an experimentellen Daten und universellen Modellformalismen sind jedoch die Mechanismen und die Kinetik der Gasaufnahme und der chemischen Transformation organischer Aerosolpartikel unzureichend erfasst. Sowohl die chemische Transformation als auch die negativen gesundheitlichen Auswirkungen von toxischen und allergenen Aerosolpartikeln, wie Ruß, polyzyklische aromatische Kohlenwasserstoffe (PAK) und Proteine, sind bislang nicht gut verstanden.rn Kinetische Fluss-Modelle für Aerosoloberflächen- und Partikelbulk-Chemie wurden auf Basis des Pöschl-Rudich-Ammann-Formalismus für Gas-Partikel-Wechselwirkungen entwickelt. Zunächst wurde das kinetische Doppelschicht-Oberflächenmodell K2-SURF entwickelt, welches den Abbau von PAK auf Aerosolpartikeln in Gegenwart von Ozon, Stickstoffdioxid, Wasserdampf, Hydroxyl- und Nitrat-Radikalen beschreibt. Kompetitive Adsorption und chemische Transformation der Oberfläche führen zu einer stark nicht-linearen Abhängigkeit der Ozon-Aufnahme bezüglich Gaszusammensetzung. Unter atmosphärischen Bedingungen reicht die chemische Lebensdauer von PAK von wenigen Minuten auf Ruß, über mehrere Stunden auf organischen und anorganischen Feststoffen bis hin zu Tagen auf flüssigen Partikeln. rn Anschließend wurde das kinetische Mehrschichtenmodell KM-SUB entwickelt um die chemische Transformation organischer Aerosolpartikel zu beschreiben. KM-SUB ist in der Lage, Transportprozesse und chemische Reaktionen an der Oberfläche und im Bulk von Aerosol-partikeln explizit aufzulösen. Es erforder im Gegensatz zu früheren Modellen keine vereinfachenden Annahmen über stationäre Zustände und radiale Durchmischung. In Kombination mit Literaturdaten und neuen experimentellen Ergebnissen wurde KM-SUB eingesetzt, um die Effekte von Grenzflächen- und Bulk-Transportprozessen auf die Ozonolyse und Nitrierung von Protein-Makromolekülen, Ölsäure, und verwandten organischen Ver¬bin-dungen aufzuklären. Die in dieser Studie entwickelten kinetischen Modelle sollen als Basis für die Entwicklung eines detaillierten Mechanismus für Aerosolchemie dienen sowie für das Herleiten von vereinfachten, jedoch realistischen Parametrisierungen für großskalige globale Atmosphären- und Klima-Modelle. rn Die in dieser Studie durchgeführten Experimente und Modellrechnungen liefern Beweise für die Bildung langlebiger reaktiver Sauerstoff-Intermediate (ROI) in der heterogenen Reaktion von Ozon mit Aerosolpartikeln. Die chemische Lebensdauer dieser Zwischenformen beträgt mehr als 100 s, deutlich länger als die Oberflächen-Verweilzeit von molekularem O3 (~10-9 s). Die ROIs erklären scheinbare Diskrepanzen zwischen früheren quantenmechanischen Berechnungen und kinetischen Experimenten. Sie spielen eine Schlüsselrolle in der chemischen Transformation sowie in den negativen Gesundheitseffekten von toxischen und allergenen Feinstaubkomponenten, wie Ruß, PAK und Proteine. ROIs sind vermutlich auch an der Zersetzung von Ozon auf mineralischem Staub und an der Bildung sowie am Wachstum von sekundären organischen Aerosolen beteiligt. Darüber hinaus bilden ROIs eine Verbindung zwischen atmosphärischen und biosphärischen Mehrphasenprozessen (chemische und biologische Alterung).rn Organische Verbindungen können als amorpher Feststoff oder in einem halbfesten Zustand vorliegen, der die Geschwindigkeit von heterogenen Reaktionenen und Mehrphasenprozessen in Aerosolen beeinflusst. Strömungsrohr-Experimente zeigen, dass die Ozonaufnahme und die oxidative Alterung von amorphen Proteinen durch Bulk-Diffusion kinetisch limitiert sind. Die reaktive Gasaufnahme zeigt eine deutliche Zunahme mit zunehmender Luftfeuchte, was durch eine Verringerung der Viskosität zu erklären ist, bedingt durch einen Phasenübergang der amorphen organischen Matrix von einem glasartigen zu einem halbfesten Zustand (feuchtigkeitsinduzierter Phasenübergang). Die chemische Lebensdauer reaktiver Verbindungen in organischen Partikeln kann von Sekunden bis zu Tagen ansteigen, da die Diffusionsrate in der halbfesten Phase bei niedriger Temperatur oder geringer Luftfeuchte um Größenordnungen absinken kann. Die Ergebnisse dieser Studie zeigen wie halbfeste Phasen die Auswirkung organischeer Aerosole auf Luftqualität, Gesundheit und Klima beeinflussen können. rn

Simulations of amphiphilic peptides at the air-water interface

Amphiphile Peptide, Pro-Glu-(Phe-Glu)n-Pro, Pro-Asp-(Phe-Asp)n-Pro, und Phe-Glu-(Phe-Glu)n-Phe, können so aus n alternierenden Sequenzen von hydrophoben und hydrophilen Aminosäuren konstruiert werden, dass sie sich in Monolagen an der Luft-Wasser Grenzfläche anordnen. In biologischen Systemen können Strukturen an der organisch-wässrigen Grenzfläche als Matrix für die Kristallisation von Hydroxyapatit dienen, ein Vorgang der für die Behandlung von Osteoporose verwendet werden kann. In der vorliegenden Arbeit wurden Computersimulationenrneingesetzt, um die Strukturen und die zugrunde liegenden Wechselwirkungen welche die Aggregation der Peptide auf mikroskopischer Ebene steuern, zu untersuchen. Atomistische Molekulardynamik-Simulationen von einzelnen Peptidsträngen zeigen, dass sie sich leicht an der Luft-Wasser Grenzfläche anordnen und die Fähigkeit haben, sich in β-Schleifen zu falten, selbst für relativ kurze Peptidlängen (n = 2). Seltene Ereignisse wie diese (i.e. Konformationsänderungen) erfordern den Einsatz fortgeschrittener Sampling-Techniken. Hier wurde “Replica Exchange” Molekulardynamik verwendet um den Einfluss der Peptidsequenzen zu untersuchen. Die Simulationsergebnisse zeigten, dass Peptide mit kürzeren azidischen Seitenketten (Asp vs. Glu) gestrecktere Konformationen aufwiesen als die mit längeren Seitenketten, die in der Lage waren die Prolin-Termini zu erreichen. Darüber hinaus zeigte sich, dass die Prolin-Termini (Pro vs. Phe) notwendig sind, um eine 2D-Ordnung innerhalb derrnAggregate zu erhalten. Das Peptid Pro-Asp-(Phe-Asp)n-Pro, das beide dieser Eigenschaften enthält, zeigt das geordnetste Verhalten, eine geringe Verdrehung der Hauptkette, und ist in der Lage die gebildeten Aggregate durch Wasserstoffbrücken zwischen den sauren Seitenketten zu stabilisieren. Somit ist dieses Peptid am besten zur Aggregation geeignet. Dies wurde auch durch die Beurteilung der Stabilität von experimentnah-aufgesetzten Peptidaggregaten, sowie der Neigung einzelner Peptide zur Selbstorganisation von anfänglich ungeordneten Konfigurationen unterstützt. Da atomistische Simulationen nur auf kleine Systemgrößen und relativ kurze Zeitskalen begrenzt sind, wird ein vergröbertes Modell entwickelt damit die Selbstorganisation auf einem größeren Maßstab studiert werden kann. Da die Selbstorganisation an der Grenzfläche vonrnInteresse ist, wurden existierenden Vergröberungsmethoden erweitert, um nicht-gebundene Potentiale für inhomogene Systeme zu bestimmen. Die entwickelte Methode ist analog zur iterativen Boltzmann Inversion, bildet aber das Update für das Interaktionspotential basierend auf der radialen Verteilungsfunktion in einer Slab-Geometrie und den Breiten des Slabs und der Grenzfläche. Somit kann ein Kompromiss zwischen der lokalen Flüssigketsstruktur und den thermodynamischen Eigenschaften der Grenzfläche erreicht werden. Die neue Methode wurde für einen Wasser- und einen Methanol-Slab im Vakuum demonstriert, sowie für ein einzelnes Benzolmolekül an der Vakuum-Wasser Grenzfläche, eine Anwendung die von besonderer Bedeutung in der Biologie ist, in der oft das thermodynamische/Grenzflächenpolymerisations-Verhalten zusätzlich der strukturellen Eigenschaften des Systems erhalten werden müssen. Daraufrnbasierend wurde ein vergröbertes Modell über einen Fragment-Ansatz parametrisiert und die Affinität des Peptids zur Vakuum-Wasser Grenzfläche getestet. Obwohl die einzelnen Fragmente sowohl die Struktur als auch die Wahrscheinlichkeitsverteilungen an der Grenzfläche reproduzierten, diffundierte das Peptid als Ganzes von der Grenzfläche weg. Jedoch führte eine Reparametrisierung der nicht-gebundenen Wechselwirkungen für eines der Fragmente der Hauptkette in einem Trimer dazu, dass das Peptid an der Grenzfläche blieb. Dies deutet darauf hin, dass die Kettenkonnektivität eine wichtige Rolle im Verhalten des Petpids an der Grenzfläche spielt.

Atomistic simulations of liquid crystals in the bulk and at their interfaces

In this thesis, we have dealt with several problems concerning liquid crystals (LC) phases, either in the bulk or at their interfaces, by the use of atomistic molecular dynamics (MD) simulations. We first focused our attention on simulating and characterizing the bulk smectic phase of 4-n-octyl-4'-cyanobiphenyl (8CB), allowing us to investigate the antiparallel molecular arrangement typical of SmAd smectic phases. A second topic of study was the characterization of the 8CB interface with vacuum by simulating freely suspended thin films, which allowed us to determine the influence of the interface on the orientational and positional order. Then we investigated the LC-water and LC-electrolyte water solution interface. This interface has recently found application in the development of sensors for several compounds, including biological molecules, and here we tried to understand the re-orientation mechanism of LC molecules at the interface which is behind the functioning of these sensors. The characterization of this peculiar interface has incidentally led us to develop a polarizable force field for the pentyl-cyanobiphenyl mesogen, whose process of parametrization and validation is reported here in detail. We have shown that this force field is a significant improvement over its previous, static charge non polarizable version in terms of density, orientational order parameter and translational diffusion.

Total OH reactivity measurements at the biosphere-atmosphere interface

The biosphere emits copiously volatile organic compounds (VOCs) into the atmosphere, which are removed again depending on the oxidative capacity of the atmosphere and physical processes such as mixing, transport and deposition. Biogenic VOCs react with the primary oxidant of the atmosphere, the hydroxyl radical (OH), and potentially lead to the formation tropospheric ozone and aerosol, which impact regional climate and air quality. The rate of OH decay in the atmosphere, the total OH reactivity is a function of the atmospheric, reactive compound's concentration and reaction velocity with OH. One way to measure the total OH reactivity, the total OH sink, is with the Comparative Reactivity Method - CRM. Basically, the reaction of OH with a reagent (here pyrrole) in clean air and in the presence of atmospheric, reactive molecules is compared. This thesis presents measurements of the total OH reactivity at the biosphere-atmosphere interface to analyze various influences and driving forces. For measurements in natural environment the instrument was automated and a direct, undisturbed sampling method developed. Additionally, an alternative detection system was tested and compared to the originally used detector (Proton Transfer Reaction-Mass Spectrometer, PTR-MS). The GC-PID (Gas Chromatographic Photo-Ionization Detector) was found as a smaller, less expensive, and robust alternative for total OH reactivity measurements. The HUMPPA-COPEC 2010 measurement campaign in the Finish forest was impacted by normal boreal forest emissions as well as prolonged heat and biomass burning emissions. The measurement of total OH reactivity was compared with a comprehensive set of monitored individual species ambient concentration levels. A significant discrepancy between those individually measured OH sinks and the total OH reactivity was observed, which was characterized in detail by the comparison of within and above the forest canopy detected OH reactivity. Direct impact of biogenic emissions on total OH reactivity was examined on Kleiner Feldberg, Germany, 2011. Trans-seasonal measurements of an enclosed Norway spruce branch were conducted via PTR-MS, for individual compound's emission rates, and CRM, for total OH reactivity emission fluxes. Especially during summertime, the individually monitored OH sink terms could not account for the measured total OH reactivity. A controlled oxidation experiment in a low NOx environment was conducted in the EUPHORE reaction chamber (CHEERS, Spain 2011). The concentration levels of the reactant isoprene and its major products were monitored and compared to total OH reactivity measurements as well as to the results of two models. The individually measured compounds could account for the total OH reactivity during this experiment as well as the traditional model-degradation scheme for isoprene (MCM 3.2). Due to previous observations of high OH levels in the isoprene-rich environment of the tropics, a novel isoprene mechanism was recently suggested. In this mechanism (MIME v4) additional OH is generated during isoprene oxidation, which could not be verified in the conditions of the CHEERS experiment.

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.

The dynamics of ultra-thin polymer films in different confinements studied by resonance enhanced dynamic light scattering

Die Kontroverse über den Glasübergang im Nanometerbereich, z. B. die Glas¬über¬gangs-temperatur Tg von dünnen Polymerfilmen, ist nicht vollständig abgeschlossen. Das dynamische Verhalten auf der Nanoskala ist stark von den einschränkenden Bedingungen abhängig, die auf die Probe wirken. Dünne Polymerfilme sind ideale Systeme um die Dynamik von Polymerketten unter der Einwirkung von Randbedingungen zu untersuchen, wie ich sie in dieser Arbeit variiert habe, um Einblick in dieses Problem zu erhalten.rnrnResonanzverstärkte dynamische Lichtstreuung ist eine Methode, frei von z.B. Fluoreszenzmarkern, die genutzt werden kann um in dünnen Polymerfilmen dynamische Phänomene

Ventilation-perfusion ratio in perflubron during partial liquid ventilation

BACKGROUND: Functional magnetic resonance imaging (fMRI) of fluorine-19 allows for the mapping of oxygen partial pressure within perfluorocarbons in the alveolar space (Pao(2)). Theoretically, fMRI-detected Pao(2) can be combined with the Fick principle approach, i.e., a mass balance of oxygen uptake by ventilation and delivery by perfusion, to quantify the ventilation-perfusion ratio (Va/Q) of a lung region: The mixed venous blood and the inspiratory oxygen fraction, which are equal for all lung regions, are measured. In addition, the local expiratory oxygen fraction and the end capillary oxygen content, both of which may differ between the lung regions, are calculated using the fMRI-detected Pao(2). We investigated this approach by numerical simulations and applied it to quantify local Va/Q in the perfluorocarbons during partial liquid ventilation. METHODS: Numerical simulations were performed to analyze the sensitivity of the Va/Q calculation and to compare this approach with another one proposed by Rizi et al. in 2004 (Magn Reson Med 2004;52:65-72). Experimentally, the method was used during partial liquid ventilation in 7 anesthetized pigs. The Pao(2) distribution in intraalveolar perflubron was measured by fluorine-19 MRI. Respiratory gas fractions together with arterial and mixed venous blood samples were taken to quantify oxygen partial pressure and content. Using the Fick principle, the local Va/Q was estimated. The impact of gravity (nondependent versus dependent) of perflubron dose (10 vs 20 mL/kg body weight) and of inspired oxygen fraction (Fio(2)) (0.4-1.0) on Va/Q was examined. RESULTS: In numerical simulations, the Fick principle proved to be appropriate over the Va/Q range from 0.02 to 2.5. Va/Q values were in acceptable agreement with the method published by Rizi et al. In the experimental setting, low mean Va/Q values were found in perflubron (confidence interval [CI] 0.08-0.29 with 20 mL/kg perflubron). At this dose, Va/Q in the nondependent lung was higher (CI 0.18-0.39) than in the dependent lung regions (CI 0.06-0.16; P = 0.006; Student t test). Differences depending on Fio(2) or perflubron dose were, however, small. CONCLUSION: The results show that derivation of Va/Q from local Po(2) measurements using fMRI in perflubron is feasible. The low detected Va/Q suggests that oxygen transport into the perflubron-filled alveolar space is significantly restrained.