Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Thermal Casimir effect at fluid interfaces: Fluctuation-induced forces between anisotropic colloids

We study the effective interaction between two ellipsoidal particles at the interface of two fluid phases which are mediated by thermal fluctuations of the interface. Within a coarse-grained picture, the properties of fluid interfaces are very well described by an effective capillary wave Hamiltonian which governs both the equilibrium interface configuration and the thermal fluctuations (capillary waves) around this equilibrium (or mean-field) position. As postulated by the Goldstone theorem the capillary waves are long-range correlated. The interface breaks the continuous translational symmetry of the system, and in the limit of vanishing external fields - like gravity - it has to be accompanied by easily excitable long wavelength (Goldstone) modes – precisely the capillary waves. In this system the restriction of the long-ranged interface fluctuations by particles gives rise to fluctuation-induced forces which are equivalent to interactions of Casimir type and which are anisotropic in the interface plane. Since the position and the orientation of the colloids with respect to the interface normal may also fluctuate, this system is an example for the Casimir effect with fluctuating boundary conditions. In the approach taken here, the Casimir interaction is rewritten as the interaction between fluctuating multipole moments of an auxiliary charge density-like field defined on the area enclosed by the contact lines. These fluctuations are coupled to fluctuations of multipole moments of the contact line position (due to the possible position and orientational fluctuations of the colloids). We obtain explicit expressions for the behavior of the Casimir interaction at large distances for arbitrary ellipsoid aspect ratios. If colloid fluctuations are suppressed, the Casimir interaction at large distances is isotropic, attractive and long ranged (double-logarithmic in the distance). If, however, colloid fluctuations are included, the Casimir interaction at large distances changes to a power law in the inverse distance and becomes anisotropic. The leading power is 4 if only vertical fluctuations of the colloid center are allowed, and it becomes 8 if also orientational fluctuations are included.

Untersuchung der Filmbildung aus Polymerdispersionen mit Hilfe der forcierten Rayleighstreuung

Der Filmbildungsprozess wässriger Polymerdispersionen wurde mit forcierter Rayleighstreuung untersucht.Diffusionskoeffizienten D von Farbstoffsonden sind damit zwischen 10^-21 und 10^-9 m^2s^-1 zugänglich. Das Diffusionsverhalten der hydrophoben FarbstoffsondeAberchrome in feuchten und trockenen Filmen aus Poly(n-butylmethacrylat-co-acrylsäure)-Dispersionen sowie bei der Wiederbefeuchtung trockener Filme wurde untersucht.Die Dynamik von Aberchrome weicht in feuchten Filmen von Fickscher Diffusion ab. Dies äußert sich in Abweichungen vomcharakteristischen q^2-Verhalten der Relaxationsrate tau^-1 (tau^-1 = Dq^2; q:Streuvektor) und im Auftreten gestrecktexponentieller Intensitätsrelaxationskurven. Diese Anomalie wurde mit einem apparenten, längenskalenabhängigen Diffusionskoeffizienten Dapp(Lambda) (Lambda= 2Pi/q) beschrieben, der für Lambda -> 0 Werte annimmt, die einem homogen hydroplastifizierten Polymermaterial entsprechen, während Dapp(Lambda) für Lambda -> Unendlich stark anwächst. Diese Anomalien verschwinden bei Wassergehalten entsprechend der Polymerwasserlöslichkeit. Weiteres Trocknen führt zum Absinken des Fickschen-Diffusionskoeffizienten auf einen Grenzwert für trockene Filme. Die Ergebnisse konnten mit einem Zwei-Zustands-Modell beschrieben werden: Die Sonde diffundiert Ficksch in einer hydrophilen Grenzflächenphase und einer langsamen Polymerphase. Austausch zwischendiesen Phasen ist ohne Einschränkung möglich. Das Modell erlaubt die Quantifizierung des Einflusses des Trocknungsprozesses auf Polymer- und Grenzflächenphaseneigenschaften.Dies wurde durch systematische Veränderungen der Grenzflächeneigenschaften demonstriert. Dies geschah durch Acrylsäuregehaltvariationen in den Poly(n-butylmethacrylat-co-acrylsäure)-Dispersionen und Emulgatorbedeckungsgradvariationen. In beiden Fällen reflektieren sich Dispersionsveränderungen in einer Veränderung der Modellparameter in der Grenzflächenphase.

