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.

Diffusion in polymer systems studied by fluorescence correlation spectroscopy

Fluorescence correlation spectroscopy (FCS) is a powerful technique to determine the diffusion of fluorescence molecules in various environments. The technique is based on detecting and analyzing the fluctuation of fluorescence light emitted by fluorescence species diffusing through a small and fixed observation volume, formed by a laser focused into the sample. Because of its great potential and high versatility in addressing the diffusion and transport properties in complex systems, FCS has been successfully applied to a great variety of systems. In my thesis, I focused on the application of FCS to study the diffusion of fluorescence molecules in organic environments, especially in polymer melts. In order to examine our FCS setup and a developed measurement protocol, I first utilized FCS to measure tracer diffusion in polystyrene (PS) solutions, for which abundance data exist in the literature. I studied molecular and polymeric tracer diffusion in polystyrene solutions over a broad range of concentrations and different tracer and matrix molecular weights (Mw). Then FCS was further established to study tracer dynamics in polymer melts. In this part I investigated the diffusion of molecular tracers in linear flexible polymer melts [polydimethylsiloxane (PDMS), polyisoprene (PI)], a miscible polymer blend [PI and poly vinyl ethylene (PVE)], and star-shaped polymer [3-arm star polyisoprene (SPI)]. The effects of tracer sizes, polymer Mw, polymer types, and temperature on the diffusion coefficients of small tracers were discussed. The distinct topology of the host polymer, i.e. star polymer melt, revealed the notably different motion of the small tracer, as compared to its linear counterpart. Finally, I emphasized the advantage of the small observation volume which allowed FCS to investigate the tracer diffusions in heterogeneous systems; a swollen cross-linked PS bead and silica inverse opals, where high spatial resolution technique was required.

Diffusion maps in the Subriemannian motion group and perceptual grouping

Il presente lavoro è motivato dal problema della constituzione di unità percettive a livello della corteccia visiva primaria V1. Si studia dettagliatamente il modello geometrico di Citti-Sarti con particolare attenzione alla modellazione di fenomeni di associazione visiva. Viene studiato nel dettaglio un modello di connettività. Il contributo originale risiede nell'adattamento del metodo delle diffusion maps, recentemente introdotto da Coifman e Lafon, alla geometria subriemanniana della corteccia visiva. Vengono utilizzati strumenti di teoria del potenziale, teoria spettrale, analisi armonica in gruppi di Lie per l'approssimazione delle autofunzioni dell'operatore del calore sul gruppo dei moti rigidi del piano. Le autofunzioni sono utilizzate per l'estrazione di unità percettive nello stimolo visivo. Sono presentate prove sperimentali e originali delle capacità performanti del metodo.

Asymmetric drift and halo flattening in disk galaxies

The main result in this work is the solution of the Jeans equations for an axisymmetric galaxy model containing a baryonic component (distributed according to a Miyamoto-Nagai profile) and a dark matter halo (described by the Binney logarithmic potential). The velocity dispersion, azimuthal velocity and some other interesting quantities such as the asymmetric drift are studied, along with the influence of the model parameters on these (observable) quantities. We also give an estimate for the velocity of the radial flow, caused by the asymmetric drift. Other than the mathematical beauty that lies in solving a model analytically, the interest of this kind of results can be mainly found in numerical simulations that study the evolution of gas flows. For example, it is important to know how certain parameters such as the shape (oblate, prolate, spherical) of a dark matter halo, or the flattening of the baryonic matter, or the mass ratio between dark and baryonic matter, have an influence on observable quantities such as the velocity dispersion. In the introductory chapter, we discuss the Jeans equations, which provide information about the velocity dispersion of a system. Next we will consider some dynamical quantities that will be useful in the rest of the work, e.g. the asymmetric drift. In Chapter 2 we discuss in some more detail the family of galaxy models we studied. In Chapter 3 we give the solution of the Jeans equations. Chapter 4 describes and illustrates the behaviour of the velocity dispersion, as a function of the several parameters, along with asymptotic expansions. In Chapter 5 we will investigate the behaviour of certain dynamical quantities for this model. We conclude with a discussion in Chapter 6.

