Magnetic structure of Li2CuO2: From ab initio calculations to macroscopic simulations

The magnetic structure of the edge-sharing cuprate compound Li2CuO2 has been investigated with highly correlated ab initio electronic structure calculations. The first- and second-neighbor in-chain magnetic interactions are calculated to be 142 and -22 K, respectively. The ratio between the two parameters is smaller than suggested previously in the literature. The interchain interactions are antiferromagnetic in nature and of the order of a few K only. Monte Carlo simulations using the ab initio parameters to define the spin model Hamiltonian result in a Nel temperature in good agreement with experiment. Spin population analysis situates the magnetic moment on the copper and oxygen ions between the completely localized picture derived from experiment and the more delocalized picture based on local-density calculations.

Analysis of Stochastic Volatility Sequences Generated by Product Autoregressive Models

The classical methods of analysing time series by Box-Jenkins approach assume that the observed series uctuates around changing levels with constant variance. That is, the time series is assumed to be of homoscedastic nature. However, the nancial time series exhibits the presence of heteroscedasticity in the sense that, it possesses non-constant conditional variance given the past observations. So, the analysis of nancial time series, requires the modelling of such variances, which may depend on some time dependent factors or its own past values. This lead to introduction of several classes of models to study the behaviour of nancial time series. See Taylor (1986), Tsay (2005), Rachev et al. (2007). The class of models, used to describe the evolution of conditional variances is referred to as stochastic volatility modelsThe stochastic models available to analyse the conditional variances, are based on either normal or log-normal distributions. One of the objectives of the present study is to explore the possibility of employing some non-Gaussian distributions to model the volatility sequences and then study the behaviour of the resulting return series. This lead us to work on the related problem of statistical inference, which is the main contribution of the thesis

Physik der Elektronen in Nanostrukturen und niedrigdimensionalen Systemen - Studie zweier Beispiele

Als Beispiele für die vielfältigen Phänomene der Physik der Elektronen in niedrigdimensionalen Systemen wurden in dieser Arbeit das Cu(110)(2x1)O-Adsorbatsystem und die violette Li0.9Mo6O17-Bronze untersucht. Das Adsorbatsystem bildet selbstorganisierte quasi-eindimensionale Nanostrukturen auf einer Kupferoberfläche. Die Li-Bronze ist ein Material, das aufgrund seiner Kristallstruktur quasi-eindimensionale elektronische Eigenschaften im Volumen aufweist. Auf der Cu(110)(2x1)O-Oberfläche kann durch Variation der Sauerstoffbedeckung die Größe der streifenartigen CuO-Domänen geändert werden und damit der Übergang von zwei Dimensionen auf eine Dimension untersucht werden. Der Einfluss der Dimensionalität wurde anhand eines unbesetzten elektronischen Oberflächenzustandes studiert. Dessen Energieposition (untere Bandkante) verschiebt mit zunehmender Einschränkung (schmalere CuO-Streifen) zu größeren Energien hin. Dies ist ein bekannter quantenmechanischer Effekt und relativ gut verstanden. Zusätzlich wurde die Lebensdauer des Zustandes auf der voll bedeckten Oberfläche (zwei Dimensionen) ermittelt und deren Veränderung mit der Breite der CuO-Streifen untersucht. Es zeigt sich, dass die Lebensdauer auf schmaleren CuO-Streifen drastisch abnimmt. Dieses Ergebnis ist neu. Es kann im Rahmen eines Fabry-Perot-Modells als Streuung in Zustände außerhalb der CuO-Streifen verstanden werden. Außer den gerade beschriebenen Effekten war es möglich die Ladungsdichte des diskutierten Zustandes orts- und energieabhängig auf den CuO-Streifen zu studieren. Die Li0.9Mo6O17-Bronze wurde im Hinblick auf das Verhalten der elektronischen Zustandsdichte an der Fermikante untersucht. Diese Fragestellung ist besonders wegen der Quasieindimensionalität des Materials interessant. Die Messungen von STS-Spektren in der Nähe der Fermienergie zeigen, dass die Elektronen in der Li0.9Mo6O17-Bronze eine sogenannte Luttingerflüssigkeit ausbilden, die anstatt einer Fermiflüssigkeit in eindimensionalen elektronischen Systemen erwartet wird. Bisher wurde Luttingerflüssigkeitsverhalten erst bei wenigen Materialien und Systemen experimentell nachgewiesen, obschon die theoretischen Voraussagen mehr als 30 Jahre zurückliegen. Ein Charakteristikum einer Luttingerflüssigkeit ist die Abnahme der Zustandsdichte an der Fermienergie mit einem Potenzgesetz. Dieses Verhalten wurde in STS-Spektren dieser Arbeit beobachtet und quantitativ im Rahmen eines Luttingerflüssigkeitsmodells beschrieben. Auch die Temperaturabhängigkeit des Phänomens im Bereich von 5K bis 55K ist konsistent mit der Beschreibung durch eine Luttingerflüssigkeit. Generell zeigen die Untersuchungen dieser Arbeit, dass die Dimensionalität, insbesondere deren Einschränkung, einen deutlichen Einfluss auf die elektronischen Eigenschaften von Systemen und Materialien haben kann.