Fourier-transform rheology on anionically synthesised polymer melts and solutions of various topology

Aim of this thesis was to further extend the applicability of the Fourier-transform (FT) rheology technique especially for non-linear mechanical characterisation of polymeric materials on the one hand and to investigated the influence of the degree of branching on the linear and non-linear relaxation behaviour of polymeric materials on the other hand. The latter was achieved by employing in particular FT-rheology and other rheological techniques to variously branched polymer melts and solutions. For these purposes, narrowly distributed linear and star-shaped polystyrene and polybutadiene homo-polymers with varying molecular weights were anionically synthesised using both high-vacuum and inert atmosphere techniques. Furthermore, differently entangled solutions of linear and star-shaped polystyrenes in di-sec-octyl phthalate (DOP) were prepared. The several linear polystyrene solutions were measured under large amplitude oscillatory shear (LAOS) conditions and the non-linear torque response was analysed in the Fourier space. Experimental results were compared with numerical predictions performed by Dr. B. Debbaut using a multi-mode differential viscoelastic fluid model obeying the Giesekus constitutive equation. Apart from the analysis of the relative intensities of the harmonics, a detailed examination of the phase information content was developed. Further on, FT-rheology allowed to distinguish polystyrene melts and solutions due to their different topologies where other rheological measurements failed. Significant differences occurred under LAOS conditions as particularly reflected in the phase difference of the third harmonic, Ħ3, which could be related to shear thinning and shear thickening behaviour.

Plasmonenresonanzen von sichelförmigen metallischen Nanoobjekten

Metallische Objekte in der Größenordnung der optischen Wellenlänge zeigen Resonanzen im optischen Spektralbereich. Mit einer Kombination aus Kolloidlithographie, Metallfilmbedampfung und reaktivem Ionenstrahl¨atzen wurden Nanosicheln aus Gold bzw. Silber mit identischer Form und Orientierung in Sichelform mit einer Größe von 60nm bis 400nm hergestellt. Der Öffnungswinkel der Nanosicheln lässt sich kontinuierlich einstellen. Durch die einheitliche Orientierung lassen sich Messungen am Ensemble direkt auf das Verhalten des Einzelobjektes übertragen, wie ein Vergleich der Extinktionsspektren einer Ensemblemessung am UV/Vis/NIR-Spektrometer mit einer Einzelpartikelmessung in einem konfokalen Mikroskop zeigt. Die optische Antwort der Nanosicheln wurde als zwei-dimensionales Modell mit einer Finite Elemente Methode berechnet. Das Ergebnis sind mehrere polarisationsabhängige Resonanzen im optischen Spektrum. Diese lassen sich durch Variation des Öffnungswinkels und der Gr¨oße der Nanosichel verschieben. Durch Beleuchten lassen sich plasmonische Schwingungen anregen, die ein stark lokalisiertes Nahfeld an den Spitzen und in der Öffnung der Nanosicheln erzeugen. Das Nahfeld der Partikelresonanz wurde mit einer Fotolackmethode nachgewiesen. Die Untersuchungen am UV/Vis/NIR-Spektrometer zeigen mehrere polarisationsabhängige Resonanzen im Spektralbereich von 300 nm bis 3200 nm. Die Resonanzen der Nanosicheln lassen sich durch den Öffnungswinkel und den Durchmesser in der Größenordnung der Halbwertbreite im optischen Spektrum verschieben. In der Anwendung als Chemo- bzw. Biosensor zeigen Gold-Nanosicheln eine ähnliche Empfindlichkeit wie vergleichbare Sensoren auf der Basis von dünnen Metallstrukturen. Das Nahfeld zeichnet sich durch eine starke Lokalisierung aus und dringt, je nach Multipolordnung, zwischen 14 nm und 70 nm in die Umgebung ein. Quantenpunkte wurden an das Nahfeld der Nanosicheln gekoppelt. Die Emission der Quantenpunkte bei einer Wellenlänge von 860nm wird durch die Resonanz der Nanosicheln verstärkt. Die Nanosicheln wurden als optische Pinzette eingesetzt. Bei einer Anregung mit einem Laser bei einer Wellenlänge von 1064 nm wurden Polystyrolkolloide mit einem Durchmesser von 40 nm von den resonanten Nanosicheln eingefangen. Die Nanosicheln zeigen außergewöhnliche optische Eigenschaften, die mithilfe der Geometrieparameter über einen großen Bereich verändert werden können. Die ersten Anwendungen haben Anknüpfungspunkte zur Verwendung in der Sensorik, Fluoreszenzspektroskopie und als optische Pinzette aufgezeigt.

