Analysis of nonlinear diffusion equations of second and fourth order

Wegen der fortschreitenden Miniaturisierung von Halbleiterbauteilen spielen Quanteneffekte eine immer wichtigere Rolle. Quantenphänomene werden gewöhnlich durch kinetische Gleichungen beschrieben, aber manchmal hat eine fluid-dynamische Beschreibung Vorteile: die bessere Nutzbarkeit für numerische Simulationen und die einfachere Vorgabe von Randbedingungen. In dieser Arbeit werden drei Diffusionsgleichungen zweiter und vierter Ordnung untersucht. Der erste Teil behandelt die implizite Zeitdiskretisierung und das Langzeitverhalten einer degenerierten Fokker-Planck-Gleichung. Der zweite Teil der Arbeit besteht aus der Untersuchung des viskosen Quantenhydrodynamischen Modells in einer Raumdimension und dessen Langzeitverhaltens. Im letzten Teil wird die Existenz von Lösungen einer parabolischen Gleichung vierter Ordnung in einer Raumdimension bewiesen, und deren Langzeitverhalten studiert.

Strukturfaktoren, Modenkopplungsgleichungen und Glasübergang molekularer Kristalle

In dieser Arbeit wird der Orientierungsglasübergang ungeordneter, molekularer Kristalle untersucht. Die theoretische Behandlung ist durch die Anisotropie der Einteilchen-Verteilungsfunktion und der Paarfunktionen erschwert. Nimmt man ein starres Gitter, wird der reziproke Raum im Gegenzug auf die 1. Brillouin-Zone eingeschränkt. Der Orientierungsglasübergang wird im Rahmen der Modenkopplungsgleichungen studiert, die dazu hergeleitet werden. Als Modell dienen harte Rotationsellipsoide auf einem starren sc Gitter. Zur Berechnung der statischen tensoriellen Strukturfaktoren wird die Ornstein-Zernike(OZ)-Gleichung molekularer Kristalle abgeleitet und selbstkonsistent zusammen mit der von molekularen Flüssigkeiten übernommenen Percus-Yevick(PY)-Näherung gelöst. Parallel dazu werden die Strukturfaktoren durch MC-Simulationen ermittelt. Die OZ-Gleichung molekularer Kristalle ähnelt der von Flüssigkeiten, direkte und totale Korrelationsfunktion kommen jedoch wegen des starren Gitters nur ohne Konstantanteile in den Winkelvariablen vor, im Gegensatz zur PY-Näherung. Die Anisotropie bringt außerdem einen nichttrivialen Zusatzfaktor. OZ/PY-Strukturfaktoren und MC-Ergebnisse stimmen gut überein. Bei den Matrixelementen der Dichte-Dichte-Korrelationsfunktion gibt es drei Hauptverläufe: oszillatorisch, monoton und unregelmäßig abfallend. Oszillationen gehören zu alternierenden Dichtefluktuationen, führen zu Maxima der Strukturfaktoren am Zonenrand und kommen bei oblaten und genügend breiten prolaten, schwächer auch bei dünnen, nicht zu langen prolaten Ellipsoiden vor. Der exponentielle monotone Abfall kommt bei allen Ellipsoiden vor und führt zu Maxima der Strukturfaktoren in der Zonenmitte, was die Tendenz zu nematischer Ordnung zeigt. Die OZ/PY-Theorie ist durch divergierende Maxima der Strukturfaktoren begrenzt. Bei den Modenkopplungsgleichungen molekularer Kristalle zeigt sich eine große Ähnlichkeit mit denen molekularer Flüssigkeiten, jedoch spielen auf einem starrem Gitter nur die Matrixelemente mit l,l' > 0 eine Rolle und es finden Umklapps von reziproken Vektoren statt. Die Anisotropie bringt auch hier nichtkonstante Zusatzfaktoren ins Spiel. Bis auf flache oblate Ellipsoide wird die Modenkopplungs-Glaslinie von der Divergenz der Strukturfaktoren bestimmt. Für sehr lange Ellipsoide müssen die Strukturfaktoren zur Divergenz hin extrapoliert werden. Daher treibt nicht der Orientierungskäfigeffekt den Glasübergang, sondern Fluktuationen an einer Phasengrenze. Nahe der Kugelform ist keine zuverlässige Glasline festlegbar. Die eingefrorenen kritischen Dichte-Dichte-Korrelatoren haben nur in wenigen Fällen die Oszillationen der statischen Korrelatoren. Der monotone Abfall bleibt dagegen für lange Zeiten meist erhalten. Folglich haben die kritischen Modenkopplungs-Nichtergodizitätsparameter abgeschwächte Maxima in der Zonenmitte, während die Maxima am Zonenrand meist verschwunden sind. Die normierten Nichtergodizitätsparameter zeigen eine Fülle von Verläufen, besonders tiefer im Glas.

Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Non-parametric estimation of the diffusion coefficient of a branching diffusion with immigration