Intrinsic channel closing in strong-field single ionization of H2

Die Ionisation von H2 in intensiven Laserpulsen wird mit Hilfe der numerischen Integration der zeitabhängigen Schrödingergleichung für ein Einelektronenmodell untersucht, das die Vibrationsbewegung berücksichtigt. Die Spektren der kinetischen Elektronenenergie hängen stark von der Vibrationsquantenzahl des erzeugten H2+ Ions ab. Für bestimmte Vibrationszustände ist die Ausbeute der Elektronen in der Mitte des Plateaus stark erhöht. Der Effekt wird "channel closings" zugeschrieben, die in Atomen durch Variation der Laserintensität beobachtet wurden. The ionization of H2 in intense laser pulses is studied by numerical integration of the time-dependent Schrödinger equation for a single-active-electron model including the vibrational motion. The electron kinetic energy spectra in high-order above-threshold ionization are strongly dependent on the vibrational quantum number of the created H2+ ion. For certain vibrational states, the electron yield in the mid-plateau region is strongly enhanced. The effect is attributed to channel closings, which were previously observed in atoms by varying the laser intensity.

Relativistic ab initio calculations for ion-atom collisions

Within the independent particle model we solve the time-dependent single-particle equation using ab initio SCF-DIRAC-FOCK-SLATER wavefunctions as a basis. To reinstate the many-particle aspect of the collision system we use the inclusive probability formalism to answer experimental questions. As an example we show an application to the case of S{^15+} on Ar where experimental data on the K-K charge transfer are available for a wide range of impact energies from 4.7 to 90 MeV. Our molecular adiabatic calculations and the evaluation using the inclusive probability formalism show good results in the low energy range from 4.7 to 16 MeV impact energy.

Inclusive probabilities for the scattering system 16 MeV-S{^16+} on Ar

We performed ab initio calculations of many particle inclusive probabilities for the scattering system 16 MeV-S{^16+} on Ar. The solution of the time-dependent DIRAC-FOCK-SLATER-equation is achieved via a set of coupled-channel equations with energy eigenvalues and matrix elements which are given by static SCF molecular many electron calculations.

Multiple vacancy inclusive probabilities for the scattering system 16 MeV S{^16+} on Ar

To evaluate single and double K-shell inclusive charge transfer probabilities in ion-atom collisions we solve the time-dependent Dirac equation. By expanding the timedependent wavefunction in a set of molecular basis states the time-dependent equation reduces to a set of coupled-channel equations. The energy eigenvalues and matrix elements are taken from self-consistent relativistic molecular many-electron Dirac-Fock-Slater calculations. We present many-electron inclusive probabilities for different final configurations as a function of impact parameter for single and double K-shell vacancy production in collisions of bare S on Ar.