Generierung leicht dispergierbarer Inhalationspulver mittels Sprühtrocknung

Das Ziel dieser Arbeit besteht darin, die Möglichkeiten der Sprühtrocknung für die Generierung von Inhalationspulvern zur Therapie von Lungenkrankheiten zu nutzen. Die Erzeugung von physikalisch stabilen und leicht dispergierbaren Partikeln steht hierbei im Vordergrund. Aufgrund von physiko-chemischen Untersuchungen (Glasübergangstemperatur, Fragilität, Relaxationsverhalten, Hygroskopizität) unterschiedlicher amorpher Hilfsstoffe (Lactose, Raffinose, Dextrane, Cyclodextrine) ist für Hydroxypropyl-β-Cyclodextrin das größte Potential für die Stabilisierung eines Wirkstoffes innerhalb einer amorphen Matrix erkennbar. Sprühgetrocknete Partikel weisen im Vergleich zu strahlgemahlenen Partikeln günstigere Dispergier- und Depositionseigenschaften auf. Dies ist vorrangig auf größere Berührungsflächen zwischen strahlgemahlenen Partikeln zurückzuführen. Kugelförmige sprühgetrocknete Partikel besitzen dagegen aufgrund einer punktförmigen Berührung geringere Haftkräfte. Versuche mit unterschiedlich stark gefalteten Partikeloberflächen weisen auf geringere Haftkräfte hin, wenn sich die Partikel an Stellen geringerer Krümmungsradien berühren. Dispergierversuche in einer definierten Rohrströmung (Deagglomerator) lassen auf einen kaskadenartigen Agglomeratzerfall schließen. Durch Sprüheinbettung unterschiedlicher Modellwirkstoffe (Salbutamolsulfat, Ipratropiumbromid, Budesonid) in Hydroxypropyl-β-Cyclodextrin konnten sowohl Einzelformulierungen als auch eine Kombinationsformulierung mit allen drei Wirkstoffen erzeugt werden. Diese weisen bei einem Wirkstoffgehalt bis max. 14% selbst nach vierwöchiger Offenlagerung bei 40°C und 75% r.F. keine bzw. nur geringfügige Veränderungen in der „Fine Particle Dose“ (FPD) auf. Die „Fine Particle Fraction“ (FPF) liegt bei diesen Formulierungen im Bereich von 40% bis 75%. In Verbindung mit einem geeigneten Pack- bzw. Trockenmittel, ist hierbei mit einer physikalischen Stabilität zu rechnen, die eine sinnvolle Produktlaufzeit eines Inhalationspulvers ermöglicht. Formulierungen mit höheren Wirkstoffkonzentrationen zeigen dagegen stärkere Veränderungen nach Stresslagerung. Als Beispiel einer kristallinen Sprühtrocknungsformulierung konnte ein Pulver bestehend aus Mannitol und Budesonid erzeugt werden.

Electronic transport in the heavy fermion superconductors UPd 2 Al 3 and UNi 2 Al 3