Wir betrachten Systeme von endlich vielen Partikeln, wobei die Partikel sich unabhängig voneinander gemäß eindimensionaler Diffusionen [dX_t = b(X_t),dt + sigma(X_t),dW_t] bewegen. Die Partikel sterben mit positionsabhängigen Raten und hinterlassen eine zufällige Anzahl an Nachkommen, die sich gemäß eines Übergangskerns im Raum verteilen. Zudem immigrieren neue Partikel mit einer konstanten Rate. Ein Prozess mit diesen Eigenschaften wird Verzweigungsprozess mit Immigration genannt. Beobachten wir einen solchen Prozess zu diskreten Zeitpunkten, so ist zunächst nicht offensichtlich, welche diskret beobachteten Punkte zu welchem Pfad gehören. Daher entwickeln wir einen Algorithmus, um den zugrundeliegenden Pfad zu rekonstruieren. Mit Hilfe dieses Algorithmus konstruieren wir einen nichtparametrischen Schätzer für den quadrierten Diffusionskoeffizienten $sigma^2(cdot),$ wobei die Konstruktion im Wesentlichen auf dem Auffüllen eines klassischen Regressionsschemas beruht. Wir beweisen Konsistenz und einen zentralen Grenzwertsatz.

Frobenius polynomials for Calabi-Yau equations

Sei $\pi:X\rightarrow S$ eine \"uber $\Z$ definierte Familie von Calabi-Yau Varietaten der Dimension drei. Es existiere ein unter dem Gauss-Manin Zusammenhang invarianter Untermodul $M\subset H^3_{DR}(X/S)$ von Rang vier, sodass der Picard-Fuchs Operator $P$ auf $M$ ein sogenannter {\em Calabi-Yau } Operator von Ordnung vier ist. Sei $k$ ein endlicher K\"orper der Charaktetristik $p$, und sei $\pi_0:X_0\rightarrow S_0$ die Reduktion von $\pi$ \uber $k$. F\ur die gew\ohnlichen (ordinary) Fasern $X_{t_0}$ der Familie leiten wir eine explizite Formel zur Berechnung des charakteristischen Polynoms des Frobeniusendomorphismus, des {\em Frobeniuspolynoms}, auf dem korrespondierenden Untermodul $M_{cris}\subset H^3_{cris}(X_{t_0})$ her. Sei nun $f_0(z)$ die Potenzreihenl\osung der Differentialgleichung $Pf=0$ in einer Umgebung der Null. Da eine reziproke Nullstelle des Frobeniuspolynoms in einem Teichm\uller-Punkt $t$ durch $f_0(z)/f_0(z^p)|_{z=t}$ gegeben ist, ist ein entscheidender Schritt in der Berechnung des Frobeniuspolynoms die Konstruktion einer $p-$adischen analytischen Fortsetzung des Quotienten $f_0(z)/f_0(z^p)$ auf den Rand des $p-$adischen Einheitskreises. Kann man die Koeffizienten von $f_0$ mithilfe der konstanten Terme in den Potenzen eines Laurent-Polynoms, dessen Newton-Polyeder den Ursprung als einzigen inneren Gitterpunkt enth\alt, ausdr\ucken,so beweisen wir gewisse Kongruenz-Eigenschaften unter den Koeffizienten von $f_0$. Diese sind entscheidend bei der Konstruktion der analytischen Fortsetzung. Enth\alt die Faser $X_{t_0}$ einen gew\ohnlichen Doppelpunkt, so erwarten wir im Grenz\ubergang, dass das Frobeniuspolynom in zwei Faktoren von Grad eins und einen Faktor von Grad zwei zerf\allt. Der Faktor von Grad zwei ist dabei durch einen Koeffizienten $a_p$ eindeutig bestimmt. Durchl\auft nun $p$ die Menge aller Primzahlen, so erwarten wir aufgrund des Modularit\atssatzes, dass es eine Modulform von Gewicht vier gibt, deren Koeffizienten durch die Koeffizienten $a_p$ gegeben sind. Diese Erwartung hat sich durch unsere umfangreichen Rechnungen best\atigt. Dar\uberhinaus leiten wir weitere Formeln zur Bestimmung des Frobeniuspolynoms her, in welchen auch die nicht-holomorphen L\osungen der Gleichung $Pf=0$ in einer Umgebung der Null eine Rolle spielen.