Tomographic imaging of molecular orbitals in length and velocity form

Recently Itatani et al. [Nature 432, 876 (2004)] introduced the new concept of molecular orbital tomography, where high harmonic generation (HHG) is used to image electronic wave functions. We describe an alternative reconstruction form, using momentum instead of dipole matrix elements for the electron recombination step in HHG. We show that using this velocity-form reconstruction, one obtains better results than using the original length-form reconstruction. We provide numerical evidence for our claim that one has to resort to extremely short pulses to perform the reconstruction for an orbital with arbitrary symmetry. The numerical evidence is based on the exact solution of the time-dependent Schrödinger equation for 2D model systems to simulate the experiment. Furthermore we show that in the case of cylindrically symmetric orbitals, such as the N2 orbital that was reconstructed in the original work, one can obtain the full 3D wave function and not only a 2D projection of it. Vor kurzem führten Itatani et al. [Nature 432, 876 (2004)] das Konzept der Molelkülorbital-Tomographie ein. Hierbei wird die Erzeugung hoher Harmonischer verwendet, um Bilder von elektronischen Wellenfunktionen zu gewinnen. Wir beschreiben eine alternative Form der Rekonstruktion, die auf Impuls- statt Dipol-Matrixelementen für den Rekombinationsschritt bei der Erzeugung der Harmonischen basiert. Wir zeigen, dass diese "Geschwindigkeitsform" der Rekonstruktion bessere Ergebnisse als die ursprüngliche "Längenform" liefert. Wir zeigen numerische Beweise für unsere Behauptung, dass man zu extrem kurzen Laserpulsen gehen muss, um Orbitale mit beliebiger Symmetrie zu rekonstruieren. Diese Ergebnisse basieren auf der exakten Lösung der zeitabhängigen Schrödingergleichung für 2D-Modellsysteme. Wir zeigen ferner, dass für zylindersymmetrische Orbitale wie das N2-Orbital, welches in der oben zitierten Arbeit rekonstruiert wurde, das volle 3D-Orbital rekonstruiert werden kann, nicht nur seine 2D-Projektion.

Realistic many electron coupled channel calculations in heavy ion scattering

The time dependent Dirac equation which describes a heavy ion-atom collision system is solved via a set of coupled channel equations with energy eigenvalues and matrix elements which are given by a selfconsistent field many electron calculation. After a brief discussion of the theoretical approximations and the connection of the many particle with the one particle interpretation we discuss first results for the systems F{^8+} - Ne and F{^6+} - Ne. The resulting P(b) curves for the creation of a Ne K-hole are in good agreement with the experimental results.

High laser field effects in multiphoton ionization of Na_2

Femtosecond pump/probe multiphoton ionization experiments on Na_2 molecules are performed. The dependence of the total Na^+_2 ion signal on the delay time and the intensity of the femtosecond laser pulses is studied in detail. It is observed that molecular vibrational wavepacket motion in different electronic states dominates the time dependence of the ion signal. For higher laser intensities the relative contributions from the A ^1 \summe^+_u and the 2 ^1 \produkt__g states change dramatically, indicating the increasing importance of a two-electron versus a one-electron process. For even stronger fields (10 ^12 W/ cm²) a vibrational wavepacket in the electronic ground state X ^1 \summe^+_g is formed and its dynamics is also observed in the transient Na^+_2 signal. Time-dependent quantum calculations are presented. The theoretical results agree well with the experiment.

Femtosecond pump-probe study of the spreading and recurrence of a vibrational wave packet in Na_2

