Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Molecular phylogeography of the European coastal plants Crithmum maritimum L., Halimione portulacoides (L.) Aellen, Salsola kali L. and Calystegia soldanella (L.) R. Br.

The intraspecific phylogeography of four European coastal plants, Crithmum maritimum, Halimione portulacoides, Salsola kali and Calystegia soldanella, was inferred from AFLP and ITS data. Only in C. maritimum, H. portulacoides and S. kali, a spatial genetic structure was revealed. The phylogeographic similarities and dissimilarities of these species include: (1) All three have distinct Black/Aegean and Adriatic Sea clusters. (2) Salsola kali and H. portulacoides show a distinct Atlantic/North Sea/Baltic Sea cluster, while Atlantic and eastern Spanish material of C. maritimum clustered together. (3) In the west Mediterranean, only S. kali forms a single cluster, while both H. portulacoides and C. maritimum display a phylogeographic break in the vicinity of the southern French coast. For S. kali, AFLP and ITS data concur in identifying separate Atlantic, east and west Mediterranean clades. All these patterns are postulated to result from both temperature changes during the last glacial and contemporary sea currents. No geographic AFLP structure was revealed in C. soldanella, both at the range-wide and population level. This was attributed to the remarkable seed dispersal ability of this species and possibly its longevity and clonal growth, preserving a random pattern of genetic variation generated by long-distance seed dispersal over long time periods.

Über das Vermessen tubulärer Strukturen in der Computertomographie