This work addresses the electronical properties of the superconductors UPd2Al3 and UNi2Al3 on the basis of thin film experiments. These isotructural compounds are ideal candiates to study the interplay of magnetism and superconductivity due to the differences of their magnetically ordered states, as well as the experimental evidence for a magnetic pairing mechanism in UPd2Al3. Epitaxial thin film samples of UPd2Al3 and UNi2Al3 were prepared using UHV Molecular Beam Epitaxy (MBE). For UPd2Al3, the change of the growth direction from the intrinsic (001) to epitaxial (100) was predicted and sucessfully demonstrated using LaAlO3 substrates cut in (110) direction. With optimized deposition process parameters for UPd2Al3 (100) on LaAlO3 (110) superconducting samples with critical temperatures up to Tc = 1.75K were obtained. UPd2Al3-AlOx-Ag mesa junctions with superconducting base electrode were prepared and shown to be in the tunneling regime. However, no signatures of a superconducting density of states were observed in the tunneling spectra. The resistive superconducting transition was probed for a possible dependence on the current direction. In contrast to UNi2Al3, the existence of such feature was excluded in UPd2Al3 (100) thin films. The second focus of this work is the dependence of the resisitive transition in UNi2Al3 (100) thin films on the current direction. The experimental fact that the resisitive transition occurs at slightly higher temperatures for I║a than for I║c can be explained within a model of two weakly coupled superconducting bands. Evidence is presented for the key assumption of the two-band model, namely that transport in and out of the ab-plane is generated on different, weakly coupled parts of the Fermi surface. Main indications are the angle dependence of the superconducting transition and the dependence of the upper critical field Bc2 on current and field orientation. Additionally, several possible alternative explanations for the directional splitting of the transition are excluded in this work. An origin due to scattering on crystal defects or impurities is ruled out, likewise a relation to ohmic heating or vortex dynamics. The shift of the transition temperature as function of the current density was found to behave as predicted by the Ginzburg-Landau theory for critical current depairing, which plays a significant role in the two-band model. In conclusion, the directional splitting of the resisitive transition has to be regarded an intrinsic and unique property of UNi2Al3 up to now. Therefore, UNi2Al3 is proposed as a role model for weakly coupled multiband superconductivity. Magnetoresistance in the normalconducting state was measured for UPd2Al3 and UNi2Al3. For UNi2Al3, a negative contribution was observed close to the antiferromagnetic ordering temperature TN only for I║a, which can be associated to reduced spin-disorder scattering. In agreement with previous results it is concluded that the magnetic moments have to be attributed to the same part of the Fermi surface which generates transport in the ab-plane.

Prozessevaluation einer kognitiv-verhaltenstherapeutischen Gruppentherapie bei Diabetes und Depression : Entwicklung und Validierung der Patienten- und Therapeuten- Gruppentherapiestundenbögen (GTS-P, GTS-T, GTS-TP) zur Vorhersage des Therapieerfolgs

Untersucht werden Prozess-Ergebnis-Zusammenhänge einer kognitiv-verhaltenstherapeutischen Gruppentherapie für Diabetes und Depression im Rahmen der DAD-Studie. rnAufgrund des Mangels an geeigneten Erhebungsinstrumenten der validen, ökonomischen und komplementären Sitzungsbewertung von Gruppenpatienten und -therapeuten wurden angelehnt an einen Patienten- (GTS-P) zwei Therapeutenstundenbögen entwickelt: der GTS-T zur Bewertung der Gesamtgruppe und der GTS-TP zur Bewertung einzelner Patienten. Die GTS-Bögen zeigen bei der Überprüfung der Testgüte insgesamt gute Itemparameter und Reliabilitäten. Das in den exploratorischen Faktorenanaylsen des GTS-P identifizierte zweifaktorielle Modell (1. wahrgenommene Zuversicht hinsichtlich der Gruppentherapie, 2. wahrgenommene persönliche Beteiligung) kann in den konfirmatorischen Faktorenanalysen bestätigt werden. Dazu wurden GTS-P-Daten aus einer Untersuchung mit Patienten mit somatoformen Störungen (Schulte, 2001) einbezogen. Den Ergebnissen der Item- und Faktorenanalysen folgend, wurden zwei Items des GTS-P und zwei weitere Items des GTS-T aus den Instrumenten ausgeschlossen. Für den GTS-T zeigt sich eine einfaktorielle, für den GTS-TP eine zum GTS-P parallele zweifaktorielle Struktur. rnIn den Mehrebenenanalysen zur Vorhersage des Therapieergebnisses (Post-Depressionssymptomatik) zeigt sich die Skala Zuversicht des GTS-P zu Therapiebeginn (1.-4. Sitzung) kontrolliert an der Skala Beteiligung und der Prä-Symptomatik, als valider Prädiktor. Das Item 5 „Anregungen“ (Skala Zuversicht) und Item 2 „Aktive Mitwirkung“ (Skala Beteiligung) sind am stärksten an diesem Effekt beteiligt, da diese Itemkombination das Therapieergebnis ebenfalls valide vorhersagen kann. Die Prognose ist schon durch die Werte der ersten Gruppentherapiesitzungen in der Remoralisierungsphase (Howard et al., 1993) möglich und verbessert sich nicht bei Berücksichtigung aller 10 Gruppensitzungen. Die Therapeutenbögen zeigen keine prädiktive Validität. Bedeutsame Zusammenhänge der Patienten- und Therapeutenbewertungen finden sich lediglich für den GTS-P und GTS-TP. Weitere Prädiktoren, wie der Diabetestyp, Diabeteskomplikationen und die Adhärenz, konnten nicht zur Verbesserung der Vorhersage beitragen. Für sekundär überprüfte Kriterien gelang die Prognose lediglich für ein weiteres Maß der Depressionssymptomatik und für eine Gesamtbewertung der Gruppentherapie durch die Patienten zu Therapieende. Bei der deskriptiven Betrachtung der Prozessqualität der DAD-Gruppentherapien zeigen sich positive, über den Verlauf der Gruppe zunehmende und nach Therapiephasen differenzierbare Bewertungsverläufe. rnDie Ergebnisse der Studie sprechen für die Relevanz von unspezifischen Wirkfaktoren für das Therapieergebnis von kognitiv-behavioralen Gruppentherapien. Die von den Gruppenpatienten wahrgenommene Zuversicht und Beteiligung als Zeichen der Ansprechbarkeit auf die Therapie sollte mit Hilfe von Stundenbögen, wie den GTS-Bögen, von Gruppentherapeuten zur Prozessoptimierung und Prävention von Therapieabbrüchen und Misserfolgen beachtet werden. rn