The motion of a vibrational wave packet in the bound A(^1 \summe^+_u) electronic state of the sodium dimer is detected in a femtosecond pump/probe molecular beam experiment. For short times harmonic motion is seen in the total ion yield of Na^+_2 as a function of delay time between the two laser pulses. The spreading of the wave packet results in the loss of the periodic variation of the ion signal. For longer delay times (47 ps) the wave packet regains its initial form which is reflected in the revival structure of the Na^+_2 signal. Time-dependent quantum calculations reproduce the measured effects.

Dynamic characterization of olfactory receptors in the fruit fly Drosophila melanogaster

Temporal changes in odor concentration are vitally important to many animals orienting and navigating in their environment. How are such temporal changes detected? Within the scope of the present work an accurate stimulation and analysis system was developed to examine the dynamics of physiological properties of Drosophila melanogaster olfactory receptor organs. Subsequently a new method for delivering odor stimuli was tested and used to present the first dynamic characterization of olfactory receptors at the level of single neurons. Initially, recordings of the whole antenna were conducted while stimulating with different odors. The odor delivery system allowed the dynamic characterization of the whole fly antenna, including its sensilla and receptor neurons. Based on the obtained electroantennogram data a new odor delivery method called digital sequence method was developed. In addition the degree of accuracy was enhanced, initially using electroantennograms, and later recordings of odorant receptor cells at the single sensilla level. This work shows for the first time that different odors evoked different responses within one neuron depending on the chemical structure of the odor. The present work offers new insights into the dynamic properties of olfactory transduction in Drosophila melanogaster and describes time dependent parameters underlying these properties.

Das Einzugsgebiet der Fulda - Untersuchung der Methoden der Hochwasserstatistik im Zeichen des Klimawandels

Numerous studies have proven an effect of a probable climate change on the hydrosphere’s different subsystems. In the 21st century global and regional redistribution of water has to be expected and it is very likely that extreme weather phenomenon will occur more frequently. From a global view the flood situation will exacerbate. In contrast to these discoveries the classical approach of flood frequency analysis provides terms like “mean flood recurrence interval”. But for this analysis to be valid there is a need for the precondition of stationary distribution parameters which implies that the flood frequencies are constant in time. Newer approaches take into account extreme value distributions with time-dependent parameters. But the latter implies a discard of the mentioned old terminology that has been used up-to-date in engineering hydrology. On the regional scale climate change affects the hydrosphere in various ways. So, the question appears to be whether in central Europe the classical approach of flood frequency analysis is not usable anymore and whether the traditional terminology should be renewed. In the present case study hydro-meteorological time series of the Fulda catchment area (6930 km²), upstream of the gauging station Bonaforth, are analyzed for the time period 1960 to 2100. At first a distributed catchment area model (SWAT2005) is build up, calibrated and finally validated. The Edertal reservoir is regulated as well by a feedback control of the catchments output in case of low water. Due to this intricacy a special modeling strategy has been necessary: The study area is divided into three SWAT basin models and an additional physically-based reservoir model is developed. To further improve the streamflow predictions of the SWAT model, a correction by an artificial neural network (ANN) has been tested successfully which opens a new way to improve hydrological models. With this extension the calibration and validation of the SWAT model for the Fulda catchment area is improved significantly. After calibration of the model for the past 20th century observed streamflow, the SWAT model is driven by high resolution climate data of the regional model REMO using the IPCC scenarios A1B, A2, and B1, to generate future runoff time series for the 21th century for the various sub-basins in the study area. In a second step flood time series HQ(a) are derived from the 21st century runoff time series (scenarios A1B, A2, and B1). Then these flood projections are extensively tested with regard to stationarity, homogeneity and statistical independence. All these tests indicate that the SWAT-predicted 21st-century trends in the flood regime are not significant. Within the projected time the members of the flood time series are proven to be stationary and independent events. Hence, the classical stationary approach of flood frequency analysis can still be used within the Fulda catchment area, notwithstanding the fact that some regional climate change has been predicted using the IPCC scenarios. It should be noted, however, that the present results are not transferable to other catchment areas. Finally a new method is presented that enables the calculation of extreme flood statistics, even if the flood time series is non-stationary and also if the latter exhibits short- and longterm persistence. This method, which is called Flood Series Maximum Analysis here, enables the calculation of maximum design floods for a given risk- or safety level and time period.