Quadraturformeln für das Qualokationsverfahren

Im Mittelpunkt dieser Arbeit steht Beweis der Existenz- und Eindeutigkeit von Quadraturformeln, die für das Qualokationsverfahren geeignet sind. Letzteres ist ein von Sloan, Wendland und Chandler entwickeltes Verfahren zur numerischen Behandlung von Randintegralgleichungen auf glatten Kurven (allgemeiner: periodische Pseudodifferentialgleichungen). Es erreicht die gleichen Konvergenzordnungen wie das Petrov-Galerkin-Verfahren, wenn man durch den Operator bestimmte Quadraturformeln verwendet. Zunächst werden die hier behandelten Pseudodifferentialoperatoren und das Qualokationsverfahren vorgestellt. Anschließend wird eine Theorie zur Existenz und Eindeutigkeit von Quadraturformeln entwickelt. Ein wesentliches Hilfsmittel hierzu ist die hier bewiesene Verallgemeinerung eines Satzes von Nürnberger über die Existenz und Eindeutigkeit von Quadraturformeln mit positiven Gewichten, die exakt für Tschebyscheff-Räume sind. Es wird schließlich gezeigt, dass es stets eindeutig bestimmte Quadraturformeln gibt, welche die in den Arbeiten von Sloan und Wendland formulierten Bedingungen erfüllen. Desweiteren werden 2-Punkt-Quadraturformeln für so genannte einfache Operatoren bestimmt, mit welchen das Qualokationsverfahren mit einem Testraum von stückweise konstanten Funktionen eine höhere Konvergenzordnung hat. Außerdem wird gezeigt, dass es für nicht-einfache Operatoren im Allgemeinen keine Quadraturformel gibt, mit der die Konvergenzordnung höher als beim Petrov-Galerkin-Verfahren ist. Das letzte Kapitel beinhaltet schließlich numerische Tests mit Operatoren mit konstanten und variablen Koeffizienten, welche die theoretischen Ergebnisse der vorangehenden Kapitel bestätigen.

Partial reconstruction of the trajectories of a discretely observed branching diffusion with immigration and an application to inference

In this treatise we consider finite systems of branching particles where the particles move independently of each other according to d-dimensional diffusions. Particles are killed at a position dependent rate, leaving at their death position a random number of descendants according to a position dependent reproduction law. In addition particles immigrate at constant rate (one immigrant per immigration time). A process with above properties is called a branching diffusion withimmigration (BDI). In the first part we present the model in detail and discuss the properties of the BDI under our basic assumptions. In the second part we consider the problem of reconstruction of the trajectory of a BDI from discrete observations. We observe positions of the particles at discrete times; in particular we assume that we have no information about the pedigree of the particles. A natural question arises if we want to apply statistical procedures on the discrete observations: How can we find couples of particle positions which belong to the same particle? We give an easy to implement 'reconstruction scheme' which allows us to redraw or 'reconstruct' parts of the trajectory of the BDI with high accuracy. Moreover asymptotically the whole path can be reconstructed. Further we present simulations which show that our partial reconstruction rule is tractable in practice. In the third part we study how the partial reconstruction rule fits into statistical applications. As an extensive example we present a nonparametric estimator for the diffusion coefficient of a BDI where the particles move according to one-dimensional diffusions. This estimator is based on the Nadaraya-Watson estimator for the diffusion coefficient of one-dimensional diffusions and it uses the partial reconstruction rule developed in the second part above. We are able to prove a rate of convergence of this estimator and finally we present simulations which show that the estimator works well even if we leave our set of assumptions.

Theoretical investigation of the interaction of a polymer film with a nanoparticle