Tracer diffusion at water-oil interfaces studied by fluorescence correlation spectroscopy

The adsorption of particles and surfactants at water-oil interfaces has attracted continuous attention because of its emulsion stabilizing effect and the possibility to form two-dimensional materials. Herein, I studied the interfacial diffusion of single molecules and nanoparticles at water-oil interfaces using fluorescence correlation spectroscopy. rnrnFluorescence correlation spectroscopy (FCS) is a promising technique to study diffusion of fluorescent tracers in diverse conditions. This technique monitors and analyzes the fluorescence fluctuation caused by single fluorescent tracers coming in and out of a diffraction-limited observation volume “one at a time”. Thus, this technique allows a combination of high precision, high spatial resolution and low tracer concentration. rnrnIn chapter 1, I discussed some controversial questions regarding the properties of water-hydrophobic interfaces and also introduced the current progress on the stability and dynamic of single nanoparticles at water-oil interfaces. The materials and setups I used in this thesis were summarized in chapter 2. rnrnIn chapter 3, I presented a new strategy to study the properties of water-oil interfaces. The two-dimensional diffusion of isolated molecular tracers at water/n-alkane interfaces was measured using fluorescence correlation spectroscopy. The diffusion coefficients of larger tracers with a hydrodynamic radius of 4.0 nm agreed well with the values calculated from the macroscopic viscosities of the two bulk phases. However, for small molecule tracers with hydrodynamic radii of only 1.0 and 0.6 nm, notable deviations were observed, indicating the existence of an interfacial region with a reduced effective viscosity. rnrnIn chapter 4, the interfacial diffusion of nanoparticles at water-oil interfaces was investigated using FCS. In stark contrast to the interfacial diffusion of molecular tracers, that of nanoparticles at any conditions is slower than the values calculated in accordance to the surrounding viscosity. The diffusion of nanoparticles at water-oil interfaces depended on the interfacial tension of liquid-liquid interfaces, the surface properties of nanoparticles, the particle sizes and the viscosities of surrounding liquid phases. In addition, the interfacial diffusion of nanoparticles with Janus motif is even slower than that of their symmetric counterparts. Based on the experimental results I obtained, I drew some possibilities to describe the origin of nanoparticle slowdown at water-oil interfaces.

The thickness of the gaseous layer in disk galaxies

il profilo verticale del gas dipende sia dal potenziale galattico, e quindi dalla distribuzione di materia, sia dalla dispersione di velocit`a del gas La nostra analisi teorica ha mostrato che l’altezza scala dell’idrogeno gassoso ad un certo raggio `e fortemente legata alla componente di massa dominate: nelle parti centrali, dove domina la materia barionica, il disco gassoso incrementa il suo spessore esponenzialmente con una lunghezza scala pari al doppio di quella del disco stellare; dalle zone intermedie fino alla fine del disco, il profilo verticale del gas `e influenzato dalla distribuzione di DM. In queste zone lo spessore del disco gassoso cresce con il raggio in maniera circa lineare. Tale fenomeno viene definito gas flaring. Lo scopo principale di questa tesi `e la ricerca di una nuova tecnica per poter stimare il profilo vertical del gas in galassie con inclinazione intermedia.

Preconditioning strategies for the numerical solution of convection-diffusion partial differential equations

Il trattamento numerico dell'equazione di convezione-diffusione con le relative condizioni al bordo, comporta la risoluzione di sistemi lineari algebrici di grandi dimensioni in cui la matrice dei coefficienti è non simmetrica. Risolutori iterativi basati sul sottospazio di Krylov sono ampiamente utilizzati per questi sistemi lineari la cui risoluzione risulta particolarmente impegnativa nel caso di convezione dominante. In questa tesi vengono analizzate alcune strategie di precondizionamento, atte ad accelerare la convergenza di questi metodi iterativi. Vengono confrontati sperimentalmente precondizionatori molto noti come ILU e iterazioni di tipo inner-outer flessibile. Nel caso in cui i coefficienti del termine di convezione siano a variabili separabili, proponiamo una nuova strategia di precondizionamento basata sull'approssimazione, mediante equazione matriciale, dell'operatore differenziale di convezione-diffusione. L'azione di questo nuovo precondizionatore sfrutta in modo opportuno recenti risolutori efficienti per equazioni matriciali lineari. Vengono riportati numerosi esperimenti numerici per studiare la dipendenza della performance dei diversi risolutori dalla scelta del termine di convezione, e dai parametri di discretizzazione.