Analyse selbsterregter Strömungsfelder im Verdichter-Leitrad durch Kopplung von Particle Image Velocimetry und Hitzdrahtanemometrie mittels Fast Fourier Transformation

Der Einsatz der Particle Image Velocimetry (PIV) zur Analyse selbsterregter Strömungsphänomene und das dafür notwendige Auswerteverfahren werden in dieser Arbeit beschrieben. Zur Untersuchung von solchen Mechanismen, die in Turbo-Verdichtern als Rotierende Instabilitäten in Erscheinung treten, wird auf Datensätze zurückgegriffen, die anhand experimenteller Untersuchungen an einem ringförmigen Verdichter-Leitrad gewonnen wurden. Die Rotierenden Instabilitäten sind zeitabhängige Strömungsphänomene, die bei hohen aerodynamischen Belastungen in Verdichtergittern auftreten können. Aufgrund der fehlenden Phaseninformation kann diese instationäre Strömung mit konventionellen PIV-Systemen nicht erfasst werden. Die Kármánsche Wirbelstraße und Rotierende Instabilitäten stellen beide selbsterregte Strömungsvorgänge dar. Die Ähnlichkeit wird genutzt um die Funktionalität des Verfahrens anhand der Kármánschen Wirbelstraße nachzuweisen. Der mittels PIV zu visualisierende Wirbeltransport erfordert ein besonderes Verfahren, da ein externes Signal zur Festlegung des Phasenwinkels dieser selbsterregten Strömung nicht zur Verfügung steht. Die Methodik basiert auf der Kopplung der PIV-Technik mit der Hitzdrahtanemometrie. Die gleichzeitige Messung mittels einer zeitlich hochaufgelösten Hitzdraht-Messung ermöglicht den Zeitpunkten der PIV-Bilder einen Phasenwinkel zuzuordnen. Hierzu wird das Hitzdrahtsignal mit einem FFT-Verfahren analysiert, um die PIV-Bilder entsprechend ihrer Phasenwinkel zu gruppieren. Dafür werden die aufgenommenen Bilder auf der Zeitachse der Hitzdrahtmessungen markiert. Eine systematische Analyse des Hitzdrahtsignals in der Umgebung der PIV-Messung liefert Daten zur Festlegung der Grundfrequenz und erlaubt es, der markierten PIV-Position einen Phasenwinkel zuzuordnen. Die sich aus den PIV-Bildern einer Klasse ergebenden Geschwindigkeitskomponenten werden anschließend gemittelt. Aus den resultierenden Bildern jeder Klasse ergibt sich das zweidimensionale zeitabhängige Geschwindigkeitsfeld, in dem die Wirbelwanderung der Kármánschen Wirbelstraße ersichtlich wird. In hierauf aufbauenden Untersuchungen werden Zeitsignale aus Messungen in einem Verdichterringgitter analysiert. Dabei zeigt sich, dass zusätzlich Filterfunktionen erforderlich sind. Im Ergebnis wird schließlich deutlich, dass die Übertragung der anhand der Kármánschen Wirbelstraße entwickelten Methode nur teilweise gelingt und weitere Forschungsarbeiten erforderlich sind.

Numerical Methods for the Unsteady Compressible Navier-Stokes Equations

We consider numerical methods for the compressible time dependent Navier-Stokes equations, discussing the spatial discretization by Finite Volume and Discontinuous Galerkin methods, the time integration by time adaptive implicit Runge-Kutta and Rosenbrock methods and the solution of the appearing nonlinear and linear equations systems by preconditioned Jacobian-Free Newton-Krylov, as well as Multigrid methods. As applications, thermal Fluid structure interaction and other unsteady flow problems are considered. The text is aimed at both mathematicians and engineers.