The fundamental aim in our investigation of the interaction of a polymer film with a nanoparticle is the extraction of information on the dynamics of the liquid using a single tracking particle. In this work two theoretical methods were used: one passive, where the motion of the particle measures the dynamics of the liquid, one active, where perturbations in the system are introduced through the particle. In the first part of this investigation a thin polymeric film on a substrate is studied using molecular dynamics simulations. The polymer is modeled via a 'bead spring' model. The particle is spheric and non structured and is able to interact with the monomers via a Lennard Jones potential. The system is micro-canonical and simulations were performed for average temperatures between the glass transition temperature of the film and its dewetting temperature. It is shown that the stability of the nanoparticle on the polymer film in the absence of gravity depends strongly on the form of the chosen interaction potential between nanoparticle and polymer. The relative position of the tracking particle to the liquid vapor interface of the polymer film shows the glass transition of the latter. The velocity correlation function and the mean square displacement of the particle has shown that it is caged when the temperature is close to the glass transition temperature. The analysis of the dynamics at long times shows the coupling of the nanoparticle to the center of mass of the polymer chains. The use of the Stokes-Einstein formula, which relates the diffusion coefficient to the viscosity, permits to use the nanoparticle as a probe for the determination of the bulk viscosity of the melt, the so called 'microrheology'. It is shown that for low frequencies the result obtained using microrheology coincides with the results of the Rouse model applied to the polymer dynamics. In the second part of this investigation the equations of Linear Hydrodynamics are solved for a nanoparticle oscillating above the film. It is shown that compressible liquids have mechanical response to external perturbations induced with the nanoparticle. These solutions show strong velocity and pressure profiles of the liquid near the interface, as well as a mechanical response of the liquid-vapor interface. The results obtained with this calculations can be employed for the interpretation of experimental results of non contact AFM microscopy

Structuring of polymer surface by evaporation of sessile microdrops