Speziation von Plutonium bei der Sorption und Diffusion in Opalinuston

Die vorliegende Arbeit wurde im Rahmen des BMWi-Verbundprojektes Wechselwirkung und Transport von Aktiniden im natürlichen Tongestein unter Berücksichtigung von Huminstoffen und Tonorganika – Wechselwirkung von Neptunium und Plutonium mit natürlichem Tongestein“ durchgeführt. Um die langfristige Sicherheit der nuklearen Endlager beurteilen zu können, muss eine mögliche Migration der radiotoxischen Abfälle in die Umwelt betrachtet werden. Wegen seiner langen Halbwertszeit (24000 a) leistet Pu-239 einen wesentlichen Beitrag zur Radiotoxizität abgebrannter Kernbrennstoffe in einem Endlager. Das redox-sensitive Pu tritt in Lösung unter umweltrelevanten Bedingungen in den Oxidationsstufen +III bis +VI auf und kann nebeneinander in bis zu vier Oxidationsstufen vorliegen. Tonsteinformationen werden als mögliches Wirtsgestein für Endlager hoch-radioaktiver Abfälle betrachtet. Deshalb sind ausführliche Informationen zur Mobilisierung und Immobilisierung des Pu durch/in das Grundwasser aus einem Endlager von besonderer Bedeutung. In dieser Arbeit wurden neue Erkenntnisse über die Wechselwirkung zwischen Pu und dem natürlichen Tongestein Opalinuston (OPA, Mont Terri, Schweiz) mit Hinblick auf die Endlagerung wärmeentwickelnder radioaktiver Abfälle in einem geologischen Tiefenlager gewonnen.rnDer Fokus der Arbeit lag dabei auf der Bestimmung der Speziation von Pu an der Mineraloberfläche nach Sorptions- und Diffusionsprozessen mittels verschiedener synchrotronbasierter Methoden (µ-XRF, µ-XANES/EXAFS, µ-XRD, XANES/EXAFS). rnDie Wechselwirkung zwischen Pu und OPA wurde zunächst in Batch- und Diffusionsexperimenten in Abhängigkeit verschiedener experimenteller Parameter (u.a. pH, Pu-Oxidationsstufe) untersucht. In Sorptionsexperimenten konnte gezeigt werden, dass einige Parameter (z.B. Temperatur, Huminsäure) einen deutlichen Einfluss auf die Sorption von Pu haben.rnDie Speziationsuntersuchungen wurden zum einen an Pulverproben aus Batchexperimenten und zum anderen an OPA-Dünnschliffen bzw. Diffusionsproben in Abhängigkeit verschiedener experimenteller Parameter durchgeführt. Die EXAFS-Messungen an der Pu LIII-Kante der Pulverproben ergaben, dass eine innersphäriche Sorption von Pu(IV) an Tongestein unabhängig von dem Ausgangsoxidationszustand des Plutoniums in Lösung stattgefunden hat. Durch die Kombination der ortsaufgelösten Methoden wurde erstmalig mittels μ-XRF die Verteilung von Pu und anderen in OPA enthaltenen Elementen bestimmt. µ-XANES-Spektren an Pu-Anreicherungen auf OPA-Dünnschliffen und in Diffusionsproben bestätigen, dass das weniger mobile Pu(IV) die dominierende Spezies nach den Sorptions- und Diffusionsprozessen ist. Darüber hinaus wurde zum ersten Mal ein Diffusionsprofil von Pu in OPA mittels µ-XRF gemessen. Die Speziationsuntersuchungen mittels μ-XANES zeigten, dass das eingesetzte Pu(V) entlang seines Diffusionspfades zunehmend zu Pu(IV) reduziert wird. Mit µ-XRD wurde Illit als dominierende Umgebung, in der Pu angereichert wurde, identifiziert und Siderit als eine redoxaktive Phase auftreten kann. Die Ergebnisse dieser Arbeit zeigen, dass die Sicherheit von OPA als Wirtsgestein eines Endlagers hoch-radioaktiver Abfälle positiv zu bewerten ist. rn

