Dynamik und Klimatologie von diabatischen Rossby-Wellen

Der Wintersturm Lothar zog am 26. Dezember 1999 über Europa und richtete in Frankreich, in Deutschland, in der Schweiz und in Österreich ungewöhnlich hohe Schäden an. Lothar entstand aus einer diabatischen Rossby Welle (DRW) und erreichte erst wenige Stunden vor dem europäischen Kontinent Orkanstärke. DRWs weisen ein interessantes atmosphärisches Strömungsmuster auf. Sie bestehen aus einer positiven PV-Anomalie in der unteren Troposphäre, die sich in einer Region mit starkem meridionalen Temperaturgradient befindet. Die positive PV-Anomalie löst eine zyklonale Strömung aus, dadurch wird östlich der PV-Anomalie warme Luft aus dem Süden herantransportiert. Während des Aufstieg der warmen Luft finden diabatische Prozesse statt, die zur Bildung einer neuen positiven PV-Anomalie in der unteren Troposphäre (PVA) führen. DRWs entstehen unabhängig von PV-Anomalien an der Tropopause. Falls sie jedoch mit ihnen in Wechselwirkung treten, kann - wie im Falle von Lothar - eine explosive Zyklogenese daraus resultieren. Im ersten Teil wird die Dynamik einer DRW am Beispiel des Wintersturms Lothar untersucht. Es wird insbesondere auf das Potential einer DRW zur explosiven Zyklogenese eingegangen. Im zweiten Teil wird das Aufretreten von DRWs in ECMWF-Vorhersagen untersucht. Es werden Unterschiede zwischen DRWs und anderen PV-Anomalien in der unteren Troposphäre hervorgehoben. Die Dynamik von DRWs wird mit Hilfe eines ECMWF-"Ensemble Prediction System" (EPS) des Wintersturms Lothar untersucht. Die 50 Modellläufe des EPS starten am 24. Dezember 1999 um 12 UTC und reichen bis zum 26. Dezember 1999 um 12 UTC. Nur 16 der 50 Modellläufe sagen einen ähnlich starken Sturm wie Lothar vorher. 10 Modellläufen sagen am 26. Dezember keine Zyklone mehr vorher. Die Ausprägung der baroklinen Zone, in der sich die DRW befindet, ist ausschlaggebend für die Intensität der DRW. Weitere wichtige Parameter sind der Feuchtegehalt der unteren Troposphäre und der latente Wärmefluss über dem Ozean. Diejenigen DRWs, die sich zu am 25. Dezember um 12 UTC näher als 400 km am Tropopausenjet befinden, entwickeln sich zu einer starken Zyklone. Alle anderen lösen sich auf oder bleiben schwache Zyklonen. Es ist schwierig, diabatische Prozesse in Wettervorhersagemodellen abzubilden, dementsprechend treten Schwierigkeiten bei der Vorhersage von PVAs auf. In den operationellen ECMWF-Vorhersagen von Juni 2004 bis Mai 2005 werden mit Hilfe eines Tracking- Algorithmus PVAs im Nordpazifik und Nordatlantik bestimmt und in fünf Kategorien eingeteilt. Die fünf Kategorien unterscheiden sich in ihrer Häufigkeit, ihrer Zugbahn und ihrer Gestalt. Im Nordpazifik entstehen doppelt so viele PVAs wie im Nordatlantik. Durchschnittlich werden im Winter weniger PVAs gefunden als im Sommer. Die Baroklinität und die Geschwindigkeit des Tropopausenjets ist in der Nähe von DRWs besonders hoch. Verglichen mit anderen PVAs weisen DRWs eine ähnliche Verteilung des reduzierten Bodendrucks auf. DRWs können in etwa gleich gut vorhergesagt werden wie andere PVAs.

Diskontinuierliche Galerkin-Verfahren für die operationelle Wettervorhersage