In this thesis, we investigated the evaporation of sessile microdroplets on different solid substrates. Three major aspects were studied: the influence of surface hydrophilicity and heterogeneity on the evaporation dynamics for an insoluble solid substrate, the influence of external process parameters and intrinsic material properties on microstructuring of soluble polymer substrates and the influence of an increased area to volume ratio in a microfluidic capillary, when evaporation is hindered. In the first part, the evaporation dynamics of pure sessile water drops on smooth self-assembled monolayers (SAMs) of thiols or disulfides on gold on mica was studied. With increasing surface hydrophilicity the drop stayed pinned longer. Thus, the total evaporation time of a given initial drop volume was shorter, since the drop surface, through which the evaporation occurs, stays longer large. Usually, for a single drop the volume decreased linearly with t1.5, t being the evaporation time, for a diffusion-controlled evaporation process. However, when we measured the total evaporation time, ttot, for multiple droplets with different initial volumes, V0, we found a scaling of the form V0 = attotb. The more hydrophilic the substrate was, the more showed the scaling exponent a tendency to an increased value up to 1.6. This can be attributed to an increasing evaporation rate through a thin water layer in the vicinity of the drop. Under the assumption of a constant temperature at the substrate surface a cooling of the droplet and thus a decreased evaporation rate could be excluded as a reason for the different scaling exponent by simulations performed by F. Schönfeld at the IMM, Mainz. In contrast, for a hairy surface, made of dialkyldisulfide SAMs with different chain lengths and a 1:1 mixture of hydrophilic and hydrophobic end groups (hydroxy versus methyl group), the scaling exponent was found to be ~ 1.4. It increased to ~ 1.5 with increasing hydrophilicity. A reason for this observation can only be speculated: in the case of longer hydrophobic alkyl chains the formation of an air layer between substrate and surface might be favorable. Thus, the heat transport to the substrate might be reduced, leading to a stronger cooling and thus decreased evaporation rate. In the second part, the microstructuring of polystyrene surfaces by drops of toluene, a good solvent, was investigated. For this a novel deposition technique was developed, with which the drop can be deposited with a syringe. The polymer substrate is lying on a motorized table, which picks up the pendant drop by an upward motion until a liquid bridge is formed. A consecutive downward motion of the table after a variable delay, i.e. the contact time between drop and polymer, leads to the deposition of the droplet, which can evaporate. The resulting microstructure is investigated in dependence of the processes parameters, i.e. the approach and the retraction speed of the substrate and the delay between them, and in dependence of the intrinsic material properties, i.e. the molar mass and the type of the polymer/solvent system. The principal equivalence with the microstructuring by the ink-jet technique was demonstrated. For a high approach and retraction speed of 9 mm/s and no delay between them, a concave microtopology was observed. In agreement with the literature, this can be explained by a flow of solvent and the dissolved polymer to the rim of the pinned droplet, where polymer is accumulated. This effect is analogue to the well-known formation of ring-like stains after the evaporation of coffee drops (coffee-stain effect). With decreasing retraction speed down to 10 µm/s the resulting surface topology changes from concave to convex. This can be explained with the increasing dissolution of polymer into the solvent drop prior to the evaporation. If the polymer concentration is high enough, gelation occurs instead of a flow to the rim and the shape of the convex droplet is received. With increasing delay time from below 0 ms to 1s the depth of the concave microwells decreases from 4.6 µm to 3.2 µm. However, a convex surface topology could not be obtained, since for longer delay times the polymer sticks to the tip of the syringe. Thus, by changing the delay time a fine-tuning of the concave structure is accomplished, while by changing the retraction speed a principal change of the microtopolgy can be achieved. We attribute this to an additional flow inside the liquid bridge, which enhanced polymer dissolution. Even if the pendant drop is evaporating about 30 µm above the polymer surface without any contact (non-contact mode), concave structures were observed. Rim heights as high as 33 µm could be generated for exposure times of 20 min. The concave structure exclusively lay above the flat polymer surface outside the structure even after drying. This shows that toluene is taken up permanently. The increasing rim height, rh, with increasing exposure time to the solvent vapor obeys a diffusion law of rh = rh0  tn, with n in the range of 0.46 ~ 0.65. This hints at a non-Fickian swelling process. A detailed analysis showed that the rim height of the concave structure is modulated, unlike for the drop deposition. This is due to the local stress relaxation, which was initiated by the increasing toluene concentration in the extruded polymer surface. By altering the intrinsic material parameters i.e. the polymer molar mass and the polymer/solvent combination, several types of microstructures could be formed. With increasing molar mass from 20.9 kDa to 1.44 MDa the resulting microstructure changed from convex, to a structure with a dimple in the center, to concave, to finally an irregular structure. This observation can be explained if one assumes that the microstructuring is dominated by two opposing effects, a decreasing solubility with increasing polymer molar mass, but an increasing surface tension gradient leading to instabilities of Marangoni-type. Thus, a polymer with a low molar mass close or below the entanglement limit is subject to a high dissolution rate, which leads to fast gelation compared to the evaporation rate. This way a coffee-rim like effect is eliminated early and a convex structure results. For high molar masses the low dissolution rate and the low polymer diffusion might lead to increased surface tension gradients and a typical local pile-up of polymer is found. For intermediate polymer masses around 200 kDa, the dissolution and evaporation rate are comparable and the typical concave microtopology is found. This interpretation was supported by a quantitative estimation of the diffusion coefficient and the evaporation rate. For a different polymer/solvent system, polyethylmethacrylate (PEMA)/ethylacetate (EA), exclusively concave structures were found. Following the statements above this can be interpreted with a lower dissolution rate. At low molar masses the concentration of PEMA in EA most likely never reaches the gelation point. Thus, a concave instead of a convex structure occurs. At the end of this section, the optically properties of such microstructures for a potential application as microlenses are studied with laser scanning confocal microscopy. In the third part, the droplet was confined into a glass microcapillary to avoid evaporation. Since here, due to an increased area to volume ratio, the surface properties of the liquid and the solid walls became important, the influence of the surface hydrophilicity of the wall on the interfacial tension between two immiscible liquid slugs was investigated. For this a novel method for measuring the interfacial tension between the two liquids within the capillary was developed. This technique was demonstrated by measuring the interfacial tensions between slugs of pure water and standard solvents. For toluene, n-hexane and chloroform 36.2, 50.9 and 34.2 mN/m were measured at 20°C, which is in a good agreement with data from the literature. For a slug of hexane in contact with a slug of pure water containing ethanol in a concentration range between 0 and 70 (v/v %), a difference of up to 6 mN/m was found, when compared to commercial ring tensiometry. This discrepancy is still under debate.

Influence of spatial constraints in non-crystalline regions on the polymer dynamics in semi-crystalline polyethylene: a solid state NMR study