Exact quantum Monte Carlo methods for correlated electron systems

Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.

Unraveling efficiency-limiting processes in organic solar cells by ultrafast spectroscopy - impact of chemical structure and morphology on photophysics and efficiency

Diese Arbeit widmet sich der Untersuchung der photophysikalischen Prozesse, die in Mischungen von Elektronendonoren mit Elektronenakzeptoren zur Anwendung in organischen Solarzellen auftreten. Als Elektronendonoren werden das Copolymer PBDTTT-C, das aus Benzodithiophen- und Thienothiophene-Einheiten besteht, und das kleine Molekül p-DTS(FBTTh2)2, welches Silizium-überbrücktes Dithiophen, sowie fluoriertes Benzothiadiazol und Dithiophen beinhaltet, verwendet. Als Elektronenakzeptor finden ein planares 3,4:9,10-Perylentetracarbonsäurediimid-(PDI)-Derivat und verschiedene Fullerenderivate Anwendung. PDI-Derivate gelten als vielversprechende Alternativen zu Fullerenen aufgrund der durch chemische Synthese abstimmbaren strukturellen, optischen und elektronischen Eigenschaften. Das gewichtigste Argument für PDI-Derivate ist deren Absorption im sichtbaren Bereich des Sonnenspektrums was den Photostrom verbessern kann. Fulleren-basierte Mischungen übertreffen jedoch für gewöhnlich die Effizienz von Donor-PDI-Mischungen.rnUm den Nachteil der PDI-basierten Mischungen im Vergleich zu den entsprechenden Fulleren-basierten Mischungen zu identifizieren, werden die verschiedenen Donor-Akzeptor-Kombinationen auf ihre optischen, elektronischen und strukturellen Eigenschaften untersucht. Zeitaufgelöste Spektroskopie, vor allem transiente Absorptionsspektroskopie (TA), wird zur Analyse der Ladungsgeneration angewendet und der Vergleich der Donor-PDI Mischfilme mit den Donor-Fulleren Mischfilmen zeigt, dass die Bildung von Ladungstransferzuständen einen der Hauptverlustkanäle darstellt.rnWeiterhin werden Mischungen aus PBDTTT-C und [6,6]-Phenyl-C61-buttersäuremethylesther (PC61BM) mittels TA-Spektroskopie auf einer Zeitskala von ps bis µs untersucht und es kann gezeigt werden, dass der Triplettzustand des Polymers über die nicht-geminale Rekombination freier Ladungen auf einer sub-ns Zeitskala bevölkert wird. Hochentwickelte Methoden zur Datenanalyse, wie multivariate curve resolution (MCR), werden angewendet um überlagernde Datensignale zu trennen. Zusätzlich kann die Regeneration von Ladungsträgern durch Triplett-Triplett-Annihilation auf einer ns-µs Zeitskala gezeigt werden. Darüber hinaus wird der Einfluss des Lösungsmitteladditivs 1,8-Diiodooctan (DIO) auf die Leistungsfähigkeit von p-DTS(FBTTh2)2:PDI Solarzellen untersucht. Die Erkenntnisse von morphologischen und photophysikalischen Experimenten werden kombiniert, um die strukturellen Eigenschaften und die Photophysik mit den relevanten Kenngrößen des Bauteils in Verbindung zu setzen. Zeitaufgelöste Photolumineszenzmessungen (time-resolved photoluminescence, TRPL) zeigen, dass der Einsatz von DIO zu einer geringeren Reduzierung der Photolumineszenz führt, was auf eine größere Phasentrennung zurückgeführt werden kann. Außerdem kann mittels TA Spektroskopie gezeigt werden, dass die Verwendung von DIO zu einer verbesserten Kristallinität der aktiven Schicht führt und die Generation freier Ladungen fördert. Zur genauen Analyse des Signalzerfalls wird ein Modell angewendet, das den gleichzeitigen Zerfall gebundener CT-Zustände und freier Ladungen berücksichtigt und optimierte Donor-Akzeptor-Mischungen zeigen einen größeren Anteil an nicht-geminaler Rekombination freier Ladungsträger.rnIn einer weiteren Fallstudie wird der Einfluss des Fullerenderivats, namentlich IC60BA und PC71BM, auf die Leistungsfähigkeit und Photophysik der Solarzellen untersucht. Eine Kombination aus einer Untersuchung der Struktur des Dünnfilms sowie zeitaufgelöster Spektroskopie ergibt, dass Mischungen, die ICBA als Elektronenakzeptor verwenden, eine schlechtere Trennung von Ladungstransferzuständen zeigen und unter einer stärkeren geminalen Rekombination im Vergleich zu PCBM-basierten Mischungen leiden. Dies kann auf die kleinere Triebkraft zur Ladungstrennung sowie auf die höhere Unordnung der ICBA-basierten Mischungen, die die Ladungstrennung hemmen, zurückgeführt werden. Außerdem wird der Einfluss reiner Fullerendomänen auf die Funktionsfähigkeit organischer Solarzellen, die aus Mischungen des Thienothienophen-basierenden Polymers pBTTT-C14 und PC61BM bestehen, untersucht. Aus diesem Grund wird die Photophysik von Filmen mit einem Donor-Akzeptor-Mischungsverhältnis von 1:1 sowie 1:4 verglichen. Während 1:1-Mischungen lediglich eine co-kristalline Phase, in der Fullerene zwischen den Seitenketten von pBTTT interkalieren, zeigen, resultiert der Überschuss an Fulleren in den 1:4-Proben in der Ausbildung reiner Fullerendomänen zusätzlich zu der co kristallinen Phase. Transiente Absorptionsspektroskopie verdeutlicht, dass Ladungstransferzustände in 1:1-Mischungen hauptsächlich über geminale Rekombination zerfallen, während in 1:4 Mischungen ein beträchtlicher Anteil an Ladungen ihre wechselseitige Coulombanziehung überwinden und freie Ladungsträger bilden kann, die schließlich nicht-geminal rekombinieren.