In dieser Arbeit wird ein neuer Dynamikkern entwickelt und in das bestehendernnumerische Wettervorhersagesystem COSMO integriert. Für die räumlichernDiskretisierung werden diskontinuierliche Galerkin-Verfahren (DG-Verfahren)rnverwendet, für die zeitliche Runge-Kutta-Verfahren. Hierdurch ist ein Verfahrenrnhoher Ordnung einfach zu realisieren und es sind lokale Erhaltungseigenschaftenrnder prognostischen Variablen gegeben. Der hier entwickelte Dynamikkern verwendetrngeländefolgende Koordinaten in Erhaltungsform für die Orographiemodellierung undrnkoppelt das DG-Verfahren mit einem Kessler-Schema für warmen Niederschlag. Dabeirnwird die Fallgeschwindigkeit des Regens, nicht wie üblich implizit imrnKessler-Schema diskretisiert, sondern explizit im Dynamikkern. Hierdurch sindrndie Zeitschritte der Parametrisierung für die Phasenumwandlung des Wassers undrnfür die Dynamik vollständig entkoppelt, wodurch auch sehr große Zeitschritte fürrndie Parametrisierung verwendet werden können. Die Kopplung ist sowohl fürrnOperatoraufteilung, als auch für Prozessaufteilung realisiert.rnrnAnhand idealisierter Testfälle werden die Konvergenz und die globalenrnErhaltungseigenschaften des neu entwickelten Dynamikkerns validiert. Die Massernwird bis auf Maschinengenauigkeit global erhalten. Mittels Bergüberströmungenrnwird die Orographiemodellierung validiert. Die verwendete Kombination ausrnDG-Verfahren und geländefolgenden Koordinaten ermöglicht die Behandlung vonrnsteileren Bergen, als dies mit dem auf Finite-Differenzenverfahren-basierendenrnDynamikkern von COSMO möglich ist. Es wird gezeigt, wann die vollernTensorproduktbasis und wann die Minimalbasis vorteilhaft ist. Die Größe desrnEinflusses auf das Simulationsergebnis der Verfahrensordnung, desrnParametrisierungszeitschritts und der Aufteilungsstrategie wirdrnuntersucht. Zuletzt wird gezeigt dass bei gleichem Zeitschritt die DG-Verfahrenrnaufgrund der besseren Skalierbarkeit in der Laufzeit konkurrenzfähig zurnFinite-Differenzenverfahren sind.

Optical interface based on a nanofiber atom-trap

In this thesis, I present the realization of a fiber-optical interface using optically trapped cesium atoms, which is an efficient tool for coupling light and atoms. The basic principle of the presented scheme relies on the trapping of neutral cesium atoms in a two-color evanescent field surrounding a nanofiber. The strong confinement of the fiber guided light, which also protrudes outside the nanofiber, provides strong confinement of the atoms as well as efficient coupling to near-resonant light propagating through the fiber. In chapter 1, the necessary physical and mathematical background describing the propagation of light in an optical fiber is presented. The exact solution of Maxwell’s equations allows us to model fiber-guided light fields which give rise to the trapping potentials and the atom-light coupling in the close vicinity of a nanofiber. Chapter 2 gives the theoretical background of light-atom interaction. A quantum mechanical model of the light-induced shifts of the relevant atomic levels is reviewed, which allows us to quantify the perturbation of the atomic states due to the presence of the trapping light-fields. The experimental realization of the fiber-based atom trap is the focus of chapter 3. Here, I analyze the properties of the fiber-based trap in terms of the confinement of the atoms and the impact of several heating mechanisms. Furthermore, I demonstrate the transportation of the trapped atoms, as a first step towards a deterministic delivery of individual atoms. In chapter 4, I present the successful interfacing of the trapped atomic ensemble and fiber-guided light. Three different approaches are discussed, i.e., those involving the measurement of either near-resonant scattering in absorption or the emission into the guided mode of the nanofiber. In the analysis of the spectroscopic properties of the trapped ensemble we find good agreement with the prediction of theoretical model discussed in chapter 2. In addition, I introduce a non-destructive scheme for the interrogation of the atoms states, which is sensitive to phase shifts of far-detuned fiber-guided light interacting with the trapped atoms. The inherent birefringence in our system, induced by the atoms, changes the state of polarization of the probe light and can be thus detected via a Stokes vector measurement.

Prediction of new crystal structures under extreme conditions

One of the most important challenges in chemistry and material science is the connection between the contents of a compound and its chemical and physical properties. In solids, these are greatly influenced by the crystal structure.rnrnThe prediction of hitherto unknown crystal structures with regard to external conditions like pressure and temperature is therefore one of the most important goals to achieve in theoretical chemistry. The stable structure of a compound is the global minimum of the potential energy surface, which is the high dimensional representation of the enthalpy of the investigated system with respect to its structural parameters. The fact that the complexity of the problem grows exponentially with the system size is the reason why it can only be solved via heuristic strategies.rnrnImprovements to the artificial bee colony method, where the local exploration of the potential energy surface is done by a high number of independent walkers, are developed and implemented. This results in an improved communication scheme between these walkers. This directs the search towards the most promising areas of the potential energy surface.rnrnThe minima hopping method uses short molecular dynamics simulations at elevated temperatures to direct the structure search from one local minimum of the potential energy surface to the next. A modification, where the local information around each minimum is extracted and used in an optimization of the search direction, is developed and implemented. Our method uses this local information to increase the probability of finding new, lower local minima. This leads to an enhanced performance in the global optimization algorithm.rnrnHydrogen is a highly relevant system, due to the possibility of finding a metallic phase and even superconductor with a high critical temperature. An application of a structure prediction method on SiH12 finds stable crystal structures in this material. Additionally, it becomes metallic at relatively low pressures.