A broad variety of solid state NMR techniques were used to investigate the chain dynamics in several polyethylene (PE) samples, including ultrahigh molecular weight PEs (UHMW-PEs) and low molecular weight PEs (LMW-PEs). Via changing the processing history, i.e. melt/solution crystallization and drawing processes, these samples gain different morphologies, leading to different molecular dynamics. Due to the long chain nature, the molecular dynamics of polyethylene can be distinguished in local fluctuation and long range motion. With the help of NMR these different kinds of molecular dynamics can be monitored separately. In this work the local chain dynamics in non-crystalline regions of polyethylene samples was investigated via measuring 1H-13C heteronuclear dipolar coupling and 13C chemical shift anisotropy (CSA). By analyzing the motionally averaged 1H-13C heteronuclear dipolar coupling and 13C CSA, the information about the local anisotropy and geometry of motion was obtained. Taking advantage of the big difference of the 13C T1 relaxation time in crystalline and non-crystalline regions of PEs, the 1D 13C MAS exchange experiment was used to investigate the cooperative chain motion between these regions. The different chain organizations in non-crystalline regions were used to explain the relationship between the local fluctuation and the long range motion of the samples. In a simple manner the cooperative chain motion between crystalline and non-crystalline regions of PE results in the experimentally observed diffusive behavior of PE chain. The morphological influences on the diffusion motion have been discussed. The morphological factors include lamellar thickness, chain organization in non-crystalline regions and chain entanglements. Thermodynamics of the diffusion motion in melt and solution crystallized UHMW-PEs is discussed, revealing entropy-controlled features of the chain diffusion in PE. This thermodynamic consideration explains the counterintuitive relationship between the local fluctuation and the long range motion of the samples. Using the chain diffusion coefficient, the rates of jump motion in crystals of the melt crystallized PE have been calculated. A concept of "effective" jump motion has been proposed to explain the difference between the values derived from the chain diffusion coefficients and those in literatures. The observations of this thesis give a clear demonstration of the strong relationship between the sample morphology and chain dynamics. The sample morphologies governed by the processing history lead to different spatial constraints for the molecular chains, leading to different features of the local and long range chain dynamics. The knowledge of the morphological influence on the microscopic chain motion has many implications in our understanding of the alpha-relaxation process in PE and the related phenomena such as crystal thickening, drawability of PE, the easy creep of PE fiber, etc.

Diffusion of laser polarized gases in MRI

The use of Magnetic Resonance Imaging (MRI) as a diagnostic tool is increasingly employing functional contrast agents to study or contrast entire mechanisms. Contrast agents in MRI can be classified in two categories. One type of contrast agents alters the NMR signal of the protons in its surrounding, e.g. lowers the T1 relaxation time. The other type enhances the Nuclear Magnetic Resonance (NMR) signal of specific nuclei. For hyperpolarized gases the NMR signal is improved up to several orders of magnitude. However, gases have a high diffusivity which strongly influences the NMR signal strength, hence the resolution and appearance of the images. The most interesting question in spatially resolved experiments is of course the achievable resolution and contrast by controlling the diffusivity of the gas. The influence of such diffusive processes scales with the diffusion coefficient, the strength of the magnetic field gradients and the timings used in the experiment. Diffusion may not only limit the MRI resolution, but also distort the line shape of MR images for samples, which contain boundaries or diffusion barriers within the sampled space. In addition, due to the large polarization in gaseous 3He and 129Xe, spin diffusion (different from particle diffusion) could play a role in MRI experiments. It is demonstrated that for low temperatures some corrections to the NMR measured diffusion coefficient have to be done, which depend on quantum exchange effects for indistinguishable particles. Physically, if these effects can not change the spin current, they can do it indirectly by modifying the velocity distribution of the different spin states separately, so that the subsequent collisions between atoms and therefore the diffusion coefficient can eventually be affected. A detailed study of the hyperpolarized gas diffusion coefficient is presented, demonstrating the absence of spin diffusion (different from particle diffusion) influence in MRI at clinical conditions. A novel procedure is proposed to control the diffusion coefficient of gases in MRI by admixture of inert buffer gases. The experimental measured diffusion agrees with theoretical simulations. Therefore, the molecular mass and concentration enter as additional parameters into the equations that describe structural contrast. This allows for setting a structural threshold up to which structures contribute to the image. For MRI of the lung this allows for images of very small structural elements (alveoli) only, or in the other extreme, all airways can be displayed with minimal signal loss due to diffusion.

Pore structure and column efficiency of silica monoliths in high performance liquid chromatography (HPLC)