Die chronisch obstruktive Lungenerkrankung (engl. chronic obstructive pulmonary disease, COPD) ist ein Überbegriff für Erkrankungen, die zu Husten, Auswurf und Dyspnoe (Atemnot) in Ruhe oder Belastung führen - zu diesen werden die chronische Bronchitis und das Lungenemphysem gezählt. Das Fortschreiten der COPD ist eng verknüpft mit der Zunahme des Volumens der Wände kleiner Luftwege (Bronchien). Die hochauflösende Computertomographie (CT) gilt bei der Untersuchung der Morphologie der Lunge als Goldstandard (beste und zuverlässigste Methode in der Diagnostik). Möchte man Bronchien, eine in Annäherung tubuläre Struktur, in CT-Bildern vermessen, so stellt die geringe Größe der Bronchien im Vergleich zum Auflösungsvermögen eines klinischen Computertomographen ein großes Problem dar. In dieser Arbeit wird gezeigt wie aus konventionellen Röntgenaufnahmen CT-Bilder berechnet werden, wo die mathematischen und physikalischen Fehlerquellen im Bildentstehungsprozess liegen und wie man ein CT-System mittels Interpretation als lineares verschiebungsinvariantes System (engl. linear shift invariant systems, LSI System) mathematisch greifbar macht. Basierend auf der linearen Systemtheorie werden Möglichkeiten zur Beschreibung des Auflösungsvermögens bildgebender Verfahren hergeleitet. Es wird gezeigt wie man den Tracheobronchialbaum aus einem CT-Datensatz stabil segmentiert und mittels eines topologieerhaltenden 3-dimensionalen Skelettierungsalgorithmus in eine Skelettdarstellung und anschließend in einen kreisfreien Graphen überführt. Basierend auf der linearen System Theorie wird eine neue, vielversprechende, integral-basierte Methodik (IBM) zum Vermessen kleiner Strukturen in CT-Bildern vorgestellt. Zum Validieren der IBM-Resultate wurden verschiedene Messungen an einem Phantom, bestehend aus 10 unterschiedlichen Silikon Schläuchen, durchgeführt. Mit Hilfe der Skelett- und Graphendarstellung ist ein Vermessen des kompletten segmentierten Tracheobronchialbaums im 3-dimensionalen Raum möglich. Für 8 zweifach gescannte Schweine konnte eine gute Reproduzierbarkeit der IBM-Resultate nachgewiesen werden. In einer weiteren, mit IBM durchgeführten Studie konnte gezeigt werden, dass die durchschnittliche prozentuale Bronchialwandstärke in CT-Datensätzen von 16 Rauchern signifikant höher ist, als in Datensätzen von 15 Nichtrauchern. IBM läßt sich möglicherweise auch für Wanddickenbestimmungen bei Problemstellungen aus anderen Arbeitsgebieten benutzen - kann zumindest als Ideengeber dienen. Ein Artikel mit der Beschreibung der entwickelten Methodik und der damit erzielten Studienergebnisse wurde zur Publikation im Journal IEEE Transactions on Medical Imaging angenommen.

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions?

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions? The collapse of polyions in a poor solvent is a complex system and is an active research subject in the theoretical polyelectrolyte community. The complexity is due to the subtle interplay between hydrophobic effects, electrostatic interactions, entropy elasticity, intrinsic excluded volume as well as specific counter-ion and co-ion properties. Long range Coulomb forces can obscure single molecule properties. The here presented approach is to use just a small amount of screening salt in combination with a very high sample dilution in order to screen intermolecular interaction whereas keeping intramolecular interaction as much as possible (polyelectrolyte concentration cp ≤ 12 mg/L, salt concentration; Cs = 10^-5 mol/L). This is so far not described in literature. During collapse, the polyion is subject to a drastic change in size along with strong reduction of free counterions in solution. Therefore light scattering was utilized to obtain the size of the polyion whereas a conductivity setup was developed to monitor the proceeding of counterion collection by the polyion. Partially quaternized PVP’s below and above the Manning limit were investigated and compared to the collapse of their uncharged precursor. The collapses were induced by an isorefractive solvent/non-solvent mixture consisting of 1-propanol and 2-pentanone, with nearly constant dielectric constant. The solvent quality for the uncharged polyion could be quantified which, for the first time, allowed the experimental investigation of the effect of electrostatic interaction prior and during polyion collapse. Given that the Manning parameter M for QPVP4.3 is as low as lB / c = 0.6 (lB the Bjerrum length and c the mean contour distance between two charges), no counterion binding should occur. However the Walden product reduces with first addition of non solvent and accelerates when the structural collapse sets in. Since the dielectric constant of the solvent remains virtually constant during the chain collapse, the counterion binding is entirely caused by the reduction in the polyion chain dimension. The collapse is shifted to lower wns with higher degrees of quaternization as the samples QPVP20 and QPVP35 show (M = 2.8 respectively 4.9). The combination of light scattering and conductivity measurement revealed for the first time that polyion chains already collect their counter ions well above the theta-dimension when the dimensions start to shrink. Due to only small amounts of screening salt, strong electrostatic interactions bias dynamic as well as static light scattering measurements. An extended Zimm formula was derived to account for this interaction and to obtain the real chain dimensions. The effective degree of dissociation g could be obtained semi quantitatively using this extrapolated static in combination with conductivity measurements. One can conclude the expansion factor a and the effective degree of ionization of the polyion to be mutually dependent. In the good solvent regime g of QPVP4.3, QPVP20 and QPVP35 exhibited a decreasing value in the order 1 > g4.3 > g20 > g35. The low values of g for QPVP20 and QPVP35 are assumed to be responsible for the prior collapse of the higher quaternized samples. Collapse theory predicts dipole-dipole attraction to increase accordingly and even predicts a collapse in the good solvent regime. This could be exactly observed for the QPVP35 sample. The experimental results were compared to a newly developed theory of uniform spherical collapse induced by concomitant counterion binding developed by M. Muthukumar and A. Kundagrami. The theory agrees qualitatively with the location of the phase boundary as well as the trend of an increasing expansion with an increase of the degree of quaternization. However experimental determined g for the samples QPVP4.3, QPVP20 and QPVP35 decreases linearly with the degree of quaternization whereas this theory predicts an almost constant value.

Climate change in Libya and desertification of Jifara Plain