Diffusion in heterogeneous systems studied by laser scanning confocal microscopy and fluorescence correlation spectroscopy

Understanding and controlling the mechanism of the diffusion of small molecules, macromolecules and nanoparticles in heterogeneous environments is of paramount fundamental and technological importance. The aim of the thesis is to show, how by studying the tracer diffusion in complex systems, one can obtain information about the tracer itself, and the system where the tracer is diffusing. rnIn the first part of my thesis I will introduce the Fluorescence Correlation Spectroscopy (FCS) which is a powerful tool to investigate the diffusion of fluorescent species in various environments. By using the main advantage of FCS namely the very small probing volume (<1µm3) I was able to track the kinetics of phase separation in polymer blends at late stages by looking on the molecular tracer diffusion in individual domains of the heterogeneous structure of the blend. The phase separation process at intermediate stages was monitored with laser scanning confocal microscopy (LSCM) in real time providing images of droplet coalescence and growth. rnIn a further project described in my thesis I will show that even when the length scale of the heterogeneities becomes smaller than the FCS probing volume one can still obtain important microscopic information by studying small tracer diffusion. To do so, I will introduce a system of star shaped polymer solutions and will demonstrate that the mobility of small molecular tracers on microscopic level is nearly not affected by the transition of the polymer system to a “glassy” macroscopic state. rnIn the last part of the thesis I will introduce and describe a new stimuli responsive system which I have developed, that combines two levels of nanoporosity. The system is based on poly-N-isopropylacrylamide (PNIPAM) and silica inverse opals (iOpals), and allows controlling the diffusion of tracer molecules. rn

Orthogonal gamma-based expansions for volatility option prices under jump-diffusion dynamics

In my work I derive closed-form pricing formulas for volatility based options by suitably approximating the volatility process risk-neutral density function. I exploit and adapt the idea, which stands behind popular techniques already employed in the context of equity options such as Edgeworth and Gram-Charlier expansions, of approximating the underlying process as a sum of some particular polynomials weighted by a kernel, which is typically a Gaussian distribution. I propose instead a Gamma kernel to adapt the methodology to the context of volatility options. VIX vanilla options closed-form pricing formulas are derived and their accuracy is tested for the Heston model (1993) as well as for the jump-diffusion SVJJ model proposed by Duffie et al. (2000).

Plasmonic nanoparticles as sensors : a novel approach to quantify adsorption and diffusion kinetics

Biosensors find wide application in clinical diagnostics, bioprocess control and environmental monitoring. They should not only show high specificity and reproducibility but also a high sensitivity and stability of the signal. Therefore, I introduce a novel sensor technology based on plasmonic nanoparticles which overcomes both of these limitations. Plasmonic nanoparticles exhibit strong absorption and scattering in the visible and near-infrared spectral range. The plasmon resonance, the collective coherent oscillation mode of the conduction band electrons against the positively charged ionic lattice, is sensitive to the local environment of the particle. I monitor these changes in the resonance wavelength by a new dark-field spectroscopy technique. Due to a strong light source and a highly sensitive detector a temporal resolution in the microsecond regime is possible in combination with a high spectral stability. This opens a window to investigate dynamics on the molecular level and to gain knowledge about fundamental biological processes.rnFirst, I investigate adsorption at the non-equilibrium as well as at the equilibrium state. I show the temporal evolution of single adsorption events of fibrinogen on the surface of the sensor on a millisecond timescale. Fibrinogen is a blood plasma protein with a unique shape that plays a central role in blood coagulation and is always involved in cell-biomaterial interactions. Further, I monitor equilibrium coverage fluctuations of sodium dodecyl sulfate and demonstrate a new approach to quantify the characteristic rate constants which is independent of mass transfer interference and long term drifts of the measured signal. This method has been investigated theoretically by Monte-Carlo simulations but so far there has been no sensor technology with a sufficient signal-to-noise ratio.rnSecond, I apply plasmonic nanoparticles as sensors for the determination of diffusion coefficients. Thereby, the sensing volume of a single, immobilized nanorod is used as detection volume. When a diffusing particle enters the detection volume a shift in the resonance wavelength is introduced. As no labeling of the analyte is necessary the hydrodynamic radius and thus the diffusion properties are not altered and can be studied in their natural form. In comparison to the conventional Fluorescence Correlation Spectroscopy technique a volume reduction by a factor of 5000-10000 is reached.