Five different methods were critically examined to characterize the pore structure of the silica monoliths. The mesopore characterization was performed using: a) the classical BJH method of nitrogen sorption data, which showed overestimated values in the mesopore distribution and was improved by using the NLDFT method, b) the ISEC method implementing the PPM and PNM models, which were especially developed for monolithic silicas, that contrary to the particulate supports, demonstrate the two inflection points in the ISEC curve, enabling the calculation of pore connectivity, a measure for the mass transfer kinetics in the mesopore network, c) the mercury porosimetry using a new recommended mercury contact angle values. rnThe results of the characterization of mesopores of monolithic silica columns by the three methods indicated that all methods were useful with respect to the pore size distribution by volume, but only the ISEC method with implemented PPM and PNM models gave the average pore size and distribution based on the number average and the pore connectivity values.rnThe characterization of the flow-through pore was performed by two different methods: a) the mercury porosimetry, which was used not only for average flow-through pore value estimation, but also the assessment of entrapment. It was found that the mass transfer from the flow-through pores to mesopores was not hindered in case of small sized flow-through pores with a narrow distribution, b) the liquid penetration where the average flow-through pore values were obtained via existing equations and improved by the additional methods developed according to Hagen-Poiseuille rules. The result was that not the flow-through pore size influences the column bock pressure, but the surface area to volume ratio of silica skeleton is most decisive. Thus the monolith with lowest ratio values will be the most permeable. rnThe flow-through pore characterization results obtained by mercury porosimetry and liquid permeability were compared with the ones from imaging and image analysis. All named methods enable a reliable characterization of the flow-through pore diameters for the monolithic silica columns, but special care should be taken about the chosen theoretical model.rnThe measured pore characterization parameters were then linked with the mass transfer properties of monolithic silica columns. As indicated by the ISEC results, no restrictions in mass transfer resistance were noticed in mesopores due to their high connectivity. The mercury porosimetry results also gave evidence that no restrictions occur for mass transfer from flow-through pores to mesopores in the small scaled silica monoliths with narrow distribution. rnThe prediction of the optimum regimes of the pore structural parameters for the given target parameters in HPLC separations was performed. It was found that a low mass transfer resistance in the mesopore volume is achieved when the nominal diameter of the number average size distribution of the mesopores is appr. an order of magnitude larger that the molecular radius of the analyte. The effective diffusion coefficient of an analyte molecule in the mesopore volume is strongly dependent on the value of the nominal pore diameter of the number averaged pore size distribution. The mesopore size has to be adapted to the molecular size of the analyte, in particular for peptides and proteins. rnThe study on flow-through pores of silica monoliths demonstrated that the surface to volume of the skeletons ratio and external porosity are decisive for the column efficiency. The latter is independent from the flow-through pore diameter. The flow-through pore characteristics by direct and indirect approaches were assessed and theoretical column efficiency curves were derived. The study showed that next to the surface to volume ratio, the total porosity and its distribution of the flow-through pores and mesopores have a substantial effect on the column plate number, especially as the extent of adsorption increases. The column efficiency is increasing with decreasing flow through pore diameter, decreasing with external porosity, and increasing with total porosity. Though this tendency has a limit due to heterogeneity of the studied monolithic samples. We found that the maximum efficiency of the studied monolithic research columns could be reached at a skeleton diameter of ~ 0.5 µm. Furthermore when the intention is to maximize the column efficiency, more homogeneous monoliths should be prepared.rn

Numerical studies of colloidal systems

The interplay of hydrodynamic and electrostatic forces is of great importance for the understanding of colloidal dispersions. Theoretical descriptions are often based on the so called standard electrokinetic model. This Mean Field approach combines the Stokes equation for the hydrodynamic flow field, the Poisson equation for electrostatics and a continuity equation describing the evolution of the ion concentration fields. In the first part of this thesis a new lattice method is presented in order to efficiently solve the set of non-linear equations for a charge-stabilized colloidal dispersion in the presence of an external electric field. Within this framework, the research is mainly focused on the calculation of the electrophoretic mobility. Since this transport coefficient is independent of the electric field only for small driving, the algorithm is based upon a linearization of the governing equations. The zeroth order is the well known Poisson-Boltzmann theory and the first order is a coupled set of linear equations. Furthermore, this set of equations is divided into several subproblems. A specialized solver for each subproblem is developed, and various tests and applications are discussed for every particular method. Finally, all solvers are combined in an iterative procedure and applied to several interesting questions, for example, the effect of the screening mechanism on the electrophoretic mobility or the charge dependence of the field-induced dipole moment and ion clouds surrounding a weakly charged sphere. In the second part a quantitative data analysis method is developed for a new experimental approach, known as "Total Internal Reflection Fluorescence Cross-Correlation Spectroscopy" (TIR-FCCS). The TIR-FCCS setup is an optical method using fluorescent colloidal particles to analyze the flow field close to a solid-fluid interface. The interpretation of the experimental results requires a theoretical model, which is usually the solution of a convection-diffusion equation. Since an analytic solution is not available due to the form of the flow field and the boundary conditions, an alternative numerical approach is presented. It is based on stochastic methods, i. e. a combination of a Brownian Dynamics algorithm and Monte Carlo techniques. Finally, experimental measurements for a hydrophilic surface are analyzed using this new numerical approach.