The study was arranged to manifest its objectives through preceding it with an intro-duction. Particular attention was paid in the second part to detect the physical settings of the study area, together with an attempt to show the climatic characteristics in Libya. In the third part, observed temporal and spatial climate change in Libya was investigated through the trends of temperature, precipitation, relative humidity and cloud amount over the peri-ods (1946-2000), (1946-1975), and (1976-2000), comparing the results with the global scales. The forth part detected the natural and human causes of climate change concentrat-ing on the greenhouse effect. The potential impacts of climate change on Libya were ex-amined in the fifth chapter. As a case study, desertification of Jifara Plain was studied in the sixth part. In the seventh chapter, projections and mitigations of climate change and desertification were discussed. Ultimately, the main results and recommendations of the study were summarized. In order to carry through the objectives outlined above, the following methods and approaches were used: a simple linear regression analysis was computed to detect the trends of climatic parameters over time; a trend test based on a trend-to-noise-ratio was applied for detecting linear or non-linear trends; the non-parametric Mann-Kendall test for trend was used to reveal the behavior of the trends and their significance; PCA was applied to construct the all-Libya climatic parameters trends; aridity index after Walter-Lieth was shown for computing humid respectively arid months in Libya; correlation coefficient, (after Pearson) for detecting the teleconnection between sun spot numbers, NAOI, SOI, GHGs, and global warming, climate changes in Libya; aridity index, after De Martonne, to elaborate the trends of aridity in Jifara Plain; Geographical Information System and Re-mote Sensing techniques were applied to clarify the illustrations and to monitor desertifi-cation of Jifara Plain using the available satellite images MSS, TM, ETM+ and Shuttle Radar Topography Mission (SRTM). The results are explained by 88 tables, 96 figures and 10 photos. Temporal and spatial temperature changes in Libya indicated remarkably different an-nual and seasonal trends over the long observation period 1946-2000 and the short obser-vation periods 1946-1975 and 1976-2000. Trends of mean annual temperature were posi-tive at all study stations except at one from 1946-2000, negative trends prevailed at most stations from 1946-1975, while strongly positive trends were computed at all study stations from 1976-2000 corresponding with the global warming trend. Positive trends of mean minimum temperatures were observed at all reference stations from 1946-2000 and 1976-2000, while negative trends prevailed at most stations over the period 1946-1975. For mean maximum temperature, positive trends were shown from 1946-2000 and from 1976-2000 at most stations, while most trends were negative from 1946-1975. Minimum tem-peratures increased at nearly more than twice the rate of maximum temperatures at most stations. In respect of seasonal temperature, warming mostly occurred in summer and au-tumn in contrast to the global observations identifying warming mostly in winter and spring in both study periods. Precipitation across Libya is characterized by scanty and sporadically totals, as well as high intensities and very high spatial and temporal variabilities. From 1946-2000, large inter-annual and intra-annual variabilities were observed. Positive trends of annual precipi-tation totals have been observed from 1946-2000, negative trends from 1976-2000 at most stations. Variabilities of seasonal precipitation over Libya are more strikingly experienced from 1976-2000 than from 1951-1975 indicating a growing magnitude of climate change in more recent times. Negative trends of mean annual relative humidity were computed at eight stations, while positive trends prevailed at seven stations from 1946-2000. For the short observation period 1976-2000, positive trends were computed at most stations. Annual cloud amount totals decreased at most study stations in Libya over both long and short periods. Re-markably large spatial variations of climate changes were observed from north to south over Libya. Causes of climate change were discussed showing high correlation between tempera-ture increasing over Libya and CO2 emissions; weakly positive correlation between pre-cipitation and North Atlantic Oscillation index; negative correlation between temperature and sunspot numbers; negative correlation between precipitation over Libya and Southern Oscillation Index. The years 1992 and 1993 were shown as the coldest in the 1990s result-ing from the eruption of Mount Pinatubo, 1991. Libya is affected by climate change in many ways, in particular, crop production and food security, water resources, human health, population settlement and biodiversity. But the effects of climate change depend on its magnitude and the rate with which it occurs. Jifara Plain, located in northwestern Libya, has been seriously exposed to desertifica-tion as a result of climate change, landforms, overgrazing, over-cultivation and population growth. Soils have been degraded, vegetation cover disappeared and the groundwater wells were getting dry in many parts. The effect of desertification on Jifara Plain appears through reducing soil fertility and crop productivity, leading to long-term declines in agri-cultural yields, livestock yields, plant standing biomass, and plant biodiversity. Desertifi-cation has also significant implications on livestock industry and the national economy. Desertification accelerates migration from rural and nomadic areas to urban areas as the land cannot support the original inhabitants. In the absence of major shifts in policy, economic growth, energy prices, and con-sumer trends, climate change in Libya and desertification of Jifara Plain are expected to continue in the future. Libya cooperated with United Nations and other international organizations. It has signed and ratified a number of international and regional agreements which effectively established a policy framework for actions to mitigate climate change and combat deserti-fication. Libya has implemented several laws and legislative acts, with a number of ancil-lary and supplementary rules to regulate. Despite the current efforts and ongoing projects being undertaken in Libya in the field of climate change and desertification, urgent actions and projects are needed to mitigate climate change and combat desertification in the near future.