Computer simulations of two-dimensional colloidal crystals in confinement

Monte Carlo simulations are used to study the effect of confinement on a crystal of point particles interacting with an inverse power law potential in d=2 dimensions. This system can describe colloidal particles at the air-water interface, a model system for experimental study of two-dimensional melting. It is shown that the state of the system (a strip of width D) depends very sensitively on the precise boundary conditions at the two ``walls'' providing the confinement. If one uses a corrugated boundary commensurate with the order of the bulk triangular crystalline structure, both orientational order and positional order is enhanced, and such surface-induced order persists near the boundaries also at temperatures where the system in the bulk is in its fluid state. However, using smooth repulsive boundaries as walls providing the confinement, only the orientational order is enhanced, but positional (quasi-) long range order is destroyed: The mean-square displacement of two particles n lattice parameters apart in the y-direction along the walls then crosses over from the logarithmic increase (characteristic for $d=2$) to a linear increase (characteristic for d=1). The strip then exhibits a vanishing shear modulus. These results are interpreted in terms of a phenomenological harmonic theory. Also the effect of incommensurability of the strip width D with the triangular lattice structure is discussed, and a comparison with surface effects on phase transitions in simple Ising- and XY-models is made

Theoretical studies of fluid membrane mechanics

Biologische Membranen sind Fettmolekül-Doppelschichten, die sich wie zweidimensionale Flüssigkeiten verhalten. Die Energie einer solchen fluiden Oberfläche kann häufig mit Hilfe eines Hamiltonians beschrieben werden, der invariant unter Reparametrisierungen der Oberfläche ist und nur von ihrer Geometrie abhängt. Beiträge innerer Freiheitsgrade und der Umgebung können in den Formalismus mit einbezogen werden. Dieser Ansatz wird in der vorliegenden Arbeit dazu verwendet, die Mechanik fluider Membranen und ähnlicher Oberflächen zu untersuchen. Spannungen und Drehmomente in der Oberfläche lassen sich durch kovariante Tensoren ausdrücken. Diese können dann z. B. dazu verwendet werden, die Gleichgewichtsposition der Kontaktlinie zu bestimmen, an der sich zwei aneinander haftende Oberflächen voneinander trennen. Mit Ausnahme von Kapillarphänomenen ist die Oberflächenenergie nicht nur abhängig von Translationen der Kontaktlinie, sondern auch von Änderungen in der Steigung oder sogar Krümmung. Die sich ergebenden Randbedingungen entsprechen den Gleichgewichtsbedingungen an Kräfte und Drehmomente, falls sich die Kontaktlinie frei bewegen kann. Wenn eine der Oberflächen starr ist, muss die Variation lokal dieser Fläche folgen. Spannungen und Drehmomente tragen dann zu einer einzigen Gleichgewichtsbedingung bei; ihre Beiträge können nicht mehr einzeln identifiziert werden. Um quantitative Aussagen über das Verhalten einer fluiden Oberfläche zu machen, müssen ihre elastischen Eigenschaften bekannt sein. Der "Nanotrommel"-Versuchsaufbau ermöglicht es, Membraneigenschaften lokal zu untersuchen: Er besteht aus einer porenüberspannenden Membran, die während des Experiments durch die Spitze eines Rasterkraftmikroskops in die Pore gedrückt wird. Der lineare Verlauf der resultierenden Kraft-Abstands-Kurven kann mit Hilfe der in dieser Arbeit entwickelten Theorie reproduziert werden, wenn der Einfluss von Adhäsion zwischen Spitze und Membran vernachlässigt wird. Bezieht man diesen Effekt in die Rechnungen mit ein, ändert sich das Resultat erheblich: Kraft-Abstands-Kurven sind nicht länger linear, Hysterese und nichtverschwindende Trennkräfte treten auf. Die Voraussagen der Rechnungen könnten in zukünftigen Experimenten dazu verwendet werden, Parameter wie die Biegesteifigkeit der Membran mit einer Auflösung im Nanometerbereich zu bestimmen. Wenn die Materialeigenschaften bekannt sind, können Probleme der Membranmechanik genauer betrachtet werden. Oberflächenvermittelte Wechselwirkungen sind in diesem Zusammenhang ein interessantes Beispiel. Mit Hilfe des oben erwähnten Spannungstensors können analytische Ausdrücke für die krümmungsvermittelte Kraft zwischen zwei Teilchen, die z. B. Proteine repräsentieren, hergeleitet werden. Zusätzlich wird das Gleichgewicht der Kräfte und Drehmomente genutzt, um mehrere Bedingungen an die Geometrie der Membran abzuleiten. Für den Fall zweier unendlich langer Zylinder auf der Membran werden diese Bedingungen zusammen mit Profilberechnungen kombiniert, um quantitative Aussagen über die Wechselwirkung zu treffen. Theorie und Experiment stoßen an ihre Grenzen, wenn es darum geht, die Relevanz von krümmungsvermittelten Wechselwirkungen in der biologischen Zelle korrekt zu beurteilen. In einem solchen Fall bieten Computersimulationen einen alternativen Ansatz: Die hier präsentierten Simulationen sagen voraus, dass Proteine zusammenfinden und Membranbläschen (Vesikel) bilden können, sobald jedes der Proteine eine Mindestkrümmung in der Membran induziert. Der Radius der Vesikel hängt dabei stark von der lokal aufgeprägten Krümmung ab. Das Resultat der Simulationen wird in dieser Arbeit durch ein approximatives theoretisches Modell qualitativ bestätigt.

Dynamics, rheology and critical properties of colloidal fluid mixtures

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

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.