Black-scholes type equations

In this work we are concerned with the analysis and numerical solution of Black-Scholes type equations arising in the modeling of incomplete financial markets and an inverse problem of determining the local volatility function in a generalized Black-Scholes model from observed option prices. In the first chapter a fully nonlinear Black-Scholes equation which models transaction costs arising in option pricing is discretized by a new high order compact scheme. The compact scheme is proved to be unconditionally stable and non-oscillatory and is very efficient compared to classical schemes. Moreover, it is shown that the finite difference solution converges locally uniformly to the unique viscosity solution of the continuous equation. In the next chapter we turn to the calibration problem of computing local volatility functions from market data in a generalized Black-Scholes setting. We follow an optimal control approach in a Lagrangian framework. We show the existence of a global solution and study first- and second-order optimality conditions. Furthermore, we propose an algorithm that is based on a globalized sequential quadratic programming method and a primal-dual active set strategy, and present numerical results. In the last chapter we consider a quasilinear parabolic equation with quadratic gradient terms, which arises in the modeling of an optimal portfolio in incomplete markets. The existence of weak solutions is shown by considering a sequence of approximate solutions. The main difficulty of the proof is to infer the strong convergence of the sequence. Furthermore, we prove the uniqueness of weak solutions under a smallness condition on the derivatives of the covariance matrices with respect to the solution, but without additional regularity assumptions on the solution. The results are illustrated by a numerical example.

Thermodynamics of aqueous biopolymer solutions and fractionation of Dextran

Flory-Huggins interaction parameters and thermal diffusion coefficients were measured for aqueous biopolymer solutions. Dextran (a water soluble polysaccharide) and bovine serum albumin (BSA, a water soluble protein) were used for this study. The former polymer is representative for chain macromolecules and the latter is for globular macromolecules. The interaction parameters for the systems water/dextran and water/BSA were determined as a function of composition by means of vapor pressure measurements, using a combination of headspace sampling and gas chromatography (HS-GC). A new theoretical approach, accounting for chain connectivity and conformational variability, describes the observed dependencies quantitatively for the system water/dextran and qualitatively for the system water/BSA. The phase diagrams of the ternary systems water/methanol/dextran and water/dextran/BSA were determined via cloud point measurements and modeled by means of the direct minimization of the Gibbs energy using the information on the binary subsystems as input parameters. The thermal diffusion of dextran was studied for aqueous solutions in the temperature range 15 < T < 55 oC. The effects of the addition of urea were also studied. In the absence of urea, the Soret coefficient ST changes its sign as T is varied; it is positive for T > 45.0 oC, but negative for T < 45.0 oC. The positive sign of ST means that the dextran molecules migrate towards the cold side of the fluid; this behavior is typical for polymer solutions. While a negative sign indicates the macromolecules move toward the hot side; this behavior has so far not been observed with any other binary aqueous polymer solutions. The addition of urea to the aqueous solution of dextran increases ST and reduces the inversion temperature. For 2 M urea, the change in the sign of ST is observed at T = 29.7 oC. At higher temperature ST is always positive in the studied temperature range. To rationalize these observations it is assumed that the addition of urea opens hydrogen bonds, similar to that induced by an increase in temperature. For a future extension of the thermodynamic studies to the effects of poly-dispersity, dextran was fractionated by means of a recently developed technique called Continuous Spin Fractionation (CSF). The solvent/precipitant/polymer system used for the thermodynamic studies served as the basis for the fractionation of dextran The starting polymer had a weight average molar mass Mw = 11.1 kg/mol and a molecular non-uniformity U= Mw / Mn -1= 1.0. Seventy grams of dextran were fractionated using water as the solvent and methanol as the precipitant. Five fractionation steps yielded four samples with Mw values between 4.36 and 18.2 kg/mol and U values ranging from 0.28 to 0.48.