Laplace operator on finite graphs and a network diffusion model for the progression of the Alzheimer disease

Nella tesi viene descritto il Network Diffusion Model, ovvero il modello di A. Ray, A. Kuceyeski, M. Weiner inerente i meccanismi di progressione della demenza senile. In tale modello si approssima l'encefalo sano con una rete cerebrale (ovvero un grafo pesato), si identifica un generale fattore di malattia e se ne analizza la propagazione che avviene secondo meccanismi analoghi a quelli di un'infezione da prioni. La progressione del fattore di malattia e le conseguenze macroscopiche di tale processo(tra cui principalmente l'atrofia corticale) vengono, poi, descritte mediante approccio matematico. I risultati teoretici vengono confrontati con quanto osservato sperimentalmente in pazienti affetti da demenza senile. Nella tesi, inoltre, si fornisce una panoramica sui recenti studi inerenti i processi neurodegenerativi e si costruisce il contesto matematico di riferimento del modello preso in esame. Si presenta una panoramica sui grafi finiti, si introduce l'operatore di Laplace sui grafi e si forniscono stime dall'alto e dal basso per gli autovalori. Al fine di costruire una cornice matematica completa si analizza la relazione tra caso discreto e continuo: viene descritto l'operatore di Laplace-Beltrami sulle varietà riemanniane compatte e vengono fornite stime dall'alto per gli autovalori dell'operatore di Laplace-Beltrami associato a tali varietà a partire dalle stime dall'alto per gli autovalori del laplaciano sui grafi finiti.

LoneStar : Design and evaluation of an energy-efficient, disk-based archival storage system

Bandlaufwerke waren bisher die vorherrschende Technologie, um die anfallenden Datenmengen in Archivsystemen zu speichern. Mit Zugriffsmustern, die immer aktiver werden, und Speichermedien wie Festplatten die kostenmäßig aufholen, muss die Architektur vor Speichersystemen zur Archivierung neu überdacht werden. Zuverlässigkeit, Integrität und Haltbarkeit sind die Haupteigenschaften der digitalen Archivierung. Allerdings nimmt auch die Zugriffsgeschwindigkeit einen erhöhten Stellenwert ein, wenn aktive Archive ihre gesamten Inhalte für den direkten Zugriff bereitstellen. Ein band-basiertes System kann die hierfür benötigte Parallelität, Latenz und Durchsatz nicht liefern, was in der Regel durch festplattenbasierte Systeme als Zwischenspeicher kompensiert wird.rnIn dieser Arbeit untersuchen wir die Herausforderungen und Möglichkeiten ein festplattenbasiertes Speichersystem zu entwickeln, das auf eine hohe Zuverlässigkeit und Energieeffizienz zielt und das sich sowohl für aktive als auch für kalte Archivumgebungen eignet. Zuerst analysieren wir die Speichersysteme und Zugriffsmuster eines großen digitalen Archivs und präsentieren damit ein mögliches Einsatzgebiet für unsere Architektur. Daraufhin stellen wir Mechanismen vor um die Zuverlässigkeit einer einzelnen Festplatte zu verbessern und präsentieren sowie evaluieren einen neuen, energieeffizienten, zwei- dimensionalen RAID Ansatz der für „Schreibe ein Mal, lese mehrfach“ Zugriffe optimiert ist. Letztlich stellen wir Protokollierungs- und Zwischenspeichermechanismen vor, die die zugrundeliegenden Ziele unterstützen und evaluieren das RAID System in einer Dateisystemumgebung.

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.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.